Predicting Language Function Post-Stroke: A Model-Based Structural Connectivity Approach.

BACKGROUND: The prediction of post-stroke language function is essential for the development of individualized treatment plans based on the personal recovery potential of aphasic stroke patients. OBJECTIVE: To establish a framework for integrating information on connectivity disruption of the language network based on routinely collected clinical magnetic resonance (MR) images into Random Forest modeling to predict post-stroke language function. METHODS: Language function was assessed in 76 stroke patients from the Non-Invasive Repeated Therapeutic Stimulation for Aphasia Recovery trial, using the Token Test (TT), Boston Naming Test (BNT), and Semantic Verbal Fluency (sVF) Test as primary outcome measures. Individual infarct masks were superimposed onto a diffusion tensor imaging tractogram reference set to calculate Change in Connectivity scores of language-relevant gray matter regions as estimates of structural connectivity disruption. Multivariable Random Forest models were derived to predict language function. RESULTS: Random Forest models explained moderate to high amount of variance at baseline and follow-up for the TT (62.7% and 76.2%), BNT (47.0% and 84.3%), and sVF (52.2% and 61.1%). Initial language function and non-verbal cognitive ability were the most important variables to predict language function. Connectivity disruption explained additional variance, resulting in a prediction error increase of up to 12.8% with variable omission. Left middle temporal gyrus (12.8%) and supramarginal gyrus (9.8%) were identified as among the most important network nodes. CONCLUSION: Connectivity disruption of the language network adds predictive value beyond lesion volume, initial language function, and non-verbal cognitive ability. Obtaining information on connectivity disruption based on routine clinical MR images constitutes a significant advancement toward practical clinical application.

Release the Krakencoder: A unified brain connectome translation and fusion tool.

Brain connectivity can be estimated in many ways, depending on modality and processing strategy. Here we present the Krakencoder, a joint connectome mapping tool that simultaneously, bidirectionally translates between structural (SC) and functional connectivity (FC), and across different atlases and processing choices via a common latent representation. These mappings demonstrate unprecedented accuracy and individual-level identifiability; the mapping between SC and FC has identifiability 42-54% higher than existing models. The Krakencoder combines all connectome flavors via a shared low-dimensional latent space. This "fusion" representation i) better reflects familial relatedness, ii) preserves age- and sex-relevant information and iii) enhances cognition-relevant information. The Krakencoder can be applied without retraining to new, out-of-age-distribution data while still preserving inter-individual differences in the connectome predictions and familial relationships in the latent representations. The Krakencoder is a significant leap forward in capturing the relationship between multi-modal brain connectomes in an individualized, behaviorally- and demographically-relevant way.

Diffusion Tensor Imaging Along Perivascular Spaces (DTI-ALPS) to Assess Effects of Age, Sex, and Head Size on Interstitial Fluid Dynamics in Healthy Subjects.

Diffusion tensor imaging along perivascular spaces (DTI-ALPS) is a novel MRI method for assessing brain interstitial fluid dynamics, potentially indexing glymphatic function. Failed glymphatic clearance is implicated in Alzheimer's disease (AD) pathophysiology. We assessed the contribution of age and female sex (strong AD risk factors) to DTI-ALPS index in healthy subjects. We also for the first time assessed the effect of head size. In accord with prior studies, we show reduced DTI-ALPS index with aging, and in men compared to women. However, head size may be a major contributing factor to this counterintuitive sex difference.

Altered structural connectivity and functional brain dynamics in individuals with heavy alcohol use.

Heavy alcohol use and its associated conditions, such as alcohol use disorder (AUD), impact millions of individuals worldwide. While our understanding of the neurobiological correlates of AUD has evolved substantially, we still lack models incorporating whole-brain neuroanatomical, functional, and pharmacological information under one framework. Here, we utilize diffusion and functional magnetic resonance imaging to investigate alterations to brain dynamics in N = 130 individuals with a high amount of current alcohol use. We compared these alcohol using individuals to N = 308 individuals with minimal use of any substances. We find that individuals with heavy alcohol use had less dynamic and complex brain activity, and through leveraging network control theory, had increased control energy to complete transitions between activation states. Further, using separately acquired positron emission tomography (PET) data, we deploy an in silico evaluation demonstrating that decreased D2 receptor levels, as found previously in individuals with AUD, may relate to our observed findings. This work demonstrates that whole-brain, multimodal imaging information can be combined under a network control framework to identify and evaluate neurobiological correlates and mechanisms of AUD.

