Discovering structure in the space of fMRI selectivity profiles

Distributed representations of behaviour-derived object dimensions in the human visual system

Object vision is commonly thought to involve a hierarchy of brain regions processing increasingly complex image features, with high-level visual cortex supporting object recognition and categorization. However, object vision supports diverse behavioural goals, suggesting basic limitations of this category-centric framework. To address these limitations, we mapped a series of dimensions derived from a large-scale analysis of human similarity judgements directly onto the brain. Our results reveal broadly distributed representations of behaviourally relevant information, demonstrating selectivity to a wide variety of novel dimensions while capturing known selectivities for visual features and categories. Behaviour-derived dimensions were superior to categories at predicting brain responses, yielding mixed selectivity in much of visual cortex and sparse selectivity in category-selective clusters. This framework reconciles seemingly disparate findings regarding regional specialization, explaining category selectivity as a special case of sparse response profiles among representational dimensions, suggesting a more expansive view on visual processing in the human brain.

A review of ADHD detection studies with machine learning methods using rsfMRI data

Attention deficit hyperactivity disorder (ADHD) is a common mental health condition that significantly affects school-age children, causing difficulties with learning and daily functioning. Early identification is crucial, and reliable and objective diagnostic tools are necessary. However, current clinical evaluations of behavioral symptoms can be inconsistent and subjective. Functional magnetic resonance imaging (fMRI) is a non-invasive technique that has proven effective in detecting brain abnormalities in individuals with ADHD. Recent studies have shown promising outcomes in using resting state fMRI (rsfMRI)-based brain functional networks to diagnose various brain disorders, including ADHD. Several review papers have examined the detection of other diseases using fMRI data and machine learning or deep learning methods. However, no review paper has specifically addressed ADHD. Therefore, this study aims to contribute to the literature by reviewing the use of rsfMRI data and machine learning methods for detection of ADHD. The study provides general information about fMRI databases and detailed knowledge of the ADHD-200 database, which is commonly used for ADHD detection. It also emphasizes the importance of examining all stages of the process, including network and atlas selection, feature extraction, and feature selection, before the classification stage. The study compares the performance, advantages, and disadvantages of previous studies in detail. This comprehensive approach may be a useful starting point for new researchers in this area.

Distributed representations of behavior-derived object dimensions in the human visual system

Object vision is commonly thought to involve a hierarchy of brain regions processing increasingly complex image features, with high-level visual cortex supporting object recognition and categorization. However, object vision supports diverse behavioral goals, suggesting basic limitations of this category-centric framework. To address these limitations, we mapped a series of dimensions derived from a large-scale analysis of human similarity judgments directly onto the brain. Our results reveal broadly distributed representations of behaviorally-relevant information, demonstrating selectivity to a wide variety of novel dimensions while capturing known selectivities for visual features and categories. Behavior-derived dimensions were superior to categories at predicting brain responses, yielding mixed selectivity in much of visual cortex and sparse selectivity in category-selective clusters. This framework reconciles seemingly disparate findings regarding regional specialization, explaining category selectivity as a special case of sparse response profiles among representational dimensions, suggesting a more expansive view on visual processing in the human brain.

Aberrant cortico-thalamo-cerebellar network interactions and their association with impaired cognitive functioning in patients with schizophrenia

Evidence indicating abnormal functional connectivity (FC) among the cortex, thalamus, and cerebellum in schizophrenia patients has increased. However, the role of the thalamus and cerebellum when integrated into intrinsic networks and how those integrated networks interact in schizophrenia patients are largely unknown. We generated an integrative network map by merging thalamic and cerebellar network maps, which were parcellated using a winner-take-all approach, onto a cortical network map. Using cognitive networks, the default mode network (DMN), the dorsal attention network (DAN), the salience network (SAL), and the central executive network (CEN) as regions of interest, the FC of 48 schizophrenia patients was compared with that of 57 healthy controls (HCs). The association between abnormal FC and cognitive impairment was also investigated in patients. FC was lower between the SAL-CEN, SAL-DMN, and DMN-CEN and within-CEN in schizophrenia patients than in HCs. Hypoconnectivity between the DMN-CEN was correlated with impaired cognition in schizophrenia patients. Our findings broadly suggest the plausible role of the thalamus and cerebellum in integrative intrinsic networks in patients, which may contribute to the disrupted triple network and cognitive dysmetria in schizophrenia.

Fine-grained functional parcellation maps of the infant cerebral cortex

Resting-state functional MRI (rs-fMRI) is widely used to examine the dynamic brain functional development of infants, but these studies typically require precise cortical parcellation maps, which cannot be directly borrowed from adult-based functional parcellation maps due to the substantial differences in functional brain organization between infants and adults. Creating infant-specific cortical parcellation maps is thus highly desired but remains challenging due to difficulties in acquiring and processing infant brain MRIs. In this study, we leveraged 1064 high-resolution longitudinal rs-fMRIs from 197 typically developing infants and toddlers from birth to 24 months who participated in the Baby Connectome Project to develop the first set of infant-specific, fine-grained, surface-based cortical functional parcellation maps. To establish meaningful cortical functional correspondence across individuals, we performed cortical co-registration using both the cortical folding geometric features and the local gradient of functional connectivity (FC). Then we generated both age-related and age-independent cortical parcellation maps with over 800 fine-grained parcels during infancy based on aligned and averaged local gradient maps of FC across individuals. These parcellation maps reveal complex functional developmental patterns, such as changes in local gradient, network size, and local efficiency, especially during the first 9 postnatal months. Our generated fine-grained infant cortical functional parcellation maps are publicly available at https://www.nitrc.org/projects/infantsurfatlas/ for advancing the pediatric neuroimaging field.

Semantic novelty modulates neural responses to visual change across the human brain

Our continuous visual experience in daily life is dominated by change. Previous research has focused on visual change due to stimulus motion, eye movements or unfolding events, but not their combined impact across the brain, or their interactions with semantic novelty. We investigate the neural responses to these sources of novelty during film viewing. We analyzed intracranial recordings in humans across 6328 electrodes from 23 individuals. Responses associated with saccades and film cuts were dominant across the entire brain. Film cuts at semantic event boundaries were particularly effective in the temporal and medial temporal lobe. Saccades to visual targets with high visual novelty were also associated with strong neural responses. Specific locations in higher-order association areas showed selectivity to either high or low-novelty saccades. We conclude that neural activity associated with film cuts and eye movements is widespread across the brain and is modulated by semantic novelty.

Frequency-specific brain network architecture in resting-state fMRI

