Visual inspection of independent components: defining a procedure for artifact removal from fMRI data

Functional integrity in children with anoxic brain injury from drowning

Drowning is a leading cause of accidental injury and death in young children. Anoxic brain injury (ABI) is a common consequence of drowning and can cause severe neurological morbidity in survivors. Assessment of functional status and prognostication in drowning victims can be extremely challenging, both acutely and chronically. Structural neuroimaging modalities (CT and MRI) have been of limited clinical value. Here, we tested the utility of resting-state functional MRI (rs-fMRI) for assessing brain functional integrity in this population. Eleven children with chronic, spastic quadriplegia due to drowning-induced ABI were investigated. All were comatose immediately after the injury and gradually regained consciousness, but with varying ability to communicate their cognitive state. Eleven neurotypical children matched for age and gender formed the control group. Resting-state fMRI and co-registered T1-weighted anatomical MRI were acquired at night during drug-aided sleep. Network integrity was quantified by independent components analysis (ICA), at both group- and per-subject levels. Functional-status assessments based on in-home observations were provided by families and caregivers. Motor ICNs were grossly compromised in ABI patients both group-wise and individually, concordant with their prominent motor deficits. Striking preservations of perceptual and cognitive ICNs were observed, and the degree of network preservation correlated (ρ = 0.74) with the per-subject functional status assessments. Collectively, our findings indicate that rs-fMRI has promise for assessing brain functional integrity in ABI and, potentially, in other disorders. Furthermore, our observations suggest that the severe motor deficits observed in this population can mask relatively intact perceptual and cognitive capabilities. Hum Brain Mapp 38:4813-4831, 2017. © 2017 Wiley Periodicals, Inc.

Somatic and Reinforcement-Based Plasticity in the Initial Stages of Human Motor Learning

As one learns to dance or play tennis, the desired somatosensory state is typically unknown. Trial and error is important as motor behavior is shaped by successful and unsuccessful movements. As an experimental model, we designed a task in which human participants make reaching movements to a hidden target and receive positive reinforcement when successful. We identified somatic and reinforcement-based sources of plasticity on the basis of changes in functional connectivity using resting-state fMRI before and after learning. The neuroimaging data revealed reinforcement-related changes in both motor and somatosensory brain areas in which a strengthening of connectivity was related to the amount of positive reinforcement during learning. Areas of prefrontal cortex were similarly altered in relation to reinforcement, with connectivity between sensorimotor areas of putamen and the reward-related ventromedial prefrontal cortex strengthened in relation to the amount of successful feedback received. In other analyses, we assessed connectivity related to changes in movement direction between trials, a type of variability that presumably reflects exploratory strategies during learning. We found that connectivity in a network linking motor and somatosensory cortices increased with trial-to-trial changes in direction. Connectivity varied as well with the change in movement direction following incorrect movements. Here the changes were observed in a somatic memory and decision making network involving ventrolateral prefrontal cortex and second somatosensory cortex. Our results point to the idea that the initial stages of motor learning are not wholly motor but rather involve plasticity in somatic and prefrontal networks related both to reward and exploration.

SIGNIFICANCE STATEMENT

In the initial stages of motor learning, the placement of the limbs is learned primarily through trial and error. In an experimental analog, participants make reaching movements to a hidden target and receive positive feedback when successful. We identified sources of plasticity based on changes in functional connectivity using resting-state fMRI. The main finding is that there is a strengthening of connectivity between reward-related prefrontal areas and sensorimotor areas in the basal ganglia and frontal cortex. There is also a strengthening of connectivity related to movement exploration in sensorimotor circuits involved in somatic memory and decision making. The results indicate that initial stages of motor learning depend on plasticity in somatic and prefrontal networks related to reward and exploration.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Impact of automated ICA-based denoising of fMRI data in acute stroke patients

Different strategies have been developed using Independent Component Analysis (ICA) to automatically de-noise fMRI data, either focusing on removing only certain components (e.g. motion-ICA-AROMA, Pruim et al., 2015a) or using more complex classifiers to remove multiple types of noise components (e.g. FIX, Salimi-Khorshidi et al., 2014 Griffanti et al., 2014). However, denoising data obtained in an acute setting might prove challenging: the presence of multiple noise sources may not allow focused strategies to clean the data enough and the heterogeneity in the data may be so great to critically undermine complex approaches. The purpose of this study was to explore what automated ICA based approach would better cope with these limitations when cleaning fMRI data obtained from acute stroke patients. The performance of a focused classifier (ICA-AROMA) and a complex classifier (FIX) approaches were compared using data obtained from twenty consecutive acute lacunar stroke patients using metrics determining RSN identification, RSN reproducibility, changes in the BOLD variance, differences in the estimation of functional connectivity and loss of temporal degrees of freedom. The use of generic-trained FIX resulted in misclassification of components and significant loss of signal (< 80%), and was not explored further. Both ICA-AROMA and patient-trained FIX based denoising approaches resulted in significantly improved RSN reproducibility (p < 0.001), localized reduction in BOLD variance consistent with noise removal, and significant changes in functional connectivity (p < 0.001). Patient-trained FIX resulted in higher RSN identifiability (p < 0.001) and wider changes both in the BOLD variance and in functional connectivity compared to ICA-AROMA. The success of ICA-AROMA suggests that by focusing on selected components the full automation can deliver meaningful data for analysis even in population with multiple sources of noise. However, the time invested to train FIX with appropriate patient data proved valuable, particularly in improving the signal-to-noise ratio.

Intrinsic Connectivity Network-Based Classification and Detection of Psychotic Symptoms in Youth With 22q11.2 Deletions

22q11.2 Deletion syndrome (22q11DS) is a genetic disorder associated with numerous phenotypic consequences and is one of the greatest known risk factors for psychosis. We investigated intrinsic-connectivity-networks (ICNs) as potential biomarkers for patient and psychosis-risk status in 2 independent cohorts, UCLA (33 22q11DS-participants, 33 demographically matched controls), and Syracuse (28 22q11DS, 28 controls). After assessing group connectivity differences, ICNs from the UCLA cohort were used to train classifiers to distinguish cases from controls, and to predict psychosis risk status within 22q11DS; classifiers were subsequently tested on the Syracuse cohort. In both cohorts we observed significant hypoconnectivity in 22q11DS relative to controls within anterior cingulate (ACC)/precuneus, executive, default mode (DMN), posterior DMN, and salience networks. Of 12 ICN-derived classifiers tested in the Syracuse replication-cohort, the ACC/precuneus, DMN, and posterior DMN classifiers accurately distinguished between 22q11DS and controls. Within 22q11DS subjects, connectivity alterations within 4 networks predicted psychosis risk status for a given individual in both cohorts: the ACC/precuneus, DMN, left executive, and salience networks. Widespread within-network-hypoconnectivity in large-scale networks implicated in higher-order cognition may be a defining characteristic of 22q11DS during adolescence and early adulthood; furthermore, loss of coherence within these networks may be a valuable biomarker for individual prediction of psychosis-risk in 22q11DS.

Comparison of large-scale human brain functional and anatomical networks in schizophrenia

Schizophrenia is a disease with disruptions in thought, emotion, and behavior. The dysconnectivity hypothesis suggests these disruptions are due to aberrant brain connectivity. Many studies have identified connectivity differences but few have been able to unify gray and white matter findings into one model. Here we develop an extension of the Network-Based Statistic (NBS) called NBSm (Multimodal Network-based statistic) to compare functional and anatomical networks in schizophrenia. Structural, resting functional, and diffusion magnetic resonance imaging data were collected from 29 chronic patients with schizophrenia and 29 healthy controls. Images were preprocessed, and average time courses were extracted for 90 regions of interest (ROI). Functional connectivity matrices were estimated by pairwise correlations between wavelet coefficients of ROI time series. Following diffusion tractography, anatomical connectivity matrices were estimated by white matter streamline counts between each pair of ROIs. Global and regional strength were calculated for each modality. NBSm was used to find significant overlap between functional and anatomical components that distinguished health from schizophrenia. Global strength was decreased in patients in both functional and anatomical networks. Regional strength was decreased in all regions in functional networks and only one region in anatomical networks. NBSm identified a distinguishing functional component consisting of 46 nodes with 113 links (p < 0.001), a distinguishing anatomical component with 47 nodes and 50 links (p = 0.002), and a distinguishing intermodal component with 26 nodes (p < 0.001). NBSm is a powerful technique for understanding network-based group differences present in both anatomical and functional data. In light of the dysconnectivity hypothesis, these results provide compelling evidence for the presence of significant overlapping anatomical and functional disruption in people with schizophrenia.

