Advancing motion denoising of multiband resting-state functional connectivity fMRI data

Characterization and Mitigation of a Simultaneous Multi-Slice fMRI Artifact: Multiband Artifact Regression in Simultaneous Slices

Simultaneous multi-slice (multiband) acceleration in fMRI has become widespread, but may be affected by novel forms of signal artifact. Here, we demonstrate a previously unreported artifact manifesting as a shared signal between simultaneously acquired slices in all resting-state and task-based multiband fMRI datasets we investigated, including publicly available consortium data from the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) Study. We propose Multiband Artifact Regression in Simultaneous Slices (MARSS), a regression-based detection and correction technique that successfully mitigates this shared signal in unprocessed data. We demonstrate that the signal isolated by MARSS correction is likely nonneural, appearing stronger in neurovasculature than gray matter. Additionally, we evaluate MARSS both against and in tandem with sICA+FIX denoising, which is implemented in HCP resting-state data, to show that MARSS mitigates residual artifact signal that is not modeled by sICA+FIX. MARSS correction leads to study-wide increases in signal-to-noise ratio, decreases in cortical coefficient of variation, and mitigation of systematic artefactual spatial patterns in participant-level task betas. Finally, MARSS correction has substantive effects on second-level t-statistics in analyses of task-evoked activation. We recommend that investigators apply MARSS to multiband fMRI datasets with moderate or higher acceleration factors, in combination with established denoising methods.

Functional Localization of the Human Auditory and Visual Thalamus Using a Thalamic Localizer Functional Magnetic Resonance Imaging Task

Functional magnetic resonance imaging (fMRI) of the auditory and visual sensory systems of the human brain is an active area of investigation in the study of human health and disease. The medial geniculate nucleus (MGN) and lateral geniculate nucleus (LGN) are key thalamic nuclei involved in the processing and relay of auditory and visual information, respectively, and are the subject of blood-oxygen-level-dependent (BOLD) fMRI studies of neural activation and functional connectivity in human participants. However, localization of BOLD fMRI signal originating from neural activity in MGN and LGN remains a technical challenge, due in part to the poor definition of boundaries of these thalamic nuclei in standard T1-weighted and T2-weighted magnetic resonance imaging sequences. Here, we report the development and evaluation of an auditory and visual sensory thalamic localizer (TL) fMRI task that produces participant-specific functionally-defined regions of interest (fROIs) of both MGN and LGN, using 3 Tesla multiband fMRI and a clustered-sparse temporal acquisition sequence, in less than 16 minutes of scan time. We demonstrate the use of MGN and LGN fROIs obtained from the TL fMRI task in standard resting-state functional connectivity (RSFC) fMRI analyses in the same participants. In RSFC analyses, we validated the specificity of MGN and LGN fROIs for signals obtained from primary auditory and visual cortex, respectively, and benchmark their performance against alternative atlas- and segmentation-based localization methods. The TL fMRI task and analysis code (written in Presentation and MATLAB, respectively) have been made freely available to the wider research community.

Brain markers of subtraction and multiplication skills in childhood: task-based functional connectivity and individualized structural similarity

Arithmetic, a high-order cognitive ability, show marked individual difference over development. Despite recent advancements in neuroimaging techniques have enabled the identification of brain markers for individual differences in high-order cognitive abilities, it remains largely unknown about the brain markers for arithmetic. This study used a data-driven connectome-based prediction model to identify brain markers of arithmetic skills from arithmetic-state functional connectivity and individualized structural similarity in 132 children aged 8 to 15 years. We found that both subtraction-state functional connectivity and individualized SS successfully predicted subtraction and multiplication skills but multiplication-state functional connectivity failed to predict either skill. Among the four successful prediction models, most predictive connections were located in frontal-parietal, default-mode, and secondary visual networks. Further computational lesion analyses revealed the essential structural role of frontal-parietal network in predicting subtraction and the essential functional roles of secondary visual, language, and ventral multimodal networks in predicting multiplication. Finally, a few shared nodes but largely nonoverlapping functional and structural connections were found to predict subtraction and multiplication skills. Altogether, our findings provide new insights into the brain markers of arithmetic skills in children and highlight the importance of studying different connectivity modalities and different arithmetic domains to advance our understanding of children's arithmetic skills.

Characterization and Mitigation of a Simultaneous Multi-Slice fMRI Artifact: Multiband Artifact Regression in Simultaneous Slices

Simultaneous multi-slice (multiband) acceleration in fMRI has become widespread, but may be affected by novel forms of signal artifact. Here, we demonstrate a previously unreported artifact manifesting as a shared signal between simultaneously acquired slices in all resting-state and task-based multiband fMRI datasets we investigated, including publicly available consortium data from the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) Study. We propose Multiband Artifact Regression in Simultaneous Slices (MARSS), a regression-based detection and correction technique that successfully mitigates this shared signal in unprocessed data. We demonstrate that the signal isolated by MARSS correction is likely non-neural, appearing stronger in neurovasculature than grey matter. Additionally, we evaluate MARSS both against and in tandem with sICA+FIX denoising, which is implemented in HCP resting-state data, to show that MARSS mitigates residual artifact signal that is not modeled by sICA+FIX. MARSS correction leads to study-wide increases in signal-to-noise ratio, decreases in cortical coefficient of variation, and mitigation of systematic artefactual spatial patterns in participant-level task betas. Finally, MARSS correction has substantive effects on second-level t-statistics in analyses of task-evoked activation. We recommend that investigators apply MARSS to multiband fMRI datasets with moderate or higher acceleration factors, in combination with established denoising methods.

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.

