Continuous evaluation of denoising strategies in resting-state fMRI connectivity using fMRIPrep and Nilearn.

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark prototypes an implementation of a reproducible framework, where the provided Jupyter Book enables readers to reproduce or modify the figures on the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep. Most of the benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing was generally effective, but is incompatible with statistical analyses requiring the continuous sampling of brain signal, for which a simpler strategy, using motion parameters, average activity in select brain compartments, and global signal regression, is preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods.

Prefrontal-Limbic Circuitry Is Associated With Reward Sensitivity in Nonhuman Primates.

BACKGROUND: Abnormal reward sensitivity is a risk factor for psychiatric disorders, including eating disorders such as overeating and binge-eating disorder, but the brain structural mechanisms that underlie it are not completely understood. Here, we sought to investigate the relationship between multimodal whole-brain structural features and reward sensitivity in nonhuman primates. METHODS: Reward sensitivity was evaluated through behavioral economic analysis in which monkeys (adult rhesus macaques; 7 female, 5 male) responded for sweetened condensed milk (10%, 30%, 56%), Gatorade, or water using an operant procedure in which the response requirement increased incrementally across sessions (i.e., fixed ratio 1, 3, 10). Animals were divided into high (n = 6) or low (n = 6) reward sensitivity groups based on essential value for 30% milk. Multimodal magnetic resonance imaging was used to measure gray matter volume and white matter microstructure. Brain structural features were compared between groups, and their correlations with reward sensitivity for various stimuli was investigated. RESULTS: Animals in the high sensitivity group had greater dorsolateral prefrontal cortex, centromedial amygdaloid complex, and middle cingulate cortex volumes than animals in the low sensitivity group. Furthermore, compared with monkeys in the low sensitivity group, high sensitivity monkeys had lower fractional anisotropy in the left dorsal cingulate bundle connecting the centromedial amygdaloid complex and middle cingulate cortex to the dorsolateral prefrontal cortex, and in the left superior longitudinal fasciculus 1 connecting the middle cingulate cortex to the dorsolateral prefrontal cortex. CONCLUSIONS: These results suggest that neuroanatomical variation in prefrontal-limbic circuitry is associated with reward sensitivity. These brain structural features may serve as predictive biomarkers for vulnerability to food-based and other reward-related disorders.

Striatal and Behavioral Responses to Reward Vary by Socioeconomic Status in Adolescents.

Disparities in socioeconomic status (SES) lead to unequal access to financial and social support. These disparities are believed to influence reward sensitivity, which in turn are hypothesized to shape how individuals respond to and pursue rewarding experiences. However, surprisingly little is known about how SES shapes reward sensitivity in adolescence. Here, we investigated how SES influenced adolescent responses to reward, both in behavior and the striatum-a brain region that is highly sensitive to reward. We examined responses to both immediate reward (tracked by phasic dopamine) and average reward rate fluctuations (tracked by tonic dopamine) as these distinct signals independently shape learning and motivation. Adolescents (n = 114; 12-14 years; 58 female) performed a gambling task during functional magnetic resonance imaging. We manipulated trial-by-trial reward and loss outcomes, leading to fluctuations between periods of reward scarcity and abundance. We found that a higher reward rate hastened behavioral responses, and increased guess switching, consistent with the idea that reward abundance increases response vigor and exploration. Moreover, immediate reward reinforced previously rewarding decisions (win-stay, lose-switch) and slowed responses (postreward pausing), particularly when rewards were scarce. Notably, lower-SES adolescents slowed down less after rare rewards than higher-SES adolescents. In the brain, striatal activations covaried with the average reward rate across time and showed greater activations during rewarding blocks. However, these striatal effects were diminished in lower-SES adolescents. These findings show that the striatum tracks reward rate fluctuations, which shape decisions and motivation. Moreover, lower SES appears to attenuate reward-driven behavioral and brain responses.

Neurocognitive mechanisms of co-occurring math difficulties in dyslexia: Differences in executive function and visuospatial processing.

Children with dyslexia frequently also struggle with math. However, studies of reading disability (RD) rarely assess math skill, and the neurocognitive mechanisms underlying co-occurring reading and math disability (RD+MD) are not clear. The current study aimed to identify behavioral and neurocognitive factors associated with co-occurring MD among 86 children with RD. Within this sample, 43% had co-occurring RD+MD and 22% demonstrated a possible vulnerability in math, while 35% had no math difficulties (RD-Only). We investigated whether RD-Only and RD+MD students differed behaviorally in their phonological awareness, reading skills, or executive functions, as well as in the brain mechanisms underlying word reading and visuospatial working memory using functional magnetic resonance imaging (fMRI). The RD+MD group did not differ from RD-Only on behavioral or brain measures of phonological awareness related to speech or print. However, the RD+MD group demonstrated significantly worse working memory and processing speed performance than the RD-Only group. The RD+MD group also exhibited reduced brain activations for visuospatial working memory relative to RD-Only. Exploratory brain-behavior correlations along a broad spectrum of math ability revealed that stronger math skills were associated with greater activation in bilateral visual cortex. These converging neuro-behavioral findings suggest that poor executive functions in general, including differences in visuospatial working memory, are specifically associated with co-occurring MD in the context of RD. RESEARCH HIGHLIGHTS: Children with reading disabilities (RD) frequently have a co-occurring math disability (MD), but the mechanisms behind this high comorbidity are not well understood. We examined differences in phonological awareness, reading skills, and executive function between children with RD only versus co-occurring RD+MD using behavioral and fMRI measures. Children with RD only versus RD+MD did not differ in their phonological processing, either behaviorally or in the brain. RD+MD was associated with additional behavioral difficulties in working memory, and reduced visual cortex activation during a visuospatial working memory task.