•

A novel method for the unified analysis of functional and anatomical networks in schizophrenia is proposed.

•

The methodology extends existing networkbased analyses to test brain network differences across imaging modalities.

•

A fundamental disrupted network was found in those with schizophrenia when compared to healthy, agematched controls.

Adverse Effects of Cannabis on Adolescent Brain Development: A Longitudinal Study

Cannabis is widely perceived as a safe recreational drug and its use is increasing in youth. It is important to understand the implications of cannabis use during childhood and adolescence on brain development. This is the first longitudinal study that compared resting functional connectivity of frontally mediated networks between 43 healthy controls (HCs; 20 females; age M = 16.5 ± 2.7) and 22 treatment-seeking adolescents with cannabis use disorder (CUD; 8 females; age M = 17.6 ± 2.4). Increases in resting functional connectivity between caudal anterior cingulate cortex (ACC) and superior frontal gyrus across time were found in HC, but not in CUD. CUD showed a decrease in functional connectivity between caudal ACC and dorsolateral and orbitofrontal cortices across time. Lower functional connectivity between caudal ACC cortex and orbitofrontal cortex at baseline predicted higher amounts of cannabis use during the following 18 months. Finally, high amounts of cannabis use during the 18-month interval predicted lower intelligence quotient and slower cognitive function measured at follow-up. These data provide compelling longitudinal evidence suggesting that repeated exposure to cannabis during adolescence may have detrimental effects on brain resting functional connectivity, intelligence, and cognitive function.

Functional Connectivity Between the Temporoparietal Cortex and Cerebellum in Autism Spectrum Disorder

The neural basis of autism spectrum disorder (ASD) is not yet understood. ASD is marked by social deficits and is strongly associated with cerebellar abnormalities. We studied the organization and cerebellar connectivity of the temporoparietal junction (TPJ), an area that plays a crucial role in social cognition. We applied localized independent component analysis to resting-state fMRI data from autistic and neurotypical adolescents to yield an unbiased parcellation of the bilateral TPJ into 11 independent components (ICs). A comparison between neurotypical and autistic adolescents showed that the organization of the TPJ was not significantly altered in ASD. Second, we used the time courses of the TPJ ICs as spatially unbiased "seeds" for a functional connectivity analysis applied to voxels within the cerebellum. We found that the cerebellum contained a fine-grained, lateralized map of the TPJ. The connectivity of the TPJ subdivisions with cerebellar zones showed one striking difference in ASD. The right dorsal TPJ showed markedly less connectivity with the left Crus II. Disturbed cerebellar input to this key region for cognition and multimodal integration may contribute to social deficits in ASD. The findings might also suggest that the right TPJ and/or left Crus II are potential targets for noninvasive brain stimulation therapies.

Relationship Between Changes in the Temporal Dynamics of the Blood-Oxygen-Level-Dependent Signal and Hypoperfusion in Acute Ischemic Stroke

Background and Purpose—

Changes in the blood-oxygen-level-dependent (BOLD) signal provide a noninvasive measure of blood flow, but a detailed comparison with established perfusion parameters in acute stroke is lacking. We investigated the relationship between BOLD signal temporal delay and dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) in stroke patients.

Methods—

In 30 patients with acute (<24 hours) ischemic stroke, we performed Pearson correlation and multiple linear regression between DSC-MRI parameters (time to maximum [Tmax], mean transit time, cerebral blood flow, and cerebral blood volume) and BOLD-based parameters (BOLD delay and coefficient of BOLD variation). Prediction of severe hypoperfusion (Tmax >6 seconds) was assessed using receiver–operator characteristic (ROC) analysis.

Results—

Correlation was highest between Tmax and BOLD delay (venous sinus reference; time shift range 7; median r =0.60; interquartile range=0.49–0.71). Coefficient of BOLD variation correlated with cerebral blood volume (median r = 0.37; interquartile range=0.24–0.51). Mean R 2 for predicting BOLD delay by DSC-MRI was 0.54 (SD=0.2) and for predicting coefficient of BOLD variation was 0.37 (SD=0.17). BOLD delay (whole-brain reference, time shift range 3) had an area under the curve of 0.76 for predicting severe hypoperfusion (sensitivity=69.2%; specificity=80%), whereas BOLD delay (venous sinus reference, time shift range 3) had an area under the curve of 0.76 (sensitivity=67.3%; specificity=83.5%).

Conclusions—

BOLD delay is related to macrovascular delay and microvascular hypoperfusion, can identify severely hypoperfused tissue in acute stroke, and is a promising alternative to gadolinium contrast agent–based perfusion assessment in acute stroke.

Clinical Trial Registration—

URL: http://www.clinicaltrials.gov . Unique identifier: NCT00715533 and NCT02077582.

Default mode functional connectivity is associated with social functioning in schizophrenia

