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.

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.

The Bouma law accounts for crowding in 50 observers.

Crowding is the failure to recognize an object due to surrounding clutter. Our visual crowding survey measured 13 crowding distances (or "critical spacings") twice in each of 50 observers. The survey includes three eccentricities (0, 5, and 10 deg), four cardinal meridians, two orientations (radial and tangential), and two fonts (Sloan and Pelli). The survey also tested foveal acuity, twice. Remarkably, fitting a two-parameter model-the well-known Bouma law, where crowding distance grows linearly with eccentricity-explains 82% of the variance for all 13 × 50 measured log crowding distances, cross-validated. An enhanced Bouma law, with factors for meridian, crowding orientation, target kind, and observer, explains 94% of the variance, again cross-validated. These additional factors reveal several asymmetries, consistent with previous reports, which can be expressed as crowding-distance ratios: 0.62 horizontal:vertical, 0.79 lower:upper, 0.78 right:left, 0.55 tangential:radial, and 0.78 Sloan-font:Pelli-font. Across our observers, peripheral crowding is independent of foveal crowding and acuity. Evaluation of the Bouma factor, b (the slope of the Bouma law), as a biomarker of visual health would be easier if there were a way to compare results across crowding studies that use different methods. We define a standardized Bouma factor b' that corrects for differences from Bouma's 25 choice alternatives, 75% threshold criterion, and linearly symmetric flanker placement. For radial crowding on the right meridian, the standardized Bouma factor b' is 0.24 for this study, 0.35 for Bouma (1970), and 0.30 for the geometric mean across five representative modern studies, including this one, showing good agreement across labs, including Bouma's. Simulations, confirmed by data, show that peeking can skew estimates of crowding (e.g., greatly decreasing the mean or doubling the SD of log b). Using gaze tracking to prevent peeking, individual differences are robust, as evidenced by the much larger 0.08 SD of log b across observers than the mere 0.03 test-retest SD of log b measured in half an hour. The ease of measurement of crowding enhances its promise as a biomarker for dyslexia and visual health.

Normal Retinotopy in Primary Visual Cortex in a Congenital Complete Unilateral Lesion of Lateral Geniculate Nucleus in Human: A Case Study.

Impairment of the geniculostriate pathway results in scotomas in the corresponding part of the visual field. Here, we present a case of patient IB with left eye microphthalmia and with lesions in most of the left geniculostriate pathway, including the Lateral Geniculate Nucleus (LGN). Despite the severe lesions, the patient has a very narrow scotoma in the peripheral part of the lower-right-hemifield only (beyond 15° of eccentricity) and complete visual field representation in the primary visual cortex. Population receptive field mapping (pRF) of the patient's visual field reveals orderly eccentricity maps together with contralateral activation in both hemispheres. With diffusion tractography, we revealed connections between superior colliculus (SC) and cortical structures in the hemisphere affected by the lesions, which could mediate the retinotopic reorganization at the cortical level. Our results indicate an astonishing case for the flexibility of the developing retinotopic maps where the contralateral thalamus receives fibers from both the nasal and temporal retinae.

Cross-dataset reproducibility of human retinotopic maps.

Population receptive field (pRF) models fit to fMRI data are used to non-invasively measure retinotopic maps in human visual cortex, and these maps are a fundamental component of visual neuroscience experiments. Here, we examined the reproducibility of retinotopic maps across two datasets: a newly acquired retinotopy dataset from New York University (NYU) (n = 44) and a public dataset from the Human Connectome Project (HCP) (n = 181). Our goal was to assess the degree to which pRF properties are similar across datasets, despite substantial differences in their experimental protocols. The two datasets simultaneously differ in their stimulus apertures, participant pool, fMRI protocol, MRI field strength, and preprocessing pipeline. We assessed the cross-dataset reproducibility of the two datasets in terms of the similarity of vertex-wise pRF estimates and in terms of large-scale polar angle asymmetries in cortical magnification. Within V1, V2, V3, and hV4, the group-median NYU and HCP vertex-wise polar angle estimates were nearly identical. Both eccentricity and pRF size estimates were also strongly correlated between the two datasets, but with a slope different from 1; the eccentricity and pRF size estimates were systematically greater in the NYU data. Next, to compare large-scale map properties, we quantified two polar angle asymmetries in V1 cortical magnification previously identified in the HCP data. The NYU dataset confirms earlier reports that more cortical surface area represents horizontal than vertical visual field meridian, and lower than upper vertical visual field meridian. Together, our findings show that the retinotopic properties of V1, V2, V3, and hV4 can be reliably measured across two datasets, despite numerous differences in their experimental design. fMRI-derived retinotopic maps are reproducible because they rely on an explicit computational model of the fMRI response. In the case of pRF mapping, the model is grounded in physiological evidence of how visual receptive fields are organized, allowing one to quantitatively characterize the BOLD signal in terms of stimulus properties (i.e., location and size). The new NYU Retinotopy Dataset will serve as a useful benchmark for testing hypotheses about the organization of visual areas and for comparison to the HCP 7T Retinotopy Dataset.

