MSM: a new flexible framework for Multimodal Surface Matching

A spatiotemporal complexity architecture of human brain activity

The human brain operates in large-scale functional networks. These networks are an expression of temporally correlated activity across brain regions, but how global network properties relate to the neural dynamics of individual regions remains incompletely understood. Here, we show that the brain's network architecture is tightly linked to critical episodes of neural regularity, visible as spontaneous "complexity drops" in functional magnetic resonance imaging signals. These episodes closely explain functional connectivity strength between regions, subserve the propagation of neural activity patterns, and reflect interindividual differences in age and behavior. Furthermore, complexity drops define neural activity states that dynamically shape the connectivity strength, topological configuration, and hierarchy of brain networks and comprehensively explain known structure-function relationships within the brain. These findings delineate a principled complexity architecture of neural activity-a human "complexome" that underpins the brain's functional network organization.

Individual-level functional connectomes predict the motor symptoms of Parkinson's disease

Abnormalities in functional connectivity networks are associated with sensorimotor networks in Parkinson's disease (PD) based on group-level mapping studies, but these results are controversial. Using individual-level cortical segmentation to construct individual brain atlases can supplement the individual information covered by group-level cortical segmentation. Functional connectivity analyses at the individual level are helpful for obtaining clinically useful markers and predicting treatment response. Based on the functional connectivity of individualized regions of interest, a support vector regression model was trained to estimate the severity of motor symptoms for each subject, and a correlation analysis between the estimated scores and clinical symptom scores was performed. Forty-six PD patients aged 50-75 years were included from the Parkinson's Progression Markers Initiative database, and 63 PD patients were included from the Beijing Rehabilitation Hospital database. Only patients below Hoehn and Yahr stage III were included. The analysis showed that the severity of motor symptoms could be estimated by the individualized functional connectivity between the visual network and sensorimotor network in early-stage disease. The results reveal individual-level connectivity biomarkers related to motor symptoms and emphasize the importance of individual differences in the prediction of the treatment response of PD.

Functionally and structurally distinct fusiform face area(s) in over 1000 participants

The fusiform face area (FFA) is a widely studied region causally involved in face perception. Even though cognitive neuroscientists have been studying the FFA for over two decades, answers to foundational questions regarding the function, architecture, and connectivity of the FFA from a large (N>1000) group of participants are still lacking. To fill this gap in knowledge, we quantified these multimodal features of fusiform face-selective regions in 1053 participants in the Human Connectome Project. After manually defining over 4,000 fusiform face-selective regions, we report five main findings. First, 68.76% of hemispheres have two cortically separate regions (pFus-faces/FFA-1 and mFus-faces/FFA-2). Second, in 26.69% of hemispheres, pFus-faces/FFA-1 and mFus-faces/FFA-2 are spatially contiguous, yet are distinct based on functional, architectural, and connectivity metrics. Third, pFus-faces/FFA-1 is more face-selective than mFus-faces/FFA-2, and the two regions have distinct functional connectivity fingerprints. Fourth, pFus-faces/FFA-1 is cortically thinner and more heavily myelinated than mFus-faces/FFA-2. Fifth, face-selective patterns and functional connectivity fingerprints of each region are more similar in monozygotic than dizygotic twins and more so than architectural gradients. As we share our areal definitions with the field, future studies can explore how structural and functional features of these regions will inform theories regarding how visual categories are represented in the brain.

Hierarchical particle optimization for cortical shape correspondence in temporal lobe resection

BACKGROUND

Anterior temporal lobe resection is an effective treatment for temporal lobe epilepsy. The post-surgical structural changes could influence the follow-up treatment. Capturing post-surgical changes necessitates a well-established cortical shape correspondence between pre- and post-surgical surfaces. Yet, most cortical surface registration methods are designed for normal neuroanatomy. Surgical changes can introduce wide ranging artifacts in correspondence, for which conventional surface registration methods may not work as intended.

METHODS

In this paper, we propose a novel particle method for one-to-one dense shape correspondence between pre- and post-surgical surfaces with temporal lobe resection. The proposed method can handle partial structural abnormality involving non-rigid changes. Unlike existing particle methods using implicit particle adjacency, we consider explicit particle adjacency to establish a smooth correspondence. Moreover, we propose hierarchical optimization of particles rather than full optimization of all particles at once to avoid trappings of locally optimal particle update.

RESULTS

We evaluate the proposed method on 25 pairs of T1-MRI with pre- and post-simulated resection on the anterior temporal lobe and 25 pairs of patients with actual resection. We show improved accuracy over several cortical regions in terms of ROI boundary Hausdorff distance with 4.29 mm and Dice similarity coefficients with average value 0.841, compared to existing surface registration methods on simulated data. In 25 patients with actual resection of the anterior temporal lobe, our method shows an improved shape correspondence in qualitative and quantitative evaluation on parcellation-off ratio with average value 0.061 and cortical thickness changes. We also show better smoothness of the correspondence without self-intersection, compared with point-wise matching methods which show various degrees of self-intersection.

CONCLUSION

The proposed method establishes a promising one-to-one dense shape correspondence for temporal lobe resection. The resulting correspondence is smooth without self-intersection. The proposed hierarchical optimization strategy could accelerate optimization and improve the optimization accuracy. According to the results on the paired surfaces with temporal lobe resection, the proposed method outperforms the compared methods and is more reliable to capture cortical thickness changes.

Functionally and structurally distinct fusiform face area(s) in over 1000 participants

The fusiform face area (FFA) is a widely studied region causally involved in face perception. Even though cognitive neuroscientists have been studying the FFA for over two decades, answers to foundational questions regarding the function, architecture, and connectivity of the FFA from a large (N>1000) group of participants are still lacking. To fill this gap in knowledge, we quantified these multimodal features of fusiform face-selective regions in 1053 participants in the Human Connectome Project. After manually defining over 4,000 fusiform face-selective regions, we report five main findings. First, 68.76% of hemispheres have two cortically separate regions (pFus-faces/FFA-1 and mFus-faces/FFA-2). Second, in 26.69% of hemispheres, pFus-faces/FFA-1 and mFus-faces/FFA-2 are spatially contiguous, yet are distinct based on functional, architectural, and connectivity metrics. Third, pFus-faces/FFA-1 is more face-selective than mFus-faces/FFA-2, and the two regions have distinct functional connectivity fingerprints. Fourth, pFus-faces/FFA-1 is cortically thinner and more heavily myelinated than mFus-faces/FFA-2. Fifth, face-selective patterns and functional connectivity fingerprints of each region are more similar in monozygotic than dizygotic twins and more so than architectural gradients. As we share our areal definitions with the field, future studies can explore how structural and functional features of these regions will inform theories regarding how visual categories are represented in the brain.

The subcortical and neurochemical organization of the ventral and dorsal attention networks

