brainlife.io: a decentralized and open-source cloud platform to support neuroscience research.

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

Reprint of: Mapping human brain lesions and their functional consequences.

Neuroscience has a long history of inferring brain function by examining the relationship between brain injury and subsequent behavioral impairments. The primary advantage of this method over correlative methods is that it can tell us if a certain brain region is necessary for a given cognitive function. In addition, lesion-based analyses provide unique insights into clinical deficits. In the last decade, statistical voxel-based lesion behavior mapping (VLBM) emerged as a powerful method for understanding the architecture of the human brain. This review illustrates how VLBM improves our knowledge of functional brain architecture, as well as how it is inherently limited by its mass-univariate approach. A wide array of recently developed methods appear to supplement traditional VLBM. This paper provides an overview of these new methods, including the use of specialized imaging modalities, the combination of structural imaging with normative connectome data, as well as multivariate analyses of structural imaging data. We see these new methods as complementing rather than replacing traditional VLBM, providing synergistic tools to answer related questions. Finally, we discuss the potential for these methods to become established in cognitive neuroscience and in clinical applications.

Using machine learning-based lesion behavior mapping to identify anatomical networks of cognitive dysfunction: Spatial neglect and attention.

Previous lesion behavior studies primarily used univariate lesion behavior mapping techniques to map the anatomical basis of spatial neglect after right brain damage. These studies led to inconsistent results and lively controversies. Given these inconsistencies, the idea of a wide-spread network that might underlie spatial orientation and neglect has been pushed forward. In such case, univariate lesion behavior mapping methods might have been inherently limited in detecting the presumed network due to limited statistical power. By comparing various univariate analyses with multivariate lesion-mapping based on support vector regression, we aimed to validate the network hypothesis directly in a large sample of 203 newly recruited right brain damaged patients. If the exact same correction factors and parameter combinations (FDR correction and dTLVC for lesion size control) were used, both univariate as well as multivariate approaches uncovered the same complex network pattern underlying spatial neglect. At the cortical level, lesion location dominantly affected the temporal cortex and its borders into inferior parietal and occipital cortices. Beyond, frontal and subcortical gray matter regions as well as white matter tracts connecting these regions were affected. Our findings underline the importance of a right network in spatial exploration and attention and specifically in the emergence of the core symptoms of spatial neglect.

Lead-DBS v2: Towards a comprehensive pipeline for deep brain stimulation imaging.

Deep brain stimulation (DBS) is a highly efficacious treatment option for movement disorders and a growing number of other indications are investigated in clinical trials. To ensure optimal treatment outcome, exact electrode placement is required. Moreover, to analyze the relationship between electrode location and clinical results, a precise reconstruction of electrode placement is required, posing specific challenges to the field of neuroimaging. Since 2014 the open source toolbox Lead-DBS is available, which aims at facilitating this process. The tool has since become a popular platform for DBS imaging. With support of a broad community of researchers worldwide, methods have been continuously updated and complemented by new tools for tasks such as multispectral nonlinear registration, structural/functional connectivity analyses, brain shift correction, reconstruction of microelectrode recordings and orientation detection of segmented DBS leads. The rapid development and emergence of these methods in DBS data analysis require us to revisit and revise the pipelines introduced in the original methods publication. Here we demonstrate the updated DBS and connectome pipelines of Lead-DBS using a single patient example with state-of-the-art high-field imaging as well as a retrospective cohort of patients scanned in a typical clinical setting at 1.5T. Imaging data of the 3T example patient is co-registered using five algorithms and nonlinearly warped into template space using ten approaches for comparative purposes. After reconstruction of DBS electrodes (which is possible using three methods and a specific refinement tool), the volume of tissue activated is calculated for two DBS settings using four distinct models and various parameters. Finally, four whole-brain tractography algorithms are applied to the patient's preoperative diffusion MRI data and structural as well as functional connectivity between the stimulation volume and other brain areas are estimated using a total of eight approaches and datasets. In addition, we demonstrate impact of selected preprocessing strategies on the retrospective sample of 51 PD patients. We compare the amount of variance in clinical improvement that can be explained by the computer model depending on the preprocessing method of choice. This work represents a multi-institutional collaborative effort to develop a comprehensive, open source pipeline for DBS imaging and connectomics, which has already empowered several studies, and may facilitate a variety of future studies in the field.

Revealing the dual streams of speech processing.