Distortion correction of functional MRI without reverse phase encoding scans or field maps.

Functional magnetic resonance images (fMRI) acquired using echo planar sequences typically suffer from spatial distortions due to susceptibility induced off-resonance fields, which may cause geometric mismatch with structural images and affect subsequent quantification and localization of brain function. State-of-the art distortion correction methods (for example, using FSL's topup or AFNI's 3dQwarp algorithms) require the collection of additional scans - either field maps or images with reverse phase encoding directions (i.e., blip-up/blip-down acquisitions) - to estimate and correct distortions. However, not all imaging protocols acquire these additional data and thus cannot take advantage of these post-acquisition corrections. In this study, we aim to enable state-of-the art processing of historical or limited datasets that do not include specific sequences for distortion correction by using only the acquired functional data and a single commonly acquired structural image. To achieve this, we synthesize an undistorted image with contrast similar to the fMRI data and use the non-distorted synthetic image as an anatomical target for distortion correction. We evaluate the efficacy of this approach, named SynBOLD-DisCo (Synthetic BOLD contrast for Distortion Correction), and show that this distortion correction process yields fMRI data that are geometrically similar to non-distorted structural images, with distortion correction virtually equivalent to acquisitions that do contain both blip-up/blip-down images. Our method is available as a Singularity container, source code, and an executable trained model to facilitate evaluation and integration into existing fMRI preprocessing pipelines.

Continuous Evaluation of Denoising Strategies in Resting-State fMRI Connectivity Using fMRIPrep and Nilearn.

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark is implemented in a fully reproducible framework, where the provided research objects enable readers to reproduce or modify core computations, as well as the figures of the article using the Jupyter Book project and the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep software package. The majority of benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing however disrupts the continuous sampling of brain images and is incompatible with some statistical analyses, e.g. auto-regressive modeling. In this case, a simple strategy using motion parameters, average activity in select brain compartments, and global signal regression should be preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods. Our reproducible benchmark infrastructure will facilitate such continuous evaluation in the future, and may also be applied broadly to different tools or even research fields.

ModelArray: An R package for statistical analysis of fixel-wise data.

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Diminished Repetition Suppression Reveals Selective and Systems-Level Face Processing Differences in ASD.

Repeated exposure to a stimulus results in reduced neural response, or repetition suppression, in brain regions responsible for processing that stimulus. This rapid accommodation to repetition is thought to underlie learning, stimulus selectivity, and strengthening of perceptual expectations. Importantly, reduced sensitivity to repetition has been identified in several neurodevelopmental, learning, and psychiatric disorders, including autism spectrum disorder (ASD), a neurodevelopmental disorder characterized by challenges in social communication and repetitive behaviors and restricted interests. Reduced ability to exploit or learn from repetition in ASD is hypothesized to contribute to sensory hypersensitivities, and parallels several theoretical frameworks claiming that ASD individuals show difficulty using regularities in the environment to facilitate behavior. Using fMRI in autistic and neurotypical human adults (females and males), we assessed the status of repetition suppression across two modalities (vision, audition) and with four stimulus categories (faces, objects, printed words, and spoken words). ASD individuals showed domain-specific reductions in repetition suppression for face stimuli only, but not for objects, printed words, or spoken words. Reduced repetition suppression for faces was associated with greater challenges in social communication in ASD. We also found altered functional connectivity between atypically adapting cortical regions and higher-order face recognition regions, and microstructural differences in related white matter tracts in ASD. These results suggest that fundamental neural mechanisms and system-wide circuits are selectively altered for face processing in ASD and enhance our understanding of how disruptions in the formation of stable face representations may relate to higher-order social communication processes.SIGNIFICANCE STATEMENT A common finding in neuroscience is that repetition results in plasticity in stimulus-specific processing regions, reflecting selectivity and adaptation (repetition suppression [RS]). RS is reduced in several neurodevelopmental and psychiatric conditions including autism spectrum disorder (ASD). Theoretical frameworks of ASD posit that reduced adaptation may contribute to associated challenges in social communication and sensory processing. However, the scope of RS differences in ASD is unknown. We examined RS for multiple categories across visual and auditory domains (faces, objects, printed words, spoken words) in autistic and neurotypical individuals. We found reduced RS in ASD for face stimuli only and altered functional connectivity and white matter microstructure between cortical face-recognition areas. RS magnitude correlated with social communication challenges among autistic individuals.

A large-scale investigation of white matter microstructural associations with reading ability.

