The effect of model order selection in group PICA

Network functional connectivity and whole-brain functional connectomics to investigate cognitive decline in neurodegenerative conditions

Non-invasive mapping of brain functional connectivity (FC) has played a fundamental role in neuroscience, and numerous scientists have been fascinated by its ability to reveal the brain's intricate morphology and functional properties. In recent years, two different techniques have been developed that are able to explore FC in pathophysiological conditions and to provide simple and non-invasive biomarkers for the detection of disease onset, severity and progression. These techniques are independent component analysis, which allows a network-based functional exploration of the brain, and graph theory, which provides a quantitative characterization of the whole-brain FC. In this paper we provide an overview of these two techniques and some examples of their clinical applications in the most common neurodegenerative disorders associated with cognitive decline, including mild cognitive impairment, Alzheimer's disease, Parkinson's disease, dementia with Lewy Bodies and behavioral variant frontotemporal dementia.

On the detection of high frequency correlations in resting state fMRI

Current studies of resting-state connectivity rely on coherent signal fluctuations at frequencies below 0.1 Hz, however, recent studies using high-speed fMRI have shown that fluctuations above 0.5 Hz may exist. This study replicates the feasibility of measuring high frequency (HF) correlations in six healthy controls and a patient with a brain tumor while analyzing non-physiological signal sources via simulation. Resting-state data were acquired using a high-speed multi-slab echo-volumar imaging pulse sequence with 136 ms temporal resolution. Bandpass frequency filtering in combination with sliding window seed-based connectivity analysis using running mean of the correlation maps was employed to map HF correlations up to 3.7 Hz. Computer simulations of Rician noise and the underlying point spread function were analyzed to estimate baseline spatial autocorrelation levels in four major networks (auditory, sensorimotor, visual, and default-mode). Using seed regions based on Brodmann areas, the auditory and default-mode networks were observed to have significant frequency band dependent HF correlations above baseline spatial autocorrelation levels. Correlations in the sensorimotor network were at trend level. The auditory network was still observed using a unilateral single voxel seed. In the patient, HF auditory correlations showed a spatial displacement near the tumor consistent with the displacement seen at low frequencies. In conclusion, our data suggest that HF connectivity in the human brain may be observable with high-speed fMRI, however, the detection sensitivity may depend on the network observed, data acquisition technique, and analysis method.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

Regional and source-based patterns of [11C]-(+)-PHNO binding potential reveal concurrent alterations in dopamine D2 and D3 receptor availability in cocaine-use disorder

Dopamine type 2 and type 3 receptors (D₂R/D₃R) appear critical to addictive disorders. Cocaine-use disorder (CUD) is associated with lower D₂R availability and greater D₃R availability in regions primarily expressing D₂R or D₃R concentrations, respectively. However, these CUD-related alterations in D₂R and D₃R have not been concurrently detected using available dopaminergic radioligands. Furthermore, receptor availability in regions of mixed D₂R/D₃R concentration in CUD remains unclear. The current study aimed to extend investigations of CUD-related alterations in D₂R and D₃R availability using regional and source-based analyses of [¹¹C]-(+)-PHNO positron emission tomography (PET) of 26 individuals with CUD and 26 matched healthy comparison (HC) participants. Regional analysis detected greater binding potential (BP_ND) in CUD in the midbrain, consistent with prior [¹¹C]-(+)-PHNO research, and lower BP_ND in CUD in the dorsal striatum, consistent with research using non-selective D₂R/D₃R radiotracers. Exploratory independent component analysis (ICA) identified three sources of BP_ND (striatopallidal, pallidonigral, and mesoaccumbens sources) that represent systems of brain regions displaying coherent variation in receptor availability. The striatopallidal source was associated with estimates of regional D₂R-related proportions of BP_ND (calculated using independent reports of [¹¹C]-(+)-PHNO receptor binding fractions), was lower in intensity in CUD and negatively associated with years of cocaine use. By comparison, the pallidonigral source was associated with estimates of regional D₃R distribution, was greater in intensity in CUD and positively associated with years of cocaine use. The current study extends previous D₂R/D₃R research in CUD, demonstrating both lower BP_ND in the D₂R-rich dorsal striatum and greater BP_ND in the D₃R-rich midbrain using a single radiotracer. In addition, exploratory ICA identified sources of [¹¹C]-(+)-PHNO BP_ND that were correlated with regional estimates of D₂R-related and D₃R-related proportions of BP_ND, were consistent with regional differences in CUD, and suggest receptor alterations in CUD may also be present in regions of mixed D₂R/D₃R concentration.

Intra- and inter-network functional alterations in Parkinson's disease with mild cognitive impairment

Mild cognitive impairment (MCI) is prevalent in 15%-40% of Parkinson's disease (PD) patients at diagnosis. In this investigation, we study brain intra- and inter-network alterations in resting state functional magnetic resonance imaging (rs-fMRI) in recently diagnosed PD patients and characterise them as either cognitive normal (PD-NC) or with MCI (PD-MCI). Patients were divided into two groups, PD-NC (N = 62) and PD-MCI (N = 37) and for comparison, healthy controls (HC, N = 30) were also included. Intra- and inter-network connectivity were investigated from participants' rs-fMRIs in 26 resting state networks (RSNs). Intra-network differences were found between both patient groups and HCs for networks associated with motor control (motor cortex), spatial attention and visual perception. When comparing both PD-NC and PD-MCI, intra-network alterations were found in RSNs related to attention, executive function and motor control (cerebellum). The inter-network analysis revealed a hyper-synchronisation between the basal ganglia network and the motor cortex in PD-NC compared with HCs. When both patient groups were compared, intra-network alterations in RSNs related to attention, motor control, visual perception and executive function were found. We also detected disease-driven negative synchronisations and synchronisation shifts from positive to negative and vice versa in both patient groups compared with HCs. The hyper-synchronisation between basal ganglia and motor cortical RSNs in PD and its synchronisation shift from negative to positive compared with HCs, suggest a compensatory response to basal dysfunction and altered basal-cortical motor control in the resting state brain of PD patients. Hum Brain Mapp 38:1702-1715, 2017. © 2016 Wiley Periodicals, Inc.

Common and distinct brain networks underlying verbal and visual creativity

Creativity is imperative to the progression of human civilization, prosperity, and well-being. Past creative researches tends to emphasize the default mode network (DMN) or the frontoparietal network (FPN) somewhat exclusively. However, little is known about how these networks interact to contribute to creativity and whether common or distinct brain networks are responsible for visual and verbal creativity. Here, we use functional connectivity analysis of resting-state functional magnetic resonance imaging data to investigate visual and verbal creativity-related regions and networks in 282 healthy subjects. We found that functional connectivity within the bilateral superior parietal cortex of the FPN was negatively associated with visual and verbal creativity. The strength of connectivity between the DMN and FPN was positively related to both creative domains. Visual creativity was negatively correlated with functional connectivity within the precuneus of the pDMN and right middle frontal gyrus of the FPN, and verbal creativity was negatively correlated with functional connectivity within the medial prefrontal cortex of the aDMN. Critically, the FPN mediated the relationship between the aDMN and verbal creativity, and it also mediated the relationship between the pDMN and visual creativity. Taken together, decreased within-network connectivity of the FPN and DMN may allow for flexible between-network coupling in the highly creative brain. These findings provide indirect evidence for the cooperative role of the default and executive control networks in creativity, extending past research by revealing common and distinct brain systems underlying verbal and visual creative cognition. Hum Brain Mapp 38:2094-2111, 2017. © 2017 Wiley Periodicals, Inc.

The Effect of Gray Matter ICA and Coefficient of Variation Mapping of BOLD Data on the Detection of Functional Connectivity Changes in Alzheimer's Disease and bvFTD

Resting-state fMRI results in neurodegenerative diseases have been somewhat conflicting. This may be due to complex partial volume effects of CSF in BOLD signal in patients with brain atrophy. To encounter this problem, we used a coefficient of variation (CV) map to highlight artifacts in the data, followed by analysis of gray matter voxels in order to minimize brain volume effects between groups. The effects of these measures were compared to whole brain ICA dual regression results in Alzheimer’s disease (AD) and behavioral variant frontotemporal dementia (bvFTD). 23 AD patients, 21 bvFTD patients and 25 healthy controls were included. The quality of the data was controlled by CV mapping. For detecting functional connectivity (FC) differences whole brain ICA (wbICA) and also segmented gray matter ICA (gmICA) followed by dual regression were conducted, both of which were performed both before and after data quality control. Decreased FC was detected in posterior DMN in the AD group and in the Salience network in the bvFTD group after combining CV quality control with gmICA. Before CV quality control, the decreased connectivity finding was not detectable in gmICA in neither of the groups. Same finding recurred when exclusion was based on randomization. The subjects excluded due to artifacts noticed in the CV maps had significantly lower temporal signal-to-noise ratio than the included subjects. Data quality measure CV is an effective tool in detecting artifacts from resting state analysis. CV reflects temporal dispersion of the BOLD signal stability and may thus be most helpful for spatial ICA, which has a blind spot in spatially correlating widespread artifacts. CV mapping in conjunction with gmICA yields results suiting previous findings both in AD and bvFTD.

