Fine-grained topographic organization within somatosensory cortex during resting-state and emotional face-matching task and its association with ASD traits.

Sensory atypicalities are particularly common in autism spectrum disorders (ASD). Nevertheless, our knowledge about the divergent functioning of the underlying somatosensory region and its association with ASD phenotype features is limited. We applied a data-driven approach to map the fine-grained variations in functional connectivity of the primary somatosensory cortex (S1) to the rest of the brain in 240 autistic and 164 neurotypical individuals from the EU-AIMS LEAP dataset, aged between 7 and 30. We estimated the S1 connection topography ('connectopy') at rest and during the emotional face-matching (Hariri) task, an established measure of emotion reactivity, and accessed its association with a set of clinical and behavioral variables. We first demonstrated that the S1 connectopy is organized along a dorsoventral axis, mapping onto the S1 somatotopic organization. We then found that its spatial characteristics were linked to the individuals' adaptive functioning skills, as measured by the Vineland Adaptive Behavior Scales, across the whole sample. Higher functional differentiation characterized the S1 connectopies of individuals with higher daily life adaptive skills. Notably, we detected significant differences between rest and the Hariri task in the S1 connectopies, as well as their projection maps onto the rest of the brain suggesting a task-modulating effect on S1 due to emotion processing. All in all, variation of adaptive skills appears to be reflected in the brain's mesoscale neural circuitry, as shown by the S1 connectivity profile, which is also differentially modulated during rest and emotional processing.

Striatal connectopic maps link to functional domains across psychiatric disorders.

Transdiagnostic approaches to psychiatry have significant potential in overcoming the limitations of conventional diagnostic paradigms. However, while frameworks such as the Research Domain Criteria have garnered significant enthusiasm among researchers and clinicians from a theoretical angle, examples of how such an approach might translate in practice to understand the biological mechanisms underlying complex patterns of behaviors in realistic and heterogeneous populations have been sparse. In a richly phenotyped clinical sample (n = 186) specifically designed to capture the complex nature of heterogeneity and comorbidity within- and between stress- and neurodevelopmental disorders, we use exploratory factor analysis on a wide range of clinical questionnaires to identify four stable functional domains that transcend diagnosis and relate to negative valence, cognition, social functioning and inhibition/arousal before replicating them in an independent dataset (n = 188). We then use connectopic mapping to map inter-individual variation in fine-grained topographical organization of functional connectivity in the striatum-a central hub in motor, cognitive, affective and reward-related brain circuits-and use multivariate machine learning (canonical correlation analysis) to show that these individualized topographic representations predict transdiagnostic functional domains out of sample (r = 0.20, p = 0.026). We propose that investigating psychiatric symptoms across disorders is a promising path to linking them to underlying biology, and can help bridge the gap between neuroscience and clinical psychiatry.

Bayesian connective field modeling using a Markov Chain Monte Carlo approach.

The majority of neurons in the human brain process signals from neurons elsewhere in the brain. Connective Field (CF) modelling is a biologically-grounded method to describe this essential aspect of the brain's circuitry. It allows characterizing the response of a population of neurons in terms of the activity in another part of the brain. CF modelling translates the concept of the receptive field (RF) into the domain of connectivity by assessing, at the voxel level, the spatial dependency between signals in distinct cortical visual field areas. Thus, the approach enables to characterize the functional cortical circuitry of the human cortex. While already very useful, the present CF modelling approach has some intrinsic limitations due to the fact that it only estimates the model's explained variance and not the probability distribution associated with the estimated parameters. If we could resolve this, CF modelling would lend itself much better for statistical comparisons at the level of single voxels and individuals. This is important when trying to gain a detailed understanding of the neurobiology and pathophysiology of the visual cortex, notably in rare cases. To enable this, we present a Bayesian approach to CF modeling (bCF). Using a Markov Chain Monte Carlo (MCMC) procedure, it estimates the posterior probability distribution underlying the CF parameters. Based on this, bCF quantifies, at the voxel level, the uncertainty associated with each parameter estimate. This information can be used in various ways to increase confidence in the CF model predictions. We applied bCF to BOLD responses recorded in the early human visual cortex using 3T fMRI. We estimated both the CF parameters and their associated uncertainties and show they are only weakly correlated. Moreover, we show how bCF facilitates the use of effect size (beta) as a data-driven parameter that can be used to select the most reliable voxels for further analysis. Finally, to further illustrate the functionality gained by bCF, we apply it to perform a voxel-level comparison of a single, circular symmetric, Gaussian versus a Difference-of-Gaussian model. We conclude that our bCF framework provides a comprehensive tool to study human functional cortical circuitry in health and disease.