The analysis of brain function in resting-state network (RSN) models, ascertained through the functional connectivity pattern of resting-state functional magnetic resonance imaging (rs-fMRI), is sufficiently powerful for studying large-scale functional integration of the brain. However, in RSN-based research, the network architecture has been regarded as the same through different frequency bands. Thus, here, we aimed to examined whether the network architecture changes with frequency. The blood oxygen level-dependent (BOLD) signal was decomposed into four frequency bands-ranging from 0.007 to 0.438 Hz-and the clustering algorithm was applied to each of them. The best clustering number was selected for each frequency band based on the overlap ratio with task activation maps. The results demonstrated that resting-state BOLD signals exhibited frequency-specific network architecture; that is, the networks finely subdivided in the lower frequency bands were integrated into fewer networks in higher frequency bands rather than reconfigured, and the default mode network and networks related to perception had sufficiently strong architecture to survive in an environment with a lower signal-to-noise ratio. These findings provide a novel framework to enable improved understanding of brain function through the multiband frequency analysis of ultra-slow rs-fMRI data.

Global diversity in individualized cortical network topography

Individualized cortical network topography (ICNT) varies between people and exhibits great variability in the association networks in the human brain. However, these findings were mainly discovered in Western populations. It remains unclear whether and how ICNT is shaped by the non-Western populations. Here, we leveraged a multisession hierarchical Bayesian model to define individualized functional networks in White American and Han Chinese populations with data from both US and Chinese Human Connectome Projects. We found that both the size and spatial topography of individualized functional networks differed between White American and Han Chinese groups, especially in the heteromodal association cortex (including the ventral attention, control, language, dorsal attention, and default mode networks). Employing a support vector machine, we then demonstrated that ethnicity-related ICNT diversity can be used to identify an individual's ethnicity with high accuracy (74%, pperm < 0.0001), with heteromodal networks contributing most to the classification. This finding was further validated through mass-univariate analyses with generalized additive models. Moreover, we reveal that the spatial heterogeneity of ethnic diversity in ICNT correlated with fundamental properties of cortical organization, including evolutionary cortical expansion, brain myelination, and cerebral blood flow. Altogether, this case study highlights a need for more globally diverse and publicly available neuroimaging datasets.

Increasing diversity in connectomics with the Chinese Human Connectome Project

Cultural differences and biological diversity play important roles in shaping human brain structure and function. To date, most large-scale multimodal neuroimaging datasets have been obtained primarily from people living in Western countries, omitting the crucial contrast with populations living in other regions. The Chinese Human Connectome Project (CHCP) aims to address these resource and knowledge gaps by acquiring imaging, genetic and behavioral data from a large sample of participants living in an Eastern culture. The CHCP collected multimodal neuroimaging data from healthy Chinese adults using a protocol comparable to that of the Human Connectome Project. Comparisons between the CHCP and Human Connectome Project revealed both commonalities and distinctions in brain structure, function and connectivity. The corresponding large-scale brain parcellations were highly reproducible across the two datasets, with the language processing task showing the largest differences. The CHCP dataset is publicly available in an effort to facilitate transcultural and cross-ethnic brain-mind studies.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

Spatial Stability of Functional Networks: A Measure to Assess the Robustness of Graph-Theoretical Metrics to Spatial Errors Related to Brain Parcellation

There is growing interest in studying human brain connectivity and in modelling the brain functional structure as a network. Brain network creation requires parcellation of the cerebral cortex to define nodes. Parcellation might be affected by possible errors due to inter- and intra-subject variability as a consequence of brain structural and physiological characteristics and shape variations related to ageing and diseases, acquisition noise, and misregistration. These errors could induce a knock-on effect on network measure variability. The aim of this study was to investigate spatial stability, a measure of functional connectivity variations induced by parcellation errors. We simulated parcellation variability with random small spatial changes and evaluated its effects on twenty-seven graph-theoretical measures. The study included subjects from three public online datasets. Two brain parcellations were performed using FreeSurfer with geometric atlases. Starting from these, 100 new parcellations were created by increasing the area of 30% of parcels, reducing the area of neighbour parcels, with a rearrangement of vertices. fMRI data were filtered with linear regression, CompCor, and motion correction. Adjacency matrices were constructed with 0.1, 0.2, 0.3, and 0.4 thresholds. Differences in spatial stability between datasets, atlases, and threshold were evaluated. The higher spatial stability resulted for Characteristic-path-length, Density, Transitivity, and Closeness-centrality, and the lower spatial stability resulted for Bonacich and Katz. Multivariate analysis showed a significant effect of atlas, datasets, and thresholds. Katz and Bonacich centrality, which was subject to larger variations, can be considered an unconventional graph measure, poorly implemented in the clinical field and not yet investigated for reliability assessment. Spatial stability (SS) is affected by threshold, and it decreases with increasing threshold for several measures. Moreover, SS seems to depend on atlas choice and scanning parameters. Our study highlights the importance of paying close attention to possible parcellation-related spatial errors, which may affect the reliability of functional connectivity measures.

Functional brain network community structure in childhood: Unfinished territories and fuzzy boundaries

Adult cortex is organized into distributed functional communities. Yet, little is known about community architecture of children's brains. Here, we uncovered the community structure of cortex in childhood using fMRI data from 670 children aged 9-11 years (48% female, replication sample n=544, 56% female) from the Adolescent Brain and Cognitive Development study. We first applied a data-driven community detection approach to cluster cortical regions into communities, then employed a generative model-based approach called the weighted stochastic block model to further probe community interactions. Children showed similar community structure to adults, as defined by Yeo and colleagues in 2011, in early-developing sensory and motor communities, but differences emerged in transmodal areas. Children have more cortical territory in the limbic community, which is involved in emotion processing, than adults. Regions in association cortex interact more flexibly across communities, creating uncertainty for the model-based assignment algorithm, and perhaps reflecting cortical boundaries that are not yet solidified. Uncertainty was highest for cingulo-opercular areas involved in flexible deployment of cognitive control. Activation and deactivation patterns during a working memory task showed that both the data-driven approach and a set of adult communities statistically capture functional organization in middle childhood. Collectively, our findings suggest that community boundaries are not solidified by middle childhood.

Autism Spectrum Disorder Studies Using fMRI Data and Machine Learning: A Review

Machine learning methods have been frequently applied in the field of cognitive neuroscience in the last decade. A great deal of attention has been attracted to introduce machine learning methods to study the autism spectrum disorder (ASD) in order to find out its neurophysiological underpinnings. In this paper, we presented a comprehensive review about the previous studies since 2011, which applied machine learning methods to analyze the functional magnetic resonance imaging (fMRI) data of autistic individuals and the typical controls (TCs). The all-round process was covered, including feature construction from raw fMRI data, feature selection methods, machine learning methods, factors for high classification accuracy, and critical conclusions. Applying different machine learning methods and fMRI data acquired from different sites, classification accuracies were obtained ranging from 48.3% up to 97%, and informative brain regions and networks were located. Through thorough analysis, high classification accuracies were found to usually occur in the studies which involved task-based fMRI data, single dataset for some selection principle, effective feature selection methods, or advanced machine learning methods. Advanced deep learning together with the multi-site Autism Brain Imaging Data Exchange (ABIDE) dataset became research trends especially in the recent 4 years. In the future, advanced feature selection and machine learning methods combined with multi-site dataset or easily operated task-based fMRI data may appear to have the potentiality to serve as a promising diagnostic tool for ASD.