Attention is a core cognitive function that filters and selects behaviourally relevant information in the environment. The cortical mapping of attentional systems identified two segregated networks that mediate stimulus-driven and goal-driven processes, the Ventral and the Dorsal Attention Networks (VAN, DAN). Deep brain electrophysiological recordings, behavioral data from phylogenetic distant species, and observations from human brain pathologies challenge purely corticocentric models. Here, we used advanced methods of functional alignment applied to resting-state functional connectivity analyses to map the subcortical architecture of the Ventral and Dorsal Attention Networks. Our investigations revealed the involvement of the pulvinar, the superior colliculi, the head of caudate nuclei, and a cluster of brainstem nuclei relevant to both networks. These nuclei are densely connected structural network hubs, as revealed by diffusion-weighted imaging tractography. Their projections establish interrelations with the acetylcholine nicotinic receptor as well as dopamine and serotonin transporters, as demonstrated in a spatial correlation analysis with a normative atlas of neurotransmitter systems. This convergence of functional, structural, and neurochemical evidence provides a comprehensive framework to understand the neural basis of attention across different species and brain diseases.

Hypothalamic interaction with reward-related regions during subjective evaluation of foods

The reward system implemented in the midbrain, ventral striatum, orbitofrontal cortex, and ventromedial prefrontal cortex evaluates and compares various types of rewards given to the organisms. It has been suggested that autonomic factors influence reward-related processing via the hypothalamus, but how the hypothalamus modulates the reward system remains elusive. In this functional magnetic resonance imaging study, the hypothalamus was parcellated into individual hypothalamic nuclei performing different autonomic functions using boundary mapping parcellation analyses. The effective interaction during subjective evaluation of foods in a reward task was then investigated between the human hypothalamic nuclei and the reward-related regions. We found significant brain activity decrease in the paraventricular nucleus (PVH) and lateral nucleus in the hypothalamus in food evaluation compared with monetary evaluation. A psychophysiological interaction analysis revealed dual interactions between the PVH and (1) midbrain region and (2) ventromedial prefrontal cortex, with the former correlated with the stronger tendency of participants toward food-seeking. A dynamic causal modeling analysis further revealed unidirectional interactions from the PVH to the midbrain and ventromedial prefrontal cortex. These results suggest that the PVH in the human hypothalamus interacts with the reward-related regions in the cerebral cortex via multiple pathways (i.e., the midbrain pathway and ventromedial prefrontal pathway) to evaluate rewards for subsequent decision-making.

Neural function underlying reward expectancy and attainment in adolescents with diverse psychiatric symptoms

Reward dysfunction has been hypothesized to play a key role in the development of psychiatric conditions during adolescence. To help capture the complexity of reward function in youth, we used the Reward Flanker fMRI Task, which enabled us to examine neural activity during expectancy and attainment of both certain and uncertain rewards. Participants were 84 psychotropic-medication-free adolescents, including 67 with diverse psychiatric conditions and 17 healthy controls. Functional MRI used high-resolution acquisition and high-fidelity processing techniques modeled after the Human Connectome Project. Analyses examined neural activation during reward expectancy and attainment, and their associations with clinical measures of depression, anxiety, and anhedonia severity, with results controlled for family-wise errors using non-parametric permutation tests. As anticipated, reward expectancy activated regions within the fronto-striatal reward network, thalamus, occipital lobe, superior parietal lobule, temporoparietal junction, and cerebellum. Unexpectedly, however, reward attainment was marked by widespread deactivation in many of these same regions, which we further explored using cosine similarity analysis. Across all subjects, striatum and thalamus activation during reward expectancy negatively correlated with anxiety severity, while activation in numerous cortical and subcortical regions during reward attainment positively correlated with both anxiety and depression severity. These findings highlight the complexity and dynamic nature of neural reward processing in youth.

neuromaps: structural and functional interpretation of brain maps

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Comparing experimentally generated maps to these reference maps facilitates cross-disciplinary scientific discovery. Although recent data sharing initiatives increase the accessibility of brain maps, data are often shared in disparate coordinate systems, precluding systematic and accurate comparisons. Here we introduce neuromaps, a toolbox for accessing, transforming and analyzing structural and functional brain annotations. We implement functionalities for generating high-quality transformations between four standard coordinate systems. The toolbox includes curated reference maps and biological ontologies of the human brain, such as molecular, microstructural, electrophysiological, developmental and functional ontologies. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. neuromaps combines open-access data with transparent functionality for standardizing and comparing brain maps, providing a systematic workflow for comprehensive structural and functional annotation enrichment analysis of the human brain.

Cortical gradient of a human functional similarity network captured by the geometry of cytoarchitectonic organization

Mapping the functional topology from a multifaceted perspective and relating it to underlying cross-scale structural principles is crucial for understanding the structural-functional relationships of the cerebral cortex. Previous works have described a sensory-association gradient axis in terms of coupling relationships between structure and function, but largely based on single specific feature, and the mesoscopic underpinnings are rarely determined. Here we show a gradient pattern encoded in a functional similarity network based on data from Human Connectome Project and further link it to cytoarchitectonic organizing principles. The spatial distribution of the primary gradient follows an inferior-anterior to superior-posterior axis. The primary gradient demonstrates converging relationships with layer-specific microscopic gene expression and mesoscopic cortical layer thickness, and is captured by the geometric representation of a myelo- and cyto-architecture based laminar differentiation theorem, involving a dual origin theory. Together, these findings provide a gradient, which describes the functional topology, and more importantly, linking the macroscale functional landscape with mesoscale laminar differentiation principles.

Evidence of familial confounding of the association between cannabis use and cerebellar-cortical functional connectivity using a twin study

Cerebellar-cortical resting-state functional connectivity (rsFC) has been reported to be altered in cannabis users. However, this association may be due to genetic and environmental confounding rather than a causal relationship between cannabis use and changes in rsFC. In this co-twin control study, linear mixed models were used to assess relationships between the number of lifetime cannabis uses (NLCU) and age of cannabis onset (ACO) with cerebellar-cortical rsFC. The rsFC with seven functional networks was evaluated in 147 monozygotic and 82 dizygotic twin pairs. Importantly, the use of genetically informed models in this twin sample facilitated examining whether shared genetic or environmental effects underlie crude associations between cannabis measures and connectivity. Individual-level phenotypic analyses (i.e., accounting for twin-pair non-independence) showed that individuals in the full sample with earlier ACO and higher NLCU had lower cerebellar rsFC within the VA, DA, and FP networks. Yet, there were no significant differences in cerebellar-cortical rsFC between monozygotic twins who were discordant for cannabis measures. These findings suggest shared genetic or environmental confounds contribute to associations between cannabis use and altered cerebellar-cortical rsFC, rather than unique causal impacts of cannabis use on cerebellar-cortical rsFC.

Surface-based tracking for short association fibre tractography