The sequence of regional structural disconnectivity due to multiple sclerosis lesions.

Prediction of disease progression is challenging in multiple sclerosis as the sequence of lesion development and retention of inflammation within a subset of chronic lesions is heterogeneous among patients. We investigated the sequence of lesion-related regional structural disconnectivity across the spectrum of disability and cognitive impairment in multiple sclerosis. In a full cohort of 482 multiple sclerosis patients (age: 41.83 ± 11.63 years, 71.57% females), the Expanded Disability Status Scale was used to classify patients into (i) no or mild (Expanded Disability Status Scale <3) versus (ii) moderate or severe disability groups (Expanded Disability Status Scale ≥3). In 363 out of 482 patients, quantitative susceptibility mapping was used to identify paramagnetic rim lesions, which are maintained by a rim of iron-laden innate immune cells. In 171 out of 482 patients, Brief International Cognitive Assessment was used to identify subjects as being cognitively preserved or impaired. Network Modification Tool was used to estimate the regional structural disconnectivity due to multiple sclerosis lesions. Discriminative event-based modelling was applied to investigate the sequence of regional structural disconnectivity due to (i) all representative T2 fluid-attenuated inversion recovery lesions, (ii) paramagnetic rim lesions versus non-paramagnetic rim lesions separately across disability groups ('no to mild disability' to 'moderate to severe disability'), (iii) all representative T2 fluid-attenuated inversion recovery lesions and (iv) paramagnetic rim lesions versus non-paramagnetic rim lesions separately across cognitive status ('cognitively preserved' to 'cognitively impaired'). In the full cohort, structural disconnection in the ventral attention and subcortical networks, particularly in the supramarginal and putamen regions, was an early biomarker of moderate or severe disability. The earliest biomarkers of disability progression were structural disconnections due to paramagnetic rim lesions in the motor-related regions. Subcortical structural disconnection, particularly in the ventral diencephalon and thalamus regions, was an early biomarker of cognitive impairment. Our data-driven model revealed that the structural disconnection in the subcortical regions, particularly in the thalamus, is an early biomarker for both disability and cognitive impairment in multiple sclerosis. Paramagnetic rim lesions-related structural disconnection in the motor cortex may identify the patients at risk for moderate or severe disability in multiple sclerosis. Such information might be used to identify people with multiple sclerosis who have an increased risk of disability progression or cognitive decline in order to provide personalized treatment plans.

Human brain responses are modulated when exposed to optimized natural images or synthetically generated images.

Understanding how human brains interpret and process information is important. Here, we investigated the selectivity and inter-individual differences in human brain responses to images via functional MRI. In our first experiment, we found that images predicted to achieve maximal activations using a group level encoding model evoke higher responses than images predicted to achieve average activations, and the activation gain is positively associated with the encoding model accuracy. Furthermore, anterior temporal lobe face area (aTLfaces) and fusiform body area 1 had higher activation in response to maximal synthetic images compared to maximal natural images. In our second experiment, we found that synthetic images derived using a personalized encoding model elicited higher responses compared to synthetic images from group-level or other subjects' encoding models. The finding of aTLfaces favoring synthetic images than natural images was also replicated. Our results indicate the possibility of using data-driven and generative approaches to modulate macro-scale brain region responses and probe inter-individual differences in and functional specialization of the human visual system.

Data-driven biomarkers outperform theory-based biomarkers in predicting stroke motor outcomes.