Individuals with schizophrenia display notable deficits in social functioning. Research indicates that neural connectivity within the default mode network (DMN) is related to social cognition and social functioning in healthy and clinical populations. However, the association between DMN connectivity, social cognition, and social functioning has not been studied in schizophrenia. For the present study, the authors used resting-state neuroimaging data to evaluate connectivity between the main DMN hubs (i.e., the medial prefrontal cortex [mPFC] and the posterior cingulate cortex-anterior precuneus [PPC]) in individuals with schizophrenia (n = 28) and controls (n = 32). The authors also examined whether DMN connectivity was associated with social functioning via social attainment (measured by the Specific Levels of Functioning Scale) and social competence (measured by the Social Skills Performance Assessment), and if social cognition mediates the association between DMN connectivity and these measures of social functioning. Results revealed that DMN connectivity did not differ between individuals with schizophrenia and controls. However, connectivity between the mPFC and PCC hubs was significantly associated with social competence and social attainment in individuals with schizophrenia but not in controls as reflected by a significant group-by-connectivity interaction. Social cognition did not mediate the association between DMN connectivity and social functioning in individuals with schizophrenia. The findings suggest that fronto-parietal DMN connectivity in particular may be differentially associated with social functioning in schizophrenia and controls. As a result, DMN connectivity may be used as a neuroimaging marker to monitor treatment response or as a potential target for interventions that aim to enhance social functioning in schizophrenia. (PsycINFO Database Record

Resting state synchrony in long-term abstinent alcoholics: Effects of a current major depressive disorder diagnosis

Alcoholism is characterized by a lack of control over an impulsive and compulsive drive toward excessive alcohol consumption despite significant negative consequences; our previous work demonstrated that successful abstinence is characterized by decreased resting-state synchrony (RSS) as measured with functional magnetic resonance imaging (fMRI), within appetitive drive networks and increased RSS in emotion regulation and inhibitory executive control networks. Our hypothesis is that LTAA (Long-Term Abstinent Alcoholics) with a current major depressive disorder (MDD) drank primarily to deal with the negative affect associated with their MDD and not because of a heightened externalizing diathesis (including heightened appetitive drive), and consequently, in achieving and maintaining abstinence, such individuals would not exhibit the RSS adaptations characteristic of pure alcoholics. We studied 69 NSAC (Non Substance Abusing Controls) and 40 LTAA (8 with current MDD, 32 without a current MDD) using resting-state fMRI and seed based connectivity analyses. In the inhibitory executive control network (nucleus accumbens vs. left dorsolateral prefrontal cortex), LTAA with a current MDD showed increased synchrony compared to NSAC. In the emotion regulation executive control network (subgenual anterior cingulate cortex vs. right dorsolateral prefrontal cortex), LTAA with current MDD did not show increased RSS. In the appetitive drive networks (nucleus accumbens vs, aspects of the caudate nucleus and thalamus), LTAA with a current MDD did not show a reduction of RSS compared to NSAC, but LTAA without a current MDD did. These results suggest different pathways to their alcohol dependence in LTAA with vs. without a current MDD, and different patterns of brain activity in long-term abstinence, suggesting different treatment needs.

The cortical connectivity of the periaqueductal gray and the conditioned response to the threat of breathlessness

Previously we observed differential activation in individual columns of the periaqueductal grey (PAG) during breathlessness and its conditioned anticipation (Faull et al., 2016b). Here, we have extended this work by determining how the individual columns of the PAG interact with higher cortical centres, both at rest and in the context of breathlessness threat. Activation was observed in ventrolateral PAG (vlPAG) and lateral PAG (lPAG), where activity scaled with breathlessness intensity ratings, revealing a potential interface between sensation and cognition during breathlessness. At rest the lPAG was functionally correlated with cortical sensorimotor areas, conducive to facilitating fight/flight responses, and demonstrated increased synchronicity with the amygdala during breathlessness. The vlPAG showed fronto-limbic correlations at rest, whereas during breathlessness anticipation, reduced functional synchronicity was seen to both lPAG and motor structures, conducive to freezing behaviours. These results move us towards understanding how the PAG might be intricately involved in human responses to threat.

Frequency-specific electrophysiologic correlates of resting state fMRI networks

Resting state functional MRI (R-fMRI) studies have shown that slow (< 0.1 Hz), intrinsic fluctuations of the blood oxygen level dependent (BOLD) signal are temporally correlated within hierarchically organized functional systems known as resting state networks (RSNs) (Doucet et al., 2011). Most broadly, this hierarchy exhibits a dichotomy between two opposed systems (Fox et al., 2005). One system engages with the environment and includes the visual, auditory, and sensorimotor (SMN) networks as well as the dorsal attention network (DAN), which controls spatial attention. The other system includes the default mode network (DMN) and the fronto-parietal control system (FPC), RSNs that instantiate episodic memory and executive control, respectively. Here, we test the hypothesis, based on the spectral specificity of electrophysiologic responses to perceptual vs. memory tasks (Klimesch, 1999; Pfurtscheller and Lopes da Silva, 1999), that these two large-scale neural systems also manifest frequency specificity in the resting state. We measured the spatial correspondence between electrocorticographic (ECoG) band-limited power (BLP) and R-fMRI correlation patterns in awake, resting, human subjects. Our results show that, while gamma BLP correspondence was common throughout the brain, theta (4–8 Hz) BLP correspondence was stronger in the DMN and FPC, whereas alpha (8–12 Hz) correspondence was stronger in the SMN and DAN. Thus, the human brain, at rest, exhibits frequency specific electrophysiology, respecting both the spectral structure of task responses and the hierarchical organization of RSNs.

Increased sensitivity to age-related differences in brain functional connectivity during continuous multiple object tracking compared to resting-state

Age-related differences in cognitive agility vary greatly between individuals and cognitive functions. This heterogeneity is partly mirrored in individual differences in brain network connectivity as revealed using resting-state functional magnetic resonance imaging (fMRI), suggesting potential imaging biomarkers for age-related cognitive decline. However, although convenient in its simplicity, the resting state is essentially an unconstrained paradigm with minimal experimental control. Here, based on the conception that the magnitude and characteristics of age-related differences in brain connectivity is dependent on cognitive context and effort, we tested the hypothesis that experimentally increasing cognitive load boosts the sensitivity to age and changes the discriminative network configurations. To this end, we obtained fMRI data from younger (n=25, mean age 24.16±5.11) and older (n=22, mean age 65.09±7.53) healthy adults during rest and two load levels of continuous multiple object tracking (MOT). Brain network nodes and their time-series were estimated using independent component analysis (ICA) and dual regression, and the edges in the brain networks were defined as the regularized partial temporal correlations between each of the node pairs at the individual level. Using machine learning based on a cross-validated regularized linear discriminant analysis (rLDA) we attempted to classify groups and cognitive load from the full set of edge-wise functional connectivity indices. While group classification using resting-state data was highly above chance (approx. 70% accuracy), functional connectivity (FC) obtained during MOT strongly increased classification performance, with 82% accuracy for the young and 95% accuracy for the old group at the highest load level. Further, machine learning revealed stronger differentiation between rest and task in young compared to older individuals, supporting the notion of network dedifferentiation in cognitive aging. Task-modulation in edgewise FC was primarily observed between attention- and sensorimotor networks; with decreased negative correlations between attention- and default mode networks in older adults. These results demonstrate that the magnitude and configuration of age-related differences in brain functional connectivity are partly dependent on cognitive context and load, which emphasizes the importance of assessing brain connectivity differences across a range of cognitive contexts beyond the resting-state.

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.

Disrupted global metastability and static and dynamic brain connectivity across individuals in the Alzheimer's disease continuum

As findings on the neuropathological and behavioral components of Alzheimer's disease (AD) continue to accrue, converging evidence suggests that macroscale brain functional disruptions may mediate their association. Recent developments on theoretical neuroscience indicate that instantaneous patterns of brain connectivity and metastability may be a key mechanism in neural communication underlying cognitive performance. However, the potential significance of these patterns across the AD spectrum remains virtually unexplored. We assessed the clinical sensitivity of static and dynamic functional brain disruptions across the AD spectrum using resting-state fMRI in a sample consisting of AD patients (n = 80) and subjects with either mild (n = 44) or subjective (n = 26) cognitive impairment (MCI, SCI). Spatial maps constituting the nodes in the functional brain network and their associated time-series were estimated using spatial group independent component analysis and dual regression, and whole-brain oscillatory activity was analyzed both globally (metastability) and locally (static and dynamic connectivity). Instantaneous phase metrics showed functional coupling alterations in AD compared to MCI and SCI, both static (putamen, dorsal and default-mode) and dynamic (temporal, frontal-superior and default-mode), along with decreased global metastability. The results suggest that brains of AD patients display altered oscillatory patterns, in agreement with theoretical premises on cognitive dynamics.

Role of testosterone and Y chromosome genes for the masculinization of the human brain

Women with complete androgen insensitivity syndrome (CAIS) have a male (46,XY) karyotype but no functional androgen receptors. Their condition, therefore, offers a unique model for studying testosterone effects on cerebral sex dimorphism. We present MRI data from 16 women with CAIS and 32 male (46,XY) and 32 female (46,XX) controls.

METHODS

FreeSurfer software was employed to measure cortical thickness and subcortical structural volumes. Axonal connections, indexed by fractional anisotropy, (FA) were measured with diffusion tensor imaging, and functional connectivity with resting state fMRI.

RESULTS

Compared to men, CAIS women displayed a "female" pattern by having thicker parietal and occipital cortices, lower FA values in the right corticospinal, superior and inferior longitudinal tracts, and corpus callosum. Their functional connectivity from the amygdala to the medial prefrontal cortex, was stronger and amygdala-connections to the motor cortex weaker than in control men. CAIS and control women also showed stronger posterior cingulate and precuneus connections in the default mode network. Thickness of the motor cortex, the caudate volume, and the FA in the callosal body followed, however, a "male" pattern.

CONCLUSION

Altogether, these data suggest that testosterone modulates the microstructure of somatosensory and visual cortices and their axonal connections to the frontal cortex. Testosterone also influenced functional connections from the amygdala, whereas the motor cortex could, in agreement with our previous reports, be moderated by processes linked to X-chromosome gene dosage. These data raise the question about other genetic factors masculinizing the human brain than the SRY gene and testosterone. Hum Brain Mapp 38:1801-1814, 2017. © 2017 Wiley Periodicals, Inc.

Excitation and inhibition in anterior cingulate predict use of past experiences

Dorsal anterior cingulate cortex (dACC) mediates updating and maintenance of cognitive models of the world used to drive adaptive reward-guided behavior. We investigated the neurochemical underpinnings of this process. We used magnetic resonance spectroscopy in humans, to measure levels of glutamate and GABA in dACC. We examined their relationship to neural signals in dACC, measured with fMRI, and cognitive task performance. Both inhibitory and excitatory neurotransmitters in dACC were predictive of the strength of neural signals in dACC and behavioral adaptation. Glutamate levels were correlated, first, with stronger neural activity representing information to be learnt about the tasks' costs and benefits and, second, greater use of this information in the guidance of behavior. GABA levels were negatively correlated with the same neural signals and the same indices of behavioral influence. Our results suggest that glutamate and GABA in dACC affect the encoding and use of past experiences to guide behavior.

Opioid suppression of conditioned anticipatory brain responses to breathlessness

Opioid painkillers are a promising treatment for chronic breathlessness, but are associated with potentially fatal side effects. In the treatment of breathlessness, their mechanisms of action are unclear. A better understanding might help to identify safer alternatives. Learned associations between previously neutral stimuli (e.g. stairs) and repeated breathlessness induce an anticipatory threat response that may worsen breathlessness, contributing to the downward spiral of decline seen in clinical populations. As opioids are known to influence associative learning, we hypothesized that they may interfere with the brain processes underlying a conditioned anticipatory response to breathlessness in relevant brain areas, including the amygdala and the hippocampus.

Healthy volunteers viewed visual cues (neutral stimuli) immediately before induction of experimental breathlessness with inspiratory resistive loading. Thus, an association was formed between the cue and breathlessness. Subsequently, this paradigm was repeated in two identical neuroimaging sessions with intravenous infusions of either low-dose remifentanil (0.7 ng/ml target-controlled infusion) or saline (randomised).

During saline infusion, breathlessness anticipation activated the right anterior insula and the adjacent operculum. Breathlessness was associated with activity in a network including the insula, operculum, dorsolateral prefrontal cortex, anterior cingulate cortex and the primary sensory and motor cortices.

Remifentanil reduced breathlessness unpleasantness but not breathlessness intensity. Remifentanil depressed anticipatory activity in the amygdala and the hippocampus that correlated with reductions in breathlessness unpleasantness. During breathlessness, remifentanil decreased activity in the anterior insula, anterior cingulate cortex and sensory motor cortices. Remifentanil-induced reduction in breathlessness unpleasantness was associated with increased activity in the rostral anterior cingulate cortex and nucleus accumbens, components of the endogenous opioid system known to decrease the perception of aversive stimuli.

These findings suggest that in addition to effects on brainstem respiratory control, opioids palliate breathlessness through an interplay of altered associative learning mechanisms. These mechanisms provide potential targets for novel ways to develop and assess treatments for chronic breathlessness.

•

The mechanisms of how low-dose opioids relieve breathlessness are unknown.

•

We tested whether low-dose opioids affect conditioned anticipation and perception of breathlessness.

•

Low-dose opioids reduced unpleasantness, but not intensity of breathlessness.

•

Reduced breathlessness unpleasantness was associated with activation of the endogenous opioid system.

•

Breathlessness relief was predicted by decreased anticipatory activity in amygdala/hippocampus.

Fusing multiple neuroimaging modalities to assess group differences in perception-action coupling

In the last few decades, non-invasive neuroimaging has revealed macro-scale brain dynamics that underlie perception, cognition and action. Advances in non-invasive neuroimaging target two capabilities; 1) increased spatial and temporal resolution of measured neural activity, and 2) innovative methodologies to extract brain-behavior relationships from evolving neuroimaging technology. We target the second. Our novel methodology integrated three neuroimaging methodologies and elucidated expertise-dependent differences in functional (fused EEG-fMRI) and structural (dMRI) brain networks for a perception-action coupling task. A set of baseball players and controls performed a Go/No-Go task designed to mimic the situation of hitting a baseball. In the functional analysis, our novel fusion methodology identifies 50ms windows with predictive EEG neural correlates of expertise and fuses these temporal windows with fMRI activity in a whole-brain 2mm voxel analysis, revealing time-localized correlations of expertise at a spatial scale of millimeters. The spatiotemporal cascade of brain activity reflecting expertise differences begins as early as 200ms after the pitch starts and lasting up to 700ms afterwards. Network differences are spatially localized to include motor and visual processing areas, providing evidence for differences in perception-action coupling between the groups. Furthermore, an analysis of structural connectivity revealed that the players have significantly more connections between cerebellar and left frontal/motor regions, and many of the functional activation differences between the groups are located within structurally defined network modules that differentiate expertise. In short, our novel method illustrates how multimodal neuroimaging can provide specific macro-scale insights into the functional and structural correlates of expertise development.