Several dual route models of human speech processing have been proposed suggesting a large-scale anatomical division between cortical regions that support motor-phonological aspects vs. lexical-semantic aspects of speech processing. However, to date, there is no complete agreement on what areas subserve each route or the nature of interactions across these routes that enables human speech processing. Relying on an extensive behavioral and neuroimaging assessment of a large sample of stroke survivors, we used a data-driven approach using principal components analysis of lesion-symptom mapping to identify brain regions crucial for performance on clusters of behavioral tasks without a priori separation into task types. Distinct anatomical boundaries were revealed between a dorsal frontoparietal stream and a ventral temporal-frontal stream associated with separate components. Collapsing over the tasks primarily supported by these streams, we characterize the dorsal stream as a form-to-articulation pathway and the ventral stream as a form-to-meaning pathway. This characterization of the division in the data reflects both the overlap between tasks supported by the two streams as well as the observation that there is a bias for phonological production tasks supported by the dorsal stream and lexical-semantic comprehension tasks supported by the ventral stream. As such, our findings show a division between two processing routes that underlie human speech processing and provide an empirical foundation for studying potential computational differences that distinguish between the two routes.

Mapping human brain lesions and their functional consequences.

Neuroscience has a long history of inferring brain function by examining the relationship between brain injury and subsequent behavioral impairments. The primary advantage of this method over correlative methods is that it can tell us if a certain brain region is necessary for a given cognitive function. In addition, lesion-based analyses provide unique insights into clinical deficits. In the last decade, statistical voxel-based lesion behavior mapping (VLBM) emerged as a powerful method for understanding the architecture of the human brain. This review illustrates how VLBM improves our knowledge of functional brain architecture, as well as how it is inherently limited by its mass-univariate approach. A wide array of recently developed methods appear to supplement traditional VLBM. This paper provides an overview of these new methods, including the use of specialized imaging modalities, the combination of structural imaging with normative connectome data, as well as multivariate analyses of structural imaging data. We see these new methods as complementing rather than replacing traditional VLBM, providing synergistic tools to answer related questions. Finally, we discuss the potential for these methods to become established in cognitive neuroscience and in clinical applications.

Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank.

UK Biobank is a large-scale prospective epidemiological study with all data accessible to researchers worldwide. It is currently in the process of bringing back 100,000 of the original participants for brain, heart and body MRI, carotid ultrasound and low-dose bone/fat x-ray. The brain imaging component covers 6 modalities (T1, T2 FLAIR, susceptibility weighted MRI, Resting fMRI, Task fMRI and Diffusion MRI). Raw and processed data from the first 10,000 imaged subjects has recently been released for general research access. To help convert this data into useful summary information we have developed an automated processing and QC (Quality Control) pipeline that is available for use by other researchers. In this paper we describe the pipeline in detail, following a brief overview of UK Biobank brain imaging and the acquisition protocol. We also describe several quantitative investigations carried out as part of the development of both the imaging protocol and the processing pipeline.

"Nonspatial" Attentional Deficits Interact with Spatial Position in Neglect.

A widely debated question concerns whether or not spatial and nonspatial components of visual attention interact in attentional performance. Spatial neglect is a common consequence of brain injury where individuals fail to respond to stimuli presented on their contralesional side. It has been argued that, beyond the spatial bias, these individuals also tend to exhibit nonspatial perceptual deficits. Here we demonstrate that the "nonspatial" deficits affecting the temporal dynamics of attentional deployment are in fact modulated by spatial position. Specifically, we observed that the pathological attentional blink of chronic neglect is enhanced when stimuli are presented on the contralesional side of the trunk while keeping retinal and head-centered coordinates constant. We did not find this pattern in right brain-damaged patients without neglect or in patients who had recovered from neglect. Our work suggests that the nonspatial attentional deficits observed in neglect are heavily modulated by egocentric spatial position. This provides strong evidence against models that suggest independent modules for spatial and nonspatial attentional functions while also providing strong evidence that trunk position plays an important role in neglect.

Decoding the anatomical network of spatial attention.

The study of stroke patients with modern lesion-symptom analysis techniques has yielded valuable insights into the representation of spatial attention in the human brain. Here we introduce an approach--multivariate pattern analysis--that no longer assumes independent contributions of brain regions but rather quantifies the joint contribution of multiple brain regions in determining behavior. In a large sample of stroke patients, we found patterns of damage more predictive of spatial neglect than the best-performing single voxel. In addition, modeling multiple brain regions--those that are frequently damaged and, importantly, spared--provided more predictive information than modeling single regions. Interestingly, we also found that the superior temporal gyrus demonstrated a consistent ability to improve classifier performance when added to other regions, implying uniquely predictive information. In sharp contrast, classifier performance for both the angular gyrus and insular cortex was reliably enhanced by the addition of other brain regions, suggesting these regions lack independent predictive information for spatial neglect. Our findings highlight the utility of multivariate pattern analysis in lesion mapping, furnishing neuroscience with a modern approach for using lesion data to study human brain function.