It is estimated that in the human brain, short association fibres (SAF) represent more than half of the total white matter volume and their involvement has been implicated in a range of neurological and psychiatric conditions. This population of fibres, however, remains relatively understudied in the neuroimaging literature. Some of the challenges pertinent to the mapping of SAF include their variable anatomical course and proximity to the cortical mantle, leading to partial volume effects and potentially affecting streamline trajectory estimation. This work considers the impact of seeding and filtering strategies and choice of scanner, acquisition, data resampling to propose a whole-brain, surface-based short (≤30-40 mm) SAF tractography approach. The framework is shown to produce longer streamlines with a predilection for connecting gyri as well as high cortical coverage. We further demonstrate that certain areas of subcortical white matter become disproportionally underrepresented in diffusion-weighted MRI data with lower angular and spatial resolution and weaker diffusion weighting; however, collecting data with stronger gradients than are usually available clinically has minimal impact, making our framework translatable to data collected on commonly available hardware. Finally, the tractograms are examined using voxel- and surface-based measures of consistency, demonstrating moderate reliability, low repeatability and high between-subject variability, urging caution when streamline count-based analyses of SAF are performed.

Altered brain responses to noxious dentoalveolar stimuli in high-impact temporomandibular disorder pain patients

High-impact temporomandibular disorder (TMD) pain may involve brain mechanisms related to maladaptive central pain modulation. We investigated brain responses to stimulation of trigeminal sites not typically associated with TMD pain by applying noxious dentoalveolar pressure to high- and low-impact TMD pain cases and pain-free controls during functional magnetic resonance imaging (fMRI). Fifty female participants were recruited and assigned to one of three groups based on the Diagnostic Criteria for Temporomandibular Disorders (DC/TMD) and Graded Chronic Pain Scale: controls (n = 17), low-impact (n = 17) and high-impact TMD (n = 16). Multimodal whole-brain MRI was acquired following the Human Connectome Project Lifespan protocol, including stimulus-evoked fMRI scans during which painful dentoalveolar pressure was applied to the buccal gingiva of participants. Group analyses were performed using non-parametric permutation tests for parcellated cortical and subcortical neuroimaging data. There were no significant between-group differences for brain activations/deactivations evoked by the noxious dentoalveolar pressure. For individual group mean activations/deactivations, a gradient in the number of parcels surviving thresholding was found according to the TMD pain grade, with the highest number seen in the high-impact group. Among the brain regions activated in chronic TMD pain groups were those previously implicated in sensory-discriminative and motivational-affective pain processing. These results suggest that dentoalveolar pressure pain evokes abnormal brain responses to sensory processing of noxious stimuli in high-impact TMD pain participants, which supports the presence of maladaptive brain plasticity in chronic TMD pain.

Modulation of the functional connectome in major depressive disorder by ketamine therapy

BACKGROUND

Subanesthetic ketamine infusion therapy can produce fast-acting antidepressant effects in patients with major depression. How single and repeated ketamine treatment modulates the whole-brain functional connectome to affect clinical outcomes remains uncharacterized.

METHODS

Data-driven whole brain functional connectivity (FC) analysis was used to identify the functional connections modified by ketamine treatment in patients with major depressive disorder (MDD). MDD patients (N = 61, mean age = 38, 19 women) completed baseline resting-state (RS) functional magnetic resonance imaging and depression symptom scales. Of these patients, n = 48 and n = 51, completed the same assessments 24 h after receiving one and four 0.5 mg/kg intravenous ketamine infusions. Healthy controls (HC) (n = 40, 24 women) completed baseline assessments with no intervention. Analysis of RS FC addressed effects of diagnosis, time, and remitter status.

RESULTS

Significant differences (p < 0.05, corrected) in RS FC were observed between HC and MDD at baseline in the somatomotor network and between association and default mode networks. These disruptions in FC in MDD patients trended toward control patterns with ketamine treatment. Furthermore, following serial ketamine infusions, significant decreases in FC were observed between the cerebellum and salience network (SN) (p < 0.05, corrected). Patient remitters showed increased FC between the cerebellum and the striatum prior to treatment that decreased following treatment, whereas non-remitters showed the opposite pattern.

CONCLUSION

Results support that ketamine treatment leads to neurofunctional plasticity between distinct neural networks that are shown as disrupted in MDD patients. Cortico-striatal-cerebellar loops that encompass the SN could be a potential biomarker for ketamine treatment.

Developmental trajectories of cortical thickness by functional brain network: The roles of pubertal timing and socioeconomic status

The human cerebral cortex undergoes considerable changes during development, with cortical maturation patterns reflecting regional heterogeneity that generally progresses in a posterior-to-anterior fashion. However, the organizing principles that govern cortical development remain unclear. In the current study, we characterized age-related differences in cortical thickness (CT) as a function of sex, pubertal timing, and two dissociable indices of socioeconomic status (i.e., income-to-needs and maternal education) in the context of functional brain network organization, using a cross-sectional sample (n = 789) diverse in race, ethnicity, and socioeconomic status from the Lifespan Human Connectome Project in Development (HCP-D). We found that CT generally followed a linear decline from 5 to 21 years of age, except for three functional networks that displayed nonlinear trajectories. We found no main effect of sex or age by sex interaction for any network. Earlier pubertal timing was associated with reduced mean CT and CT in seven networks. We also found a significant age by maternal education interaction for mean CT across cortex and CT in the dorsal attention network, where higher levels of maternal education were associated with steeper age-related decreases in CT. Taken together, our results suggest that these biological and environmental variations may impact the emerging functional connectome.

Accurate predictions of individual differences in task-evoked brain activity from resting-state fMRI using a sparse ensemble learner

Modelling and predicting individual differences in task-fMRI activity can have a wide range of applications from basic to clinical neuroscience. It has been shown that models based on resting-state activity can have high predictive accuracy. Here we propose several improvements to such models. Using a sparse ensemble learner, we show that (i) features extracted using Stochastic Probabilistic Functional Modes (sPROFUMO) outperform the previously proposed dual-regression approach, (ii) that the shape and overall intensity of individualised task activations can be modelled separately and explicitly, (iii) training the model on predicting residual differences in brain activity further boosts individualised predictions. These results hold for both surface-based analyses of the Human Connectome Project data as well as volumetric analyses of UK-biobank data. Overall, our model achieves state of the art prediction accuracy on par with the test-retest reliability of task-fMRI scans, suggesting that it has potential to supplement traditional task localisers.

An evaluation of how connectopic mapping reveals visual field maps in V1

Abstract