Resting-State Functional Connectivity in the Infant Brain: Methods, Pitfalls, and Potentiality

Early brain development is characterized by rapid growth and perpetual reconfiguration, driven by a dynamic milieu of heterogeneous processes. Postnatal brain plasticity is associated with increased vulnerability to environmental stimuli. However, little is known regarding the ontogeny and temporal manifestations of inter- and intra-regional functional connectivity that comprise functional brain networks. Resting-state functional magnetic resonance imaging (rs-fMRI) has emerged as a promising non-invasive neuroinvestigative tool, measuring spontaneous fluctuations in blood oxygen level dependent (BOLD) signal at rest that reflect baseline neuronal activity. Over the past decade, its application has expanded to infant populations providing unprecedented insight into functional organization of the developing brain, as well as early biomarkers of abnormal states. However, many methodological issues of rs-fMRI analysis need to be resolved prior to standardization of the technique to infant populations. As a primary goal, this methodological manuscript will (1) present a robust methodological protocol to extract and assess resting-state networks in early infancy using independent component analysis (ICA), such that investigators without previous knowledge in the field can implement the analysis and reliably obtain viable results consistent with previous literature; (2) review the current methodological challenges and ethical considerations associated with emerging field of infant rs-fMRI analysis; and (3) discuss the significance of rs-fMRI application in infants for future investigations of neurodevelopment in the context of early life stressors and pathological processes. The overarching goal is to catalyze efforts toward development of robust, infant-specific acquisition, and preprocessing pipelines, as well as promote greater transparency by researchers regarding methods used.

Hand classification of fMRI ICA noise components

We present a practical "how-to" guide to help determine whether single-subject fMRI independent components (ICs) characterise structured noise or not. Manual identification of signal and noise after ICA decomposition is required for efficient data denoising: to train supervised algorithms, to check the results of unsupervised ones or to manually clean the data. In this paper we describe the main spatial and temporal features of ICs and provide general guidelines on how to evaluate these. Examples of signal and noise components are provided from a wide range of datasets (3T data, including examples from the UK Biobank and the Human Connectome Project, and 7T data), together with practical guidelines for their identification. Finally, we discuss how the data quality, data type and preprocessing can influence the characteristics of the ICs and present examples of particularly challenging datasets.

Monitoring cardiac and respiratory physiology during FMRI

This article will consider how physiological monitoring can be used both as an intrinsic part of an experiment, or for removing unwanted physiological signals from the FMRI time series. As functional MRI is used for a wide variety of applications beyond the identification of regions involved in a task, different sources of noise in the time series become important. The use of arterial spin labelling sequences, either in isolation or combined with BOLD imaging, means that temporal noise must be dealt with differently. Moreover, when these are combined with global cerebrovascular stimuli, such as respiratory challenges, the standard analysis tools must be employed with great care so as not to detrimentally distort the data. Acquiring and analysing physiological data is sometimes more art than science, and this article attempts to provide some insight into common techniques as well as advice on identifying and correcting some of the problems that may be encountered.

Cortical networks involved in visual awareness independent of visual attention

It is now well established that visual attention, as measured with standard spatial attention tasks, and visual awareness, as measured by report, can be dissociated. It is possible to attend to a stimulus with no reported awareness of the stimulus. We used a behavioral paradigm in which people were aware of a stimulus in one condition and unaware of it in another condition, but the stimulus drew a similar amount of spatial attention in both conditions. The paradigm allowed us to test for brain regions active in association with awareness independent of level of attention. Participants performed the task in an MRI scanner. We looked for brain regions that were more active in the aware than the unaware trials. The largest cluster of activity was obtained in the temporoparietal junction (TPJ) bilaterally. Local independent component analysis (ICA) revealed that this activity contained three distinct, but overlapping, components: a bilateral, anterior component; a left dorsal component; and a right dorsal component. These components had brain-wide functional connectivity that partially overlapped the ventral attention network and the frontoparietal control network. In contrast, no significant activity in association with awareness was found in the banks of the intraparietal sulcus, a region connected to the dorsal attention network and traditionally associated with attention control. These results show the importance of separating awareness and attention when testing for cortical substrates. They are also consistent with a recent proposal that awareness is associated with ventral attention areas, especially in the TPJ.