Aberrant neural networks for the recognition memory of socially relevant information in patients with schizophrenia

INTRODUCTION

Patients with schizophrenia exhibit several cognitive deficits, including memory impairment. Problems with recognition memory can hinder socially adaptive behavior. Previous investigations have suggested that altered activation of the frontotemporal area plays an important role in recognition memory impairment. However, the cerebral networks related to these deficits are not known. The aim of this study was to elucidate the brain networks required for recognizing socially relevant information in patients with schizophrenia performing an old-new recognition task.

METHODS

Sixteen patients with schizophrenia and 16 controls participated in this study. First, the subjects performed the theme-identification task during functional magnetic resonance imaging. In this task, pictures depicting social situations were presented with three words, and the subjects were asked to select the best theme word for each picture. The subjects then performed an old-new recognition task in which they were asked to discriminate whether the presented words were old or new. Task performance and neural responses in the old-new recognition task were compared between the subject groups. An independent component analysis of the functional connectivity was performed.

RESULTS

The patients with schizophrenia exhibited decreased discriminability and increased activation of the right superior temporal gyrus compared with the controls during correct responses. Furthermore, aberrant network activities were found in the frontopolar and language comprehension networks in the patients.

CONCLUSIONS

The functional connectivity analysis showed aberrant connectivity in the frontopolar and language comprehension networks in the patients with schizophrenia, and these aberrations possibly contribute to their low recognition performance and social dysfunction. These results suggest that the frontopolar and language comprehension networks are potential therapeutic targets in patients with schizophrenia.

Segregation between the parietal memory network and the default mode network: effects of spatial smoothing and model order in ICA

A brain network consisting of two key parietal nodes, the precuneus and the posterior cingulate cortex, has emerged from recent fMRI studies. Though it is anatomically adjacent to and spatially overlaps with the default mode network (DMN), its function has been associated with memory processing, and it has been referred to as the parietal memory network (PMN). Independent component analysis (ICA) is the most common data-driven method used to extract PMN and DMN simultaneously. However, the effects of data preprocessing and parameter determination in ICA on PMN–DMN segregation are completely unknown. Here, we employ three typical algorithms of group ICA to assess how spatial smoothing and model order influence the degree of PMN–DMN segregation. Our findings indicate that PMN and DMN can only be stably separated using a combination of low-level spatial smoothing and high model order across the three ICA algorithms. We thus argue for more considerations on parametric settings for interpreting DMN data.

Chronnectomic patterns and neural flexibility underlie executive function

Despite extensive research into executive function (EF), the precise relationship between brain dynamics and flexible cognition remains unknown. Using a large, publicly available dataset (189 participants), we find that functional connections measured throughout 56min of resting state fMRI data comprise five distinct connectivity states. Elevated EF performance as measured outside of the scanner was associated with greater episodes of more frequently occurring connectivity states, and fewer episodes of less frequently occurring connectivity states. Frequently occurring states displayed metastable properties, where cognitive flexibility may be facilitated by attenuated correlations and greater functional connection variability. Less frequently occurring states displayed properties consistent with low arousal and low vigilance. These findings suggest that elevated EF performance may be associated with the propensity to occupy more frequently occurring brain configurations that enable cognitive flexibility, while avoiding less frequently occurring brain configurations related to low arousal/vigilance states. The current findings offer a novel framework for identifying neural processes related to individual differences in executive function.

Multimodal Classification of Schizophrenia Patients with MEG and fMRI Data Using Static and Dynamic Connectivity Measures

Mental disorders like schizophrenia are currently diagnosed by physicians/psychiatrists through clinical assessment and their evaluation of patient's self-reported experiences as the illness emerges. There is great interest in identifying biological markers of prognosis at the onset of illness, rather than relying on the evolution of symptoms across time. Functional network connectivity, which indicates a subject's overall level of "synchronicity" of activity between brain regions, demonstrates promise in providing individual subject predictive power. Many previous studies reported functional connectivity changes during resting-state using only functional magnetic resonance imaging (fMRI). Nevertheless, exclusive reliance on fMRI to generate such networks may limit the inference of the underlying dysfunctional connectivity, which is hypothesized to be a factor in patient symptoms, as fMRI measures connectivity via hemodynamics. Therefore, combination of connectivity assessments using fMRI and magnetoencephalography (MEG), which more directly measures neuronal activity, may provide improved classification of schizophrenia than either modality alone. Moreover, recent evidence indicates that metrics of dynamic connectivity may also be critical for understanding pathology in schizophrenia. In this work, we propose a new framework for extraction of important disease related features and classification of patients with schizophrenia based on using both fMRI and MEG to investigate functional network components in the resting state. Results of this study show that the integration of fMRI and MEG provides important information that captures fundamental characteristics of functional network connectivity in schizophrenia and is helpful for prediction of schizophrenia patient group membership. Combined fMRI/MEG methods, using static functional network connectivity analyses, improved classification accuracy relative to use of fMRI or MEG methods alone (by 15 and 12.45%, respectively), while combined fMRI/MEG methods using dynamic functional network connectivity analyses improved classification up to 5.12% relative to use of fMRI alone and up to 17.21% relative to use of MEG alone.

PreSurgMapp: a MATLAB Toolbox for Presurgical Mapping of Eloquent Functional Areas Based on Task-Related and Resting-State Functional MRI

The main goal of brain tumor surgery is to maximize tumor resection while minimizing the risk of irreversible postoperative functional sequelae. Eloquent functional areas should be delineated preoperatively, particularly for patients with tumors near eloquent areas. Functional magnetic resonance imaging (fMRI) is a noninvasive technique that demonstrates great promise for presurgical planning. However, specialized data processing toolkits for presurgical planning remain lacking. Based on several functions in open-source software such as Statistical Parametric Mapping (SPM), Resting-State fMRI Data Analysis Toolkit (REST), Data Processing Assistant for Resting-State fMRI (DPARSF) and Multiple Independent Component Analysis (MICA), here, we introduce an open-source MATLAB toolbox named PreSurgMapp. This toolbox can reveal eloquent areas using comprehensive methods and various complementary fMRI modalities. For example, PreSurgMapp supports both model-based (general linear model, GLM, and seed correlation) and data-driven (independent component analysis, ICA) methods and processes both task-based and resting-state fMRI data. PreSurgMapp is designed for highly automatic and individualized functional mapping with a user-friendly graphical user interface (GUI) for time-saving pipeline processing. For example, sensorimotor and language-related components can be automatically identified without human input interference using an effective, accurate component identification algorithm using discriminability index. All the results generated can be further evaluated and compared by neuro-radiologists or neurosurgeons. This software has substantial value for clinical neuro-radiology and neuro-oncology, including application to patients with low- and high-grade brain tumors and those with epilepsy foci in the dominant language hemisphere who are planning to undergo a temporal lobectomy.

Building an EEG-fMRI Multi-Modal Brain Graph: A Concurrent EEG-fMRI Study

The topological architecture of brain connectivity has been well-characterized by graph theory based analysis. However, previous studies have primarily built brain graphs based on a single modality of brain imaging data. Here we develop a framework to construct multi-modal brain graphs using concurrent EEG-fMRI data which are simultaneously collected during eyes open (EO) and eyes closed (EC) resting states. FMRI data are decomposed into independent components with associated time courses by group independent component analysis (ICA). EEG time series are segmented, and then spectral power time courses are computed and averaged within 5 frequency bands (delta; theta; alpha; beta; low gamma). EEG-fMRI brain graphs, with EEG electrodes and fMRI brain components serving as nodes, are built by computing correlations within and between fMRI ICA time courses and EEG spectral power time courses. Dynamic EEG-fMRI graphs are built using a sliding window method, versus static ones treating the entire time course as stationary. In global level, static graph measures and properties of dynamic graph measures are different across frequency bands and are mainly showing higher values in eyes closed than eyes open. Nodal level graph measures of a few brain components are also showing higher values during eyes closed in specific frequency bands. Overall, these findings incorporate fMRI spatial localization and EEG frequency information which could not be obtained by examining only one modality. This work provides a new approach to examine EEG-fMRI associations within a graph theoretic framework with potential application to many topics.

Large-scale functional network overlap is a general property of brain functional organization: Reconciling inconsistent fMRI findings from general-linear-model-based analyses