Mapping dopaminergic projections in the human brain with resting-state fMRI.

The striatum receives dense dopaminergic projections, making it a key region of the dopaminergic system. Its dysfunction has been implicated in various conditions including Parkinson's disease (PD) and substance use disorder. However, the investigation of dopamine-specific functioning in humans is problematic as current MRI approaches are unable to differentiate between dopaminergic and other projections. Here, we demonstrate that 'connectopic mapping' - a novel approach for characterizing fine-grained, overlapping modes of functional connectivity - can be used to map dopaminergic projections in striatum. We applied connectopic mapping to resting-state functional MRI data of the Human Connectome Project (population cohort; N = 839) and selected the second-order striatal connectivity mode for further analyses. We first validated its specificity to dopaminergic projections by demonstrating a high spatial correlation (r = 0.884) with dopamine transporter availability - a marker of dopaminergic projections - derived from DaT SPECT scans of 209 healthy controls. Next, we obtained the subject-specific second-order modes from 20 controls and 39 PD patients scanned under placebo and under dopamine replacement therapy (L-DOPA), and show that our proposed dopaminergic marker tracks PD diagnosis, symptom severity, and sensitivity to L-DOPA. Finally, across 30 daily alcohol users and 38 daily smokers, we establish strong associations with self-reported alcohol and nicotine use. Our findings provide evidence that the second-order mode of functional connectivity in striatum maps onto dopaminergic projections, tracks inter-individual differences in PD symptom severity and L-DOPA sensitivity, and exhibits strong associations with levels of nicotine and alcohol use, thereby offering a new biomarker for dopamine-related (dys)function in the human brain.

Connectivity gradients on tractography data: Pipeline and example applications.

Gray matter connectivity can be described in terms of its topographical organization, but the differential role of white matter connections underlying that organization is often unknown. In this study, we propose a method for unveiling principles of organization of both gray and white matter based on white matter connectivity as assessed using diffusion magnetic ressonance imaging (MRI) tractography with spectral embedding gradient mapping. A key feature of the proposed approach is its capacity to project the individual connectivity gradients it reveals back onto its input data in the form of projection images, allowing one to assess the contributions of specific white matter tracts to the observed gradients. We demonstrate the ability of our proposed pipeline to identify connectivity gradients in prefrontal and occipital gray matter. Finally, leveraging the use of tractography, we demonstrate that it is possible to observe gradients within the white matter bundles themselves. Together, the proposed framework presents a generalized way to assess both the topographical organization of structural brain connectivity and the anatomical features driving it.

Atypical Integration of Sensory-to-Transmodal Functional Systems Mediates Symptom Severity in Autism.