Functional gradients, in which response properties change gradually across the cortical surface, have been proposed as a key organising principle of the brain. However, the presence of these gradients remains undetermined in many brain regions. Resting-state neuroimaging studies have suggested these gradients can be reconstructed from patterns of functional connectivity. Here we investigate the accuracy of these reconstructions and establish whether it is connectivity or the functional properties within a region that determine these “connectopic maps”. Different manifold learning techniques were used to recover visual field maps while participants were at rest or engaged in natural viewing. We benchmarked these reconstructions against maps measured by traditional visual field mapping. We report an initial exploratory experiment of a publicly available naturalistic imaging dataset, followed by a preregistered replication using larger resting-state and naturalistic imaging datasets from the Human Connectome Project. Connectopic mapping accurately predicted visual field maps in primary visual cortex, with better predictions for eccentricity than polar angle maps. Non-linear manifold learning methods outperformed simpler linear embeddings. We also found more accurate predictions during natural viewing compared to resting-state. Varying the source of the connectivity estimates had minimal impact on the connectopic maps, suggesting the key factor is the functional topography within a brain region. The application of these standardised methods for connectopic mapping will allow the discovery of functional gradients across the brain.

Protocol registration

The stage 1 protocol for this Registered Report was accepted in principle on 19 April 2022. The protocol, as accepted by the journal, can be found at 10.6084/m9.figshare.19771717.

Attentional M100 gain modulation localizes to auditory sensory cortex and is deficient in first-episode psychosis

Selective attention is impaired in first-episode psychosis (FEP). Selective attention effects can be detected during auditory tasks as increased sensory activity. We previously reported electroencephalography scalp-measured N100 enhancement is reduced in FEP. Here, we localized magnetoencephalography (MEG) M100 source activity within the auditory cortex, making novel use of the Human Connectome Project multimodal parcellation (HCP-MMP) to identify precise auditory cortical areas involved in attention modulation and its impairment in FEP. MEG was recorded from 27 FEP and 31 matched healthy controls (HC) while individuals either ignored frequent standard and rare oddball tones while watching a silent movie or attended tones by pressing a button to oddballs. Because M100 arises mainly in the auditory cortices, MEG activity during the M100 interval was projected to the auditory sensory cortices defined by the HCP-MMP (A1, lateral belt, and parabelt parcels). FEP had less auditory sensory cortex M100 activity in both conditions. In addition, there was a significant interaction between group and attention. HC enhanced source activity with attention, but FEP did not. These results demonstrate deficits in both sensory processing and attentional modulation of the M100 in FEP. Novel use of the HCP-MMP revealed the precise cortical areas underlying attention modulation of auditory sensory activity in healthy individuals and impairments in FEP. The sensory reduction and attention modulation impairment indicate local and systems-level pathophysiology proximal to disease onset that may be critical for etiology. Further, M100 and N100 enhancement may serve as outcome variables for targeted intervention to improve attention in early psychosis.

Physical distance to sensory-motor landmarks predicts language function

Auditory language comprehension recruits cortical regions that are both close to sensory-motor landmarks (supporting auditory and motor features) and far from these landmarks (supporting word meaning). We investigated whether the responsiveness of these regions in task-based functional MRI is related to individual differences in their physical distance to primary sensorimotor landmarks. Parcels in the auditory network, that were equally responsive across story and math tasks, showed stronger activation in individuals who had less distance between these parcels and transverse temporal sulcus, in line with the predictions of the "tethering hypothesis," which suggests that greater proximity to input regions might increase the fidelity of sensory processing. Conversely, language and default mode parcels, which were more active for the story task, showed positive correlations between individual differences in activation and sensory-motor distance from primary sensory-motor landmarks, consistent with the view that physical separation from sensory-motor inputs supports aspects of cognition that draw on semantic memory. These results demonstrate that distance from sensorimotor regions provides an organizing principle of functional differentiation within the cortex. The relationship between activation and geodesic distance to sensory-motor landmarks is in opposite directions for cortical regions that are proximal to the heteromodal (DMN and language network) and unimodal ends of the principal gradient of intrinsic connectivity.

Anterior default mode network and posterior insular connectivity is predictive of depressive symptom reduction following serial ketamine infusion

BACKGROUND

Ketamine is a rapidly-acting antidepressant treatment with robust response rates. Previous studies have reported that serial ketamine therapy modulates resting state functional connectivity in several large-scale networks, though it remains unknown whether variations in brain structure, function, and connectivity impact subsequent treatment success. We used a data-driven approach to determine whether pretreatment multimodal neuroimaging measures predict changes along symptom dimensions of depression following serial ketamine infusion.

METHODS

Patients with depression (n = 60) received structural, resting state functional, and diffusion MRI scans before treatment. Depressive symptoms were assessed using the 17-item Hamilton Depression Rating Scale (HDRS-17), the Inventory of Depressive Symptomatology (IDS-C), and the Rumination Response Scale (RRS) before and 24 h after patients received four (0.5 mg/kg) infusions of racemic ketamine over 2 weeks. Nineteen unaffected controls were assessed at similar timepoints. Random forest regression models predicted symptom changes using pretreatment multimodal neuroimaging and demographic measures.

RESULTS

Two HDRS-17 subscales, the HDRS-6 and core mood and anhedonia (CMA) symptoms, and the RRS: reflection (RRSR) scale were predicted significantly with 19, 27, and 1% variance explained, respectively. Increased right medial prefrontal cortex/anterior cingulate and posterior insula (PoI) and lower kurtosis of the superior longitudinal fasciculus predicted reduced HDRS-6 and CMA symptoms following treatment. RRSR change was predicted by global connectivity of the left posterior cingulate, left insula, and right superior parietal lobule.

CONCLUSIONS

Our findings support that connectivity of the anterior default mode network and PoI may serve as potential biomarkers of antidepressant outcomes for core depressive symptoms.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Person-specific and precision neuroimaging: Current methods and future directions

Most neuroimaging studies of brain function analyze data in normalized space to identify regions of common activation across participants. These studies treat interindividual differences in brain organization as noise, but this approach can obscure important information about the brain's functional architecture. Recently, a number of studies have adopted a person-specific approach that aims to characterize these individual differences and explore their reliability and implications for behavior. A subset of these studies has taken a precision imaging approach that collects multiple hours of data from each participant to map brain function on a finer scale. In this review, we provide a broad overview of how person-specific and precision imaging techniques have used resting-state measures to examine individual differences in the brain's organization and their impact on behavior, followed by how task-based activity continues to add detail to these discoveries. We argue that person-specific and precision approaches demonstrate substantial promise in uncovering new details of the brain's functional organization and its relationship to behavior in many areas of cognitive neuroscience. We also discuss some current limitations in this new field and some new directions it may take.

Relationship of neurite architecture to brain activity during task-based fMRI