Chronic motor impairments are a leading cause of disability after stroke. Previous studies have predicted motor outcomes based on the degree of damage to predefined structures in the motor system, such as the corticospinal tract. However, such theory-based approaches may not take full advantage of the information contained in clinical imaging data. The present study uses data-driven approaches to predict chronic motor outcomes after stroke and compares the accuracy of these predictions to previously-identified theory-based biomarkers. Using a cross-validation framework, regression models were trained using lesion masks and motor outcomes data from 789 stroke patients (293 female/496 male) from the ENIGMA Stroke Recovery Working Group (age 64.9±18.0 years; time since stroke 12.2±0.2 months; normalised motor score 0.7±0.5 (range [0,1]). The out-of-sample prediction accuracy of two theory-based biomarkers was assessed: lesion load of the corticospinal tract, and lesion load of multiple descending motor tracts. These theory-based prediction accuracies were compared to the prediction accuracy from three data-driven biomarkers: lesion load of lesion-behaviour maps, lesion load of structural networks associated with lesion-behaviour maps, and measures of regional structural disconnection. In general, data-driven biomarkers had better prediction accuracy - as measured by higher explained variance in chronic motor outcomes - than theory-based biomarkers. Data-driven models of regional structural disconnection performed the best of all models tested (R2 = 0.210, p < 0.001), performing significantly better than predictions using the theory-based biomarkers of lesion load of the corticospinal tract (R2 = 0.132, p< 0.001) and of multiple descending motor tracts (R2 = 0.180, p < 0.001). They also performed slightly, but significantly, better than other data-driven biomarkers including lesion load of lesion-behaviour maps (R2 =0.200, p < 0.001) and lesion load of structural networks associated with lesion-behaviour maps (R2 =0.167, p < 0.001). Ensemble models - combining basic demographic variables like age, sex, and time since stroke - improved prediction accuracy for theory-based and data-driven biomarkers. Finally, combining both theory-based and data-driven biomarkers with demographic variables improved predictions, and the best ensemble model achieved R2 = 0.241, p < 0.001. Overall, these results demonstrate that models that predict chronic motor outcomes using data-driven features, particularly when lesion data is represented in terms of structural disconnection, perform better than models that predict chronic motor outcomes using theory-based features from the motor system. However, combining both theory-based and data-driven models provides the best predictions.

Association Between the Attention Network Test, Neuropsychological Measures, and Disability in Post-Acute Traumatic Brain Injury.

Cognitive impairment after traumatic brain injury (TBI) is persistent and disabling. Assessing cognitive function in a reliable and valid manner, using measures that are sensitive to the integrity of underlying neural substrates, is crucial in clinical research. The Attention Network Test (ANT) is one such assessment measure that has demonstrated associations with neural regions involved in attention; however, clinical utility of the ANT is limited because its relationship with neuropsychological measures of cognitive function (i.e., its construct validity) has not yet been established in TBI. We evaluated the association between the ANT and 1) a neuropsychological battery assessing executive function and memory and 2) global function assessed by the Glasgow Outcome Scale-Extended (GOSE). Forty-eight adults with complicated mild-severe TBI were evaluated ∼5 months post-injury. Using principal component analysis and multi-variate linear regression adjusted for age, gender, education, and cause of injury, we found that ANT reaction time and executive network scores predicted a principal component assessing processing speed and executive function. Conversely, the ANT did not predict a principal component assessing memory. The ANT was weakly associated with the GOSE. Among persons with TBI during the post-acute phase of recovery, the ANT has good construct validity as evidenced by its associations with neuropsychological measures of processing speed and executive function, but not memory. Given that ANT networks are known to relate to specific neuroanatomical regions, the ANT may be a useful outcome measure for evaluating novel therapeutics targeting attention and executive functions after TBI.

Larger lesion volume in people with multiple sclerosis is associated with increased transition energies between brain states and decreased entropy of brain activity.

Quantifying the relationship between the brain's functional activity patterns and its structural backbone is crucial when relating the severity of brain pathology to disability in multiple sclerosis (MS). Network control theory (NCT) characterizes the brain's energetic landscape using the structural connectome and patterns of brain activity over time. We applied NCT to investigate brain-state dynamics and energy landscapes in controls and people with MS (pwMS). We also computed entropy of brain activity and investigated its association with the dynamic landscape's transition energy and lesion volume. Brain states were identified by clustering regional brain activity vectors, and NCT was applied to compute the energy required to transition between these brain states. We found that entropy was negatively correlated with lesion volume and transition energy, and that larger transition energies were associated with pwMS with disability. This work supports the notion that shifts in the pattern of brain activity in pwMS without disability results in decreased transition energies compared to controls, but, as this shift evolves over the disease, transition energies increase beyond controls and disability occurs. Our results provide the first evidence in pwMS that larger lesion volumes result in greater transition energy between brain states and decreased entropy of brain activity.