Automated lesion detection on MRI scans using combined unsupervised and supervised methods.

BACKGROUND: Accurate and precise detection of brain lesions on MR images (MRI) is paramount for accurately relating lesion location to impaired behavior. In this paper, we present a novel method to automatically detect brain lesions from a T1-weighted 3D MRI. The proposed method combines the advantages of both unsupervised and supervised methods. METHODS: First, unsupervised methods perform a unified segmentation normalization to warp images from the native space into a standard space and to generate probability maps for different tissue types, e.g., gray matter, white matter and fluid. This allows us to construct an initial lesion probability map by comparing the normalized MRI to healthy control subjects. Then, we perform non-rigid and reversible atlas-based registration to refine the probability maps of gray matter, white matter, external CSF, ventricle, and lesions. These probability maps are combined with the normalized MRI to construct three types of features, with which we use supervised methods to train three support vector machine (SVM) classifiers for a combined classifier. Finally, the combined classifier is used to accomplish lesion detection. RESULTS: We tested this method using T1-weighted MRIs from 60 in-house stroke patients. Using leave-one-out cross validation, the proposed method can achieve an average Dice coefficient of 73.1% when compared to lesion maps hand-delineated by trained neurologists. Furthermore, we tested the proposed method on the T1-weighted MRIs in the MICCAI BRATS 2012 dataset. The proposed method can achieve an average Dice coefficient of 66.5% in comparison to the expert annotated tumor maps provided in MICCAI BRATS 2012 dataset. In addition, on these two test datasets, the proposed method shows competitive performance to three state-of-the-art methods, including Stamatakis et al., Seghier et al., and Sanjuan et al. CONCLUSIONS: In this paper, we introduced a novel automated procedure for lesion detection from T1-weighted MRIs by combining both an unsupervised and a supervised component. In the unsupervised component, we proposed a method to identify lesioned hemisphere to help normalize the patient MRI with lesions and initialize/refine a lesion probability map. In the supervised component, we extracted three different-order statistical features from both the tissue/lesion probability maps obtained from the unsupervised component and the original MRI intensity. Three support vector machine classifiers are then trained for the three features respectively and combined for final voxel-based lesion classification.

Distinct brain morphometry patterns revealed by deep learning improve prediction of post-stroke aphasia severity.

BACKGROUND: Emerging evidence suggests that post-stroke aphasia severity depends on the integrity of the brain beyond the lesion. While measures of lesion anatomy and brain integrity combine synergistically to explain aphasic symptoms, substantial interindividual variability remains unaccounted. One explanatory factor may be the spatial distribution of morphometry beyond the lesion (e.g., atrophy), including not just specific brain areas, but distinct three-dimensional patterns. METHODS: Here, we test whether deep learning with Convolutional Neural Networks (CNNs) on whole brain morphometry (i.e., segmented tissue volumes) and lesion anatomy better predicts chronic stroke individuals with severe aphasia (N = 231) than classical machine learning (Support Vector Machines; SVMs), evaluating whether encoding spatial dependencies identifies uniquely predictive patterns. RESULTS: CNNs achieve higher balanced accuracy and F1 scores, even when SVMs are nonlinear or integrate linear or nonlinear dimensionality reduction. Parity only occurs when SVMs access features learned by CNNs. Saliency maps demonstrate that CNNs leverage distributed morphometry patterns, whereas SVMs focus on the area around the lesion. Ensemble clustering of CNN saliencies reveals distinct morphometry patterns unrelated to lesion size, consistent across individuals, and which implicate unique networks associated with different cognitive processes as measured by the wider neuroimaging literature. Individualized predictions depend on both ipsilateral and contralateral features outside the lesion. CONCLUSIONS: Three-dimensional network distributions of morphometry are directly associated with aphasia severity, underscoring the potential for CNNs to improve outcome prognostication from neuroimaging data, and highlighting the prospective benefits of interrogating spatial dependence at different scales in multivariate feature space.

Some stroke survivors experience difficulties understanding and producing language. We performed brain imaging to capture information about brain structure in stroke survivors and used it to predict which survivors have more severe language problems. We found that a type of artificial intelligence (AI) specifically designed to find patterns in spatial data was more accurate at this task than more traditional methods. AI found more complex patterns of brain structure that distinguish stroke survivors with severe language problems by analyzing the brain's spatial properties. Our findings demonstrate that AI tools can provide new information about brain structure and function following stroke. With further developments, these models may be able to help clinicians understand the extent to which language problems can be improved after a stroke.