Brain parcellation selection: An overlooked decision point with meaningful effects on individual differences in resting-state functional connectivity

Over the past decade extensive research has examined the segregation of the human brain into large-scale functional networks. The resulting network maps, i.e. parcellations, are now commonly used for the a priori identification of functional networks. However, the use of these parcellations, particularly in developmental and clinical samples, hinges on four fundamental assumptions: (1) the various parcellations are equally able to recover the networks of interest; (2) adult-derived parcellations well represent the networks in children’s brains; (3) network properties, such as within-network connectivity, are reliably measured across parcellations; and (4) parcellation selection does not impact the results with regard to individual differences in given network properties. In the present study we examined these assumptions using eight common parcellation schemes in two independent developmental samples. We found that the parcellations are equally able to capture networks of interest in both children and adults. However, networks bearing the same name across parcellations (e.g., default network) do not produce reliable within-network measures of functional connectivity. Critically, parcellation selection significantly impacted the magnitude of associations of functional connectivity with age, poverty, and cognitive ability, producing meaningful differences in interpretation of individual differences in functional connectivity based on parcellation choice. Our findings suggest that work employing parcellations may benefit from the use of multiple schemes to confirm the robustness and generalizability of results. Furthermore, researchers looking to gain insight into functional networks may benefit from employing more nuanced network identification approaches such as using densely-sampled data to produce individual-derived network parcellations. A transition towards precision neuroscience will provide new avenues in the characterization of functional brain organization across development and within clinical populations.

Functional connectivity of EEG is subject-specific, associated with phenotype, and different from fMRI

A variety of psychiatric, behavioral and cognitive phenotypes have been linked to brain ''functional connectivity'' -- the pattern of correlation observed between different brain regions. Most commonly assessed using functional magnetic resonance imaging (fMRI), here, we investigate the connectivity-phenotype associations with functional connectivity measured with electroencephalography (EEG), using phase-coupling. We analyzed data from the publicly available Healthy Brain Network Biobank. This database compiles a growing sample of children and adolescents, currently encompassing 1657 individuals. Among a variety of assessment instruments we focus on ten phenotypic and additional demographic measures that capture most of the variance in this sample. The largest effect sizes are found for age and sex for both fMRI and EEG. We replicate previous findings of an association of Intelligence Quotient (IQ) and Attention Deficit Hyperactivity Disorder (ADHD) with the pattern of fMRI functional connectivity. We also find an association with socioeconomic status, anxiety and the Child Behavior Checklist Score. For EEG we find a significant connectivity-phenotype relationship with IQ. The actual spatial patterns of functional connectivity are quite different between fMRI and source-space EEG. However, within EEG we observe clusters of functional connectivity that are consistent across frequency bands. Additionally we analyzed reproducibility of functional connectivity. We compare connectivity obtained with different tasks, including resting state, a video and a visual flicker task. For both EEG and fMRI the variation between tasks was smaller than the variability observed between subjects. We also found an increase of reliability with increasing frequency of the EEG, and increased sampling duration. We conclude that, while the patterns of functional connectivity are distinct between fMRI and phase-coupling of EEG, they are nonetheless similar in their robustness to the task, and similar in that idiosyncratic patterns of connectivity predict individual phenotypes.

Sociality and interaction envelope organize visual action representations

Humans observe a wide range of actions in their surroundings. How is the visual cortex organized to process this diverse input? Using functional neuroimaging, we measured brain responses while participants viewed short videos of everyday actions, then probed the structure in these responses using voxel-wise encoding modeling. Responses are well fit by feature spaces that capture the body parts involved in an action and the action’s targets (i.e. whether the action was directed at an object, another person, the actor, and space). Clustering analyses reveal five large-scale networks that summarize the voxel tuning: one related to social aspects of an action, and four related to the scale of the interaction envelope, ranging from fine-scale manipulations directed at objects, to large-scale whole-body movements directed at distant locations. We propose that these networks reveal the major representational joints in how actions are processed by visual regions of the brain.

How is action perception organized in the brain? Here, the authors report evidence for five networks tuned to actions’ social content and the scale of their effect on the world and propose that sociality and interaction envelope are organizing dimensions of visual action representation.

A connectome-based mechanistic model of focal cortical dysplasia

Neuroimaging studies have consistently shown distributed brain anomalies in epilepsy syndromes associated with a focal structural lesion, particularly mesiotemporal sclerosis. Conversely, a system-level approach to focal cortical dysplasia has been rarely considered, likely due to methodological difficulties in addressing variable location and topography. Given the known heterogeneity in focal cortical dysplasia histopathology, we hypothesized that lesional connectivity consists of subtypes with distinct structural signatures. Furthermore, in light of mounting evidence for focal anomalies impacting whole-brain systems, we postulated that patterns of focal cortical dysplasia connectivity may exert differential downstream effects on global network topology. We studied a cohort of patients with histologically verified focal cortical dysplasia type II (n = 27), and age- and sex-matched healthy controls (n = 34). We subdivided each lesion into similarly sized parcels and computed their connectivity to large-scale canonical functional networks (or communities). We then dichotomized connectivity profiles of lesional parcels into those belonging to the same functional community as the focal cortical dysplasia (intra-community) and those adhering to other communities (inter-community). Applying hierarchical clustering to community-reconfigured connectome profiles identified three lesional classes with distinct patterns of functional connectivity: decreased intra- and inter-community connectivity, a selective decrease in intra-community connectivity, and increased intra- as well as inter-community connectivity. Hypo-connectivity classes were mainly composed of focal cortical dysplasia type IIB, while the hyperconnected lesions were type IIA. With respect to whole-brain networks, patients with hypoconnected focal cortical dysplasia and marked structural damage showed only mild imbalances, while those with hyperconnected subtle lesions had more pronounced topological alterations. Correcting for interictal epileptic discharges did not impact connectivity patterns. Multivariate structural equation analysis provided a mechanistic model of such complex, diverging interactions, whereby the focal cortical dysplasia structural makeup shapes its functional connectivity, which in turn modulates whole-brain network topology.

The parietal operculum preferentially encodes heat pain and not salience