Functional MRI (fMRI) has been widely used to examine changes in neuronal activity during cognitive tasks. Commonly used measures of gray matter macrostructure (e.g., cortical thickness, surface area, volume) do not consistently appear to serve as structural correlates of brain function. In contrast, gray matter microstructure, measured using neurite orientation dispersion and density imaging (NODDI), enables the estimation of indices of neurite density (neurite density index; NDI) and organization (orientation dispersion index; ODI) in gray matter. Our study explored the relationship among neurite architecture, BOLD (blood-oxygen-level-dependent) fMRI, and cognition, using a large sample (n = 750) of young adults of the human connectome project (HCP) and two tasks that index more cortical (working memory) and more subcortical (emotion processing) targeting of brain functions. Using NODDI, fMRI, structural MRI and task performance data, hierarchical regression analyses revealed that higher working memory- and emotion processing-evoked BOLD activity was related to lower ODI in the right DLPFC, and lower ODI and NDI values in the right and left amygdala, respectively. Common measures of brain macrostructure (i.e., DLPFC thickness/surface area and amygdala volume) did not explain any additional variance (beyond neurite architecture) in BOLD activity. A moderating effect of neurite architecture on the relationship between emotion processing task-evoked BOLD response and performance was observed. Our findings provide evidence that neuro-/social-affective cognition-related BOLD activity is partially driven by the local neurite organization and density with direct impact on emotion processing. In vivo gray matter microstructure represents a new target of investigation providing strong potential for clinical translation.

A causal role of anterior prefrontal-putamen circuit for response inhibition revealed by transcranial ultrasound stimulation in humans

Stopping an inappropriate response requires the involvement of the prefrontal-subthalamic hyperdirect pathway. However, how the prefrontal-striatal indirect pathway contributes to stopping is poorly understood. In this study, transcranial ultrasound stimulation is used to perform interventions in a task-related region in the striatum. Functional magnetic resonance imaging (MRI) reveals activation in the right anterior part of the putamen during response inhibition, and ultrasound stimulation to the anterior putamen, as well as the subthalamic nucleus, results in significant impairments in stopping performance. Diffusion imaging further reveals prominent structural connections between the anterior putamen and the right anterior part of the inferior frontal cortex (IFC), and ultrasound stimulation to the anterior IFC also shows significant impaired stopping performance. These results demonstrate that the right anterior putamen and right anterior IFC causally contribute to stopping and suggest that the anterior IFC-anterior putamen circuit in the indirect pathway serves as an essential route for stopping.

Heritability and cross-species comparisons of human cortical functional organization asymmetry

The human cerebral cortex is symmetrically organized along large-scale axes but also presents inter-hemispheric differences in structure and function. The quantified contralateral homologous difference, that is asymmetry, is a key feature of the human brain left-right axis supporting functional processes, such as language. Here, we assessed whether the asymmetry of cortical functional organization is heritable and phylogenetically conserved between humans and macaques. Our findings indicate asymmetric organization along an axis describing a functional trajectory from perceptual/action to abstract cognition. Whereas language network showed leftward asymmetric organization, frontoparietal network showed rightward asymmetric organization in humans. These asymmetries were heritable in humans and showed a similar spatial distribution with macaques, in the case of intra-hemispheric asymmetry of functional hierarchy. This suggests (phylo)genetic conservation. However, both language and frontoparietal networks showed a qualitatively larger asymmetry in humans relative to macaques. Overall, our findings suggest a genetic basis for asymmetry in intrinsic functional organization, linked to higher order cognitive functions uniquely developed in humans.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

The resting-state causal human connectome is characterized by hub connectivity of executive and attentional networks

We demonstrate a data-driven approach for calculating a "causal connectome" of directed connectivity from resting-state fMRI data using a greedy adjacency search and pairwise non-Gaussian edge orientations. We used this approach to construct n = 442 causal connectomes. These connectomes were very sparse in comparison to typical Pearson correlation-based graphs (roughly 2.25% edge density) yet were fully connected in nearly all cases. Prominent highly connected hubs of the causal connectome were situated in attentional (dorsal attention) and executive (frontoparietal and cingulo-opercular) networks. These hub networks had distinctly different connectivity profiles: attentional networks shared incoming connections with sensory regions and outgoing connections with higher cognitive networks, while executive networks primarily connected to other higher cognitive networks and had a high degree of bidirected connectivity. Virtual lesion analyses accentuated these findings, demonstrating that attentional and executive hub networks are points of critical vulnerability in the human causal connectome. These data highlight the central role of attention and executive control networks in the human cortical connectome and set the stage for future applications of data-driven causal connectivity analysis in psychiatry.

Individualized Functional Subnetworks Connect Human Striatum and Frontal Cortex

The striatum and cerebral cortex are interconnected via multiple recurrent loops that play a major role in many neuropsychiatric conditions. Primate corticostriatal connections can be precisely mapped using invasive tract-tracing. However, noninvasive human research has not mapped these connections with anatomical precision, limited in part by the practice of averaging neuroimaging data across individuals. Here we utilized highly sampled resting-state functional connectivity MRI for individual-specific precision functional mapping (PFM) of corticostriatal connections. We identified ten individual-specific subnetworks linking cortex-predominately frontal cortex-to striatum, most of which converged with nonhuman primate tract-tracing work. These included separable connections between nucleus accumbens core/shell and orbitofrontal/medial frontal gyrus; between anterior striatum and dorsomedial prefrontal cortex; between dorsal caudate and lateral prefrontal cortex; and between middle/posterior putamen and supplementary motor/primary motor cortex. Two subnetworks that did not converge with nonhuman primates were connected to cortical regions associated with human language function. Thus, precision subnetworks identify detailed, individual-specific, neurobiologically plausible corticostriatal connectivity that includes human-specific language networks.

Dynamical differential covariance recovers directional network structure in multiscale neural systems

Investigating neural interactions is essential to understanding the neural basis of behavior. Many statistical methods have been used for analyzing neural activity, but estimating the direction of network interactions correctly and efficiently remains a difficult problem. Here, we derive dynamical differential covariance (DDC), a method based on dynamical network models that detects directional interactions with low bias and high noise tolerance under nonstationarity conditions. Moreover, DDC scales well with the number of recording sites and the computation required is comparable to that needed for covariance. DDC was validated and compared favorably with other methods on networks with false positive motifs and multiscale neural simulations where the ground-truth connectivity was known. When applied to recordings of resting-state functional magnetic resonance imaging (rs-fMRI), DDC consistently detected regional interactions with strong structural connectivity in over 1,000 individual subjects obtained by diffusion MRI (dMRI). DDC is a promising family of methods for estimating connectivity that can be generalized to a wide range of dynamical models and recording techniques and to other applications where system identification is needed.

Resting Hypoconnectivity of Theoretically Defined Addiction Networks during Early Abstinence Predicts Subsequent Relapse in Alcohol Use Disorder

Theoretical models of addiction suggest that alterations in addiction domains including incentive salience, negative emotionality, and executive control lead to relapse in alcohol use disorder (AUD). To determine whether the functional organization of neural networks underlying these domains predict subsequent relapse, we generated theoretically defined addiction networks. We collected resting functional magnetic resonance imaging data from 45 individuals with AUD during early abstinence (number of days abstinent M = 25.40, SD = 16.51) and calculated the degree of resting-state functional connectivity (RSFC) within these networks. Regression analyses determined whether the RSFC strength in domain-defined addiction networks measured during early abstinence predicted subsequent relapse (dichotomous or continuous relapse metrics). RSFC within each addiction network measured during early abstinence was significantly lower in those that relapsed (vs. abstained) and predicted subsequent time to relapse. Lower incentive salience RSFC during early abstinence increased the odds of relapsing. Neither RSFC in a control network nor clinical self-report measures predicted relapse. The association between low incentive salience RSFC and faster relapse highlights the need to design timely interventions that enhance RSFC in AUD individuals at risk of relapsing faster.