Medial Prefrontal Cortex Dysfunction Mediates Working Memory Deficits in Patients With Schizophrenia

Background

Schizophrenia (SCZ) is marked by working memory (WM) deficits, which predict poor functional outcome. While most functional magnetic resonance imaging studies of WM in SCZ have focused on the dorsolateral prefrontal cortex (PFC), some recent work suggests that the medial PFC (mPFC) may play a role. We investigated whether task-evoked mPFC deactivation is associated with WM performance and whether it mediates deficits in SCZ. In addition, we investigated associations between mPFC deactivation and cortical dopamine release.

Methods

Patients with SCZ (n = 41) and healthy control participants (HCs) (n = 40) performed a visual object n-back task during functional magnetic resonance imaging. Dopamine release capacity in mPFC was quantified with [11C]FLB457 in a subset of participants (9 SCZ, 14 HCs) using an amphetamine challenge. Correlations between task-evoked deactivation and performance were assessed in mPFC and dorsolateral PFC masks and were further examined for relationships with diagnosis and dopamine release.

Results

mPFC deactivation was associated with WM task performance, but dorsolateral PFC activation was not. Deactivation in the mPFC was reduced in patients with SCZ relative to HCs and mediated the relationship between diagnosis and WM performance. In addition, mPFC deactivation was significantly and inversely associated with dopamine release capacity across groups and in HCs alone, but not in patients.

Conclusions

Reduced WM task-evoked mPFC deactivation is a mediator of, and potential substrate for, WM impairment in SCZ, although our study design does not rule out the possibility that these findings could relate to cognition in general rather than WM specifically. We further present preliminary evidence of an inverse association between deactivation during WM tasks and dopamine release capacity in the mPFC.

Less is more: balancing noise reduction and data retention in fMRI with data-driven scrubbing

Functional MRI (fMRI) data may be contaminated by artifacts arising from a myriad of sources, including subject head motion, respiration, heartbeat, scanner drift, and thermal noise. These artifacts cause deviations from common distributional assumptions, introduce spatial and temporal outliers, and reduce the signal-to-noise ratio of the data-all of which can have negative consequences for the accuracy and power of downstream statistical analysis. Scrubbing is a technique for excluding fMRI volumes thought to be contaminated by artifacts and generally comes in two flavors. Motion scrubbing based on subject head motion-derived measures is popular but suffers from a number of drawbacks, among them the need to choose a threshold, a lack of generalizability to multiband acquisitions, and high rates of censoring of individual volumes and entire subjects. Alternatively, data-driven scrubbing methods like DVARS are based on observed noise in the processed fMRI timeseries and may avoid some of these issues. Here we propose "projection scrubbing", a novel data-driven scrubbing method based on a statistical outlier detection framework and strategic dimension reduction, including independent component analysis (ICA), to isolate artifactual variation. We undertake a comprehensive comparison of motion scrubbing with data-driven projection scrubbing and DVARS. We argue that an appropriate metric for the success of scrubbing is maximal data retention subject to reasonable performance on typical benchmarks such as the validity, reliability, and identifiability of functional connectivity. We find that stringent motion scrubbing yields worsened validity, worsened reliability, and produced small improvements to fingerprinting. Meanwhile, data-driven scrubbing methods tend to yield greater improvements to fingerprinting while not generally worsening validity or reliability. Importantly, however, data-driven scrubbing excludes a fraction of the number of volumes or entire sessions compared to motion scrubbing. The ability of data-driven fMRI scrubbing to improve data retention without negatively impacting the quality of downstream analysis has major implications for sample sizes in population neuroscience research.

Can this data be saved? Techniques for high motion in resting state scans of first grade children

Motion remains a significant technical hurdle in fMRI studies of young children. Our aim was to develop a straightforward and effective method for obtaining and preprocessing resting state data from a high-motion pediatric cohort. This approach combines real-time monitoring of head motion with a preprocessing pipeline that uses volume censoring and concatenation alongside independent component analysis based denoising. We evaluated this method using a sample of 108 first grade children (age 6-8) enrolled in a longitudinal study of math development. Data quality was assessed by analyzing the correlation between participant head motion and two key metrics for resting state data, temporal signal-to-noise and functional connectivity. These correlations should be minimal in the absence of noise-related artifacts. We compared these data quality indicators using several censoring thresholds to determine the necessary degree of censoring. Volume censoring was highly effective at removing motion-corrupted volumes and ICA denoising removed much of the remaining motion artifact. With the censoring threshold set to exclude volumes that exceeded a framewise displacement of 0.3 mm, preprocessed data met rigorous standards for data quality while retaining a large majority of subjects (83 % of participants). Overall, results show it is possible to obtain usable resting-state data despite extreme motion in a group of young, untrained subjects.