Substantial controversy exists as to which part of brain activity is genuinely attributable to pain-related percepts and which activity is due to general aspects of sensory stimulation, such as its salience, or the accompanying arousal. The challenge posed by this question rests largely in the fact that pain per se exhibits highly intense but unspecific characteristics. These therefore should be matched by potential control conditions. Here, we used a unique combination of functional magnetic resonance imaging (fMRI) and behavioral and autonomic measures to address this longstanding debate in pain research. Subjects rated perceived intensity in a sequence alternating between heat and sound stimuli. Neuronal activity was monitored using fMRI. Either modality was presented in 6 different intensities, 3 of which lay above the pain threshold (for heat) or the unpleasantness threshold (for sound). We performed our analysis on 26 volunteers in which psychophysiological responses (as per skin conductance responses [SCRs]) did not differ between the 2 stimulus modalities. Having thus ascertained a comparable amount of stimulation-related but unspecific activation, we analyzed stimulus-response functions (SRFs) after painful stimulation and contrasted them with those of the matched acoustic control condition. Furthermore, analysis of fMRI data was performed on the brain surface to circumvent blurring issues stemming from the close proximity of several regions of interest located in heavily folded cortical areas. We focused our analyses on insular and peri-insular regions that are strongly involved in processing of painful stimuli. We employed an axiomatic approach to determine areas showing higher activation in painful compared to nonpainful heat and, at the same time, showing a steeper SRF for painful heat compared to unpleasant sound. Intriguingly, an area in the posterior parietal operculum emerged, whose response showed a pain preference after satisfying all axiomatic constraints. This result has important implications for the interpretation of functional imaging findings in pain research, because it clearly demonstrates that there are areas where activity following painful stimulation is not due to general attributes or results of sensory stimulation, such as salience or arousal. Conversely, several areas did not conform to the formulated axioms to rule out general factors as explanations.

The brain activity detected during pain could be due merely to the fact that pain is arousing and attention-grabbing, rather than being directly attributable to the pain itself. This study identifies an area of the brain — the parietal operculum — whose activity can only be explained by the painfulness of pain.

Functional Brain Parcellations of the Infant Brain and the Associated Developmental Trends

Resting-state functional connectivity studies have dramatically improved our understanding of the early human brain functional development during the past decade. However, one emerging problem that could potentially impede future progresses in the field is the definition of regions of interest (ROI), since it is well known that functional connectivity estimation can be seriously contaminated by within-ROI signal heterogeneity. In this study, based on a large-scale rsfMRI data set in human infants (230 neonates, 143 1-year olds, and 107 2-year olds), we aimed to derive a set of anatomically constrained, infant-specific functional brain parcellations using functional connectivity-based clustering. Our results revealed significantly higher levels of signal homogeneity within the newly defined functional parcellations compared with other schemes. Importantly, the global functional connectivity patterns associated with the newly defined functional subunits demonstrated significantly increasing levels of differentiation with age, confirming increasing levels of local specialization. Subsequent whole brain connectivity analysis revealed intriguing patterns of regional-level functional connectivity developments and system-level hub redistribution during infancy. Overall, the newly derived infant-specific functional brain parcellations and the associated novel developmental patterns will likely prove valuable for future early developmental studies using the functional connectivity technique.

Spatial Topography of Individual-Specific Cortical Networks Predicts Human Cognition, Personality, and Emotion

Resting-state functional magnetic resonance imaging (rs-fMRI) offers the opportunity to delineate individual-specific brain networks. A major question is whether individual-specific network topography (i.e., location and spatial arrangement) is behaviorally relevant. Here, we propose a multi-session hierarchical Bayesian model (MS-HBM) for estimating individual-specific cortical networks and investigate whether individual-specific network topography can predict human behavior. The multiple layers of the MS-HBM explicitly differentiate intra-subject (within-subject) from inter-subject (between-subject) network variability. By ignoring intra-subject variability, previous network mappings might confuse intra-subject variability for inter-subject differences. Compared with other approaches, MS-HBM parcellations generalized better to new rs-fMRI and task-fMRI data from the same subjects. More specifically, MS-HBM parcellations estimated from a single rs-fMRI session (10 min) showed comparable generalizability as parcellations estimated by 2 state-of-the-art methods using 5 sessions (50 min). We also showed that behavioral phenotypes across cognition, personality, and emotion could be predicted by individual-specific network topography with modest accuracy, comparable to previous reports predicting phenotypes based on connectivity strength. Network topography estimated by MS-HBM was more effective for behavioral prediction than network size, as well as network topography estimated by other parcellation approaches. Thus, similar to connectivity strength, individual-specific network topography might also serve as a fingerprint of human behavior.

Hierarchy of speech-driven spectrotemporal receptive fields in human auditory cortex

Existing data indicate that cortical speech processing is hierarchically organized. Numerous studies have shown that early auditory areas encode fine acoustic details while later areas encode abstracted speech patterns. However, it remains unclear precisely what speech information is encoded across these hierarchical levels. Estimation of speech-driven spectrotemporal receptive fields (STRFs) provides a means to explore cortical speech processing in terms of acoustic or linguistic information associated with characteristic spectrotemporal patterns. Here, we estimate STRFs from cortical responses to continuous speech in fMRI. Using a novel approach based on filtering randomly-selected spectrotemporal modulations (STMs) from aurally-presented sentences, STRFs were estimated for a group of listeners and categorized using a data-driven clustering algorithm. 'Behavioral STRFs' highlighting STMs crucial for speech recognition were derived from intelligibility judgments. Clustering revealed that STRFs in the supratemporal plane represented a broad range of STMs, while STRFs in the lateral temporal lobe represented circumscribed STM patterns important to intelligibility. Detailed analysis recovered a bilateral organization with posterior-lateral regions preferentially processing STMs associated with phonological information and anterior-lateral regions preferentially processing STMs associated with word- and phrase-level information. Regions in lateral Heschl's gyrus preferentially processed STMs associated with vocalic information (pitch).

Multiple Frequency Bands Analysis of Large Scale Intrinsic Brain Networks and Its Application in Schizotypal Personality Disorder

The human brain is a complex system composed by several large scale intrinsic networks with distinct functions. The low frequency oscillation (LFO) signal of blood oxygen level dependent (BOLD), measured through resting-state fMRI, reflects the spontaneous neural activity of these networks. We propose to characterize these networks by applying the multiple frequency bands analysis (MFBA) to the LFO time courses (TCs) resulted from the group independent component analysis (ICA). Specifically, seven networks, including the default model network (DMN), dorsal attention network (DAN), control executive network (CEN), salience network, sensorimotor network, visual network and limbic network, are identified. After the power spectral density (PSD) analysis, the amplitude of low frequency fluctuation (ALFF) and the fractional amplitude of low frequency fluctuation (fALFF) is determined in three bands: <0.1 Hz; slow-5; and slow-4. Moreover, the MFBA method is applied to reveal the frequency-dependent alternations of fALFF for seven networks in schizotypal personality disorder (SPD). It is found that seven networks can be divided into three categories: the advanced cognitive networks, primary sensorimotor networks and limbic networks, and their fALFF successively decreases in both slow-4 and slow-5 bands. Comparing to normal control group, the fALFF of DMN, DAN and CEN in SPD tends to be higher in slow-5 band, but lower in slow-4. Higher fALFF of sensorimotor and visual networks in slow-5, higher fALFF of limbic network in both bands have been observed for SPD group. The results of ALFF are consistent with those of fALFF. The proposed MFBA method may help distinguish networks or oscillators in the human brain, reveal subtle alternations of networks through locating their dominant frequency band, and present potential to interpret the neuropathology disruptions.

Network Analysis in Disorders of Consciousness: Four Problems and One Proposed Solution (Exponential Random Graph Models)