Functional magnetic resonance imaging (fMRI) studies regularly use univariate general-linear-model-based analyses (GLM). Their findings are often inconsistent across different studies, perhaps because of several fundamental brain properties including functional heterogeneity, balanced excitation and inhibition (E/I), and sparseness of neuronal activities. These properties stipulate heterogeneous neuronal activities in the same voxels and likely limit the sensitivity and specificity of GLM. This paper selectively reviews findings of histological and electrophysiological studies and fMRI spatial independent component analysis (sICA) and reports new findings by applying sICA to two existing datasets. The extant and new findings consistently demonstrate several novel features of brain functional organization not revealed by GLM. They include overlap of large-scale functional networks (FNs) and their concurrent opposite modulations, and no significant modulations in activity of most FNs across the whole brain during any task conditions. These novel features of brain functional organization are highly consistent with the brain's properties of functional heterogeneity, balanced E/I, and sparseness of neuronal activity, and may help reconcile inconsistent GLM findings.

Amplitude variability over trials in hemodynamic responses in adolescents with ADHD: The role of the anterior default mode network and the non-specific role of the striatum

It has been suggested that intra-individual variability (IIV) in performance on attention and other cognitive tasks might be a cognitive endophenotype in individuals with ADHD. Despite robust IIV findings in behavioral data, only sparse data exist on how what type of brain dysfunction underlies variable response times. In this study, we asked whether ADHD IIV in reaction time on a commonly-used test of attention might be related to variation in hemodynamic responses (HRs) observed trial-to-trial. Based on previous studies linking IIV to regions within the “default mode” network (DMN), we predicted that adolescents with ADHD would have higher HR variability in the DMN compared with controls, and this in turn would be related to behavioral IIV. We also explored the influence of social anxiety on HR variability in ADHD as means to test whether higher arousal associated with high trait anxiety would affect the neural abnormalities. We assessed single-trial variability of HRs, estimated from fMRI event-related responses elicited during an auditory oddball paradigm in adolescents with ADHD and healthy controls (11–18 years old; N = 46). Adolescents with ADHD had higher HR variability compared with controls in anterior regions of the DMN. This effect was specific to ADHD and not associated with traits of age, IQ and anxiety. However, an ADHD effect of higher HR variability also appeared in a basal ganglia network, but for these brain regions the relationships of HR variability and social anxiety levels were more complex. Performance IIV correlated significantly with variability of HRs in both networks. These results suggest that assessment of trial-to-trial HR variability in ADHD provides information beyond that detectable through analysis of behavioral data and average brain activation levels, revealing specific neural correlates of a possible ADHD IIV endophenotype.

•

We studied if the behavioral variability in ADHD is also found on a neuronal level.

•

Independent component analysis was combined with BOLD amplitude variability.

•

Adolescents with ADHD had higher amplitude variability than healthy controls.

•

Higher amplitude variability was shown in an anterior default mode network.

•

Social anxiety in ADHD associated with high amplitude variability in the striatum

Dynamic interactions of the cortical networks during thought suppression

OBJECTIVES

Thought suppression has spurred extensive research in clinical and preclinical fields, particularly with regard to the paradoxical aspects of this behavior. However, the involvement of the brain's inhibitory system in the dynamics underlying the continuous effort to suppress thoughts has yet to be clarified. This study aims to provide a unified perspective for the volitional suppression of internal events incorporating the current understanding of the brain's inhibitory system.

MATERIALS AND METHODS

Twenty healthy volunteers underwent functional magnetic resonance imaging while they performed thought suppression blocks alternating with visual imagery blocks. The whole dataset was decomposed by group-independent component analysis into 30 components. After discarding noise components, the 20 valid components were subjected to further analysis of their temporal properties including task-relatedness and between-component residual correlation.

RESULTS

Combining a long task period and a data-driven approach, we observed a right-side-dominant, lateral frontoparietal network to be strongly suppression related. This network exhibited increased fluctuation during suppression, which is compatible with the well-known difficulty of suppression maintenance.

CONCLUSIONS

Between-network correlation provided further insight into the coordinated engagement of the executive control and dorsal attention networks, as well as the reciprocal activation of imagery-related components, thus revealing neural substrates associated with the rivalry between intrusive thoughts and the suppression process.

Is the Volume of the Caudate Nuclei Associated With Area of Secondary Hyperalgesia? - Protocol for a 3-Tesla MRI Study of Healthy Volunteers

Background

Experience and development of pain may be influenced by a number of physiological, psychological, and psychosocial factors. In a previous study we found differences in neuronal activation to noxious stimulation, and microstructural neuroanatomical differences, when comparing healthy volunteers with differences in size of the area of secondary hyperalgesia following a standardized burn injury.

Objective

We aim to investigate the degree of association between the volume of pain-relevant structures in the brain and the size of the area of secondary hyperalgesia following brief thermal sensitization.

Methods

The study consists of one experimental day, in which whole-brain magnetic resonance imaging (MRI) scans will be conducted including T1-weighed three-dimensional anatomy scan, diffusion tensor imaging, and resting state functional MRI. Before the experimental day, all included participants will undergo experimental pain testing in a parallel study (Clinicaltrials.gov Identifier: NCT02527395). Results from this experimental pain testing, as well as the size of the area of secondary hyperalgesia from the included participants, will be extracted from this parallel study.

Results

The association between the volume of pain-relevant structures in the brain and the area of secondary hyperalgesia will be investigated by linear regression of the estimated best linear unbiased predictors on the individual volumes of the pain relevant brain structures.

Conclusions

We plan to investigate the association between experimental pain testing parameters and the volume, connectivity, and resting state activity of pain-relevant structures in the brain. These results may improve our knowledge of the mechanisms responsible for the development of acute and chronic pain.

ClinicalTrial

Danish Research Ethics Committee (identifier: H-15010473). Danish Data Protection Agency (identifier: RH-2015-149). Clinicaltrials.gov NCT02567318; http://clinicaltrials.gov/ct2/show/NCT02567318 (Archived by WebCite at http://www.webcitation.org/6i4OtP0Oi)

Aberrant salience network and its functional coupling with default and executive networks in minimal hepatic encephalopathy: a resting-state fMRI study

The purposes of this study are to explore functional alterations in salience network (SN) and its functional coupling with default mode (DMN) and central executive (CEN) networks in minimal hepatic encephalopathy (MHE). Twenty cirrhotic patients with MHE, 23 cirrhotic patients without MHE (NHE), and 18 controls underwent resting-state fMRI and psychometric hepatic encephalopathy score (PHES) test. Independent component analysis was performed to obtain DMN (including three subsystems: anterior, inferior-posterior, and superior-posterior DMN [a/ip/spDMN]), SN, and CEN (including three subsystems: left-ventral, right-ventral, and dorsal CEN [lv/rv/dCEN]). The intrinsic functional connectivity (iFC) within (intra-iFC) and between (inter-iFC and time-lagged inter-iFC) networks was measured. MHE patients showed decreased intra-iFC within aDMN, SN, lvCEN, and rvCEN; and decreased inter-iFC and time-lagged inter-iFC between SN and ipDMN/spDMN/lvCEN and increased inter-iFC and time-lagged inter-iFC between SN and aDMN, compared with controls. A progressive trend in connectivity alterations was found as the disease developed from NHE to MHE. The inter-iFC between ipDMN/spDMN and SN was significantly correlated with PHES score. In conclusion, an aberrant SN and its functional interaction with the DMN/CEN are core features of MHE that are associated with disease progression and may play an important role in neurocognitive dysfunction in MHE.

A Preliminary Prospective Study of an Escalation in 'Maximum Daily Drinks', Fronto-Parietal Circuitry and Impulsivity-Related Domains in Young Adult Drinkers

Excessive alcohol use in young adults is associated with greater impulsivity and neurobiological alterations in executive control systems. The maximum number of drinks consumed during drinking occasions ('MaxDrinks') represents a phenotype linked to vulnerability of alcohol use disorders, and an increase, or 'escalation', in MaxDrinks may be indicative of greater risk for problematic drinking. Thirty-six young adult drinkers performed a Go/No-Go task during fMRI, completed impulsivity-related assessments, and provided monthly reports of alcohol use during a 12-month follow-up period. Participants were characterized by MaxDrinks at baseline and after follow-up, identifying 18 escalating drinkers and 18 constant drinkers. Independent component analysis was used to investigate functional brain networks associated with response inhibition, and relationships with principal component analysis derived impulsivity-related domains were examined. Greater baseline MaxDrinks was associated with an average reduction in the engagement of a right-lateralized fronto-parietal functional network, while an escalation in MaxDrinks was associated with a greater difference in fronto-parietal engagement between successful inhibitions and error trials. Escalating drinkers displayed greater impulsivity/compulsivity-related domain scores that were positively associated with fronto-parietal network engagement and change in MaxDrinks during follow-up. In young adults, an escalating MaxDrinks trajectory was prospectively associated with altered fronto-parietal control mechanisms and greater impulsivity/compulsivity scores. Continued longitudinal studies of MaxDrinks trajectories, functional network activity, and impulsivity/compulsivity-related features may lend further insight into an intermediate phenotype vulnerable for alcohol use and addictive disorders.