Modulating human brain responses via optimal natural image selection and synthetic image generation.

One of the main goals of neuroscience is to understand how biological brains interpret and process incoming environmental information. Building computational encoding models that map images to neural responses is one way to pursue this goal. Moreover, generating or selecting visual stimuli designed to achieve specific patterns of responses allows exploration and control of neuronal firing rates or regional brain activity responses. Here, we investigated the brain's regional activation selectivity and inter-individual differences in human brain responses to various sets of natural and synthetic (generated) images via two functional MRI (fMRI) studies. For our first fMRI study, we used a pre-trained group-level neural model for selecting or synthesizing images that are predicted to maximally activate targeted brain regions. We then presented these images to subjects while collecting their fMRI data. Our results show that optimized images indeed evoke larger magnitude responses than other images predicted to achieve average levels of activation.Furthermore, the activation gain is positively associated with the encoding model accuracy. While most regions' activations in response to maximal natural images and maximal synthetic images were not different, two regions, namely anterior temporal lobe faces (aTLfaces) and fusiform body area 1 (FBA1), had significantly higher activation in response to maximal synthetic images compared to maximal natural images. On the other hand, three regions; medial temporal lobe face area (mTLfaces), ventral word form area 1 (VWFA1) and ventral word form area 2 (VWFA2), had higher activation in response to maximal natural images compared to maximal synthetic images. In our second fMRI experiment, we focused on probing inter-individual differences in face regions' responses and found that individual-specific synthetic (and not natural) images derived using a personalized encoding model elicited significantly higher responses compared to synthetic images derived from the group-level or other subjects' encoding models. Finally, we replicated the finding showing synthetic images elicited larger activation responses in the aTLfaces region compared to natural image responses in that region. Here, for the first time, we leverage our data-driven and generative modeling framework NeuroGen to probe inter-individual differences in and functional specialization of the human visual system. Our results indicate that NeuroGen can be used to modulate macro-scale brain regions in specific individuals using synthetically generated visual stimuli.

Comparisons of electrophysiological markers of impaired executive attention after traumatic brain injury and in healthy aging.

Executive attention impairments are a persistent and debilitating consequence of traumatic brain injury (TBI). To make headway towards treating and predicting outcomes following heterogeneous TBI, cognitive impairment specific pathophysiology first needs to be characterized. In a prospective observational study, we measured EEG during the attention network test aimed at detecting alerting, orienting, executive attention and processing speed. The sample (N = 110) of subjects aged 18-86 included those with and without traumatic brain injury: n = 27, complicated mild TBI; n = 5, moderate TBI; n = 10, severe TBI; n = 63, non-brain-injured controls. Subjects with TBI had impairments in processing speed and executive attention. Electrophysiological markers of executive attention processing in the midline frontal regions reveal that, as a group, those with TBI and elderly non-brain-injured controls have reduced responses. We also note that those with TBI and elderly controls have responses that are similar for both low and high-demand trials. In subjects with moderate-severe TBI, reductions in frontal cortical activation and performance profiles are both similar to that of controls who are ∼4 to 7 years older. Our specific observations of frontal response reductions in subjects with TBI and in older adults is consistent with the suggested role of the anterior forebrain mesocircuit as underlying cognitive impairments. Our results provide novel correlative data linking specific pathophysiological mechanisms underlying domain-specific cognitive deficits following TBI and with normal aging. Collectively, our findings provide biomarkers that may serve to track therapeutic interventions and guide development of targeted therapeutics following brain injuries.

Glymphatic clearance estimated using diffusion tensor imaging along perivascular spaces is reduced after traumatic brain injury and correlates with plasma neurofilament light, a biomarker of injury severity.