Three dimensional MRF obtains highly repeatable and reproducible multi-parametric estimations in the healthy human brain at 1.5T and 3T.

Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories. Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated a high repeatability (CVs 0.7-1.3% for T1, 2.0-7.8% for T2, 1.4-2.5% for normalized PD) and reproducibility (CVs of 2.0-5.8% for T1, 7.4-10.2% for T2, 5.2-9.2% for normalized PD) in gray and white matter. Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain.

PURPOSE: To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF). METHODS: Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps. RESULTS: A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T1 and T2 estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T1 and 18 ms for T2 values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations. CONCLUSION: We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

Conservation of cortical crowding distance across individuals in human V4.

Crowding is the failure to recognize an object due to insufficient spacing, which slows daily tasks such as reading and search. Across 49 observers, we found large variations in psychophysical crowding distance and retinotopic map size. These measures covary, conserving a 1.4-mm cortical crowding distance (threshold object spacing on the cortical surface) in the human V4 map, but not V1-V3. This links the spacing limit of visual recognition to overall V4 size.

EasyEyes - Accurate fixation for online vision testing of crowding and beyond.

Online methods allow testing of larger, more diverse populations, with much less effort than in-lab testing. However, many psychophysical measurements, including visual crowding, require accurate eye fixation, which is classically achieved by testing only experienced observers who have learned to fixate reliably, or by using a gaze tracker to restrict testing to moments when fixation is accurate. Alas, both approaches are impractical online since online observers tend to be inexperienced, and online gaze tracking, using the built-in webcam, has a low precision (±4 deg, Papoutsaki et al., 2016). The EasyEyes open-source software reliably measures peripheral thresholds online with accurate fixation achieved in a novel way, without gaze tracking. EasyEyes tells observers to use the cursor to track a moving crosshair. At a random time during successful tracking, a brief target is presented in the periphery. The observer responds by identifying the target. To evaluate EasyEyes fixation accuracy and thresholds, we tested 12 naive observers in three ways in a counterbalanced order: first, in the lab, using gaze-contingent stimulus presentation (Kurzawski et al., 2023; Pelli et al., 2016); second, in the lab, using EasyEyes while independently monitoring gaze; third, online at home, using EasyEyes. We find that crowding thresholds are consistent (no significant differences in mean and variance of thresholds across ways) and individual differences are conserved. The small root mean square (RMS) fixation error (0.6 deg) during target presentation eliminates the need for gaze tracking. Thus, EasyEyes enables fixation-dependent measurements online, for easy testing of larger and more diverse populations.

EasyEyes - A new method for accurate fixation in online vision testing.

Online methods allow testing of larger, more diverse populations, with much less effort than in-lab testing. However, many psychophysical measurements, including visual crowding, require accurate eye fixation, which is classically achieved by testing only experienced observers who have learned to fixate reliably, or by using a gaze tracker to restrict testing to moments when fixation is accurate. Alas, both approaches are impractical online as online observers tend to be inexperienced, and online gaze tracking, using the built-in webcam, has a low precision (±4 deg). EasyEyes open-source software reliably measures peripheral thresholds online with accurate fixation achieved in a novel way, without gaze tracking. It tells observers to use the cursor to track a moving crosshair. At a random time during successful tracking, a brief target is presented in the periphery. The observer responds by identifying the target. To evaluate EasyEyes fixation accuracy and thresholds, we tested 12 naive observers in three ways in a counterbalanced order: first, in the laboratory, using gaze-contingent stimulus presentation; second, in the laboratory, using EasyEyes while independently monitoring gaze using EyeLink 1000; third, online at home, using EasyEyes. We find that crowding thresholds are consistent and individual differences are conserved. The small root mean square (RMS) fixation error (0.6 deg) during target presentation eliminates the need for gaze tracking. Thus, this method enables fixation-dependent measurements online, for easy testing of larger and more diverse populations.

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.