Classification of schizophrenia and bipolar patients using static and dynamic resting-state fMRI brain connectivity

Recently, functional network connectivity (FNC, defined as the temporal correlation among spatially distant brain networks) has been used to examine the functional organization of brain networks in various psychiatric illnesses. Dynamic FNC is a recent extension of the conventional FNC analysis that takes into account FNC changes over short periods of time. While such dynamic FNC measures may be more informative about various aspects of connectivity, there has been no detailed head-to-head comparison of the ability of static and dynamic FNC to perform classification in complex mental illnesses. This paper proposes a framework for automatic classification of schizophrenia, bipolar and healthy subjects based on their static and dynamic FNC features. Also, we compare cross-validated classification performance between static and dynamic FNC. Results show that the dynamic FNC significantly outperforms the static FNC in terms of predictive accuracy, indicating that features from dynamic FNC have distinct advantages over static FNC for classification purposes. Moreover, combining static and dynamic FNC features does not significantly improve the classification performance over the dynamic FNC features alone, suggesting that static FNC does not add any significant information when combined with dynamic FNC for classification purposes. A three-way classification methodology based on static and dynamic FNC features discriminates individual subjects into appropriate diagnostic groups with high accuracy. Our proposed classification framework is potentially applicable to additional mental disorders.

Test-retest reliability of the default mode network in a multi-centric fMRI study of healthy elderly: Effects of data-driven physiological noise correction techniques

Understanding how to reduce the influence of physiological noise in resting state fMRI data is important for the interpretation of functional brain connectivity. Limited data is currently available to assess the performance of physiological noise correction techniques, in particular when evaluating longitudinal changes in the default mode network (DMN) of healthy elderly participants. In this 3T harmonized multisite fMRI study, we investigated how different retrospective physiological noise correction (rPNC) methods influence the within-site test-retest reliability and the across-site reproducibility consistency of DMN-derived measurements across 13 MRI sites. Elderly participants were scanned twice at least a week apart (five participants per site). The rPNC methods were: none (NPC), Tissue-based regression, PESTICA and FSL-FIX. The DMN at the single subject level was robustly identified using ICA methods in all rPNC conditions. The methods significantly affected the mean z-scores and, albeit less markedly, the cluster-size in the DMN; in particular, FSL-FIX tended to increase the DMN z-scores compared to others. Within-site test-retest reliability was consistent across sites, with no differences across rPNC methods. The absolute percent errors were in the range of 5-11% for DMN z-scores and cluster-size reliability. DMN pattern overlap was in the range 60-65%. In particular, no rPNC method showed a significant reliability improvement relative to NPC. However, FSL-FIX and Tissue-based physiological correction methods showed both similar and significant improvements of reproducibility consistency across the consortium (ICC = 0.67) for the DMN z-scores relative to NPC. Overall these findings support the use of rPNC methods like tissue-based or FSL-FIX to characterize multisite longitudinal changes of intrinsic functional connectivity. Hum Brain Mapp 37:2114-2132, 2016. © 2016 Wiley Periodicals, Inc.

Dynamic functional network connectivity reveals unique and overlapping profiles of insula subdivisions

The human insular cortex consists of functionally diverse subdivisions that engage during tasks ranging from interoception to cognitive control. The multiplicity of functions subserved by insular subdivisions calls for a nuanced investigation of their functional connectivity profiles. Four insula subdivisions (dorsal anterior, dAI; ventral, VI; posterior, PI; middle, MI) derived using a data-driven approach were subjected to static- and dynamic functional network connectivity (s-FNC and d-FNC) analyses. Static-FNC analyses replicated previous work demonstrating a cognition-emotion-interoception division of the insula, where the dAI is functionally connected to frontal areas, the VI to limbic areas, and the PI and MI to sensorimotor areas. Dynamic-FNC analyses consisted of k-means clustering of sliding windows to identify variable insula connectivity states. The d-FNC analysis revealed that the most frequently occurring dynamic state mirrored the cognition-emotion-interoception division observed from the s-FNC analysis, with less frequently occurring states showing overlapping and unique subdivision connectivity profiles. In two of the states, all subdivisions exhibited largely overlapping profiles, consisting of subcortical, sensory, motor, and frontal connections. Two other states showed the dAI exhibited a unique connectivity profile compared with other insula subdivisions. Additionally, the dAI exhibited the most variable functional connections across the s-FNC and d-FNC analyses, and was the only subdivision to exhibit dynamic functional connections with regions of the default mode network. These results highlight how a d-FNC approach can capture functional dynamics masked by s-FNC approaches, and reveal dynamic functional connections enabling the functional flexibility of the insula across time. Hum Brain Mapp 37:1770-1787, 2016. © 2016 Wiley Periodicals, Inc.

Analysis of Residual Dependencies of Independent Components Extracted from fMRI Data

Independent component analysis (ICA) of functional magnetic resonance imaging (fMRI) data can be employed as an exploratory method. The lack in the ICA model of strong a priori assumptions about the signal or about the noise leads to difficult interpretations of the results. Moreover, the statistical independence of the components is only approximated. Residual dependencies among the components can reveal informative structure in the data. A major problem is related to model order selection, that is, the number of components to be extracted. Specifically, overestimation may lead to component splitting. In this work, a method based on hierarchical clustering of ICA applied to fMRI datasets is investigated. The clustering algorithm uses a metric based on the mutual information between the ICs. To estimate the similarity measure, a histogram-based technique and one based on kernel density estimation are tested on simulated datasets. Simulations results indicate that the method could be used to cluster components related to the same task and resulting from a splitting process occurring at different model orders. Different performances of the similarity measures were found and discussed. Preliminary results on real data are reported and show that the method can group task related and transiently task related components.

Modern Methods for Interrogating the Human Connectome

OBJECTIVES

Connectionist theories of brain function took hold with the seminal contributions of Norman Geschwind a half century ago. Modern neuroimaging techniques have expanded the scientific interest in the study of brain connectivity to include the intact as well as disordered brain.

METHODS

In this review, we describe the most common techniques used to measure functional and structural connectivity, including resting state functional MRI, diffusion MRI, and electroencephalography and magnetoencephalography coherence. We also review the most common analytical approaches used for examining brain interconnectivity associated with these various imaging methods.

RESULTS

This review presents a critical analysis of the assumptions, as well as methodological limitations, of each imaging and analysis approach.

CONCLUSIONS

The overall goal of this review is to provide the reader with an introduction to evaluating the scientific methods underlying investigations that probe the human connectome.

Resting-state functional network connectivity in prefrontal regions differs between unmedicated patients with bipolar and major depressive disorders

BACKGROUND

Differentiating bipolar disorder (BD) from major depressive disorder (MDD) often poses a major clinical challenge, and optimal clinical care can be hindered by misdiagnoses. This study investigated the differences between BD and MDD in resting-state functional network connectivity (FNC) using a data-driven image analysis method.

METHODS

In this study, fMRI data were collected from unmedicated subjects including 13 BD, 40 MDD and 33 healthy controls (HC). The FNC was calculated between functional brain networks derived from fMRI using group independent component analysis (ICA). Group comparisons were performed on connectivity strengths and other graph measures of FNC matrices.

RESULTS

Statistical tests showed that, compared to MDD, the FNC in BD was characterized by more closely connected and more efficient topological structures as assessed by graph theory. The differences were found at both the whole-brain-level and the functional-network-level in prefrontal networks located in the dorsolateral/ventrolateral prefrontal cortex (DLPFC, VLPFC) and anterior cingulate cortex (ACC). Furthermore, interconnected structures in these networks in both patient groups were negatively associated with symptom severity on depression rating scales.

LIMITATIONS

As patients were unmedicated, the sample sizes were relatively small, although they were comparable to those in previous fMRI studies comparing BD and MDD.

CONCLUSIONS

Our results suggest that the differences in FNC of the PFC reflect distinct pathophysiological mechanisms in BD and MDD. Such findings ultimately may elucidate the neural pathways in which distinct functional changes can give rise to the clinical differences observed between these syndromes.

Metabolic spatial connectivity in amyotrophic lateral sclerosis as revealed by independent component analysis

OBJECTIVES