Brain connectivity aberrations in anabolic-androgenic steroid users

Sustained anabolic-androgenic steroid (AAS) use has adverse behavioral consequences, including aggression, violence and impulsivity. Candidate mechanisms include disruptions of brain networks with high concentrations of androgen receptors and critically involved in emotional and cognitive regulation. Here, we tested the effects of AAS on resting-state functional brain connectivity in the largest sample of AAS-users to date. We collected resting-state functional magnetic resonance imaging (fMRI) data from 151 males engaged in heavy resistance strength training. 50 users tested positive for AAS based on the testosterone to epitestosterone (T/E) ratio and doping substances in urine. 16 previous users and 59 controls tested negative. We estimated brain network nodes and their time-series using ICA and dual regression and defined connectivity matrices as the between-node partial correlations. In line with the emotional and behavioral consequences of AAS, current users exhibited reduced functional connectivity between key nodes involved in emotional and cognitive regulation, in particular reduced connectivity between the amygdala and default-mode network (DMN) and between the dorsal attention network (DAN) and a frontal node encompassing the superior and inferior frontal gyri (SFG/IFG) and the anterior cingulate cortex (ACC), with further reductions as a function of dependency, lifetime exposure, and cycle state (on/off).

•

Sustained AAS use has adverse behavioral consequences, including aggression, violence and impulsivity.

•

We obtained r-fMRI data from 50 male users testing positive for AAS and 16 previous users and 59 controls testing negative.

•

We used ICA and dual regression, and defined connectivity matrices as the between-node temporal partial correlations.

•

Current users showed significantly reduced connectivity between amygdala and DMN and between DAN and a SFG/IFG/ACC node.

Reduction of Motion Artifacts and Noise Using Independent Component Analysis in Task-Based Functional MRI for Preoperative Planning in Patients with Brain Tumor

BACKGROUND AND PURPOSE

Although it is a potentially powerful presurgical tool, fMRI can be fraught with artifacts, leading to interpretive errors, many of which are not fully accounted for in routinely applied correction methods. The purpose of this investigation was to evaluate the effects of data denoising by independent component analysis in patients undergoing preoperative evaluation for glioma resection compared with more routinely applied correction methods such as realignment or motion scrubbing.

MATERIALS AND METHODS

Thirty-five functional runs (both motor and language) in 12 consecutive patients with glioma were analyzed retrospectively by double-blind review. Data were processed and compared with the following: 1) realignment alone, 2) motion scrubbing, 3) independent component analysis denoising, and 4) both independent component analysis denoising and motion scrubbing. Primary outcome measures included a change in false-positives, false-negatives, z score, and diagnostic rating.

RESULTS

Independent component analysis denoising reduced false-positives in 63% of studies versus realignment alone. There was also an increase in the z score in areas of true activation in 71.4% of studies. Areas of new expected activation (previous false-negatives) were revealed in 34.4% of cases with independent component analysis denoising versus motion scrubbing or realignment alone. Of studies deemed nondiagnostic with realignment or motion scrubbing alone, 65% were considered diagnostic after independent component analysis denoising.

CONCLUSIONS

The addition of independent component analysis denoising of fMRI data in preoperative patients with glioma has a significant impact on data quality, resulting in reduced false-positives and an increase in true-positives compared with more commonly applied motion scrubbing or simple realignment methods.

Variable disruption of a syntactic processing network in primary progressive aphasia

Syntactic processing deficits are highly variable in individuals with primary progressive aphasia. Damage to left inferior frontal cortex has been associated with syntactic deficits in primary progressive aphasia in a number of structural and functional neuroimaging studies. However, a contrasting picture of a broader syntactic network has emerged from neuropsychological studies in other aphasic cohorts, and functional imaging studies in healthy controls. To reconcile these findings, we used functional magnetic resonance imaging to investigate the functional neuroanatomy of syntactic comprehension in 51 individuals with primary progressive aphasia, composed of all clinical variants and a range of degrees of syntactic processing impairment. We used trial-by-trial reaction time as a proxy for syntactic processing load, to determine which regions were modulated by syntactic processing in each patient, and how the set of regions recruited was related to whether syntactic processing was ultimately successful or unsuccessful. Relationships between functional abnormalities and patterns of cortical atrophy were also investigated. We found that the individual degree of syntactic comprehension impairment was predicted by left frontal atrophy, but also by functional disruption of a broader syntactic processing network, comprising left posterior frontal cortex, left posterior temporal cortex, and the left intraparietal sulcus and adjacent regions. These regions were modulated by syntactic processing in healthy controls and in patients with primary progressive aphasia with relatively spared syntax, but they were modulated to a lesser extent or not at all in primary progressive aphasia patients whose syntax was relatively impaired. Our findings suggest that syntactic comprehension deficits in primary progressive aphasia reflect not only structural and functional changes in left frontal cortex, but also disruption of a wider syntactic processing network.

Tuning to Binaural Cues in Human Auditory Cortex

Interaural level and time differences (ILD and ITD), the primary binaural cues for sound localization in azimuth, are known to modulate the tuned responses of neurons in mammalian auditory cortex (AC). The majority of these neurons respond best to cue values that favor the contralateral ear, such that contralateral bias is evident in the overall population response and thereby expected in population-level functional imaging data. Human neuroimaging studies, however, have not consistently found contralaterally biased binaural response patterns. Here, we used functional magnetic resonance imaging (fMRI) to parametrically measure ILD and ITD tuning in human AC. For ILD, contralateral tuning was observed, using both univariate and multivoxel analyses, in posterior superior temporal gyrus (pSTG) in both hemispheres. Response-ILD functions were U-shaped, revealing responsiveness to both contralateral and—to a lesser degree—ipsilateral ILD values, consistent with rate coding by unequal populations of contralaterally and ipsilaterally tuned neurons. In contrast, for ITD, univariate analyses showed modest contralateral tuning only in left pSTG, characterized by a monotonic response-ITD function. A multivoxel classifier, however, revealed ITD coding in both hemispheres. Although sensitivity to ILD and ITD was distributed in similar AC regions, the differently shaped response functions and different response patterns across hemispheres suggest that basic ILD and ITD processes are not fully integrated in human AC. The results support opponent-channel theories of ILD but not necessarily ITD coding, the latter of which may involve multiple types of representation that differ across hemispheres.

Interhemispheric Functional Brain Connectivity in Neonates with Prenatal Alcohol Exposure: Preliminary Findings

BACKGROUND

Children exposed to alcohol in utero demonstrate reduced white matter microstructural integrity. While early evidence suggests altered functional brain connectivity in the lateralization of motor networks in school-age children with prenatal alcohol exposure (PAE), the specific effects of alcohol exposure on the establishment of intrinsic connectivity in early infancy have not been explored.

METHODS

Sixty subjects received functional imaging at 2 to 4 weeks of age for 6 to 8 minutes during quiet natural sleep. Thirteen alcohol-exposed (PAE) and 14 age-matched control (CTRL) participants with usable data were included in a multivariate model of connectivity between sensorimotor intrinsic functional connectivity networks. Seed-based analyses of group differences in interhemispheric connectivity of intrinsic motor networks were also conducted. The Dubowitz neurological assessment was performed at the imaging visit.

RESULTS

Alcohol exposure was associated with significant increases in connectivity between somatosensory, motor networks, brainstem/thalamic, and striatal intrinsic networks. Reductions in interhemispheric connectivity of motor and somatosensory networks did not reach significance.

CONCLUSIONS

Although results are preliminary, findings suggest PAE may disrupt the temporal coherence in blood oxygenation utilization in intrinsic networks underlying motor performance in newborn infants. Studies that employ longitudinal designs to investigate the effects of in utero alcohol exposure on the evolving resting-state networks will be key in establishing the distribution and timing of connectivity disturbances already described in older children.

Denoising spinal cord fMRI data: Approaches to acquisition and analysis