Act natural: Functional connectivity from naturalistic stimuli fMRI outperforms resting-state in predicting brain activity

The search for an 'ideal' approach to investigate the functional connections in the human brain is an ongoing challenge for the neuroscience community. While resting-state functional magnetic resonance imaging (fMRI) has been widely used to study individual functional connectivity patterns, recent work has highlighted the benefits of collecting functional connectivity data while participants are exposed to naturalistic stimuli, such as watching a movie or listening to a story. For example, functional connectivity data collected during movie-watching were shown to predict cognitive and emotional scores more accurately than resting-state-derived functional connectivity. We have previously reported a tight link between resting-state functional connectivity and task-derived neural activity, such that the former successfully predicts the latter. In the current work we use data from the Human Connectome Project to demonstrate that naturalistic-stimulus-derived functional connectivity predicts task-induced brain activation maps more accurately than resting-state-derived functional connectivity. We then show that activation maps predicted using naturalistic stimuli are better predictors of individual intelligence scores than activation maps predicted using resting-state. We additionally examine the influence of naturalistic-stimulus type on prediction accuracy. Our findings emphasize the potential of naturalistic stimuli as a promising alternative to resting-state fMRI for connectome-based predictive modelling of individual brain activity and cognitive traits.

Modeling spatial, developmental, physiological, and topological constraints on human brain connectivity

The complex connectivity of nervous systems is thought to have been shaped by competitive selection pressures to minimize wiring costs and support adaptive function. Accordingly, recent modeling work indicates that stochastic processes, shaped by putative trade-offs between the cost and value of each connection, can successfully reproduce many topological properties of macroscale human connectomes measured with diffusion magnetic resonance imaging. Here, we derive a new formalism that more accurately captures the competing pressures of wiring cost minimization and topological complexity. We further show that model performance can be improved by accounting for developmental changes in brain geometry and associated wiring costs, and by using interregional transcriptional or microstructural similarity rather than topological wiring rules. However, all models struggled to capture topographical (i.e., spatial) network properties. Our findings highlight an important role for genetics in shaping macroscale brain connectivity and indicate that stochastic models offer an incomplete account of connectome organization.

A Riemannian approach to predicting brain function from the structural connectome

Ongoing brain function is largely determined by the underlying wiring of the brain, but the specific rules governing this relationship remain unknown. Emerging literature has suggested that functional interactions between brain regions emerge from the structural connections through mono- as well as polysynaptic mechanisms. Here, we propose a novel approach based on diffusion maps and Riemannian optimization to emulate this dynamic mechanism in the form of random walks on the structural connectome and predict functional interactions as a weighted combination of these random walks. Our proposed approach was evaluated in two different cohorts of healthy adults (Human Connectome Project, HCP; Microstructure-Informed Connectomics, MICs). Our approach outperformed existing approaches and showed that performance plateaus approximately around the third random walk. At macroscale, we found that the largest number of walks was required in nodes of the default mode and frontoparietal networks, underscoring an increasing relevance of polysynaptic communication mechanisms in transmodal cortical networks compared to primary and unimodal systems.

The retrocalcarine sulcus maps different retinotopic representations in macaques and humans

Primate cerebral cortex is highly convoluted with much of the cortical surface buried in sulcal folds. The origins of cortical folding and its functional relevance have been a major focus of systems and cognitive neuroscience, especially when considering stereotyped patterns of cortical folding that are shared across individuals within a primate species and across multiple species. However, foundational questions regarding organizing principles shared across species remain unanswered. Taking a cross-species comparative approach with a careful consideration of historical observations, we investigate cortical folding relative to primary visual cortex (area V1). We identify two macroanatomical structures-the retrocalcarine and external calcarine sulci-in 24 humans and 6 macaque monkeys. We show that within species, these sulci are identifiable in all individuals, fall on a similar part of the V1 retinotopic map, and thus, serve as anatomical landmarks predictive of functional organization. Yet, across species, the underlying eccentricity representations corresponding to these macroanatomical structures differ strikingly across humans and macaques. Thus, the correspondence between retinotopic representation and cortical folding for an evolutionarily old structure like V1 is species-specific and suggests potential differences in developmental and experiential constraints across primates.

Dynamic trajectories of connectome state transitions are heritable