Positron emission tomography (PET) and volume of interest (VOI) analysis have recently shown in amyotrophic lateral sclerosis (ALS) an accuracy of 93% in differentiating patients from controls. The aim of this study was to disclose by spatial independent component analysis (ICA) the brain networks involved in ALS pathological processes and evaluate their discriminative value in separating patients from controls.

EXPERIMENTAL DESIGN

Two hundred fifty-nine ALS patients and 40 age- and sex-matched control subjects underwent brain 18F-2-fluoro-2-deoxy-D-glucose PET (FDG-PET). Spatial ICA of the preprocessed FDG-PET images was performed. Intensity values were converted to z-scores and binary masks were used as data-driven VOIs. The accuracy of this classifier was tested versus a validated system processing intensity signals in 27 brain meta-VOIs. A support vector machine was independently applied to both datasets and the 'leave-one-out' technique verified the general validity of results.

PRINCIPAL OBSERVATIONS

The 8 components selected as pathophysiologically meaningful discriminated patients from controls with 99.0% accuracy, the discriminating value of bilateral cerebellum/midbrain alone representing 96.3%. Among the meta-VOIs, right temporal lobe alone reached an accuracy of 93.7%.

CONCLUSIONS

Spatial ICA identified in a very large cohort of ALS patients distinct spatial networks showing a high discriminatory value, improving substantially on the previously obtained accuracy. The cerebellar/midbrain component accounted for the highest accuracy in separating ALS patients from controls. Spatial ICA and multivariate analysis perform better than univariate semi-quantification methods in identifying the neurodegenerative features of ALS and pave the way for inclusion of PET in clinical trials and early diagnosis.

Artifact removal in the context of group ICA: A comparison of single-subject and group approaches

Independent component analysis (ICA) has been widely applied to identify intrinsic brain networks from fMRI data. Group ICA computes group-level components from all data and subsequently estimates individual-level components to recapture intersubject variability. However, the best approach to handle artifacts, which may vary widely among subjects, is not yet clear. In this work, we study and compare two ICA approaches for artifacts removal. One approach, recommended in recent work by the Human Connectome Project, first performs ICA on individual subject data to remove artifacts, and then applies a group ICA on the cleaned data from all subjects. We refer to this approach as Individual ICA based artifacts Removal Plus Group ICA (IRPG). A second proposed approach, called Group Information Guided ICA (GIG-ICA), performs ICA on group data, then removes the group-level artifact components, and finally performs subject-specific ICAs using the group-level non-artifact components as spatial references. We used simulations to evaluate the two approaches with respect to the effects of data quality, data quantity, variable number of sources among subjects, and spatially unique artifacts. Resting-state test-retest datasets were also employed to investigate the reliability of functional networks. Results from simulations demonstrate GIG-ICA has greater performance compared with IRPG, even in the case when single-subject artifacts removal is perfect and when individual subjects have spatially unique artifacts. Experiments using test-retest data suggest that GIG-ICA provides more reliable functional networks. Based on high estimation accuracy, ease of implementation, and high reliability of functional networks, we find GIG-ICA to be a promising approach.

The nature of the language input affects brain activation during learning from a natural language

Artificial language studies have demonstrated that learners are able to segment individual word-like units from running speech using the transitional probability information. However, this skill has rarely been examined in the context of natural languages, where stimulus parameters can be quite different. In this study, two groups of English-speaking learners were exposed to Norwegian sentences over the course of three fMRI scans. One group was provided with input in which transitional probabilities predicted the presence of target words in the sentences. This group quickly learned to identify the target words and fMRI data revealed an extensive and highly dynamic learning network. These results were markedly different from activation seen for a second group of participants. This group was provided with highly similar input that was modified so that word learning based on syllable co-occurrences was not possible. These participants showed a much more restricted network. The results demonstrate that the nature of the input strongly influenced the nature of the network that learners employ to learn the properties of words in a natural language.

Connectivity-based whole brain dual parcellation by group ICA reveals tract structures and decreased connectivity in schizophrenia

Mapping brain connectivity based on neuroimaging data is a promising new tool for understanding brain structure and function. In this methods paper, we demonstrate that group independent component analysis (GICA) can be used to perform a dual parcellation of the brain based on its connectivity matrix (cmICA). This dual parcellation consists of a set of spatially independent source maps, and a corresponding set of paired dual maps that define the connectivity of each source map to the brain. These dual maps are called the connectivity profiles of the source maps. Traditional analysis of connectivity matrices has been used previously for brain parcellation, but the present method provides additional information on the connectivity of these segmented regions. In this paper, the whole brain structural connectivity matrices were calculated on a 5 mm(3) voxel scale from diffusion imaging data based on the probabilistic tractography method. The effect of the choice of the number of components (30 and 100) and their stability were examined. This method generated a set of spatially independent components that are consistent with the canonical brain tracts provided by previous anatomic descriptions, with the high order model yielding finer segmentations. The corpus-callosum example shows how this method leads to a robust parcellation of a brain structure based on its connectivity properties. We applied cmICA to study structural connectivity differences between a group of schizophrenia subjects and healthy controls. The connectivity profiles at both model orders showed similar regions with reduced connectivity in schizophrenia patients. These regions included forceps major, right inferior fronto-occipital fasciculus, uncinate fasciculus, thalamic radiation, and corticospinal tract. This paper provides a novel unsupervised data-driven framework that summarizes the information in a large global connectivity matrix and tests for brain connectivity differences. It has the potential for capturing important brain changes related to disease in connectivity-based disorders.

Neural Processes in the Human Temporoparietal Cortex Separated by Localized Independent Component Analysis

The human temporoparietal junction (TPJ) is a topic of intense research. Imaging studies have identified TPJ activation in association with many higher-order functions such as theory-of-mind, episodic memory, and attention, causing debate about the distribution of different processes. One major challenge is the lack of consensus about the anatomical location and extent of the TPJ. Here, we address this problem using data-driven analysis to test the hypothesis that the bilateral TPJ can be parcellated into subregions. We applied independent component analysis (ICA) to task-free fMRI data within a local region around the bilateral TPJ, iterating the ICA at multiple model orders and in several datasets. The localized analysis allowed finer separation of processes and the use of multiple dimensionalities provided qualitative information about lateralization. We identified four subdivisions that were bilaterally symmetrical and one that was right biased. To test whether the independent components (ICs) reflected true subdivisions, we performed functional connectivity analysis using the IC coordinates as seeds. This confirmed that the subdivisions belonged to distinct networks. The right-biased IC was connected with a network often associated with attentional processing. One bilateral subdivision was connected to sensorimotor regions and another was connected to auditory regions. One subdivision that presented as distinct left- and right-biased ICs was connected to frontoparietal regions. Another subdivision that also had left- and right-biased ICs was connected to social or default mode networks. Our results show that the TPJ in both hemispheres hosts multiple neural processes with connectivity patterns consistent with well developed specialization and lateralization.

Altered resting perfusion and functional connectivity of default mode network in youth with autism spectrum disorder

BACKGROUND

Neuroimaging studies can shed light on the neurobiological underpinnings of autism spectrum disorders (ASD). Studies of the resting brain have shown both altered baseline metabolism from PET/SPECT and altered functional connectivity (FC) of intrinsic brain networks based on resting-state fMRI. To date, however, no study has investigated these two physiological parameters of resting brain function jointly, or explored the relationship between these measures and ASD symptom severity.

METHODS

Here, we used pseudo-continuous arterial spin labeling with 3D background-suppressed GRASE to assess resting cerebral blood flow (CBF) and FC in 17 youth with ASD and 22 matched typically developing (TD) children.

RESULTS

A pattern of altered resting perfusion was found in ASD versus TD children including frontotemporal hyperperfusion and hypoperfusion in the dorsal anterior cingulate cortex. We found increased local FC in the anterior module of the default mode network (DMN) accompanied by decreased CBF in the same area. In our cohort, both alterations were associated with greater social impairments as assessed with the Social Responsiveness Scale (SRS-total T scores). While FC was correlated with CBF in TD children, this association between FC and baseline perfusion was disrupted in children with ASD. Furthermore, there was reduced long-range FC between anterior and posterior modules of the DMN in children with ASD.

CONCLUSION

Taken together, the findings of this study--the first to jointly assess resting CBF and FC in ASD--highlight new avenues for identifying novel imaging markers of ASD symptomatology.

Detection of short-term activity avalanches in human brain default mode network with ultrafast MR encephalography