A notable characteristic of autism spectrum disorder (ASD) is co-occurring deficits in low-level sensory processing and high-order social interaction. While there is evidence indicating detrimental cascading effects of sensory anomalies on the high-order cognitive functions in ASD, the exact pathological mechanism underlying their atypical functional interaction across the cortical hierarchy has not been systematically investigated. To address this gap, here we assessed the functional organisation of sensory and motor areas in ASD, and their relationship with subcortical and high-order trandmodal systems. In a resting-state fMRI data of 107 ASD and 113 neurotypical individuals, we applied advanced connectopic mapping to probe functional organization of primary sensory/motor areas, together with targeted seed-based intrinsic functional connectivity (iFC) analyses. In ASD, the connectopic mapping revealed topological anomalies (i.e., excessively more segregated iFC) in the motor and visual areas, the former of which patterns showed association with the symptom severity of restricted and repetitive behaviors. Moreover, the seed-based analysis found diverging patterns of ASD-related connectopathies: decreased iFCs within the sensory/motor areas but increased iFCs between sensory and subcortical structures. While decreased iFCs were also found within the higher-order functional systems, the overall proportion of this anomaly tends to increase along the level of cortical hierarchy, suggesting more dysconnectivity in the higher-order functional networks. Finally, we demonstrated that the association between low-level sensory/motor iFCs and clinical symptoms in ASD was mediated by the high-order transmodal systems, suggesting pathogenic functional interactions along the cortical hierarchy. Findings were largely replicated in the independent dataset. These results highlight that atypical integration of sensory-to-high-order systems contributes to the complex ASD symptomatology.

Functional co-activation of the default mode network in APOE ε4-carriers: A replication study.

Structural and functional alterations of the brain in persons genetically at-risk for Alzheimer's disease (AD) are crucial in unravelling AD development. Filippini et al. found that the default mode network (DMN) is already affected in young APOE ε4-carriers, with increased co-activation of the DMN during rest and increased hippocampal task activation. We aimed to replicate the early findings of Filippini et al, using the APOE gene, still the principal AD risk gene, and extended this with a polygenic risk score (PRS) analysis for AD, using the Human Connectome Project dataset (HCP). We included participants from the HCP S1200 dataset (age range: 22-36 years). We studied morphometric features, functional DMN co-activation and functional task activation of recollection performance. Permutation Analysis of Linear Models (PALM) was used to test for group differences between APOE ε4-carriers and non-carriers, and to test the association with PRS. PALM controls for biases induced by the family structure of the HCP sample. Results were family-wise error rate corrected at p < 0.05. Our primary analysis did not replicate the early findings of Filippini et al. (2009). However, compared with non-carriers, APOE ε4-carriers showed increased functional activation during the encoding of subsequently recollected items in areas related to facial recognition (p<0.05, t>756.11). This increased functional activation was also positively associated with PRS (APOE variants included) (p<0.05, t>647.55). Our results are supportive for none to limited genetic effects on brain structure and function in young adults. Taking the methodological considerations of replication studies into account, the true effect of APOE ε4-carriership is likely smaller than indicated in the Filippini paper. However, it still holds that we may not yet be able to detect already present measurable effects decades before a clinical expression of AD. Since the mechanistic pathway of AD is likely to encompass many different factors, further research should be focused on the interactions of genetic risk, biomarkers, aging and lifestyle factors over the life course. Sensitive functional neuroimaging as used here may help disentangling these complex interactions.

Mesoscale hierarchical organization of primary somatosensory cortex captured by resting-state-fMRI in humans.

The primary somatosensory cortex (S1) plays a key role in the processing and integration of afferent somatosensory inputs along an anterior-to-posterior axis, contributing towards necessary human function. It is believed that anatomical connectivity can be used to probe hierarchical organization, however direct characterization of this principle in-vivo within humans remains elusive. Here, we use resting-state functional connectivity as a complement to anatomical connectivity to investigate topographical principles of human S1. We employ a novel approach to examine mesoscopic variations of functional connectivity, and demonstrate a topographic organisation spanning the region's hierarchical axis that strongly correlates with underlying microstructure while tracing along architectonic Brodmann areas. Our findings characterize anatomical hierarchy of S1 as a 'continuous spectrum' with evidence supporting a functional boundary between areas 3b and 1. The identification of this topography bridges the gap between structure and connectivity, and may be used to help further current understanding of sensorimotor deficits.

Assessing Uncertainty and Reliability of Connective Field Estimations From Resting State fMRI Activity at 3T.