Functional magnetic resonance imaging (fMRI) of the human spinal cord is a difficult endeavour due to the cord's small cross-sectional diameter, signal drop-out as well as image distortion due to magnetic field inhomogeneity, and the confounding influence of physiological noise from cardiac and respiratory sources. Nevertheless, there is great interest in spinal fMRI due to the spinal cord's role as the principal sensorimotor interface between the brain and the body and its involvement in a variety of sensory and motor pathologies. In this review, we give an overview of the various methods that have been used to address the technical challenges in spinal fMRI, with a focus on reducing the impact of physiological noise. We start out by describing acquisition methods that have been tailored to the special needs of spinal fMRI and aim to increase the signal-to-noise ratio and reduce distortion in obtained images. Following this, we concentrate on image processing and analysis approaches that address the detrimental effects of noise. While these include variations of standard pre-processing methods such as motion correction and spatial filtering, the main focus lies on denoising techniques that can be applied to task-based as well as resting-state data sets. We review both model-based approaches that rely on externally acquired respiratory and cardiac signals as well as data-driven approaches that estimate and correct for noise using the data themselves. We conclude with an outlook on techniques that have been successfully applied for noise reduction in brain imaging and whose use might be beneficial for fMRI of the human spinal cord.

Bout-associated intrinsic functional network changes in cluster headache: A longitudinal resting-state functional MRI study

Background Previous imaging studies on the pathogenesis of cluster headache (CH) have implicated the hypothalamus and multiple brain networks. However, very little is known regarding dynamic bout-associated, large-scale resting state functional network changes related to CH. Methods Resting-state functional magnetic resonance imaging data were obtained from CH patients and matched controls. Data were analyzed using independent component analysis for exploratory assessment of the changes in intrinsic brain networks and their relationship between in-bout and out-of-bout periods, as well as correlations with clinical observations. Results Compared to healthy controls, CH patients had functional connectivity (FC) changes in the temporal, frontal, salience, default mode, somatosensory, dorsal attention, and visual networks, independent of bout period. Compared to out-of-bout scans, in-bout scans showed altered FC in the frontal and dorsal attention networks. Lower frontal network FC correlated with longer duration of CH. Conclusions The present findings suggest that episodic CH with dynamic bout period shifts may involve bout-associated FC changes in multiple discrete cortical areas within networks outside traditional pain processing areas. Dynamic changes in FC in frontal and dorsal attention networks between bout periods could be important for understanding episodic CH pathophysiology.

Noise contributions to the fMRI signal: An overview

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

Brain Responses during the Anticipation of Dyspnea

Dyspnea is common in many cardiorespiratory diseases. Already the anticipation of this aversive symptom elicits fear in many patients resulting in unfavorable health behaviors such as activity avoidance and sedentary lifestyle. This study investigated brain mechanisms underlying these anticipatory processes. We induced dyspnea using resistive-load breathing in healthy subjects during functional magnetic resonance imaging. Blocks of severe and mild dyspnea alternated, each preceded by anticipation periods. Severe dyspnea activated a network of sensorimotor, cerebellar, and limbic areas. The left insular, parietal opercular, and cerebellar cortices showed increased activation already during dyspnea anticipation. Left insular and parietal opercular cortex showed increased connectivity with right insular and anterior cingulate cortex when severe dyspnea was anticipated, while the cerebellum showed increased connectivity with the amygdala. Notably, insular activation during dyspnea perception was positively correlated with midbrain activation during anticipation. Moreover, anticipatory fear was positively correlated with anticipatory activation in right insular and anterior cingulate cortex. The results demonstrate that dyspnea anticipation activates brain areas involved in dyspnea perception. The involvement of emotion-related areas such as insula, anterior cingulate cortex, and amygdala during dyspnea anticipation most likely reflects anticipatory fear and might underlie the development of unfavorable health behaviors in patients suffering from dyspnea.

Brain Signal Variability Differentially Affects Cognitive Flexibility and Cognitive Stability

Recent research yielded the intriguing conclusion that, in healthy adults, higher levels of variability in neuronal processes are beneficial for cognitive functioning. Beneficial effects of variability in neuronal processing can also be inferred from neurocomputational theories of working memory, albeit this holds only for tasks requiring cognitive flexibility. However, cognitive stability, i.e., the ability to maintain a task goal in the face of irrelevant distractors, should suffer under high levels of brain signal variability. To directly test this prediction, we studied both behavioral and brain signal variability during cognitive flexibility (i.e., task switching) and cognitive stability (i.e., distractor inhibition) in a sample of healthy human subjects and developed an efficient and easy-to-implement analysis approach to assess BOLD-signal variability in event-related fMRI task paradigms. Results show a general positive effect of neural variability on task performance as assessed by accuracy measures. However, higher levels of BOLD-signal variability in the left inferior frontal junction area result in reduced error rate costs during task switching and thus facilitate cognitive flexibility. In contrast, variability in the same area has a detrimental effect on cognitive stability, as shown in a negative effect of variability on response time costs during distractor inhibition. This pattern was mirrored at the behavioral level, with higher behavioral variability predicting better task switching but worse distractor inhibition performance. Our data extend previous results on brain signal variability by showing a differential effect of brain signal variability that depends on task context, in line with predictions from computational theories.

SIGNIFICANCE STATEMENT

Recent neuroscientific research showed that the human brain signal is intrinsically variable and suggested that this variability improves performance. Computational models of prefrontal neural networks predict differential effects of variability for different behavioral situations requiring either cognitive flexibility or stability. However, this hypothesis has so far not been put to an empirical test. In this study, we assessed cognitive flexibility and cognitive stability, and, besides a generally positive effect of neural variability on accuracy measures, we show that neural variability in a prefrontal brain area at the inferior frontal junction is differentially associated with performance: higher levels of variability are beneficial for the effectiveness of task switching (cognitive flexibility) but detrimental for the efficiency of distractor inhibition (cognitive stability).

Bootstrap Enhanced Penalized Regression for Variable Selection with Neuroimaging Data

Recent advances in fMRI research highlight the use of multivariate methods for examining whole-brain connectivity. Complementary data-driven methods are needed for determining the subset of predictors related to individual differences. Although commonly used for this purpose, ordinary least squares (OLS) regression may not be ideal due to multi-collinearity and over-fitting issues. Penalized regression is a promising and underutilized alternative to OLS regression. In this paper, we propose a nonparametric bootstrap quantile (QNT) approach for variable selection with neuroimaging data. We use real and simulated data, as well as annotated R code, to demonstrate the benefits of our proposed method. Our results illustrate the practical potential of our proposed bootstrap QNT approach. Our real data example demonstrates how our method can be used to relate individual differences in neural network connectivity with an externalizing personality measure. Also, our simulation results reveal that the QNT method is effective under a variety of data conditions. Penalized regression yields more stable estimates and sparser models than OLS regression in situations with large numbers of highly correlated neural predictors. Our results demonstrate that penalized regression is a promising method for examining associations between neural predictors and clinically relevant traits or behaviors. These findings have important implications for the growing field of functional connectivity research, where multivariate methods produce numerous, highly correlated brain networks.

Temporal Non-Local Means Filtering Reveals Real-Time Whole-Brain Cortical Interactions in Resting fMRI

Intensity variations over time in resting BOLD fMRI exhibit spatial correlation patterns consistent with a set of large scale cortical networks. However, visualizations of this data on the brain surface, even after extensive preprocessing, are dominated by local intensity fluctuations that obscure larger scale behavior. Our novel adaptation of non-local means (NLM) filtering, which we refer to as temporal NLM or tNLM, reduces these local fluctuations without the spatial blurring that occurs when using standard linear filtering methods. We show examples of tNLM filtering that allow direct visualization of spatio-temporal behavior on the cortical surface. These results reveal patterns of activity consistent with known networks as well as more complex dynamic changes within and between these networks. This ability to directly visualize brain activity may facilitate new insights into spontaneous brain dynamics. Further, temporal NLM can also be used as a preprocessor for resting fMRI for exploration of dynamic brain networks. We demonstrate its utility through application to graph-based functional cortical parcellation. Simulations with known ground truth functional regions demonstrate that tNLM filtering prior to parcellation avoids the formation of false parcels that can arise when using linear filtering. Application to resting fMRI data from the Human Connectome Project shows significant improvement, in comparison to linear filtering, in quantitative agreement with functional regions identified independently using task-based experiments as well as in test-retest reliability.