The glymphatic system is a perivascular fluid clearance system, most active during sleep, considered important for clearing the brain of waste products and toxins. Glymphatic failure is hypothesized to underlie brain protein deposition in neurodegenerative disorders like Alzheimer's disease. Preclinical evidence suggests that a functioning glymphatic system is also essential for recovery from traumatic brain injury, which involves release of debris and toxic proteins that need to be cleared from the brain. In a cross-sectional observational study, we estimated glymphatic clearance using diffusion tensor imaging along perivascular spaces, an MRI-derived measure of water diffusivity surrounding veins in the periventricular region, in 13 non-injured controls and 37 subjects who had experienced traumatic brain injury ∼5 months previously. We additionally measured the volume of the perivascular space using T2-weighted MRI. We measured plasma concentrations of neurofilament light chain, a biomarker of injury severity, in a subset of subjects. Diffusion tensor imaging along perivascular spaces index was modestly though significantly lower in subjects with traumatic brain injury compared with controls when covarying for age. Diffusion tensor imaging along perivascular spaces index was significantly, negatively correlated with blood levels of neurofilament light chain. Perivascular space volume did not differ in subjects with traumatic brain injury as compared with controls and did not correlate with blood levels of neurofilament light chain, suggesting it may be a less sensitive measure for injury-related perivascular clearance changes. Glymphatic impairment after traumatic brain injury could be due to mechanisms such as mislocalization of glymphatic water channels, inflammation, proteinopathy and/or sleep disruption. Diffusion tensor imaging along perivascular spaces is a promising method for estimating glymphatic clearance, though additional work is needed to confirm results and assess associations with outcome. Understanding changes in glymphatic functioning following traumatic brain injury could inform novel therapies to improve short-term recovery and reduce later risk of neurodegeneration.

Graph-matching distance between individuals' functional connectomes varies with relatedness, age, and cognitive score.

Functional connectomes (FCs), represented by networks or graphs that summarize coactivation patterns between pairs of brain regions, have been related at a population level to age, sex, cognitive/behavioral scores, life experience, genetics, and disease/disorders. However, quantifying FC differences between individuals also provides a rich source of information with which to map to differences in those individuals' biology, experience, genetics or behavior. In this study, graph matching is used to create a novel inter-individual FC metric, called swap distance, that quantifies the distance between pairs of individuals' partial FCs, with a smaller swap distance indicating the individuals have more similar FC. We apply graph matching to align FCs between individuals from the the Human Connectome Project N = 997 and find that swap distance (i) increases with increasing familial distance, (ii) increases with subjects' ages, (iii) is smaller for pairs of females compared to pairs of males, and (iv) is larger for females with lower cognitive scores compared to females with larger cognitive scores. Regions that contributed most to individuals' swap distances were in higher-order networks, that is, default-mode and fronto-parietal, that underlie executive function and memory. These higher-order networks' regions also had swap frequencies that varied monotonically with familial relatedness of the individuals in question. We posit that the proposed graph matching technique provides a novel way to study inter-subject differences in FC and enables quantification of how FC may vary with age, relatedness, sex, and behavior.

The sequence of regional structural disconnectivity due to multiple sclerosis lesions.

Objective: Prediction of disease progression is challenging in multiple sclerosis (MS) as the sequence of lesion development and retention of inflammation within a subset of chronic lesions is heterogeneous among patients. We investigated the sequence of lesion-related regional structural disconnectivity across the spectrum of disability and cognitive impairment in MS. Methods: In a full cohort of 482 patients, the Expanded Disability Status Scale was used to classify patients into (i) no or mild vs (ii) moderate or severe disability groups. In 363 out of 482 patients, Quantitative Susceptibility Mapping was used to identify paramagnetic rim lesions (PRL), which are maintained by a rim of iron-laden innate immune cells. In 171 out of 482 patients, Brief International Cognitive Assessment was used to identify subjects with cognitive impairment. Network Modification Tool was used to estimate the regional structural disconnectivity due to MS lesions. Discriminative event-based modeling was applied to investigate the sequence of regional structural disconnectivity due to all representative lesions across the spectrum of disability and cognitive impairment. Results: Structural disconnection in the ventral attention and subcortical networks was an early biomarker of moderate or severe disability. The earliest biomarkers of disability progression were structural disconnections due to PRL in the motor-related regions. Subcortical structural disconnection was an early biomarker of cognitive impairment. Interpretation: MS lesion-related structural disconnections in the subcortex is an early biomarker for both disability and cognitive impairment in MS. PRL-related structural disconnection in the motor cortex may identify the patients at risk for moderate or severe disability in MS.

