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.

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.

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.

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.

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.

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.

Connectopic mapping with resting-state fMRI.

Brain regions are often topographically connected: nearby locations within one brain area connect with nearby locations in another area. Mapping these connection topographies, or 'connectopies' in short, is crucial for understanding how information is processed in the brain. Here, we propose principled, fully data-driven methods for mapping connectopies using functional magnetic resonance imaging (fMRI) data acquired at rest by combining spectral embedding of voxel-wise connectivity 'fingerprints' with a novel approach to spatial statistical inference. We apply the approach in human primary motor and visual cortex, and show that it can trace biologically plausible, overlapping connectopies in individual subjects that follow these regions' somatotopic and retinotopic maps. As a generic mechanism to perform inference over connectopies, the new spatial statistics approach enables rigorous statistical testing of hypotheses regarding the fine-grained spatial profile of functional connectivity and whether that profile is different between subjects or between experimental conditions. The combined framework offers a fundamental alternative to existing approaches to investigating functional connectivity in the brain, from voxel- or seed-pair wise characterizations of functional association, towards a full, multivariate characterization of spatial topography.

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.

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).

Preserved retinotopic brain connectivity in macular degeneration.

PURPOSE: The eye disease macular degeneration (MD) is a leading cause of blindness worldwide. There is no cure for MD, but several promising treatments aimed at restoring vision at the level of the retina are currently under investigation. These treatments assume that the patient's brain can still process appropriately the retinal input once it is restored, but whether this assumption is correct has yet to be determined. METHODS: We used functional magnetic resonance imaging (fMRI) and connective field modelling to determine whether the functional connectivity between the input-deprived portions of primary visual cortex (V1) and early extrastriate areas (V2/3) is still retinotopically organised. Specifically, in both patients with juvenile macular degeneration and age-matched controls with simulated retinal lesions, we assessed the extent to which the V1-referred connective fields of extrastriate voxels, as estimated on the basis of spontaneous fMRI signal fluctuations, adhered to retinotopic organisation. RESULTS: We found that functional connectivity between the input-deprived portions of visual areas V1 and extrastriate cortex is still largely retinotopically organised in MD, although on average less so than in controls. Patients with stable fixation exhibited normal retinotopic connectivity, however, suggesting that for the patients with unstable fixation, eye-movements resulted in spurious, homogeneous signal modulations across the entire input-deprived cortex, which would have hampered our ability to assess their spatial structure of connectivity. CONCLUSIONS: Despite the prolonged loss of visual input due to MD, the cortico-cortical connections of input-deprived visual cortex remain largely intact. This suggests that the restoration of sight in macular degeneration can rely on a largely unchanged retinotopic representation in early visual cortex following loss of central retinal function.

Equalization and decorrelation in primary visual cortex.

There are many theories on the purpose of neural adaptation, but evidence remains elusive. Here, we discuss the recent work by Benucci et al. (Nat Neurosci 16: 724-729, 2013), who measured for the first time the immediate effects of adaptation on the overall activity of a neuronal population. These measurements confirm two long-standing hypotheses about the purpose of adaptation, namely that adaptation counteracts biases in the statistics of the environment, and that it maintains decorrelation in neuronal stimulus selectivity.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization.

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Connective field modeling.

The traditional way to study the properties of visual neurons is to measure their responses to visually presented stimuli. A second way to understand visual neurons is to characterize their responses in terms of activity elsewhere in the brain. Understanding the relationships between responses in distinct locations in the visual system is essential to clarify this network of cortical signaling pathways. Here, we describe and validate connective field modeling, a model-based analysis for estimating the dependence between signals in distinct cortical regions using functional magnetic resonance imaging (fMRI). Just as the receptive field of a visual neuron predicts its response as a function of stimulus position, the connective field of a neuron predicts its response as a function of activity in another part of the brain. Connective field modeling opens up a wide range of research opportunities to study information processing in the visual system and other topographically organized cortices.

Linking functional and structural brain organisation with behaviour in autism: a multimodal EU-AIMS Longitudinal European Autism Project (LEAP) study.

Neuroimaging analyses of brain structure and function in autism have typically been conducted in isolation, missing the sensitivity gains of linking data across modalities. Here we focus on the integration of structural and functional organisational properties of brain regions. We aim to identify novel brain-organisation phenotypes of autism. We utilised multimodal MRI (T1-, diffusion-weighted and resting state functional), behavioural and clinical data from the EU AIMS Longitudinal European Autism Project (LEAP) from autistic (n = 206) and non-autistic (n = 196) participants. Of these, 97 had data from 2 timepoints resulting in a total scan number of 466. Grey matter density maps, probabilistic tractography connectivity matrices and connectopic maps were extracted from respective MRI modalities and were then integrated with Linked Independent Component Analysis. Linear mixed-effects models were used to evaluate the relationship between components and group while accounting for covariates and non-independence of participants with longitudinal data. Additional models were run to investigate associations with dimensional measures of behaviour. We identified one component that differed significantly between groups (coefficient = 0.33, padj = 0.02). This was driven (99%) by variance of the right fusiform gyrus connectopic map 2. While there were multiple nominal (uncorrected p < 0.05) associations with behavioural measures, none were significant following multiple comparison correction. Our analysis considered the relative contributions of both structural and functional brain phenotypes simultaneously, finding that functional phenotypes drive associations with autism. These findings expanded on previous unimodal studies by revealing the topographic organisation of functional connectivity patterns specific to autism and warrant further investigation.

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.

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.

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.

Correction: Protocol of the Healthy Brain Study: An accessible resource for understanding the human brain and how it dynamically and individually operates in its bio-social context.

[This corrects the article DOI: 10.1371/journal.pone.0260952.].

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.