The brain’s functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly, Fractional Occupancy (FO) and Transition Probability (TP), representing the duration spent in each connectivity configuration and the frequency of shifting between configurations, respectively. These effects were robust against methodological choices of number of states and global signal regression. Genetic effects explained a substantial proportion of phenotypic variance of these features (h² = 0.39, 95% CI = [.24,.54] for FO; h² = 0. 43, 95% CI = [.29,.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states’ Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the dynamic trajectory of state transitions (captured by FO and TP), and such temporal features may act as endophenotypes for cognitive abilities.

Local structure-function relationships in human brain networks across the lifespan

A growing number of studies have used stylized network models of communication to predict brain function from structure. Most have focused on a small set of models applied globally. Here, we compare a large number of models at both global and regional levels. We find that globally most predictors perform poorly. At the regional level, performance improves but heterogeneously, both in terms of variance explained and the optimal model. Next, we expose synergies among predictors by using pairs to jointly predict FC. Finally, we assess age-related differences in global and regional coupling across the human lifespan. We find global decreases in the magnitude of structure-function coupling with age. We find that these decreases are driven by reduced coupling in sensorimotor regions, while higher-order cognitive systems preserve local coupling with age. Our results describe patterns of structure-function coupling across the cortex and how this may change with age.

Cortico-subcortical interactions in overlapping communities of edge functional connectivity

Both cortical and subcortical regions can be functionally organized into networks. Regions of the basal ganglia are extensively interconnected with the cortex via reciprocal connections that relay and modulate cortical function. Here we employ an edge-centric approach, which computes co-fluctuations among region pairs in a network to investigate the role and interaction of subcortical regions with cortical systems. By clustering edges into communities, we show that cortical systems and subcortical regions couple via multiple edge communities, with hippocampus and amygdala having a distinct pattern from striatum and thalamus. We show that the edge community structure of cortical networks is highly similar to one obtained from cortical nodes when the subcortex is present in the network. Additionally, we show that the edge community profile of both cortical and subcortical nodes can be estimates solely from cortico-subcortical interactions. Finally, we used a motif analysis focusing on edge community triads where a subcortical region coupled to two cortical regions and found that two community triads where one community couples the subcortex to the cortex were overrepresented. In summary, our results show organized coupling of the subcortex to the cortex that may play a role in cortical organization of primary sensorimotor/attention and heteromodal systems and puts forth the motif analysis of edge community triads as a promising method for investigation of communication patterns in networks.

Cortical thickness systematically varies with curvature and depth in healthy human brains

Cortical thickness varies throughout the cortex in a systematic way. However, it is challenging to investigate the patterns of cortical thickness due to the intricate geometry of the cortex. The cortex has a folded nature both in radial and tangential directions which forms not only gyri and sulci but also tangential folds and intersections. In this article, cortical curvature and depth are used to characterize the spatial distribution of the cortical thickness with much higher resolution than conventional regional atlases. To do this, a computational pipeline was developed that is capable of calculating a variety of quantitative measures such as surface area, cortical thickness, curvature (mean curvature, Gaussian curvature, shape index, intrinsic curvature index, and folding index), and sulcal depth. By analyzing 501 neurotypical adult human subjects from the ABIDE-I dataset, we show that cortex has a very organized structure and cortical thickness is strongly correlated with local shape. Our results indicate that cortical thickness consistently increases along the gyral-sulcal spectrum from concave to convex shape, encompassing the saddle shape along the way. Additionally, tangential folds influence cortical thickness in a similar way as gyral and sulcal folds; outer folds are consistently thicker than inner.

Disentangling the variability of the superficial white matter organization using regional-tractogram-based population stratification

Each variation of the cortical folding pattern implies a particular rearrangement of the geometry of the fibers of the underlying white matter. While this rearrangement only impacts the ends of the long pathways, it may affect most of the trajectory of the short bundles. Therefore, mapping the short fibers of the human brain using diffusion-based tractography requires a dedicated strategy to overcome the variability of the folding patterns. In this paper, we propose a fiber-based stratification strategy splitting the population into homogeneous groups for disentangling the superficial white matter bundle organization. This strategy introduces a new refined fiber distance which includes angular considerations for inferring fine-grained atlases of the short bundles surrounding a specific sulcus and a subtractogram distance that quantifies the similitude between fiber sets of two different subjects. The stratification splits the population into groups with similar regional fiber organization using manifold learning. We first successfully test the hypothesis that the main source of variability of the regional fiber organization is the variability of the regional folding pattern. Then, in each group, we proceed with the automatic identification of the most stable bundles, at a higher granularity level than what can be achieved with the non-stratified whole population, enabling the disentanglement of the very variable configuration of the short fibers. Finally, the method searches for bundle correspondence across groups to build a population level atlas. As a proof of concept, the atlas refinement achieved by this strategy is illustrated for the fibers that surround the central sulcus and the superior temporal sulcus using the HCP dataset.

The medial occipital longitudinal tract supports early stage encoding of visuospatial information

Visuospatial learning depends on the parahippocampal place area (PPA), a functionally heterogenous area which current visuospatial processing models place downstream from parietal cortex and only from area V4 of early visual cortex (EVC). However, evidence for anatomical connections between the PPA and other EVC areas is inconsistent, and these connections are not discussed in current models. Through a data-driven analysis based on diffusion MRI tractography, we present evidence that the PPA sits at the confluence of two white matter systems. The first conveys information from the retrosplenial complex to the anterior PPA and runs within the cingulum bundle. The second system connects all peripheral EVC areas to the posterior PPA and corresponds to the medial occipital longitudinal tract (MOLT), a white matter pathway that is distinct from the cingulum and that we describe here in detail. Based on further functional connectivity analysis and meta-analytic data, we propose that the MOLT supports early stage encoding of visuospatial information by allowing direct reciprocal exchange between the PPA and EVC. Our findings may improve symptom interpretation in stroke and tumour patients with damage to the medial occipito-temporal region and call for revisiting current visuospatial processing models.

A white matter pathway (termed, MOLT) connecting the parahippocampal place area and the medial early visual cortex contributes to visuospatial learning in humans.

Control theory illustrates the energy efficiency in the dynamic reconfiguration of functional connectivity

The brain's functional connectivity fluctuates over time instead of remaining steady in a stationary mode even during the resting state. This fluctuation establishes the dynamical functional connectivity that transitions in a non-random order between multiple modes. Yet it remains unexplored how the transition facilitates the entire brain network as a dynamical system and what utility this mechanism for dynamic reconfiguration can bring over the widely used graph theoretical measurements. To address these questions, we propose to conduct an energetic analysis of functional brain networks using resting-state fMRI and behavioral measurements from the Human Connectome Project. Through comparing the state transition energy under distinct adjacent matrices, we justify that dynamic functional connectivity leads to 60% less energy cost to support the resting state dynamics than static connectivity when driving the transition through default mode network. Moreover, we demonstrate that combining graph theoretical measurements and our energy-based control measurements as the feature vector can provide complementary prediction power for the behavioral scores (Combination vs. Control: t = 9.41, p = 1.64e-13; Combination vs. Graph: t = 4.92, p = 3.81e-6). Our approach integrates statistical inference and dynamical system inspection towards understanding brain networks.

Neonatal multi-modal cortical profiles predict 18-month developmental outcomes

Developmental delays in infanthood often persist, turning into life-long difficulties, and coming at great cost for the individual and community. By examining the developing brain and its relation to developmental outcomes we can start to elucidate how the emergence of brain circuits is manifested in variability of infant motor, cognitive and behavioural capacities. In this study, we examined if cortical structural covariance at birth, indexing coordinated development, is related to later infant behaviour. We included 193 healthy term-born infants from the Developing Human Connectome Project (dHCP). An individual cortical connectivity matrix derived from morphological and microstructural features was computed for each subject (morphometric similarity networks, MSNs) and was used as input for the prediction of behavioural scores at 18 months using Connectome-Based Predictive Modeling (CPM). Neonatal MSNs successfully predicted social-emotional performance. Predictive edges were distributed between and within known functional cortical divisions with a specific important role for primary and posterior cortical regions. These results reveal that multi-modal neonatal cortical profiles showing coordinated maturation are related to developmental outcomes and that network organization at birth provides an early infrastructure for future functional skills.

Effects of virtual lesions on temporal dynamics in cortical networks based on personalized dynamic models

The human brain dynamically shifts between a predominantly integrated state and a predominantly segregated state, each with different roles in supporting cognition and behavior. However, no studies to date have investigated lesions placed in different regions of the cerebral cortex to determine the effects on the temporal balance of integration and segregation. Here, we used the integrated state occurrence rate to measure the temporal balance of integration and segregation in the resting brain. Based on dynamic mean-field models coupled through the individual's structural white matter connections, neural activity was modeled. By lesioning different individual nodes from the model, our results implied that the impact of lesions reaches far beyond focal damage and could impair cognition by affecting system-level dynamics. Lesions in a portion of the nodes in the visual, somatomotor, limbic, and default networks significantly compromised the temporal balance of integration and segregation. Crucially, the effects of lesions in different regions could be predicted based on the hierarchical axis inferred from the T1w/T2w map and specific graph measures of structural brain networks. Taken together, our findings indicate the possibility of significant contributions of anatomical heterogeneity to the dynamics of functional network topology. Although still in its infancy, computational modeling may provide an entry point for understanding brain disorders at a causal mechanistic level, possibly leading to novel, more effective therapeutic interventions.

Neural spatio-temporal patterns of information processing related to cognitive conflict and correct or false recognitions

Using a visual short-term memory task and employing a new methodological approach, we analyzed neural responses from the perspective of the conflict level and correctness/erroneous over a longer time window. Sixty-five participants performed the short-term memory task in the fMRI scanner. We explore neural spatio-temporal patterns of information processing in the context of correct or erroneous response and high or low level of cognitive conflict using classical fMRI analysis, surface-based cortical data, temporal analysis of interpolated mean activations, and machine learning classifiers. Our results provide evidence that information processing dynamics during the retrieval process vary depending on the correct or false recognition—for stimuli inducing a high level of cognitive conflict and erroneous response, information processing is prolonged. The observed phenomenon may be interpreted as the manifestation of the brain’s preparation for future goal-directed action.

Diffeomorphic registration for retinotopic maps of multiple visual regions

Retinotopic map, the mapping between visual inputs on the retina and neuronal responses on the cortical surface, is one of the central topics in vision science. Typically, human retinotopic maps are constructed by analyzing functional magnetic resonance responses to designed visual stimuli on the cortical surface. Although it is widely used in visual neuroscience, retinotopic maps are limited by the signal-to-noise ratio and spatial resolution of fMRI. One promising approach to improve the quality of retinotopic maps is to register individual subject's retinotopic maps to a retinotopic template. However, none of the existing retinotopic registration methods has explicitly quantified the diffeomorphic condition, that is, retinotopic maps shall be aligned by stretching/compressing without tearing up the cortical surface. Here, we developed Diffeomorphic Registration for Retinotopic Maps (DRRM) to simultaneously align retinotopic maps in multiple visual regions under the diffeomorphic condition. Specifically, we used the Beltrami coefficient to model the diffeomorphic condition and performed surface registration based on retinotopic coordinates. The overall framework preserves the topological condition defined in the template. We further developed a unique evaluation protocol and compared the performance of the new method with several existing registration methods on both synthetic and real datasets. The results showed that DRRM is superior to the existing methods in achieving diffeomorphic  registration in synthetic and empirical data from 3T and 7T MRI systems. DRRM may improve the interpretation of low-quality retinotopic maps and facilitate applications of retinotopic maps in clinical settings.

Fiber tracing and microstructural characterization among audiovisual integration brain regions in neonates compared with young adults

Audiovisual integration has been related with cognitive-processing and behavioral advantages, as well as with various socio-cognitive disorders. While some studies have identified brain regions instantiating this ability shortly after birth, little is known about the structural pathways connecting them. The goal of the present study was to reconstruct fiber tracts linking AVI regions in the newborn in-vivo brain and assess their adult-likeness by comparing them with analogous fiber tracts of young adults. We performed probabilistic tractography and compared connective probabilities between a sample of term-born neonates (N = 311; the Developing Human Connectome Project (dHCP, http://www.developingconnectome.org) and young adults (N = 311 The Human Connectome Project; https://www.humanconnectome.org/) by means of a classification algorithm. Furthermore, we computed Dice coefficients to assess between-group spatial similarity of the reconstructed fibers and used diffusion metrics to characterize neonates' AVI brain network in terms of microstructural properties, interhemispheric differences and the association with perinatal covariates and biological sex. Overall, our results indicate that the AVI fiber bundles were successfully reconstructed in a vast majority of neonates, similarly to adults. Connective probability distributional similarities and spatial overlaps of AVI fibers between the two groups differed across the reconstructed fibers. There was a rank-order correspondence of the fibers' connective strengths across the groups. Additionally, the study revealed patterns of diffusion metrics in line with early white matter developmental trajectories and a developmental advantage for females. Altogether, these findings deliver evidence of meaningful structural connections among AVI regions in the newborn in-vivo brain.

Browsing Multiple Subjects When the Atlas Adaptation Cannot Be Achieved via a Warping Strategy

Brain mapping studies often need to identify brain structures or functional circuits into a set of individual brains. To this end, multiple atlases have been published to represent such structures based on different modalities, subject sets, and techniques. The mainstream approach to exploit these atlases consists in spatially deforming each individual data onto a given atlas using dense deformation fields, which supposes the existence of a continuous mapping between atlases and individuals. However, this continuity is not always verified, and this "iconic" approach has limits. We present in this study an alternative, complementary, "structural" approach, which consists in extracting structures from the individual data, and comparing them without deformation. A "structural atlas" is thus a collection of annotated individual data with a common structure nomenclature. It may be used to characterize structure shape variability across individuals or species, or to train machine learning systems. This study exhibits Anatomist, a powerful structural 3D visualization software dedicated to building, exploring, and editing structural atlases involving a large number of subjects. It has been developed primarily to decipher the cortical folding variability; cortical sulci vary enormously in both size and shape, and some may be missing or have various topologies, which makes iconic approaches inefficient to study them. We, therefore, had to build structural atlases for cortical sulci, and use them to train sulci identification algorithms. Anatomist can display multiple subject data in multiple views, supports all kinds of neuroimaging data, including compound structural object graphs, handles arbitrary coordinate transformation chains between data, and has multiple display features. It is designed as a programming library in both C++ and Python languages, and may be extended or used to build dedicated custom applications. Its generic design makes all the display and structural aspects used to explore the variability of the cortical folding pattern work in other applications, for instance, to browse axonal fiber bundles, deep nuclei, functional activations, or other kinds of cortical parcellations. Multimodal, multi-individual, or inter-species display is supported, and adaptations to large scale screen walls have been developed. These very original features make it a unique viewer for structural atlas browsing.

The Social Cerebellum: A Large-Scale Investigation of Functional and Structural Specificity and Connectivity

The cerebellum has been traditionally disregarded in relation to nonmotor functions, but recent findings indicate it may be involved in language, affective processing, and social functions. Mentalizing, or Theory of Mind (ToM), is the ability to infer mental states of others and this skill relies on a distributed network of brain regions. Here, we leveraged large-scale multimodal neuroimaging data to elucidate the structural and functional role of the cerebellum in mentalizing. We used functional activations to determine whether the cerebellum has a domain-general or domain-specific functional role, and effective connectivity and probabilistic tractography to map the cerebello-cerebral mentalizing network. We found that the cerebellum is organized in a domain-specific way and that there is a left cerebellar effective and structural lateralization, with more and stronger effective connections from the left cerebellar hemisphere to the right cerebral mentalizing areas, and greater cerebello-thalamo-cortical and cortico-ponto-cerebellar streamline counts from and to the left cerebellum. Our study provides novel insights to the network organization of the cerebellum, an overlooked brain structure, and mentalizing, one of humans' most essential abilities to navigate the social world.