Connective Field (CF) modeling estimates the local spatial integration between signals in distinct cortical visual field areas. As we have shown previously using 7T data, CF can reveal the visuotopic organization of visual cortical areas even when applied to BOLD activity recorded in the absence of external stimulation. This indicates that CF modeling can be used to evaluate cortical processing in participants in which the visual input may be compromised. Furthermore, by using Bayesian CF modeling it is possible to estimate the co-variability of the parameter estimates and therefore, apply CF modeling to single cases. However, no previous studies evaluated the (Bayesian) CF model using 3T resting-state fMRI data. This is important since 3T scanners are much more abundant and more often used in clinical research compared to 7T scanners. Therefore in this study, we investigate whether it is possible to obtain meaningful CF estimates from 3T resting state (RS) fMRI data. To do so, we applied the standard and Bayesian CF modeling approaches on two RS scans, which were separated by the acquisition of visual field mapping data in 12 healthy participants. Our results show good agreement between RS- and visual field (VF)- based maps using either the standard or Bayesian CF approach. In addition to quantify the uncertainty associated with each estimate in both RS and VF data, we applied our Bayesian CF framework to provide the underlying marginal distribution of the CF parameters. Finally, we show how an additional CF parameter, beta, can be used as a data-driven threshold on the RS data to further improve CF estimates. We conclude that Bayesian CF modeling can characterize local functional connectivity between visual cortical areas from RS data at 3T. Moreover, observations obtained using 3T scanners were qualitatively similar to those reported for 7T. In particular, we expect the ability to assess parameter uncertainty in individual participants will be important for future clinical studies.

Understanding brain organisation in the face of functional heterogeneity and functional multiplicity.

Understanding the fundamental organisation of the brain in terms of functional specialisation and integration is one of the principal aims of imaging neuroscience. Many investigations into the functional organisation of the brain are predicated on parcellating the brain into patches of assumed piece-wise constant connectivity. There are, however, many brain areas where the assumption of piece-wise constant organisation is violated. Connectivity, and by extension function, often varies continuously across the grey matter according to multiple overlapping modes of change. The organisation is governed by functional heterogeneity (continuous change) as well as functional multiplicity (overlapping modes). Functional heterogeneity and multiplicity have important implications for how we can and should analyse our data and how we ought to interpret the results, both in the classical context of parcellated modes and under models that allow for overlapping modes of continuous change. The goal of this opinion paper is to raise awareness of these issues and highlight recent methodological developments toward accounting for these important fundamental features of brain organisation.

Functional connectivity of the Precuneus reflects effectiveness of visual restitution training in chronic hemianopia.

Visual field defects in chronic hemianopia can improve through visual restitution training, yet not all patients benefit equally from this long and exhaustive procedure. Here, we asked if resting-state functional connectivity prior to visual restitution could predict training success. In two training sessions of eight weeks each, 20 patients with chronic hemianopia performed a visual discrimination task by directing spatial selective attention towards stimuli presented in either hemifield, while suppressing eye movements. We examined two effects: a sensitivity change in the attended (trained) minus the unattended (control) hemifield (i.e., a training-specific improvement), and an overall improvement (i.e., a total change in sensitivity after both sessions). We then identified five visual resting-state networks and evaluated their functional connectivity in relation to both training effects. We found that the functional connectivity strength between the anterior Precuneus and the Occipital Pole Network was positively related to the attention modulated (i.e., training-specific) improvement. No such relationship was found for the overall improvement or for the other visual networks of interest. Our finding suggests that the anterior Precuneus plays a role in attention-modulated visual field improvements. The resting-state functional connectivity between the anterior Precuneus and the Occipital Pole Network may thus serve as an imaging-based biomarker that quantifies a patient's potential capacity to direct spatial attention. This may help to identify hemianopia patients that are most likely to benefit from visual restitution training.

Principles of temporal association cortex organisation as revealed by connectivity gradients.