Reading involves the functioning of a widely distributed brain network, and white matter tracts are responsible for transmitting information between constituent network nodes. Several studies have analyzed fiber bundle microstructural properties to shed insights into the neural basis of reading abilities and disabilities. Findings have been inconsistent, potentially due to small sample sizes and varying methodology. To address this, we analyzed a large data set of 686 children ages 5-18 using state-of-the-art neuroimaging acquisitions and processing techniques. We searched for associations between fractional anisotropy (FA) and single-word and single-nonword reading skills in children with diverse reading abilities across multiple tracts previously thought to contribute to reading. We also looked for group differences in tract FA between typically reading children and children with reading disabilities. FA of the white matter increased with age across all participants. There were no significant correlations between overall reading abilities and tract FAs across all children, and no significant group differences in tract FA between children with and without reading disabilities. There were associations between FA and nonword reading ability in older children (ages 9 and above). Higher FA in the right superior longitudinal fasciculus (SLF) and left inferior cerebellar peduncle (ICP) correlated with better nonword reading skills. These results suggest that letter-sound correspondence skills, as measured by nonword reading, are associated with greater white matter coherence among older children in these two tracts, as indexed by higher FA.

Personalized Connectome Mapping to Guide Targeted Therapy and Promote Recovery of Consciousness in the Intensive Care Unit.

There are currently no therapies proven to promote early recovery of consciousness in patients with severe brain injuries in the intensive care unit (ICU). For patients whose families face time-sensitive, life-or-death decisions, treatments that promote recovery of consciousness are needed to reduce the likelihood of premature withdrawal of life-sustaining therapy, facilitate autonomous self-expression, and increase access to rehabilitative care. Here, we present the Connectome-based Clinical Trial Platform (CCTP), a new paradigm for developing and testing targeted therapies that promote early recovery of consciousness in the ICU. We report the protocol for STIMPACT (Stimulant Therapy Targeted to Individualized Connectivity Maps to Promote ReACTivation of Consciousness), a CCTP-based trial in which intravenous methylphenidate will be used for targeted stimulation of dopaminergic circuits within the subcortical ascending arousal network (ClinicalTrials.gov NCT03814356). The scientific premise of the CCTP and the STIMPACT trial is that personalized brain network mapping in the ICU can identify patients whose connectomes are amenable to neuromodulation. Phase 1 of the STIMPACT trial is an open-label, safety and dose-finding study in 22 patients with disorders of consciousness caused by acute severe traumatic brain injury. Patients in Phase 1 will receive escalating daily doses (0.5-2.0 mg/kg) of intravenous methylphenidate over a 4-day period and will undergo resting-state functional magnetic resonance imaging and electroencephalography to evaluate the drug's pharmacodynamic properties. The primary outcome measure for Phase 1 relates to safety: the number of drug-related adverse events at each dose. Secondary outcome measures pertain to pharmacokinetics and pharmacodynamics: (1) time to maximal serum concentration; (2) serum half-life; (3) effect of the highest tolerated dose on resting-state functional MRI biomarkers of connectivity; and (4) effect of each dose on EEG biomarkers of cerebral cortical function. Predetermined safety and pharmacodynamic criteria must be fulfilled in Phase 1 to proceed to Phase 2A. Pharmacokinetic data from Phase 1 will also inform the study design of Phase 2A, where we will test the hypothesis that personalized connectome maps predict therapeutic responses to intravenous methylphenidate. Likewise, findings from Phase 2A will inform the design of Phase 2B, where we plan to enroll patients based on their personalized connectome maps. By selecting patients for clinical trials based on a principled, mechanistic assessment of their neuroanatomic potential for a therapeutic response, the CCTP paradigm and the STIMPACT trial have the potential to transform the therapeutic landscape in the ICU and improve outcomes for patients with severe brain injuries.

Does data cleaning improve brain state classification?

BACKGROUND: Neuroscientists routinely seek to identify and remove noisy or artifactual observations from their data. They do so with the belief that removing such data improves power to detect relations between neural activity and behavior, which are often subtle and can be overwhelmed by noise. Whereas standard methods can exclude certain well-defined noise sources (e.g., 50/60 Hz electrical noise), in many situations there is not a clear difference between noise and signals so it is not obvious how to separate the two. Here we ask whether methods routinely used to "clean" human electrophysiological recordings lead to greater power to detect brain-behavior relations. NEW METHOD: This, to the authors' knowledge, is the first large-scale simultaneous evaluation of multiple commonly used methods for removing noise from intracranial EEG recordings. RESULTS: We find that several commonly used data cleaning methods (automated methods based on statistical signal properties and manual methods based on expert review) do not increase the power to detect univariate and multivariate electrophysiological biomarkers of successful episodic memory encoding, a well-characterized broadband pattern of neural activity observed across the brain. COMPARISON WITH EXISTING METHODS: Researchers may be more likely to increase statistical power to detect physiological phenomena of interest by allocating resources away from cleaning noisy data and toward collecting more within-patient observations. CONCLUSIONS: These findings highlight the challenge of partitioning signal and noise in the analysis of brain-behavior relations, and suggest increasing sample size and numbers of observations, rather than data cleaning, as the best approach to improving statistical power.