Short-term plasticity in the human visual thalamus.

While there is evidence that the visual cortex retains a potential for plasticity in adulthood, less is known about the subcortical stages of visual processing. Here, we asked whether short-term ocular dominance plasticity affects the human visual thalamus. We addressed this question in normally sighted adult humans, using ultra-high field (7T) magnetic resonance imaging combined with the paradigm of short-term monocular deprivation. With this approach, we previously demonstrated transient shifts of perceptual eye dominance and ocular dominance in visual cortex (Binda et al., 2018). Here, we report evidence for short-term plasticity in the ventral division of the pulvinar (vPulv), where the deprived eye representation was enhanced over the nondeprived eye. This vPulv plasticity was similar as previously seen in visual cortex and it was correlated with the ocular dominance shift measured behaviorally. In contrast, there was no effect of monocular deprivation in two adjacent thalamic regions: dorsal pulvinar and Lateral Geniculate Nucleus. We conclude that the visual thalamus retains potential for short-term plasticity in adulthood; the plasticity effect differs across thalamic subregions, possibly reflecting differences in their corticofugal connectivity.

Improving the accuracy of single-trial fMRI response estimates using GLMsingle.

Advances in artificial intelligence have inspired a paradigm shift in human neuroscience, yielding large-scale functional magnetic resonance imaging (fMRI) datasets that provide high-resolution brain responses to thousands of naturalistic visual stimuli. Because such experiments necessarily involve brief stimulus durations and few repetitions of each stimulus, achieving sufficient signal-to-noise ratio can be a major challenge. We address this challenge by introducing GLMsingle, a scalable, user-friendly toolbox available in MATLAB and Python that enables accurate estimation of single-trial fMRI responses (glmsingle.org). Requiring only fMRI time-series data and a design matrix as inputs, GLMsingle integrates three techniques for improving the accuracy of trial-wise general linear model (GLM) beta estimates. First, for each voxel, a custom hemodynamic response function (HRF) is identified from a library of candidate functions. Second, cross-validation is used to derive a set of noise regressors from voxels unrelated to the experiment. Third, to improve the stability of beta estimates for closely spaced trials, betas are regularized on a voxel-wise basis using ridge regression. Applying GLMsingle to the Natural Scenes Dataset and BOLD5000, we find that GLMsingle substantially improves the reliability of beta estimates across visually-responsive cortex in all subjects. Comparable improvements in reliability are also observed in a smaller-scale auditory dataset from the StudyForrest experiment. These improvements translate into tangible benefits for higher-level analyses relevant to systems and cognitive neuroscience. We demonstrate that GLMsingle: (i) helps decorrelate response estimates between trials nearby in time; (ii) enhances representational similarity between subjects within and across datasets; and (iii) boosts one-versus-many decoding of visual stimuli. GLMsingle is a publicly available tool that can significantly improve the quality of past, present, and future neuroimaging datasets sampling brain activity across many experimental conditions.

Learning residual motion correction for fast and robust 3D multiparametric MRI.

Voluntary and involuntary patient motion is a major problem for data quality in clinical routine of Magnetic Resonance Imaging (MRI). It has been thoroughly investigated and, yet it still remains unresolved. In quantitative MRI, motion artifacts impair the entire temporal evolution of the magnetization and cause errors in parameter estimation. Here, we present a novel strategy based on residual learning for retrospective motion correction in fast 3D whole-brain multiparametric MRI. We propose a 3D multiscale convolutional neural network (CNN) that learns the non-linear relationship between the motion-affected quantitative parameter maps and the residual error to their motion-free reference. For supervised model training, despite limited data availability, we propose a physics-informed simulation to generate self-contained paired datasets from a priori motion-free data. We evaluate motion-correction performance of the proposed method for the example of 3D Quantitative Transient-state Imaging at 1.5T and 3T. We show the robustness of the motion correction for various motion regimes and demonstrate the generalization capabilities of the residual CNN in terms of real-motion in vivo data of healthy volunteers and clinical patient cases, including pediatric and adult patients with large brain lesions. Our study demonstrates that the proposed motion correction outperforms current state of the art, reliably providing a high, clinically relevant image quality for mild to pronounced patient movements. This has important implications in clinical setups where large amounts of motion affected data must be discarded as they are rendered diagnostically unusable.