To establish the link between structure and function of any large area of the neocortex, it is helpful to identify its principles of organisation. One way to establish such principles is to investigate how differences in whole-brain connectivity are structured across the area. Here, we use Laplacian eigenmaps on diffusion MRI tractography data to investigate the organisational principles of the human temporal association cortex. We identify three overlapping gradients of connectivity that are, for the most part, consistent across hemispheres. The first gradient reveals an inferior-superior organisation of predominantly longitudinal tracts and separates visual and auditory unimodal and multimodal cortices. The second gradient radiates outward from the posterior middle temporal cortex with the arcuate fascicle as a distinguishing feature; the third gradient is concentrated in the anterior temporal lobe and emanates towards its posterior end. We describe the functional relevance of each of these gradients through the meta-analysis of data from the neuroimaging literature. Together, these results unravel the overlapping dimensions of structural organization of the human temporal cortex and provide a framework underlying its functional multiplicity.

Stimulus- and Neural-Referred Visual Receptive Field Properties following Hemispherectomy: A Case Study Revisited.

Damage to the visual system can result in (a partial) loss of vision, in response to which the visual system may functionally reorganize. Yet the timing, extent, and conditions under which this occurs are not well understood. Hence, studies in individuals with diverse congenital and acquired conditions and using various methods are needed to better understand this. In the present study, we examined the visual system of a young girl who received a hemispherectomy at the age of three and who consequently suffered from hemianopia. We did so by evaluating the corticocortical and retinocortical projections in the visual system of her remaining hemisphere. For the examination of these aspects, we analyzed the characteristics of the connective fields ("neural-referred" receptive fields) based on both resting-state (RS) and retinotopy data. The evaluation of RS data, reflecting brain activity independent from visual stimulation, is of particular interest as it is not biased by the patient's atypical visual percept. We found that, primarily when the patient was at rest, the connective fields between V1 and both early and late visual areas were larger than normal. These abnormally large connective fields could be a sign either of functional reorganization or of unmasked suppressive feedback signals that are normally masked by interhemispheric signals. Furthermore, we confirmed our previous finding of abnormal retinocortical or "stimulus-referred" projections in both early and late visual areas. More specifically, we found an enlarged foveal representation and smaller population receptive fields. These differences could also be a sign of functional reorganization or rather a reflection of the interruption visual information that travels, via the remainder of the visual pathway, from the retina to the visual cortex. To conclude, while we do find indications for relatively subtle changes in visual field map properties, we found no evidence of large-scale reorganization-even though the patient could have benefitted from this. Our work suggests that at a later developmental stage, large-scale reorganization of the visual system no longer occurs, while small-scale properties may still change to facilitate adaptive processing and viewing strategies.

Plasticity versus stability across the human cortical visual connectome.

Whether and how the balance between plasticity and stability varies across the brain is an important open question. Within a processing hierarchy, it is thought that plasticity is increased at higher levels of cortical processing, but direct quantitative comparisons between low- and high-level plasticity have not been made so far. Here, we address this issue for the human cortical visual system. We quantify plasticity as the complement of the heritability of resting-state functional connectivity and thereby demonstrate a non-monotonic relationship between plasticity and hierarchical level, such that plasticity decreases from early to mid-level cortex, and then increases further of the visual hierarchy. This non-monotonic relationship argues against recent theory that the balance between plasticity and stability is governed by the costs of the "coding-catastrophe", and can be explained by a concurrent decline of short-term adaptation and rise of long-term plasticity up the visual processing hierarchy.

The functional organisation of the hippocampus along its long axis is gradual and predicts recollection.

Understanding the functional organisation of the hippocampus is crucial for understanding its role in cognition and disorders in which it is implicated. Different views have been proposed of how function is distributed along its long axis: one view suggests segregation, whereas the alternative view postulates a more gradual organisation. Here, we applied a novel 'connectopic mapping' data-analysis approach to the resting-state fMRI data of participants of the Human Connectome Project, and demonstrate that the functional organisation of the hippocampal longitudinal axis is gradual rather than segregated into parcels. In addition, we show that inter-individual variations in this gradual organisation predict variations in recollection memory better than a characterisation based on functional parcellation. These results present an important step forward in understanding the functional organisation of the human hippocampus and have important implications for translating between rodent and human research.