In recent years, the study of the neural basis of consciousness, particularly in the context of patients recovering from severe brain injury, has greatly benefited from the application of sophisticated network analysis techniques to functional brain data. Yet, current graph theoretic approaches, as employed in the neuroimaging literature, suffer from four important shortcomings. First, they require arbitrary fixing of the number of connections (i.e., density) across networks which are likely to have different "natural" (i.e., stable) density (e.g., patients vs. controls, vegetative state vs. minimally conscious state patients). Second, when describing networks, they do not control for the fact that many characteristics are interrelated, particularly some of the most popular metrics employed (e.g., nodal degree, clustering coefficient)-which can lead to spurious results. Third, in the clinical domain of disorders of consciousness, there currently are no methods for incorporating structural connectivity in the characterization of functional networks which clouds the interpretation of functional differences across groups with different underlying pathology as well as in longitudinal approaches where structural reorganization processes might be operating. Finally, current methods do not allow assessing the dynamics of network change over time. We present a different framework for network analysis, based on Exponential Random Graph Models, which overcomes the above limitations and is thus particularly well suited for clinical populations with disorders of consciousness. We demonstrate this approach in the context of the longitudinal study of recovery from coma. First, our data show that throughout recovery from coma, brain graphs vary in their natural level of connectivity (from 10.4 to 14.5%), which conflicts with the standard approach of imposing arbitrary and equal density thresholds across networks (e.g., time-points, subjects, groups). Second, we show that failure to consider the interrelation between network measures does lead to spurious characterization of both inter- and intra-regional brain connectivity. Finally, we show that Separable Temporal ERGM can be employed to describe network dynamics over time revealing the specific pattern of formation and dissolution of connectivity that accompany recovery from coma.

Reliable new measures capturing low-frequency fluctuations from resting-state functional MRI

Resting-state functional MRI (rsfMRI) is one of the most important neuroimaging modalities for investigating alterations in the resting-state networks of the human brain, given that abnormal neural activity during the resting state is associated with neurological disorders. However, neuroimaging results obtained from rsfMRI have rarely been replicated with repeated measurements. Therefore, we aimed to develop new measures to extract highly reliable and reproducible functional neuroimaging metrics from rsfMRI data. Preprocessed rsfMRI data from 30 patients with 10 sessions of rsfMRI scans taken within 1 month were obtained from the Consortium for Reliability and Reproducibility. We developed a time-domain measure to capture low-frequency fluctuation (LFF) using a general linear model with three different periodic regressors: boxcar, triangular, and sinusoidal functions. Then, test-retest reliability for the proposed methods was evaluated using the intraclass correlation (ICC). Our approaches for evaluating LFF from rsfMRI data significantly identified the default mode network areas (corrected P<0.05). The regression model with the sinusoidal basis function produced the most reliable results (ICC=0.6) compared with the boxcar (ICC=0.32) or triangular (ICC=0.34) functions. Taken together, the proposed methods successfully identified the default mode network regions. In addition, our results suggest that new functional metrics aiming to extract LFF components by modeling rsfMRI time-series data might provide a reliable biomarker to identify neurological disorders accompanying abnormal functional activity.

Idiosyncratic organization of cortical networks in autism spectrum disorder

Neuroimaging studies of Autism Spectrum Disorder (ASD) have yielded inconsistent results indicating either increases or decreases in functional connectivity, or both. Recent findings suggest that these seemingly divergent results might be underpinned by greater inter-individual variability in brain network connectivity in ASD. We tested the hypothesis that the spatial patterns of intrinsic connectivity networks (ICNs) are more idiosyncratic in ASD, and demonstrated that this increased variability is associated with symptomatology. We estimated whole brain functional connectivity based on resting state functional magnetic resonance imaging (fMRI) data obtained from the Autism Brain Imaging Data Exchange I & II (ABIDE I & II) repository: 422 (69 females) participants with ASD and 424 (59 females) typically developing (TD) participants between 6 and 30 years of age. We clustered individuals' patterns of resting state functional connectivity into seven networks, each representing an ICN, and assessed the heterogeneity of each vertex on the cortical surface across individuals in terms of its incorporation into a particular ICN. We found that the incorporation of individual anatomical locations (vertices) to a common network was less consistent across individuals in ASD, indicating a more idiosyncratic organization of ICNs in the ASD brain. This spatial shifting effect was particularly pronounced in the Sensory-Motor Network (SMN) and the Default Mode Network (DMN). We also found that this idiosyncrasy in large-scale brain network organization was correlated with ASD symptomatology (ADOS). These results support the view that idiosyncratic functional connectivity is a hallmark of the ASD brain. We provide the first evidence that the anatomical organization of ICNs is idiosyncratic in ASD, as well as providing evidence that such abnormalities in brain network organization may contribute to the symptoms of ASD.

Control networks and hubs

Executive control functions are associated with frontal, parietal, cingulate, and insular brain regions that interact through distributed large-scale networks. Here, we discuss how fMRI functional connectivity can shed light on the organization of control networks and how they interact with other parts of the brain. In the first section of our review, we present convergent evidence from fMRI functional connectivity, activation, and lesion studies that there are multiple dissociable control networks in the brain with distinct functional properties. In the second section, we discuss how graph theoretical concepts can help illuminate the mechanisms by which control networks interact with other brain regions to carry out goal-directed functions, focusing on the role of specialized hub regions for mediating cross-network interactions. Again, we use a combination of functional connectivity, lesion, and task activation studies to bolster this claim. We conclude that a large-scale network perspective provides important neurobiological constraints on the neural underpinnings of executive control, which will guide future basic and translational research into executive function and its disruption in disease.

Infinite von Mises–Fisher Mixture Modeling of Whole Brain fMRI Data

Cluster analysis of functional magnetic resonance imaging (fMRI) data is often performed using gaussian mixture models, but when the time series are standardized such that the data reside on a hypersphere, this modeling assumption is questionable. The consequences of ignoring the underlying spherical manifold are rarely analyzed, in part due to the computational challenges imposed by directional statistics. In this letter, we discuss a Bayesian von Mises–Fisher (vMF) mixture model for data on the unit hypersphere and present an efficient inference procedure based on collapsed Markov chain Monte Carlo sampling. Comparing the vMF and gaussian mixture models on synthetic data, we demonstrate that the vMF model has a slight advantage inferring the true underlying clustering when compared to gaussian-based models on data generated from both a mixture of vMFs and a mixture of gaussians subsequently normalized. Thus, when performing model selection, the two models are not in agreement. Analyzing multisubject whole brain resting-state fMRI data from healthy adult subjects, we find that the vMF mixture model is considerably more reliable than the gaussian mixture model when comparing solutions across models trained on different groups of subjects, and again we find that the two models disagree on the optimal number of components. The analysis indicates that the fMRI data support more than a thousand clusters, and we confirm this is not a result of overfitting by demonstrating better prediction on data from held-out subjects. Our results highlight the utility of using directional statistics to model standardized fMRI data and demonstrate that whole brain segmentation of fMRI data requires a very large number of functional units in order to adequately account for the discernible statistical patterns in the data.