The visual white matter connecting human area prostriata and the thalamus is retinotopically organized.

The human visual system is capable of processing visual information from fovea to the far peripheral visual field. Recent fMRI studies have shown a full and detailed retinotopic map in area prostriata, located ventro-dorsally and anterior to the calcarine sulcus along the parieto-occipital sulcus with strong preference for peripheral and wide-field stimulation. Here, we report the anatomical pattern of white matter connections between area prostriata and the thalamus encompassing the lateral geniculate nucleus (LGN). To this end, we developed and utilized an automated pipeline comprising a series of Apps that run openly on the cloud computing platform brainlife.io to analyse 139 subjects of the Human Connectome Project (HCP). We observe a continuous and extended bundle of white matter fibers from which two subcomponents can be extracted: one passing ventrally parallel to the optic radiations (OR) and another passing dorsally circumventing the lateral ventricle. Interestingly, the loop travelling dorsally connects the thalamus with the central visual field representation of prostriata located anteriorly, while the other loop travelling more ventrally connects the LGN with the more peripheral visual field representation located posteriorly. We then analyse an additional cohort of 10 HCP subjects using a manual plane extraction method outside brainlife.io to study the relationship between the two extracted white matter subcomponents and eccentricity, myelin and cortical thickness gradients within prostriata. Our results are consistent with a retinotopic segregation recently demonstrated in the OR, connecting the LGN and V1 in humans and reveal for the first time a retinotopic segregation regarding the trajectory of a fiber bundle between the thalamus and an associative visual area.

Response to short-term deprivation of the human adult visual cortex measured with 7T BOLD.

Sensory deprivation during the post-natal 'critical period' leads to structural reorganization of the developing visual cortex. In adulthood, the visual cortex retains some flexibility and adapts to sensory deprivation. Here we show that short-term (2 hr) monocular deprivation in adult humans boosts the BOLD response to the deprived eye, changing ocular dominance of V1 vertices, consistent with homeostatic plasticity. The boost is strongest in V1, present in V2, V3 and V4 but absent in V3a and hMT+. Assessment of spatial frequency tuning in V1 by a population Receptive-Field technique shows that deprivation primarily boosts high spatial frequencies, consistent with a primary involvement of the parvocellular pathway. Crucially, the V1 deprivation effect correlates across participants with the perceptual increase of the deprived eye dominance assessed with binocular rivalry, suggesting a common origin. Our results demonstrate that visual cortex, particularly the ventral pathway, retains a high potential for homeostatic plasticity in the human adult.

Area Prostriata in the Human Brain.

Area prostriata is a cortical area at the fundus of the calcarine sulcus, described anatomically in humans [1-5] and other primates [6-9]. It is lightly myelinated and lacks the clearly defined six-layer structure evident throughout the cerebral cortex, with a thinner layer 4 and thicker layer 2 [10], characteristic of limbic cortex [11]. In the marmoset and rhesus monkey, area prostriata has cortical connections with MT+ [12], the cingulate motor cortex [8], the auditory cortex [13], the orbitofrontal cortex, and the frontal polar cortices [14]. Here we use functional magnetic resonance together with a wide-field projection system to study its functional properties in humans. With population receptive field mapping [15], we show that area prostriata has a complete representation of the visual field, clearly distinct from the adjacent area V1. As in the marmoset, the caudal-dorsal border of human prostriata-abutting V1-represents the far peripheral visual field, with eccentricities decreasing toward its rostral boundary. Area prostriata responds strongly to very fast motion, greater than 500°/s. The functional properties of area prostriata suggest that it may serve to alert the brain quickly to fast visual events, particularly in the peripheral visual field.