Recent studies pinpoint visually cued networks of avalanches with MEG/EEG data. Co-activation pattern (CAP) analysis can be used to detect single brain volume activity profiles and hemodynamic fingerprints of neuronal avalanches as sudden high signal activity peaks in classical fMRI data. In this study, we aimed to detect dynamic patterns of brain activity spreads with the use of ultrafast MR encephalography (MREG). MREG achieves 10 Hz whole brain sampling, allowing the estimation of spatial spread of an avalanche, even with the inherent hemodynamic delay of the BOLD signal. We developed a novel computational method to separate avalanche type fast activity spreads from motion artifacts, vasomotor fluctuations, and cardio-respiratory noise in human brain default mode network (DMN). Reproducible and classical DMN sources were identified using spatial ICA prior to advanced noise removal in order to assure that ICA converges to reproducible networks. Brain activity peaks were identified from parts of the DMN, and normalized MREG data around each peak were extracted individually to show dynamic avalanche type spreads as video clips within the DMN. Individual activity spread video clips of specific parts of the DMN were then averaged over the group of subjects. The experiments show that the high BOLD values around the peaks are mostly spreading along the spatial pattern of the particular DMN segment detected with ICA. With also the spread size and lifetime resembling the expected power law distributions, this indicates that the detected peaks are parts of activity avalanches, starting from (or crossing) the DMN. Furthermore, the split, one-sided sub-networks of the DMN show different spread directions within the same DMN framework. The results open possibilities to follow up brain activity avalanches in the hope to understand more about the system wide properties of diseases related to DMN dysfunction.

A group ICA based framework for evaluating resting fMRI markers when disease categories are unclear: application to schizophrenia, bipolar, and schizoaffective disorders

Schizophrenia (SZ), bipolar disorder (BP) and schizoaffective disorder (SAD) share some common symptoms, and there is still a debate about whether SAD is an independent category. To the best of our knowledge, no study has been done to differentiate these three disorders or to investigate the distinction of SAD as an independent category using fMRI data. This study is aimed to explore biomarkers from resting-state fMRI networks for differentiating these disorders and investigate the relationship among these disorders based on fMRI networks with an emphasis on SAD. Firstly, a novel group ICA method, group information guided independent component analysis (GIG-ICA), was applied to extract subject-specific brain networks from fMRI data of 20 healthy controls (HC), 20 SZ patients, 20 BP patients, 20 patients suffering from SAD with manic episodes (SADM), and 13 patients suffering from SAD with depressive episodes exclusively (SADD). Then, five-level one-way analysis of covariance and multiclass support vector machine recursive feature elimination were employed to identify discriminative regions from the networks. Subsequently, the t-distributed stochastic neighbor embedding (t-SNE) projection and the hierarchical clustering were implemented to investigate the relationship among those groups. Finally, to evaluate the generalization ability, 16 new subjects were classified based on the found regions and the trained model using original 93 subjects. Results show that the discriminative regions mainly included frontal, parietal, precuneus, cingulate, supplementary motor, cerebellar, insula and supramarginal cortices, which performed well in distinguishing different groups. SADM and SADD were the most similar to each other, although SADD had greater similarity to SZ compared to other groups, which indicates that SAD may be an independent category. BP was closer to HC compared with other psychotic disorders. In summary, resting-state fMRI brain networks extracted via GIG-ICA provide a promising potential to differentiate SZ, BP, and SAD.

Automated iterative reclustering framework for determining hierarchical functional networks in resting state fMRI

To spatially cluster resting state-functional magnetic resonance imaging (rs-fMRI) data into potential networks, there are only a few general approaches that determine the number of networks/clusters, despite a wide variety of techniques proposed for clustering. For individual subjects, extraction of a large number of spatially disjoint clusters results in multiple small networks that are spatio-temporally homogeneous but irreproducible across subjects. Alternatively, extraction of a small number of clusters creates spatially large networks that are temporally heterogeneous but spatially reproducible across subjects. We propose a fully automatic, iterative reclustering framework in which a small number of spatially large, heterogeneous networks are initially extracted to maximize spatial reproducibility. Subsequently, the large networks are iteratively subdivided to create spatially reproducible subnetworks until the overall within-network homogeneity does not increase substantially. The proposed approach discovers a rich network hierarchy in the brain while simultaneously optimizing spatial reproducibility of networks across subjects and individual network homogeneity. We also propose a novel metric to measure the connectivity of brain regions, and in a simulation study show that our connectivity metric and framework perform well in the face of low signal to noise and initial segmentation errors. Experimental results generated using real fMRI data show that the proposed metric improves stability of network clusters across subjects, and generates a meaningful pattern for spatially hierarchical structure of the brain.

Dimensionality of ICA in resting-state fMRI investigated by feature optimized classification of independent components with SVM

Different machine learning algorithms have recently been used for assisting automated classification of independent component analysis (ICA) results from resting-state fMRI data. The success of this approach relies on identification of artifact components and meaningful functional networks. A limiting factor of ICA is the uncertainty of the number of independent components (NIC). We aim to develop a framework based on support vector machines (SVM) and optimized feature-selection for automated classification of independent components (ICs) and use the framework to investigate the effects of input NIC on the ICA results. Seven different resting-state fMRI datasets were studied. 18 features were devised by mimicking the empirical criteria for manual evaluation. The five most significant (p < 0.01) features were identified by general linear modeling and used to generate a classification model for the framework. This feature-optimized classification of ICs with SVM (FOCIS) framework was used to classify both group and single subject ICA results. The classification results obtained using FOCIS and previously published FSL-FIX were compared against manually evaluated results. On average the false negative rate in identifying artifact contaminated ICs for FOCIS and FSL-FIX were 98.27 and 92.34%, respectively. The number of artifact and functional network components increased almost linearly with the input NIC. Through tracking, we demonstrate that incrementing NIC affects most ICs when NIC < 33, whereas only a few limited ICs are affected by direct splitting when NIC is incremented beyond NIC > 40. For a given IC, its changes with increasing NIC are individually specific irrespective whether the component is a potential resting-state functional network or an artifact component. Using FOCIS, we investigated experimentally the ICA dimensionality of resting-state fMRI datasets and found that the input NIC can critically affect the ICA results of resting-state fMRI data.

Specific default mode subnetworks support mentalizing as revealed through opposing network recruitment by social and semantic FMRI tasks

The ability to attribute mental states to others, or "mentalizing," is posited to involve specific subnetworks within the overall default mode network (DMN), but this question needs clarification. To determine which default mode (DM) subnetworks are engaged by mentalizing processes, we assessed task-related recruitment of DM subnetworks. Spatial independent component analysis (sICA) applied to fMRI data using relatively high-order model (75 components). Healthy participants (n = 53, ages 17-60) performed two fMRI tasks: an interactive game involving mentalizing (Domino), a semantic memory task (SORT), and a resting state fMRI scan. sICA of the two tasks split the DMN into 10 subnetworks located in three core regions: medial prefrontal cortex (mPFC; five subnetworks), posterior cingulate/precuneus (PCC/PrC; three subnetworks), and bilateral temporoparietal junction (TPJ). Mentalizing events increased recruitment in five of 10 DM subnetworks, located in all three core DMN regions. In addition, three of these five DM subnetworks, one dmPFC subnetwork, one PCC/PrC subnetwork, and the right TPJ subnetwork, showed reduced recruitment by semantic memory task events. The opposing modulation by the two tasks suggests that these three DM subnetworks are specifically engaged in mentalizing. Our findings, therefore, suggest the unique involvement of mentalizing processes in only three of 10 DM subnetworks, and support the importance of the dmPFC, PCC/PrC, and right TPJ in mentalizing as described in prior studies.

Functional topography of the thalamocortical system in human

Various studies have indicated that the thalamus is involved in controlling both cortico-cortical information flow and cortical communication with the rest of the brain. Detailed anatomy and functional connectivity patterns of the thalamocortical system are essential to understanding the cortical organization and pathophysiology of a wide range of thalamus-related neurological and neuropsychiatric diseases. The current study used resting-state fMRI to investigate the topography of the human thalamocortical system from a functional perspective. The thalamus-related cortical networks were identified by performing independent component analysis on voxel-based thalamic functional connectivity maps across a large group of subjects. The resulting functional brain networks were very similar to well-established resting-state network maps. Using these brain network components in a spatial regression model with each thalamic voxel's functional connectivity map, we localized the thalamic subdivisions related to each brain network. For instance, the medial dorsal nucleus was shown to be associated with the default mode, the bilateral executive, the medial visual networks; and the pulvinar nucleus was involved in both the dorsal attention and the visual networks. These results revealed that a single nucleus may have functional connections with multiple cortical regions or even multiple functional networks, and may be potentially related to the function of mediation or modulation of multiple cortical networks. This observed organization of thalamocortical system provided a reference for studying the functions of thalamic sub-regions. The importance of intrinsic connectivity-based mapping of the thalamocortical relationship is discussed, as well as the applicability of the approach for future studies.

Enhanced disease characterization through multi network functional normalization in fMRI