Human V2A: A map of the peripheral visual hemifield with functional connections to scene-selective cortex.

Humans can recognize a scene in the blink of an eye. This gist-based visual scene perception is thought to be underpinned by specialized visual processing emphasizing the visual periphery at a cortical locus relatively low in the visual processing hierarchy. Using wide-field retinotopic mapping and population receptive field (pRF) modeling, we identified a new visual hemifield map anterior of area V2d and inferior to area V6, which we propose to call area V2A. Based on its location relative to other visual areas, V2A may correspond to area 23V described in nonhuman primates. The pRF analysis revealed unique receptive field properties for V2A: a large (FWHM ∼23°) and constant receptive field size across the central ∼70° of the visual field. Resting-state fMRI connectivity analysis further suggests that V2A is ideally suited to quickly feed the scene-processing network with information that is not biased towards the center of the visual field. Our findings not only indicate a likely cortical locus for the initial stages of gist-based visual scene perception, but also suggest a reappraisal of the organization of human dorsomedial occipital cortex with a strip of separate hemifield representations anterior to the early visual areas (V1, V2d, and V3d).

Thresholding functional connectomes by means of mixture modeling.

Functional connectivity has been shown to be a very promising tool for studying the large-scale functional architecture of the human brain. In network research in fMRI, functional connectivity is considered as a set of pair-wise interactions between the nodes of the network. These interactions are typically operationalized through the full or partial correlation between all pairs of regional time series. Estimating the structure of the latent underlying functional connectome from the set of pair-wise partial correlations remains an open research problem though. Typically, this thresholding problem is approached by proportional thresholding, or by means of parametric or non-parametric permutation testing across a cohort of subjects at each possible connection. As an alternative, we propose a data-driven thresholding approach for network matrices on the basis of mixture modeling. This approach allows for creating subject-specific sparse connectomes by modeling the full set of partial correlations as a mixture of low correlation values associated with weak or unreliable edges in the connectome and a sparse set of reliable connections. Consequently, we propose to use alternative thresholding strategy based on the model fit using pseudo-False Discovery Rates derived on the basis of the empirical null estimated as part of the mixture distribution. We evaluate the method on synthetic benchmark fMRI datasets where the underlying network structure is known, and demonstrate that it gives improved performance with respect to the alternative methods for thresholding connectomes, given the canonical thresholding levels. We also demonstrate that mixture modeling gives highly reproducible results when applied to the functional connectomes of the visual system derived from the n-back Working Memory task in the Human Connectome Project. The sparse connectomes obtained from mixture modeling are further discussed in the light of the previous knowledge of the functional architecture of the visual system in humans. We also demonstrate that with use of our method, we are able to extract similar information on the group level as can be achieved with permutation testing even though these two methods are not equivalent. We demonstrate that with both of these methods, we obtain functional decoupling between the two hemispheres in the higher order areas of the visual cortex during visual stimulation as compared to the resting state, which is in line with previous studies suggesting lateralization in the visual processing. However, as opposed to permutation testing, our approach does not require inference at the cohort level and can be used for creating sparse connectomes at the level of a single subject.

Connectivity-Based Parcellation of the Amygdala Predicts Social Skills in Adolescents with Autism Spectrum Disorder.

Amygdala dysfunction plays a role in the social impairments in autism spectrum disorders (ASD), but it is unclear which of its subregions are abnormal in ASD. This study compared the volume and functional connectivity (FC) strength of three FC-defined amygdala subregions between ASD and controls, and assessed their relation to social skills in ASD. A subregion associated with the social perception network was enlarged in ASD (F1 = 7.842, p = .008) and its volume correlated significantly with symptom severity (social skills: r = .548, p = .009). Posthoc analysis revealed that the enlargement was driven by the vmPFC amygdala network. These findings refine our understanding of abnormal amygdala connectivity in ASD and may inform future strategies for therapeutic interventions targeting the amygdalofrontal pathway.