Resting state functional MRI reveals abnormal network connectivity in orthostatic tremor

Very little is known about the pathogenesis of orthostatic tremor (OT). We have observed that OT patients might have deficits in specific aspects of neuropsychological function, particularly those thought to rely on the integrity of the prefrontal cortex, which suggests a possible involvement of frontocerebellar circuits. We examined whether resting-state functional magnetic resonance imaging (fMRI) might provide further insights into the pathogenesis on OT. Resting-state fMRI data in 13 OT patients (11 women and 2 men) and 13 matched healthy controls were analyzed using independent component analysis, in combination with a "dual-regression" technique, to identify group differences in several resting-state networks (RSNs). All participants also underwent neuropsychological testing during the same session. Relative to healthy controls, OT patients showed increased connectivity in RSNs involved in cognitive processes (default mode network [DMN] and frontoparietal networks), and decreased connectivity in the cerebellum and sensorimotor networks. Changes in network integrity were associated not only with duration (DMN and medial visual network), but also with cognitive function. Moreover, in at least 2 networks (DMN and medial visual network), increased connectivity was associated with worse performance on different cognitive domains (attention, executive function, visuospatial ability, visual memory, and language). In this exploratory study, we observed selective impairments of RSNs in OT patients. This and other future resting-state fMRI studies might provide a novel method to understand the pathophysiological mechanisms of motor and nonmotor features of OT.

Anticipation-related brain connectivity in bipolar and unipolar depression: a graph theory approach

Bipolar disorder is often misdiagnosed as major depressive disorder, which leads to inadequate treatment. Depressed individuals versus healthy control subjects, show increased expectation of negative outcomes. Due to increased impulsivity and risk for mania, however, depressed individuals with bipolar disorder may differ from those with major depressive disorder in neural mechanisms underlying anticipation processes. Graph theory methods for neuroimaging data analysis allow the identification of connectivity between multiple brain regions without prior model specification, and may help to identify neurobiological markers differentiating these disorders, thereby facilitating development of better therapeutic interventions. This study aimed to compare brain connectivity among regions involved in win/loss anticipation in depressed individuals with bipolar disorder (BDD) versus depressed individuals with major depressive disorder (MDD) versus healthy control subjects using graph theory methods. The study was conducted at the University of Pittsburgh Medical Center and included 31 BDD, 39 MDD, and 36 healthy control subjects. Participants were scanned while performing a number guessing reward task that included the periods of win and loss anticipation. We first identified the anticipatory network across all 106 participants by contrasting brain activation during all anticipation periods (win anticipation + loss anticipation) versus baseline, and win anticipation versus loss anticipation. Brain connectivity within the identified network was determined using the Independent Multiple sample Greedy Equivalence Search (IMaGES) and Linear non-Gaussian Orientation, Fixed Structure (LOFS) algorithms. Density of connections (the number of connections in the network), path length, and the global connectivity direction ('top-down' versus 'bottom-up') were compared across groups (BDD/MDD/healthy control subjects) and conditions (win/loss anticipation). These analyses showed that loss anticipation was characterized by denser top-down fronto-striatal and fronto-parietal connectivity in healthy control subjects, by bottom-up striatal-frontal connectivity in MDD, and by sparse connectivity lacking fronto-striatal connections in BDD. Win anticipation was characterized by dense connectivity of medial frontal with striatal and lateral frontal cortical regions in BDD, by sparser bottom-up striatum-medial frontal cortex connectivity in MDD, and by sparse connectivity in healthy control subjects. In summary, this is the first study to demonstrate that BDD and MDD with comparable levels of current depression differed from each other and healthy control subjects in density of connections, connectivity path length, and connectivity direction as a function of win or loss anticipation. These findings suggest that different neurobiological mechanisms may underlie aberrant anticipation processes in BDD and MDD, and that distinct therapeutic strategies may be required for these individuals to improve coping strategies during expectation of positive and negative outcomes.

Deficits in Visual System Functional Connectivity after Blast-Related Mild TBI are Associated with Injury Severity and Executive Dysfunction

INTRODUCTION

Approximately, 275,000 American service members deployed to Iraq or Afghanistan have sustained a mild traumatic brain injury (mTBI), with 75% of these incidents involving an explosive blast. Visual processing problems and cognitive dysfunction are common complaints following blast-related mTBI.

METHODS

In 127 veterans, we examined resting fMRI functional connectivity (FC) of four key nodes within the visual system: lateral geniculate nucleus (LGN), primary visual cortex (V1), lateral occipital gyrus (LO), and fusiform gyrus (FG). Regression analyses were performed (i) to obtain correlations between time-series from each seed and all voxels in the brain, and (ii) to identify brain regions in which FC variability was related to blast mTBI severity. Blast-related mTBI severity was quantified as the sum of the severity scores assigned to each of the three most significant blast-related injuries self-reported by subjects. Correlations between FC and performance on executive functioning tasks were performed across participants with available behavioral data (n = 94).

RESULTS

Greater blast mTBI severity scores were associated with lower FC between: (A) LGN seed and (i) medial frontal gyrus, (ii) lingual gyrus, and (iii) right ventral anterior nucleus of thalamus; (B) V1 seed and precuneus; (C) LO seed and middle and superior frontal gyri; (D) FG seed and (i) superior and medial frontal gyrus, and (ii) left middle frontal gyrus. Finally, lower FC between visual network regions and frontal cortical regions predicted worse performance on the WAIS digit-symbol coding task.

CONCLUSION

These are the first published results that directly illustrate the relationship between blast-related mTBI severity, visual pathway neural networks, and executive dysfunction - results that highlight the detrimental relationship between blast-related brain injury and the integration of visual sensory input and executive processes.

Building a Science of Individual Differences from fMRI

To date, fMRI research has been concerned primarily with evincing generic principles of brain function through averaging data from multiple subjects. Given rapid developments in both hardware and analysis tools, the field is now poised to study fMRI-derived measures in individual subjects, and to relate these to psychological traits or genetic variations. We discuss issues of validity, reliability and statistical assessment that arise when the focus shifts to individual subjects and that are applicable also to other imaging modalities. We emphasize that individual assessment of neural function with fMRI presents specific challenges and necessitates careful consideration of anatomical and vascular between-subject variability as well as sources of within-subject variability.

Functional connectivity associated with social networks in older adults: A resting-state fMRI study