Conventionally, structural topology is used for spatial normalization during the pre-processing of fMRI. The co-existence of multiple intrinsic networks which can be detected in the resting brain are well-studied. Also, these networks exhibit temporal and spatial modulation during cognitive task vs. rest which shows the existence of common spatial excitation patterns between these identified networks. Previous work (Khullar et al., 2011) has shown that structural and functional data may not have direct one-to-one correspondence and functional activation patterns in a well-defined structural region can vary across subjects even for a well-defined functional task. The results of this study and the existence of the neural activity patterns in multiple networks motivates us to investigate multiple resting-state networks as a single fusion template for functional normalization for multi groups of subjects. We extend the previous approach (Khullar et al., 2011) by co-registering multi group of subjects (healthy control and schizophrenia patients) and by utilizing multiple resting-state networks (instead of just one) as a single fusion template for functional normalization. In this paper we describe the initial steps toward using multiple resting-state networks as a single fusion template for functional normalization. A simple wavelet-based image fusion approach is presented in order to evaluate the feasibility of combining multiple functional networks. Our results showed improvements in both the significance of group statistics (healthy control and schizophrenia patients) and the spatial extent of activation when a multiple resting-state network applied as a single fusion template for functional normalization after the conventional structural normalization. Also, our results provided evidence that the improvement in significance of group statistics lead to better accuracy results for classification of healthy controls and schizophrenia patients.

How many separable sources? Model selection in independent components analysis

Unlike mixtures consisting solely of non-Gaussian sources, mixtures including two or more Gaussian components cannot be separated using standard independent components analysis methods that are based on higher order statistics and independent observations. The mixed Independent Components Analysis/Principal Components Analysis (mixed ICA/PCA) model described here accommodates one or more Gaussian components in the independent components analysis model and uses principal components analysis to characterize contributions from this inseparable Gaussian subspace. Information theory can then be used to select from among potential model categories with differing numbers of Gaussian components. Based on simulation studies, the assumptions and approximations underlying the Akaike Information Criterion do not hold in this setting, even with a very large number of observations. Cross-validation is a suitable, though computationally intensive alternative for model selection. Application of the algorithm is illustrated using Fisher's iris data set and Howells' craniometric data set. Mixed ICA/PCA is of potential interest in any field of scientific investigation where the authenticity of blindly separated non-Gaussian sources might otherwise be questionable. Failure of the Akaike Information Criterion in model selection also has relevance in traditional independent components analysis where all sources are assumed non-Gaussian.

Modality-spanning deficits in attention-deficit/hyperactivity disorder in functional networks, gray matter, and white matter

Previous neuroimaging investigations in attention-deficit/hyperactivity disorder (ADHD) have separately identified distributed structural and functional deficits, but interconnections between these deficits have not been explored. To unite these modalities in a common model, we used joint independent component analysis, a multivariate, multimodal method that identifies cohesive components that span modalities. Based on recent network models of ADHD, we hypothesized that altered relationships between large-scale networks, in particular, default mode network (DMN) and task-positive networks (TPNs), would co-occur with structural abnormalities in cognitive regulation regions. For 756 human participants in the ADHD-200 sample, we produced gray and white matter volume maps with voxel-based morphometry, as well as whole-brain functional connectomes. Joint independent component analysis was performed, and the resulting transmodal components were tested for differential expression in ADHD versus healthy controls. Four components showed greater expression in ADHD. Consistent with our a priori hypothesis, we observed reduced DMN-TPN segregation co-occurring with structural abnormalities in dorsolateral prefrontal cortex and anterior cingulate cortex, two important cognitive control regions. We also observed altered intranetwork connectivity in DMN, dorsal attention network, and visual network, with co-occurring distributed structural deficits. There was strong evidence of spatial correspondence across modalities: For all four components, the impact of the respective component on gray matter at a region strongly predicted the impact on functional connectivity at that region. Overall, our results demonstrate that ADHD involves multiple, cohesive modality spanning deficits, each one of which exhibits strong spatial overlap in the pattern of structural and functional alterations.

Functional network overlap as revealed by fMRI using sICA and its potential relationships with functional heterogeneity, balanced excitation and inhibition, and sparseness of neuron activity

Functional magnetic resonance imaging (fMRI) studies traditionally use general linear model-based analysis (GLM-BA) and regularly report task-related activation, deactivation, or no change in activation in separate brain regions. However, several recent fMRI studies using spatial independent component analysis (sICA) find extensive overlap of functional networks (FNs), each exhibiting different task-related modulation (e.g., activation vs. deactivation), different from the dominant findings of GLM-BA. This study used sICA to assess overlap of FNs extracted from four datasets, each related to a different cognitive task. FNs extracted from each dataset overlapped with each other extensively across most or all brain regions and showed task-related concurrent increases, decreases, or no changes in activity. These findings indicate that neural substrates showing task-related concurrent but different modulations in activity intermix with each other and distribute across most of the brain. Furthermore, spatial correlation analyses found that most FNs were highly consistent in spatial patterns across different datasets. This finding indicates that these FNs probably reflect large-scale patterns of task-related brain activity. We hypothesize that FN overlaps as revealed by sICA might relate to functional heterogeneity, balanced excitation and inhibition, and population sparseness of neuron activity, three fundamental properties of the brain. These possibilities deserve further investigation.

High-Dimensional ICA Analysis Detects Within-Network Functional Connectivity Damage of Default-Mode and Sensory-Motor Networks in Alzheimer's Disease

High-dimensional independent component analysis (ICA), compared to low-dimensional ICA, allows to conduct a detailed parcellation of the resting-state networks. The purpose of this study was to give further insight into functional connectivity (FC) in Alzheimer's disease (AD) using high-dimensional ICA. For this reason, we performed both low- and high-dimensional ICA analyses of resting-state fMRI data of 20 healthy controls and 21 patients with AD, focusing on the primarily altered default-mode network (DMN) and exploring the sensory-motor network. As expected, results obtained at low dimensionality were in line with previous literature. Moreover, high-dimensional results allowed us to observe either the presence of within-network disconnections and FC damage confined to some of the resting-state subnetworks. Due to the higher sensitivity of the high-dimensional ICA analysis, our results suggest that high-dimensional decomposition in subnetworks is very promising to better localize FC alterations in AD and that FC damage is not confined to the DMN.

Aberrant Functional Connectivity in the Default Mode and Central Executive Networks in Subjects with Schizophrenia - A Whole-Brain Resting-State ICA Study

Neurophysiological changes of schizophrenia are currently linked to disturbances in connectivity between functional brain networks. Functional magnetic resonance imaging studies on schizophrenia have focused on a few selected networks. Also previously, it has not been possible to discern whether the functional alterations in schizophrenia originate from spatial shifting or amplitude alterations of functional connectivity. In this study, we aim to discern the differences in schizophrenia patients with respect to spatial shifting vs. signal amplitude changes in functional connectivity in the whole-brain connectome. We used high model order-independent component analysis to study some 40 resting-state networks (RSN) covering the whole cortex. Group differences were analyzed with dual regression coupled with y-concat correction for multiple comparisons. We investigated the RSNs with and without variance normalization in order to discern spatial shifting from signal amplitude changes in 43 schizophrenia patients and matched controls from the Northern Finland 1966 Birth Cohort. Voxel-level correction for multiple comparisons revealed 18 RSNs with altered functional connectivity, 6 of which had both spatial and signal amplitude changes. After adding the multiple comparison, y-concat correction to the analysis for including the 40 RSNs as well, we found that four RSNs showed still changes. These robust changes actually seem encompass parcellations of the default mode network and central executive networks. These networks both have spatially shifted connectivity and abnormal signal amplitudes. Interestingly the networks seem to mix their functional representations in areas like left caudate nucleus and dorsolateral prefrontal cortex. These changes overlapped with areas that have been related to dopaminergic alterations in patients with schizophrenia compared to controls.

Increased Amplitude of Low Frequency Fluctuations but Normal Hippocampal-Default Mode Network Connectivity in Schizophrenia

BACKGROUND

Clinical and preclinical studies have established that the hippocampus is hyperactive in schizophrenia, making it a possible biomarker for drug development. Increased hippocampal connectivity, which can be studied conveniently with resting state imaging, has been proposed as a readily accessible corollary of hippocampal hyperactivity. Here, we tested the hypothesis that hippocampal activity and connectivity are increased in patients with schizophrenia.

METHODS

Sixty-three schizophrenia patients and 71 healthy control subjects completed a resting state functional magnetic resonance imaging scan. We assessed hippocampal activity with the amplitude of low frequency fluctuations. We analyzed hippocampal functional connectivity with the default mode network using three common methods: group and single subject level independent component analysis, and seed-based functional connectivity.

RESULTS

In patients with schizophrenia, we observed increased amplitude of low frequency fluctuations but normal hippocampal connectivity using independent component and seed-based analyses.