Personalized visual encoding model construction with small data.

Quantifying population heterogeneity in brain stimuli-response mapping may allow insight into variability in bottom-up neural systems that can in turn be related to individual's behavior or pathological state. Encoding models that predict brain responses to stimuli are one way to capture this relationship. However, they generally need a large amount of fMRI data to achieve optimal accuracy. Here, we propose an ensemble approach to create encoding models for novel individuals with relatively little data by modeling each subject's predicted response vector as a linear combination of the other subjects' predicted response vectors. We show that these ensemble encoding models trained with hundreds of image-response pairs, achieve accuracy not different from models trained on 20,000 image-response pairs. Importantly, the ensemble encoding models preserve patterns of inter-individual differences in the image-response relationship. We also show the proposed approach is robust against domain shift by validating on data with a different scanner and experimental setup. Additionally, we show that the ensemble encoding models are able to discover the inter-individual differences in various face areas' responses to images of animal vs human faces using a recently developed NeuroGen framework. Our approach shows the potential to use existing densely-sampled data, i.e. large amounts of data collected from a single individual, to efficiently create accurate, personalized encoding models and, subsequently, personalized optimal synthetic images for new individuals scanned under different experimental conditions.

Frontoparietal network activation is associated with motor recovery in ischemic stroke patients.

Strokes cause lesions that damage brain tissue, disrupt normal brain activity patterns and can lead to impairments in motor function. Although modulation of cortical activity is central to stimulation-based rehabilitative therapies, aberrant and adaptive patterns of brain activity after stroke have not yet been fully characterized. Here, we apply a brain dynamics analysis approach to study longitudinal brain activity patterns in individuals with ischemic pontine stroke. We first found 4 commonly occurring brain states largely characterized by high amplitude activations in the visual, frontoparietal, default mode, and motor networks. Stroke subjects spent less time in the frontoparietal state compared to controls. For individuals with dominant-hand CST damage, more time spent in the frontoparietal state from 1 week to 3-6 months post-stroke was associated with better motor recovery over the same time period, an association which was independent of baseline impairment. Furthermore, the amount of time spent in brain states was linked empirically to functional connectivity. This work suggests that when the dominant-hand CST is compromised in stroke, resting state configurations may include increased activation of the frontoparietal network, which may facilitate compensatory neural pathways that support recovery of motor function when traditional motor circuits of the dominant-hemisphere are compromised.

Early effect of thrombolysis on structural brain network organisation after anterior-circulation stroke in the randomized WAKE-UP trial.

The symptoms of acute ischemic stroke can be attributed to disruption of the brain network architecture. Systemic thrombolysis is an effective treatment that preserves structural connectivity in the first days after the event. Its effect on the evolution of global network organisation is, however, not well understood. We present a secondary analysis of 269 patients from the randomized WAKE-UP trial, comparing 127 imaging-selected patients treated with alteplase with 142 controls who received placebo. We used indirect network mapping to quantify the impact of ischemic lesions on structural brain network organisation in terms of both global parameters of segregation and integration, and local disruption of individual connections. Network damage was estimated before randomization and again 22 to 36 h after administration of either alteplase or placebo. Evolution of structural network organisation was characterised by a loss in integration and gain in segregation, and this trajectory was attenuated by the administration of alteplase. Preserved brain network organization was associated with excellent functional outcome. Furthermore, the protective effect of alteplase was spatio-topologically nonuniform, concentrating on a subnetwork of high centrality supported in the salvageable white matter surrounding the ischemic cores. This interplay between the location of the lesion, the pathophysiology of the ischemic penumbra, and the spatial embedding of the brain network explains the observed potential of thrombolysis to attenuate topological network damage early after stroke. Our findings might, in the future, lead to new brain network-informed imaging biomarkers and improved prognostication in ischemic stroke.

Non-neural factors influencing BOLD response magnitudes within individual subjects.