A flexible graphical model for multi-modal parcellation of the cortex

Advances in neuroimaging have provided a tremendous amount of in-vivo information on the brain's organisation. Its anatomy and cortical organisation can be investigated from the point of view of several imaging modalities, many of which have been studied for mapping functionally specialised cortical areas. There is strong evidence that a single modality is not sufficient to fully identify the brain's cortical organisation. Combining multiple modalities in the same parcellation task has the potential to provide more accurate and robust subdivisions of the cortex. Nonetheless, existing brain parcellation methods are typically developed and tested on single modalities using a specific type of information. In this paper, we propose Graph-based Multi-modal Parcellation (GraMPa), an iterative framework designed to handle the large variety of available input modalities to tackle the multi-modal parcellation task. At each iteration, we compute a set of parcellations from different modalities and fuse them based on their local reliabilities. The fused parcellation is used to initialise the next iteration, forcing the parcellations to converge towards a set of mutually informed modality specific parcellations, where correspondences are established. We explore two different multi-modal configurations for group-wise parcellation using resting-state fMRI, diffusion MRI tractography, myelin maps and task fMRI. Quantitative and qualitative results on the Human Connectome Project database show that integrating multi-modal information yields a stronger agreement with well established atlases and more robust connectivity networks that provide a better representation of the population.

The impact of hemodynamic variability and signal mixing on the identifiability of effective connectivity structures in BOLD fMRI

PURPOSE

Multiple computational studies have demonstrated that essentially all current analytical approaches to determine effective connectivity perform poorly when applied to synthetic functional Magnetic Resonance Imaging (fMRI) datasets. In this study, we take a theoretical approach to investigate the potential factors facilitating and hindering effective connectivity research in fMRI.

MATERIALS AND METHODS

In this work, we perform a simulation study with use of Dynamic Causal Modeling generative model in order to gain new insights on the influence of factors such as the slow hemodynamic response, mixed signals in the network and short time series, on the effective connectivity estimation in fMRI studies.

RESULTS

First, we perform a Linear Discriminant Analysis study and find that not the hemodynamics itself but mixed signals in the neuronal networks are detrimental to the signatures of distinct connectivity patterns. This result suggests that for statistical methods (which do not involve lagged signals), deconvolving the BOLD responses is not necessary, but at the same time, functional parcellation into Regions of Interest (ROIs) is essential. Second, we study the impact of hemodynamic variability on the inference with use of lagged methods. We find that the local hemodynamic variability provide with an upper bound on the success rate of the lagged methods. Furthermore, we demonstrate that upsampling the data to TRs lower than the TRs in state-of-the-art datasets does not influence the performance of the lagged methods.

CONCLUSIONS

Factors such as background scale-free noise and hemodynamic variability have a major impact on the performance of methods for effective connectivity research in functional Magnetic Resonance Imaging.

Interpreting the dimensions of neural feature representations revealed by dimensionality reduction

Recent progress in understanding the structure of neural representations in the cerebral cortex has centred around the application of multivariate classification analyses to measurements of brain activity. These analyses have proved a sensitive test of whether given brain regions provide information about specific perceptual or cognitive processes. An exciting extension of this approach is to infer the structure of this information, thereby drawing conclusions about the underlying neural representational space. These approaches rely on exploratory data-driven dimensionality reduction to extract the natural dimensions of neural spaces, including natural visual object and scene representations, semantic and conceptual knowledge, and working memory. However, the efficacy of these exploratory methods is unknown, because they have only been applied to representations in brain areas for which we have little or no secondary knowledge. One of the best-understood areas of the cerebral cortex is area MT of primate visual cortex, which is known to be important in motion analysis. To assess the effectiveness of dimensionality reduction for recovering neural representational space we applied several dimensionality reduction methods to multielectrode measurements of spiking activity obtained from area MT of marmoset monkeys, made while systematically varying the motion direction and speed of moving stimuli. Despite robust tuning at individual electrodes, and high classifier performance, dimensionality reduction rarely revealed dimensions for direction and speed. We use this example to illustrate important limitations of these analyses, and suggest a framework for how to best apply such methods to data where the structure of the neural representation is unknown.

Human brain mapping: A systematic comparison of parcellation methods for the human cerebral cortex

The macro-connectome elucidates the pathways through which brain regions are structurally connected or functionally coupled to perform a specific cognitive task. It embodies the notion of representing and understanding all connections within the brain as a network, while the subdivision of the brain into interacting functional units is inherent in its architecture. As a result, the definition of network nodes is one of the most critical steps in connectivity network analysis. Although brain atlases obtained from cytoarchitecture or anatomy have long been used for this task, connectivity-driven methods have arisen only recently, aiming to delineate more homogeneous and functionally coherent regions. This study provides a systematic comparison between anatomical, connectivity-driven and random parcellation methods proposed in the thriving field of brain parcellation. Using resting-state functional MRI data from the Human Connectome Project and a plethora of quantitative evaluation techniques investigated in the literature, we evaluate 10 subject-level and 24 groupwise parcellation methods at different resolutions. We assess the accuracy of parcellations from four different aspects: (1) reproducibility across different acquisitions and groups, (2) fidelity to the underlying connectivity data, (3) agreement with fMRI task activation, myelin maps, and cytoarchitectural areas, and (4) network analysis. This extensive evaluation of different parcellations generated at the subject and group level highlights the strengths and shortcomings of the various methods and aims to provide a guideline for the choice of parcellation technique and resolution according to the task at hand. The results obtained in this study suggest that there is no optimal method able to address all the challenges faced in this endeavour simultaneously.

sGraSP: A graph-based method for the derivation of subject-specific functional parcellations of the brain

BACKGROUND

Resting-state fMRI (rs-fMRI) has emerged as a prominent tool for the study of functional connectivity. The identification of the regions associated with the different brain functions has received significant interest. However, most of the studies conducted so far have focused on the definition of a common set of regions, valid for an entire population. The variation of the functional regions within a population has rarely been accounted for.

NEW METHOD

In this paper, we propose sGraSP, a graph-based approach for the derivation of subject-specific functional parcellations. Our method generates first a common parcellation for an entire population, which is then adapted to each subject individually.

RESULTS

Several cortical parcellations were generated for 859 children being part of the Philadelphia Neurodevelopmental Cohort. The stability of the parcellations generated by sGraSP was tested by mixing population and subject rs-fMRI signals, to generate subject-specific parcels increasingly closer to the population parcellation. We also checked if the parcels generated by our method were better capturing a development trend underlying our data than the original parcels, defined for the entire population.

COMPARISON WITH EXISTING METHODS

We compared sGraSP with a simpler and faster approach based on a Voronoi tessellation, by measuring their ability to produce functionally coherent parcels adapted to the subject data.

CONCLUSIONS

Our parcellations outperformed the Voronoi tessellations. The parcels generated by sGraSP vary consistently with respect to signal mixing, the results are highly reproducible and the neurodevelopmental trend is better captured with the subject-specific parcellation, under all the signal mixing conditions.