Subject position affects EEG magnitudes.

EEG (electroencephalography) has been used for decades in thousands of research studies and is today a routine clinical tool despite the small magnitude of measured scalp potentials. It is widely accepted that the currents originating in the brain are strongly influenced by the high resistivity of skull bone, but it is less well known that the thin layer of CSF (cerebrospinal fluid) has perhaps an even more important effect on EEG scalp magnitude by spatially blurring the signals. Here it is shown that brain shift and the resulting small changes in CSF layer thickness, induced by changing the subject's position, have a significant effect on EEG signal magnitudes in several standard visual paradigms. For spatially incoherent high-frequency activity the effect produced by switching from prone to supine can be dramatic, increasing occipital signal power by several times for some subjects (on average 80%). MRI measurements showed that the occipital CSF layer between the brain and skull decreases by approximately 30% in thickness when a subject moves from prone to supine position. A multiple dipole model demonstrated that this can indeed lead to occipital EEG signal power increases in the same direction and order of magnitude as those observed here. These results suggest that future EEG studies should control for subjects' posture, and that some studies may consider placing their subjects into the most favorable position for the experiment. These findings also imply that special consideration should be given to EEG measurements from subjects with brain atrophy due to normal aging or neurodegenerative diseases, since the resulting increase in CSF layer thickness could profoundly decrease scalp potential measurements.

Distinct brain morphometry patterns revealed by deep learning improve prediction of aphasia severity.

Emerging evidence suggests that post-stroke aphasia severity depends on the integrity of the brain beyond the stroke lesion. While measures of lesion anatomy and brain integrity combine synergistically to explain aphasic symptoms, significant interindividual variability remains unaccounted for. A possible explanatory factor may be the spatial distribution of brain atrophy beyond the lesion. This includes not just the specific brain areas showing atrophy, but also distinct three-dimensional patterns of atrophy. Here, we tested whether deep learning with Convolutional Neural Networks (CNN) on whole brain morphometry (i.e., segmented tissue volumes) and lesion anatomy can better predict which individuals with chronic stroke (N=231) have severe aphasia, and whether encoding spatial dependencies in the data might be capable of improving predictions by identifying unique individualized spatial patterns. We observed that CNN achieves significantly higher accuracy and F1 scores than Support Vector Machine (SVM), even when the SVM is nonlinear or integrates linear and nonlinear dimensionality reduction techniques. Performance parity was only achieved when the SVM was directly trained on the latent features learned by the CNN. Saliency maps demonstrated that the CNN leveraged widely distributed patterns of brain atrophy predictive of aphasia severity, whereas the SVM focused almost exclusively on the area around the lesion. Ensemble clustering of CNN saliency maps revealed distinct morphometry patterns that were unrelated to lesion size, highly consistent across individuals, and implicated unique brain networks associated with different cognitive processes as measured by the wider neuroimaging literature. Individualized predictions of severity depended on both ipsilateral and contralateral features outside of the location of stroke. Our findings illustrate that three-dimensional network distributions of atrophy in individuals with aphasia are directly associated with aphasia severity, underscoring the potential for deep learning to improve prognostication of behavioral outcomes from neuroimaging data, and highlighting the prospective benefits of interrogating spatial dependence at different scales in multivariate feature space.

Cerebral blood flow in patients recovered from mild COVID-19.

BACKGROUND AND PURPOSE: Cerebral hypoperfusion has been described in both severe and mild forms of symptomatic Coronavirus Disease 2019 (COVID-19) infection. The purpose of this study was to investigate global and regional cerebral blood flow (CBF) in asymptomatic COVID-19 patients. METHODS: Cases with mild COVID-19 infection and age-, sex-, and race-matched healthy controls were drawn from the Aging Brain Consortium at The University of South Carolina data repository. Demographics, risk factors, and data from the Montreal Cognitive Assessment were collected. Mean CBF values for gray matter (GM), white matter (WM), and the whole brain were calculated by averaging CBF values of standard space-normalized CBF image values falling within GM and WM masks. Whole brain region of interest-based analyses were used to create standardized CBF maps and explore differences between groups. RESULTS: Twenty-eight cases with prior mild COVID-19 infection were compared with 28 controls. Whole-brain CBF (46.7 ± 5.6 vs. 49.3 ± 3.7, p = .05) and WM CBF (29.3 ± 2.6 vs. 31.0 ± 1.6, p = .03) were noted to be significantly lower in COVID-19 cases as compared to controls. Predictive models based on these data predicted COVID-19 group membership with a high degree of accuracy (85.2%, p < .001), suggesting CBF patterns are an imaging marker of mild COVID-19 infection. CONCLUSION: In this study, lower WM CBF, as well as widespread regional CBF changes identified using quantitative MRI, was found in mild COVID-19 patients. Further studies are needed to determine the reliability of this newly identified COVID-19 brain imaging marker and determine what drives these CBF changes.