To what extent is the size of the blood-oxygen-level-dependent (BOLD) response influenced by factors other than neural activity? In a re-analysis of three neuroimaging datasets (male and female human participants), we find large systematic inhomogeneities in the BOLD response magnitude in primary visual cortex (V1): stimulus-evoked BOLD responses, expressed in units of percent signal change, are up to 50% larger along the representation of the horizontal meridian than the vertical meridian. To assess whether this surprising effect can be interpreted as differences in local neural activity, we quantified several factors that potentially contribute to the size of the BOLD response. We find relationships between BOLD response magnitude and cortical thickness, curvature, depth and macrovasculature. These relationships are consistently found across subjects and datasets and suggest that variation in BOLD response magnitudes across cortical locations reflects, in part, differences in anatomy and vascularization. To compensate for these factors, we implement a regression-based correction method and show that after correction, BOLD responses become more homogeneous across V1. The correction reduces the horizontal/vertical difference by about half, indicating that some of the difference is likely not due to neural activity differences. We conclude that interpretation of variation in BOLD response magnitude across cortical locations should consider the influence of the potential confounding factors of thickness, curvature, depth and vascularization.SIGNIFICANCE STATEMENTThe magnitude of the BOLD signal is often used as a surrogate of neural activity, but the exact factors that contribute to its strength have not been studied on a voxel-wise level. Here, we examined several anatomical and measurement-related factors to assess their relationship with BOLD signal magnitude. We find that BOLD magnitude correlates with cortical anatomy, depth and macrovasculature. To remove the contribution of these factors, we propose a simple, data-driven correction method that can be used in any functional magnetic resonance imaging (fMRI) experiment. After accounting for the confounding factors, BOLD magnitude becomes more spatially homogenous. Our correction method improves the ability to make more accurate inferences about local neural activity from fMRI data.

Shared functional connections within and between cortical networks predict cognitive abilities in adult males and females.

A thorough understanding of sex-independent and sex-specific neurobiological features that underlie cognitive abilities in healthy individuals is essential for the study of neurological illnesses in which males and females differentially experience and exhibit cognitive impairment. Here, we evaluate sex-independent and sex-specific relationships between functional connectivity and individual cognitive abilities in 392 healthy young adults (196 males) from the Human Connectome Project. First, we establish that sex-independent models comparably predict crystallised abilities in males and females, but only successfully predict fluid abilities in males. Second, we demonstrate sex-specific models comparably predict crystallised abilities within and between sexes, and generally fail to predict fluid abilities in either sex. Third, we reveal that largely overlapping connections between visual, dorsal attention, ventral attention, and temporal parietal networks are associated with better performance on crystallised and fluid cognitive tests in males and females, while connections within visual, somatomotor, and temporal parietal networks are associated with poorer performance. Together, our findings suggest that shared neurobiological features of the functional connectome underlie crystallised and fluid abilities across the sexes.

NeuroGen: Activation optimized image synthesis for discovery neuroscience.

Functional MRI (fMRI) is a powerful technique that has allowed us to characterize visual cortex responses to stimuli, yet such experiments are by nature constructed based on a priori hypotheses, limited to the set of images presented to the individual while they are in the scanner, are subject to noise in the observed brain responses, and may vary widely across individuals. In this work, we propose a novel computational strategy, which we call NeuroGen, to overcome these limitations and develop a powerful tool for human vision neuroscience discovery. NeuroGen combines an fMRI-trained neural encoding model of human vision with a deep generative network to synthesize images predicted to achieve a target pattern of macro-scale brain activation. We demonstrate that the reduction of noise that the encoding model provides, coupled with the generative network's ability to produce images of high fidelity, results in a robust discovery architecture for visual neuroscience. By using only a small number of synthetic images created by NeuroGen, we demonstrate that we can detect and amplify differences in regional and individual human brain response patterns to visual stimuli. We then verify that these discoveries are reflected in the several thousand observed image responses measured with fMRI. We further demonstrate that NeuroGen can create synthetic images predicted to achieve regional response patterns not achievable by the best-matching natural images. The NeuroGen framework extends the utility of brain encoding models and opens up a new avenue for exploring, and possibly precisely controlling, the human visual system.