CONCLUSION

Our results indicate that although intrinsic hippocampal activity may be increased in schizophrenia, this finding does not extend to functional connectivity. Neuroimaging methods that directly assess hippocampal activity may be more promising for the identification of a biomarker for schizophrenia.

Assessing dynamic brain graphs of time-varying connectivity in fMRI data: application to healthy controls and patients with schizophrenia

Graph theory-based analysis has been widely employed in brain imaging studies, and altered topological properties of brain connectivity have emerged as important features of mental diseases such as schizophrenia. However, most previous studies have focused on graph metrics of stationary brain graphs, ignoring that brain connectivity exhibits fluctuations over time. Here we develop a new framework for accessing dynamic graph properties of time-varying functional brain connectivity in resting-state fMRI data and apply it to healthy controls (HCs) and patients with schizophrenia (SZs). Specifically, nodes of brain graphs are defined by intrinsic connectivity networks (ICNs) identified by group independent component analysis (ICA). Dynamic graph metrics of the time-varying brain connectivity estimated by the correlation of sliding time-windowed ICA time courses of ICNs are calculated. First- and second-level connectivity states are detected based on the correlation of nodal connectivity strength between time-varying brain graphs. Our results indicate that SZs show decreased variance in the dynamic graph metrics. Consistent with prior stationary functional brain connectivity works, graph measures of identified first-level connectivity states show lower values in SZs. In addition, more first-level connectivity states are disassociated with the second-level connectivity state which resembles the stationary connectivity pattern computed by the entire scan. Collectively, the findings provide new evidence about altered dynamic brain graphs in schizophrenia, which may underscore the abnormal brain performance in this mental illness.

Functional connectivity in BOLD and CBF data: similarity and reliability of resting brain networks

Resting-state functional connectivity (FC) fMRI (rs-fcMRI) offers an appealing approach to mapping the brain's intrinsic functional organization. Blood oxygen level dependent (BOLD) and arterial spin labeling (ASL) are the two main rs-fcMRI approaches to assess alterations in brain networks associated with individual differences, behavior and psychopathology. While the BOLD signal is stronger with a higher temporal resolution, ASL provides quantitative, direct measures of the physiology and metabolism of specific networks. This study systematically investigated the similarity and reliability of resting brain networks (RBNs) in BOLD and ASL. A 2 × 2 × 2 factorial design was employed where each subject underwent repeated BOLD and ASL rs-fcMRI scans on two occasions on two MRI scanners respectively. Both independent and joint FC analyses revealed common RBNs in ASL and BOLD rs-fcMRI with a moderate to high level of spatial overlap, verified by Dice Similarity Coefficients. Test-retest analyses indicated more reliable spatial network patterns in BOLD (average modal Intraclass Correlation Coefficients: 0.905 ± 0.033 between-sessions; 0.885 ± 0.052 between-scanners) than ASL (0.545 ± 0.048; 0.575 ± 0.059). Nevertheless, ASL provided highly reproducible (0.955 ± 0.021; 0.970 ± 0.011) network-specific CBF measurements. Moreover, we observed positive correlations between regional CBF and FC in core areas of all RBNs indicating a relationship between network connectivity and its baseline metabolism. Taken together, the combination of ASL and BOLD rs-fcMRI provides a powerful tool for characterizing the spatiotemporal and quantitative properties of RBNs. These findings pave the way for future BOLD and ASL rs-fcMRI studies in clinical populations that are carried out across time and scanners.

Effects of Beta-Amyloid on Resting State Functional Connectivity Within and Between Networks Reflect Known Patterns of Regional Vulnerability

Beta-amyloid (Aβ) deposition is one of the hallmarks of Alzheimer's disease (AD). However, it is also present in some cognitively normal elderly adults and may represent a preclinical disease state. While AD patients exhibit disrupted functional connectivity (FC) both within and between resting-state networks, studies of preclinical cases have focused primarily on the default mode network (DMN). The extent to which Aβ-related effects occur outside of the DMN and between networks remains unclear. In the present study, we examine how within- and between-network FC are related to both global and regional Aβ deposition as measured by [(11)C]PIB-PET in 92 cognitively normal older people. We found that within-network FC changes occurred in multiple networks, including the DMN. Changes of between-network FC were also apparent, suggesting that regions maintaining connections to multiple networks may be particularly susceptible to Aβ-induced alterations. Cortical regions showing altered FC clustered in parietal and temporal cortex, areas known to be susceptible to AD pathology. These results likely represent a mix of local network disruption, compensatory reorganization, and impaired control network function. They indicate the presence of Aβ-related dysfunction of neural systems in cognitively normal people well before these areas become hypometabolic with the onset of cognitive decline.

Dynamic connectivity states estimated from resting fMRI Identify differences among Schizophrenia, bipolar disorder, and healthy control subjects

Schizophrenia (SZ) and bipolar disorder (BP) share significant overlap in clinical symptoms, brain characteristics, and risk genes, and both are associated with dysconnectivity among large-scale brain networks. Resting state functional magnetic resonance imaging (rsfMRI) data facilitates studying macroscopic connectivity among distant brain regions. Standard approaches to identifying such connectivity include seed-based correlation and data-driven clustering methods such as independent component analysis (ICA) but typically focus on average connectivity. In this study, we utilize ICA on rsfMRI data to obtain intrinsic connectivity networks (ICNs) in cohorts of healthy controls (HCs) and age matched SZ and BP patients. Subsequently, we investigated difference in functional network connectivity, defined as pairwise correlations among the timecourses of ICNs, between HCs and patients. We quantified differences in both static (average) and dynamic (windowed) connectivity during the entire scan duration. Disease-specific differences were identified in connectivity within different dynamic states. Notably, results suggest that patients make fewer transitions to some states (states 1, 2, and 4) compared to HCs, with most such differences confined to a single state. SZ patients showed more differences from healthy subjects than did bipolars, including both hyper and hypo connectivity in one common connectivity state (dynamic state 3). Also group differences between SZ and bipolar patients were identified in patterns (states) of connectivity involving the frontal (dynamic state 1) and frontal-parietal regions (dynamic state 3). Our results provide new information about these illnesses and strongly suggest that state-based analyses are critical to avoid averaging together important factors that can help distinguish these clinical groups.

Altered functional connectivity between emotional and cognitive resting state networks in euthymic bipolar I disorder patients

Bipolar disorder is characterized by a functional imbalance between hyperactive ventral/limbic areas and hypoactive dorsal/cognitive brain regions potentially contributing to affective and cognitive symptoms. Resting-state studies in bipolar disorder have identified abnormal functional connectivity between these brain regions. However, most of these studies used a seed-based approach, thus restricting the number of regions that were analyzed. Using data-driven approaches, researchers identified resting state networks whose spatial maps overlap with frontolimbic areas such as the default mode network, the frontoparietal networks, the salient network, and the meso/paralimbic network. These networks are specifically engaged during affective and cognitive tasks and preliminary evidence suggests that functional connectivity within and between some of these networks is impaired in bipolar disorder. The present study used independent component analysis and functional network connectivity approaches to investigate functional connectivity within and between these resting state networks in bipolar disorder. We compared 30 euthymic bipolar I disorder patients and 35 age- and gender-matched healthy controls. Inter-network connectivity analysis revealed increased functional connectivity between the meso/paralimbic and the right frontoparietal network in bipolar disorder. This abnormal connectivity pattern did not correlate with variables related to the clinical course of the disease. The present finding may reflect abnormal integration of affective and cognitive information in ventral-emotional and dorsal-cognitive networks in euthymic bipolar patients. Furthermore, the results provide novel insights into the role of the meso/paralimbic network in bipolar disorder.

The rich get richer: brain injury elicits hyperconnectivity in core subnetworks

There remains much unknown about how large-scale neural networks accommodate neurological disruption, such as moderate and severe traumatic brain injury (TBI). A primary goal in this study was to examine the alterations in network topology occurring during the first year of recovery following TBI. To do so we examined 21 individuals with moderate and severe TBI at 3 and 6 months after resolution of posttraumatic amnesia and 15 age- and education-matched healthy adults using functional MRI and graph theoretical analyses. There were two central hypotheses in this study: 1) physical disruption results in increased functional connectivity, or hyperconnectivity, and 2) hyperconnectivity occurs in regions typically observed to be the most highly connected cortical hubs, or the "rich club". The current findings generally support the hyperconnectivity hypothesis showing that during the first year of recovery after TBI, neural networks show increased connectivity, and this change is disproportionately represented in brain regions belonging to the brain's core subnetworks. The selective increases in connectivity observed here are consistent with the preferential attachment model underlying scale-free network development. This study is the largest of its kind and provides the unique opportunity to examine how neural systems adapt to significant neurological disruption during the first year after injury.