Deriving reproducible biomarkers from multi-site resting-state data: An Autism-based example

Resting-state functional Magnetic Resonance Imaging (R-fMRI) holds the promise to reveal functional biomarkers of neuropsychiatric disorders. However, extracting such biomarkers is challenging for complex multi-faceted neuropathologies, such as autism spectrum disorders. Large multi-site datasets increase sample sizes to compensate for this complexity, at the cost of uncontrolled heterogeneity. This heterogeneity raises new challenges, akin to those face in realistic diagnostic applications. Here, we demonstrate the feasibility of inter-site classification of neuropsychiatric status, with an application to the Autism Brain Imaging Data Exchange (ABIDE) database, a large (N=871) multi-site autism dataset. For this purpose, we investigate pipelines that extract the most predictive biomarkers from the data. These R-fMRI pipelines build participant-specific connectomes from functionally-defined brain areas. Connectomes are then compared across participants to learn patterns of connectivity that differentiate typical controls from individuals with autism. We predict this neuropsychiatric status for participants from the same acquisition sites or different, unseen, ones. Good choices of methods for the various steps of the pipeline lead to 67% prediction accuracy on the full ABIDE data, which is significantly better than previously reported results. We perform extensive validation on multiple subsets of the data defined by different inclusion criteria. These enables detailed analysis of the factors contributing to successful connectome-based prediction. First, prediction accuracy improves as we include more subjects, up to the maximum amount of subjects available. Second, the definition of functional brain areas is of paramount importance for biomarker discovery: brain areas extracted from large R-fMRI datasets outperform reference atlases in the classification tasks.

Directional Clustering Through Matrix Factorization

This paper deals with a clustering problem where feature vectors are clustered depending on the angle between feature vectors, that is, feature vectors are grouped together if they point roughly in the same direction. This directional distance measure arises in several applications, including document classification and human brain imaging. Using ideas from the field of constrained low-rank matrix factorization and sparse approximation, a novel approach is presented that differs from classical clustering methods, such as seminonnegative matrix factorization, K -EVD, or k -means clustering, yet combines some aspects of all these. As in nonnegative matrix factorization and K -EVD, the matrix decomposition is iteratively refined to optimize a data fidelity term; however, no positivity constraint is enforced directly nor do we need to explicitly compute eigenvectors. As in k -means and K -EVD, each optimization step is followed by a hard cluster assignment. This leads to an efficient algorithm that is shown here to outperform common competitors in terms of clustering performance and/or computation speed. In addition to a detailed theoretical analysis of some of the algorithm's main properties, the approach is empirically evaluated on a range of toy problems, several standard text clustering data sets, and a high-dimensional problem in brain imaging, where functional magnetic resonance imaging data are used to partition the human cerebral cortex into distinct functional regions.

Group-wise parcellation of the cortex through multi-scale spectral clustering

The delineation of functionally and structurally distinct regions as well as their connectivity can provide key knowledge towards understanding the brain's behaviour and function. Cytoarchitecture has long been the gold standard for such parcellation tasks, but has poor scalability and cannot be mapped in vivo. Functional and diffusion magnetic resonance imaging allow in vivo mapping of brain's connectivity and the parcellation of the brain based on local connectivity information. Several methods have been developed for single subject connectivity driven parcellation, but very few have tackled the task of group-wise parcellation, which is essential for uncovering group specific behaviours. In this paper, we propose a group-wise connectivity-driven parcellation method based on spectral clustering that captures local connectivity information at multiple scales and directly enforces correspondences between subjects. The method is applied to diffusion Magnetic Resonance Imaging driven parcellation on two independent groups of 50 subjects from the Human Connectome Project. Promising quantitative and qualitative results in terms of information loss, modality comparisons, group consistency and inter-group similarities demonstrate the potential of the method.

The role of diversity in complex ICA algorithms for fMRI analysis

BACKGROUND

The widespread use of data-driven methods, such as independent component analysis (ICA), for the analysis of functional magnetic resonance imaging data (fMRI) has enabled deeper understanding of neural function. However, most popular ICA algorithms for fMRI analysis make several simplifying assumptions, thus ignoring sources of statistical information, types of "diversity," and limiting their performance.

NEW METHOD

We propose the use of complex entropy rate bound minimization (CERBM) for the analysis of actual fMRI data in its native, complex, domain. Though CERBM achieves enhanced performance through the exploitation of the three types of diversity inherent to complex fMRI data: noncircularity, non-Gaussianity, and sample-to-sample dependence, CERBM produces results that are more variable than simpler methods. This motivates the development of a minimum spanning tree (MST)-based stability analysis that mitigates the variability of CERBM.

COMPARISON WITH EXISTING METHODS

In order to validate our method, we compare the performance of CERBM with the popular CInfomax as well as complex entropy bound minimization (CEBM).

RESULTS

We show that by leveraging CERBM and the MST-based stability analysis, we are able to consistently produce components that have a greater number of activated voxels in physically meaningful regions and can more accurately classify patients with schizophrenia than components generated using simpler models.

CONCLUSIONS

Our results demonstrate the advantages of using ICA algorithms that can exploit all inherent types of diversity for the analysis of fMRI data when coupled with appropriate stability analyses.

Transcriptional profiles of supragranular-enriched genes associate with corticocortical network architecture in the human brain

The human brain is patterned with disproportionately large, distributed cerebral networks that connect multiple association zones in the frontal, temporal, and parietal lobes. The expansion of the cortical surface, along with the emergence of long-range connectivity networks, may be reflected in changes to the underlying molecular architecture. Using the Allen Institute's human brain transcriptional atlas, we demonstrate that genes particularly enriched in supragranular layers of the human cerebral cortex relative to mouse distinguish major cortical classes. The topography of transcriptional expression reflects large-scale brain network organization consistent with estimates from functional connectivity MRI and anatomical tracing in nonhuman primates. Microarray expression data for genes preferentially expressed in human upper layers (II/III), but enriched only in lower layers (V/VI) of mouse, were cross-correlated to identify molecular profiles across the cerebral cortex of postmortem human brains (n = 6). Unimodal sensory and motor zones have similar molecular profiles, despite being distributed across the cortical mantle. Sensory/motor profiles were anticorrelated with paralimbic and certain distributed association network profiles. Tests of alternative gene sets did not consistently distinguish sensory and motor regions from paralimbic and association regions: (i) genes enriched in supragranular layers in both humans and mice, (ii) genes cortically enriched in humans relative to nonhuman primates, (iii) genes related to connectivity in rodents, (iv) genes associated with human and mouse connectivity, and (v) 1,454 gene sets curated from known gene ontologies. Molecular innovations of upper cortical layers may be an important component in the evolution of long-range corticocortical projections.

The cortical analysis of speech-specific temporal structure revealed by responses to sound quilts