brainlife.io: A decentralized and open source cloud platform to support neuroscience research.

Neuroscience research has expanded dramatically over the past 30 years by advancing standardization and tool development to support rigor and transparency. Consequently, the complexity of the data pipeline has also increased, hindering access to FAIR data analysis to portions of the worldwide research community. brainlife.io was developed to reduce these burdens and democratize modern neuroscience research across institutions and career levels. Using community software and hardware infrastructure, the platform provides open-source data standardization, management, visualization, and processing and simplifies the data pipeline. brainlife.io automatically tracks the provenance history of thousands of data objects, supporting simplicity, efficiency, and transparency in neuroscience research. Here brainlife.io's technology and data services are described and evaluated for validity, reliability, reproducibility, replicability, and scientific utility. Using data from 4 modalities and 3,200 participants, we demonstrate that brainlife.io's services produce outputs that adhere to best practices in modern neuroscience research.

Age-specific CT and MRI templates for spatial normalization.

Spatial normalization reshapes an individual's brain to match the shape and size of a template image. This is a crucial step required for group-level statistical analyses. The most popular standard templates are derived from MRI scans of young adults. We introduce specialized templates that allow normalization algorithms to be applied to stroke-aged populations. First, we developed a CT template: while this is the dominant modality for many clinical situations, there are no modern CT templates and popular algorithms fail to successfully normalize CT scans. Importantly, our template was based on healthy individuals with ages similar to what is commonly seen in stroke (mean 65 years old). This template allows studies where only CT scans are available. Second, we derived a MRI template that approximately matches the shape of our CT template as well as processing steps that aid the normalization of scans from older individuals (including lesion masking and the ability to generate high quality cortical renderings despite brain injury). The benefit of this strategy is that the resulting templates can be used in studies where mixed modalities are present. We have integrated these templates and processing algorithms into a simple SPM toolbox (http://www.mricro.com/clinical-toolbox/spm8-scripts).

Disconnectomics to unravel the network underlying deficits of spatial exploration and attention.

Spatial attention and exploration are related to a predominantly right hemispheric network structure. However, the areas of the brain involved and their exact role is still debated. Spatial neglect following right hemispheric stroke lesions has been frequently viewed as a model to study these processes in humans. Previous investigations on the anatomical basis on spatial neglect predominantly focused on focal brain damage and lesion-behaviour mapping analyses. This approach might not be suited to detect remote areas structurally spared but which might contribute to the behavioural deficit. In the present study of a sample of 203 right hemispheric stroke patients, we combined connectome lesion-symptom mapping with multivariate support vector regression to unravel the complex and disconnected network structure in spatial neglect. We delineated three central nodes that were extensively disconnected from other intrahemispheric areas, namely the right superior parietal lobule, the insula, and the temporal pole. Additionally, the analysis allocated central roles within this network to the inferior frontal gyrus (pars triangularis and opercularis), right middle temporal gyrus, right temporal pole and left and right orbitofrontal cortices, including interhemispheric disconnection. Our results suggest that these structures-although not necessarily directly damaged-might play a role within the network underlying spatial neglect in humans.

Right hemispheric white matter hyperintensities improve the prediction of spatial neglect severity in acute stroke.

White matter hyperintensities (WMH) are frequently observed in brain scans of elderly people. They are associated with an increased risk of stroke, cognitive decline, and dementia. However, it is unknown yet if measures of WMH provide information that improve the understanding of poststroke outcome compared to only state-of-the-art stereotaxic structural lesion data. We implemented high-dimensional machine learning models, based on support vector regression, to predict the severity of spatial neglect in 103 acute right hemispheric stroke patients. We found that (1) the additional information of right hemispheric or bilateral voxel-based topographic WMH extent indeed yielded a significant improvement in predicting acute neglect severity (compared to the voxel-based stroke lesion map alone). (2) Periventricular WMH appeared more relevant for prediction than deep subcortical WMH. (3) Among different measures of WMH, voxel-based maps as measures of topographic extent allowed more accurate predictions compared to the use of traditional ordinally assessed visual rating scales (Fazekas-scale, Cardiovascular Health Study-scale). In summary, topographic WMH appear to be a valuable clinical imaging biomarker for predicting the severity of cognitive deficits and bears great potential for rehabilitation guidance of acute stroke patients.