Poor social networks and decreased levels of social support are associated with worse mood, health, and cognition in younger and older adults. Yet, we know very little about the brain substrates associated with social networks and social support, particularly in older adults. This study examined functional brain substrates associated with social networks using the Social Network Index (SNI) and resting-state functional magnetic resonance imaging (fMRI). Resting-state fMRI data from 28 non-demented older adults were analyzed with independent components analyses. As expected, four established resting-state networks-previously linked to motor, vision, speech, and other language functions-correlated with the quality (SNI-1: total number of high-contact roles of a respondent) and quantity (SNI-2: total number of individuals in a respondent's social network) of social networks: a sensorimotor, a visual, a vestibular/insular, and a left frontoparietal network. Moreover, SNI-1 was associated with greater functional connectivity in the lateral prefrontal regions of the left frontoparietal network, while SNI-2 was associated with greater functional connectivity in the medial prefrontal regions of this network. Thus, lateral prefrontal regions may be particularly linked to the quality of social networks while medial prefrontal regions may be particularly linked to the quantity of social networks.

Is Traumatic Brain Injury Associated with Reduced Inter-Hemispheric Functional Connectivity? A Study of Large-Scale Resting State Networks following Traumatic Brain Injury

Traumatic brain injury (TBI) often has long-term debilitating sequelae in cognitive and behavioral domains. Understanding how TBI impacts functional integrity of brain networks that underlie these domains is key to guiding future approaches to TBI rehabilitation. In the current study, we investigated the differences in inter-hemispheric functional connectivity (FC) of resting state networks (RSNs) between chronic mild-to-severe TBI patients and normal comparisons (NC), focusing on two externally oriented networks (i.e., the fronto-parietal network [FPN] and the executive control network [ECN]), one internally oriented network (i.e., the default mode network [DMN]), and one somato-motor network (SMN). Seed voxel correlation analysis revealed that TBI patients displayed significantly less FC between lateralized seeds and both homologous and non-homologous regions in the opposite hemisphere for externally oriented networks but not for DMN or SMN; conversely, TBI patients showed increased FC within regions of the DMN, especially precuneus and parahippocampal gyrus. Region of interest correlation analyses confirmed the presence of significantly higher inter-hemispheric FC in NC for the FPN (p < 0.01), and ECN (p < 0.05), but not for the DMN (p > 0.05) or SMN (p > 0.05). Further analysis revealed that performance on a neuropsychological test measuring organizational skills and visuo-spatial abilities administered to the TBI group, the Rey-Osterrieth Complex Figure Test, positively correlated with FC between the right FPN and homologous regions. Our findings suggest that distinct RSNs display specific patterns of aberrant FC following TBI; this represents a step forward in the search for biomarkers useful for early diagnosis and treatment of TBI-related cognitive impairment.

Very Early Brain Damage Leads to Remodeling of the Working Memory System in Adulthood: A Combined fMRI/Tractography Study

The human brain can adapt to overcome injury even years after an initial insult. One hypothesis states that early brain injury survivors, by taking advantage of critical periods of high plasticity during childhood, should recover more successfully than those who suffer injury later in life. This hypothesis has been challenged by recent studies showing worse cognitive outcome in individuals with early brain injury, compared with individuals with later brain injury, with working memory particularly affected. We invited individuals who suffered perinatal brain injury (PBI) for an fMRI/diffusion MRI tractography study of working memory and hypothesized that, 30 years after the initial injury, working memory deficits in the PBI group would remain, despite compensatory activation in areas outside the typical working memory network. Furthermore we hypothesized that the amount of functional reorganization would be related to the level of injury to the dorsal cingulum tract, which connects medial frontal and parietal working memory structures. We found that adults who suffered PBI did not significantly differ from controls in working memory performance. They exhibited less activation in classic frontoparietal working memory areas and a relative overactivation of bilateral perisylvian cortex compared with controls. Structurally, the dorsal cingulum volume and hindrance-modulated orientational anisotropy was significantly reduced in the PBI group. Furthermore there was uniquely in the PBI group a significant negative correlation between the volume of this tract and activation in the bilateral perisylvian cortex and a positive correlation between this activation and task performance. This provides the first evidence of compensatory plasticity of the working memory network following PBI.

Validity and reliability of four language mapping paradigms

Language areas of the brain can be mapped in individual participants with functional MRI. We investigated the validity and reliability of four language mapping paradigms that may be appropriate for individuals with acquired aphasia: sentence completion, picture naming, naturalistic comprehension, and narrative comprehension. Five neurologically normal older adults were scanned on each of the four paradigms on four separate occasions. Validity was assessed in terms of whether activation patterns reflected the known typical organization of language regions, that is, lateralization to the left hemisphere, and involvement of the left inferior frontal gyrus and the left middle and/or superior temporal gyri. Reliability (test-retest reproducibility) was quantified in terms of the Dice coefficient of similarity, which measures overlap of activations across time points. We explored the impact of different absolute and relative voxelwise thresholds, a range of cluster size cutoffs, and limitation of analyses to a priori potential language regions. We found that the narrative comprehension and sentence completion paradigms offered the best balance of validity and reliability. However, even with optimal combinations of analysis parameters, there were many scans on which known features of typical language organization were not demonstrated, and test-retest reproducibility was only moderate for realistic parameter choices. These limitations in terms of validity and reliability may constitute significant limitations for many clinical or research applications that depend on identifying language regions in individual participants.

•

Validity and reliability were investigated for four language mapping paradigms.

•

Narrative comprehension and sentence completion paradigms performed best.

•

Lateralization to the left hemisphere was not always apparent.

•

Test-retest reproducibility was only moderate.

Conditioned respiratory threat in the subdivisions of the human periaqueductal gray

The sensation of breathlessness is the most threatening symptom of respiratory disease. The different subdivisions of the midbrain periaqueductal gray (PAG) are intricately (and differentially) involved in integrating behavioural responses to threat in animals, while the PAG has previously only been considered as a single entity in human research. Here we investigate how these individual PAG columns are differently involved with respiratory threat. Eighteen healthy subjects were conditioned to associate shapes with certain or uncertain impending respiratory load, and scanned the following day during anticipation and application of inspiratory loading using 7 T functional MRI. We showed activity in the ventrolateral PAG (vlPAG) during anticipation of resistive loading, with activity in the lateral PAG (lPAG) during resistive loading, revealing spatially and temporally distinct functions within this structure. We propose that lPAG is involved with sensorimotor responses to breathlessness, while the vlPAG operates within the threat perception network for impending breathlessness.

DOI:http://dx.doi.org/10.7554/eLife.12047.001

Many people find feeling breathless one of the most upsetting symptoms of respiratory diseases, and breathlessness often causes anxiety that makes the condition seem more threatening than it is. Studies in animals suggest that a small cluster of neurons called the periaqueductal gray is important for responding to threats. This cluster, located at the top of the brainstem, is divided into parallel columns running from top to bottom. In animals, these columns are known to have distinct roles, but human research has tended to consider the periaqueductal gray as a single, uniform entity.

Faull et al. wanted to find out whether different columns of the human periaqueductal gray have distinct roles in the perception of respiratory threat. During the study, participants breathed through a tube while watching shapes appear on a screen. This tube could be altered to make breathing more or less difficult – much like breathing through a narrow drinking straw.

A conditioning session was first conducted so that participants learned that certain shapes on the screen signalled that their breathing was about to become difficult, while other shapes signalled normal breathing. A second session was then conducted in a brain scanner, using a technique called functional magnetic resonance imaging. This allowed Faull et al. to compare brain activity during the anticipation of difficult breathing with the brain activity during the breathing challenge itself.

The results show that the column at the front of the periaqueductal gray (the ventrolateral column) was more active when participants saw the shape that signaled upcoming breathing difficulty. In contrast, difficult breathing was associated with activity in the lateral column (at the side of the periaqueductal gray).

Thus, the different columns of the human periaqueductal gray have different roles in the response to respiratory threat. Future studies could investigate how these columns interact with each other and with other brain regions. Such understanding is important for a range of conditions that may be influenced by the activity of the periaqueductal gray, including disruptions in bladder control, hypertension, chronic pain, and asthma.

DOI:http://dx.doi.org/10.7554/eLife.12047.002

Salience and Default Mode Network Coupling Predicts Cognition in Aging and Parkinson's Disease

OBJECTIVES

Cognitive impairment is common in Parkinson's disease (PD). Three neurocognitive networks support efficient cognition: the salience network, the default mode network, and the central executive network. The salience network is thought to switch between activating and deactivating the default mode and central executive networks. Anti-correlated interactions between the salience and default mode networks in particular are necessary for efficient cognition. Our previous work demonstrated altered functional coupling between the neurocognitive networks in non-demented individuals with PD compared to age-matched control participants. Here, we aim to identify associations between cognition and functional coupling between these neurocognitive networks in the same group of participants.

METHODS

We investigated the extent to which intrinsic functional coupling among these neurocognitive networks is related to cognitive performance across three neuropsychological domains: executive functioning, psychomotor speed, and verbal memory. Twenty-four non-demented individuals with mild to moderate PD and 20 control participants were scanned at rest and evaluated on three neuropsychological domains.

RESULTS

PD participants were impaired on tests from all three domains compared to control participants. Our imaging results demonstrated that successful cognition across healthy aging and Parkinson's disease participants was related to anti-correlated coupling between the salience and default mode networks. Individuals with poorer performance scores across groups demonstrated more positive salience network/default-mode network coupling.

CONCLUSIONS

Successful cognition relies on healthy coupling between the salience and default mode networks, which may become dysfunctional in PD. These results can help inform non-pharmacological interventions (repetitive transcranial magnetic stimulation) targeting these specific networks before they become vulnerable in early stages of Parkinson's disease.