Speech contains temporal structure that the brain must analyze to enable linguistic processing. To investigate the neural basis of this analysis, we used sound quilts, stimuli constructed by shuffling segments of a natural sound, approximately preserving its properties on short timescales while disrupting them on longer scales. We generated quilts from foreign speech to eliminate language cues and manipulated the extent of natural acoustic structure by varying the segment length. Using functional magnetic resonance imaging, we identified bilateral regions of the superior temporal sulcus (STS) whose responses varied with segment length. This effect was absent in primary auditory cortex and did not occur for quilts made from other natural sounds or acoustically matched synthetic sounds, suggesting tuning to speech-specific spectrotemporal structure. When examined parametrically, the STS response increased with segment length up to ∼500 ms. Our results identify a locus of speech analysis in human auditory cortex that is distinct from lexical, semantic or syntactic processes.

Reconfigurable task-dependent functional coupling modes cluster around a core functional architecture

Functional coupling across distributed brain regions varies across task contexts, yet there are stable features. To better understand the range and central tendencies of network configurations, coupling patterns were explored using functional MRI (fMRI) across 14 distinct continuously performed task states ranging from passive fixation to increasingly demanding classification tasks. Mean global correlation profiles across the cortex ranged from 0.69 to 0.82 between task states. Network configurations from both passive fixation and classification tasks similarly predicted task coactivation patterns estimated from meta-analysis of the literature. Thus, even across markedly different task states, central tendencies dominate the coupling configurations. Beyond these shared components, distinct task states displayed significant differences in coupling patterns in response to their varied demands. One possibility is that anatomical connectivity provides constraints that act as attractors pulling network configurations towards a limited number of robust states. Reconfigurable coupling modes emerge as significant modifications to a core functional architecture.

GraSP: geodesic Graph-based Segmentation with Shape Priors for the functional parcellation of the cortex

Resting-state functional MRI is a powerful technique for mapping the functional organization of the human brain. However, for many types of connectivity analysis, high-resolution voxelwise analyses are computationally infeasible and dimensionality reduction is typically used to limit the number of network nodes. Most commonly, network nodes are defined using standard anatomic atlases that do not align well with functional neuroanatomy or regions of interest covering a small portion of the cortex. Data-driven parcellation methods seek to overcome such limitations, but existing approaches are highly dependent on initialization procedures and produce spatially fragmented parcels or overly isotropic parcels that are unlikely to be biologically grounded. In this paper, we propose a novel graph-based parcellation method that relies on a discrete Markov Random Field framework. The spatial connectedness of the parcels is explicitly enforced by shape priors. The shape of the parcels is adapted to underlying data through the use of functional geodesic distances. Our method is initialization-free and rapidly segments the cortex in a single optimization. The performance of the method was assessed using a large developmental cohort of more than 850 subjects. Compared to two prevalent parcellation methods, our approach provides superior reproducibility for a similar data fit. Furthermore, compared to other methods, it avoids incoherent parcels. Finally, the method's utility is demonstrated through its ability to detect strong brain developmental effects that are only weakly observed using other methods.

Decoupling function and anatomy in atlases of functional connectivity patterns: language mapping in tumor patients

In this paper we construct an atlas that summarizes functional connectivity characteristics of a cognitive process from a population of individuals. The atlas encodes functional connectivity structure in a low-dimensional embedding space that is derived from a diffusion process on a graph that represents correlations of fMRI time courses. The functional atlas is decoupled from the anatomical space, and thus can represent functional networks with variable spatial distribution in a population. In practice the atlas is represented by a common prior distribution for the embedded fMRI signals of all subjects. We derive an algorithm for fitting this generative model to the observed data in a population. Our results in a language fMRI study demonstrate that the method identifies coherent and functionally equivalent regions across subjects. The method also successfully maps functional networks from a healthy population used as a training set to individuals whose language networks are affected by tumors.

Which fMRI clustering gives good brain parcellations?

Analysis and interpretation of neuroimaging data often require one to divide the brain into a number of regions, or parcels, with homogeneous characteristics, be these regions defined in the brain volume or on the cortical surface. While predefined brain atlases do not adapt to the signal in the individual subject images, parcellation approaches use brain activity (e.g., found in some functional contrasts of interest) and clustering techniques to define regions with some degree of signal homogeneity. In this work, we address the question of which clustering technique is appropriate and how to optimize the corresponding model. We use two principled criteria: goodness of fit (accuracy), and reproducibility of the parcellation across bootstrap samples. We study these criteria on both simulated and two task-based functional Magnetic Resonance Imaging datasets for the Ward, spectral and k-means clustering algorithms. We show that in general Ward’s clustering performs better than alternative methods with regard to reproducibility and accuracy and that the two criteria diverge regarding the preferred models (reproducibility leading to more conservative solutions), thus deferring the practical decision to a higher level alternative, namely the choice of a trade-off between accuracy and stability.

Estimates of segregation and overlap of functional connectivity networks in the human cerebral cortex

The organization of the human cerebral cortex has recently been explored using techniques for parcellating the cortex into distinct functionally coupled networks. The divergent and convergent nature of cortico-cortical anatomic connections suggests the need to consider the possibility of regions belonging to multiple networks and hierarchies among networks. Here we applied the Latent Dirichlet Allocation (LDA) model and spatial independent component analysis (ICA) to solve for functionally coupled cerebral networks without assuming that cortical regions belong to a single network. Data analyzed included 1000 subjects from the Brain Genomics Superstruct Project (GSP) and 12 high quality individual subjects from the Human Connectome Project (HCP). The organization of the cerebral cortex was similar regardless of whether a winner-take-all approach or the more relaxed constraints of LDA (or ICA) were imposed. This suggests that large-scale networks may function as partially isolated modules. Several notable interactions among networks were uncovered by the LDA analysis. Many association regions belong to at least two networks, while somatomotor and early visual cortices are especially isolated. As examples of interaction, the precuneus, lateral temporal cortex, medial prefrontal cortex and posterior parietal cortex participate in multiple paralimbic networks that together comprise subsystems of the default network. In addition, regions at or near the frontal eye field and human lateral intraparietal area homologue participate in multiple hierarchically organized networks. These observations were replicated in both datasets and could be detected (and replicated) in individual subjects from the HCP.

Timing and tuning for familiarity of cortical responses to faces

Different kinds of known faces activate brain areas to dissimilar degrees. However, the tuning to type of knowledge, and the temporal course of activation, of each area have not been well characterized. Here we measured, with functional magnetic resonance imaging, brain activity elicited by unfamiliar, visually familiar, and personally-familiar faces. We assessed response amplitude and duration using flexible hemodynamic response functions, as well as the tuning to face type, of regions within the face processing system. Core face processing areas (occipital and fusiform face areas) responded to all types of faces with only small differences in amplitude and duration. In contrast, most areas of the extended face processing system (medial orbito-frontal, anterior and posterior cingulate) had weak responses to unfamiliar and visually-familiar faces, but were highly tuned and exhibited prolonged responses to personally-familiar faces. This indicates that the neural processing of different types of familiar faces not only differs in degree, but is probably mediated by qualitatively distinct mechanisms.