Beware (Surprisingly Common) Left-Right Flips in Your MRI Data: An Efficient and Robust Method to Check MRI Dataset Consistency Using AFNI

Solving the Pervasive Problem of Protocol Non-Compliance in MRI using an Open-Source tool mrQA

Pooling data across diverse sources acquired by multisite consortia requires compliance with a predefined reference protocol i.e., ensuring different sites and scanners for a given project have used identical or compatible MR physics parameter values. Traditionally, this has been an arduous and manual process due to difficulties in working with the complicated DICOM standard and lack of resources allocated towards protocol compliance. Moreover, issues of protocol compliance is often overlooked for lack of realization that parameter values are routinely improvised/modified locally at various sites. The inconsistencies in acquisition protocols can reduce SNR, statistical power, and in the worst case, may invalidate the results altogether. An open-source tool, mrQA was developed to automatically assess protocol compliance on standard dataset formats such as DICOM and BIDS, and to study the patterns of non-compliance in over 20 open neuroimaging datasets, including the large ABCD study. The results demonstrate that the lack of compliance is rather pervasive. The frequent sources of non-compliance include but are not limited to deviations in Repetition Time, Echo Time, Flip Angle, and Phase Encoding Direction. It was also observed that GE and Philips scanners exhibited higher rates of non-compliance relative to the Siemens scanners in the ABCD dataset. Continuous monitoring for protocol compliance is strongly recommended before any pre/post-processing, ideally right after the acquisition, to avoid the silent propagation of severe/subtle issues. Although, this study focuses on neuroimaging datasets, the proposed tool mrQA can work with any DICOM-based datasets.

Design and methodology for a proof of mechanism study of individualized neuronavigated continuous Theta burst stimulation for auditory processing in adolescents with autism spectrum disorder

It has been suggested that aberrant excitation/inhibition (E/I) balance and dysfunctional structure and function of relevant brain networks may underlie the symptoms of autism spectrum disorder (ASD). However, the nomological network linking these constructs to quantifiable measures and mechanistically relating these constructs to behavioral symptoms of ASD is lacking. Herein we describe a within-subject, controlled, proof-of-mechanism study investigating the pathophysiology of auditory/language processing in adolescents with ASD. We utilize neurophysiological and neuroimaging techniques including magnetic resonance spectroscopy (MRS), diffusion-weighted imaging (DWI), functional magnetic resonance imaging (fMRI), and magnetoencephalography (MEG) metrics of language network structure and function. Additionally, we apply a single, individually targeted session of continuous theta burst stimulation (cTBS) as an experimental probe of the impact of perturbation of the system on these neurophysiological and neuroimaging outcomes. MRS, fMRI, and MEG measures are evaluated at baseline and immediately prior to and following cTBS over the posterior superior temporal cortex (pSTC), a region involved in auditory and language processing deficits in ASD. Also, behavioral measures of ASD and language processing and DWI measures of auditory/language network structures are obtained at baseline to characterize the relationship between the neuroimaging and neurophysiological measures and baseline symptom presentation. We hypothesize that local gamma-aminobutyric acid (GABA) and glutamate concentrations (measured with MRS), and structural and functional activity and network connectivity (measured with DWI and fMRI), will significantly predict MEG indices of auditory/language processing and behavioral deficits in ASD. Furthermore, a single session of cTBS over left pSTC is hypothesized to lead to significant, acute changes in local glutamate and GABA concentration, functional activity and network connectivity, and MEG indices of auditory/language processing. We have completed the pilot phase of the study (n=20 Healthy Volunteer adults) and have begun enrollment for the main phase with adolescents with ASD (n=86; age 14-17). If successful, this study will establish a nomological network linking local E/I balance measures to functional and structural connectivity within relevant brain networks, ultimately connecting them to ASD symptoms. Furthermore, this study will inform future therapeutic trials using cTBS to treat the symptoms of ASD.

Novel inductively coupled ear-bars (ICEs) to enhance restored fMRI signal from susceptibility compensation in rats

Functional magnetic resonance imaging faces inherent challenges when applied to deep-brain areas in rodents, e.g. entorhinal cortex, due to the signal loss near the ear cavities induced by susceptibility artifacts and reduced sensitivity induced by the long distance from the surface array coil. Given the pivotal roles of deep brain regions in various diseases, optimized imaging techniques are needed. To mitigate susceptibility-induced signal losses, we introduced baby cream into the middle ear. To enhance the detection sensitivity of deep brain regions, we implemented inductively coupled ear-bars, resulting in approximately a 2-fold increase in sensitivity in entorhinal cortex. Notably, the inductively coupled ear-bar can be seamlessly integrated as an add-on device, without necessitating modifications to the scanner interface. To underscore the versatility of inductively coupled ear-bars, we conducted echo-planner imaging-based task functional magnetic resonance imaging in rats modeling Alzheimer's disease. As a proof of concept, we also demonstrated resting-state-functional magnetic resonance imaging connectivity maps originating from the left entorhinal cortex-a central hub for memory and navigation networks-to amygdala hippocampal area, Insular Cortex, Prelimbic Systems, Cingulate Cortex, Secondary Visual Cortex, and Motor Cortex. This work demonstrates an optimized procedure for acquiring large-scale networks emanating from a previously challenging seed region by conventional magnetic resonance imaging detectors, thereby facilitating improved observation of functional magnetic resonance imaging outcomes.

BOLD Response is more than just magnitude: Improving detection sensitivity through capturing hemodynamic profiles

Typical fMRI analyses often assume a canonical hemodynamic response function (HRF) that primarily focuses on the peak height of the overshoot, neglecting other morphological aspects. Consequently, reported analyses often reduce the overall response curve to a single scalar value. In this study, we take a data-driven approach to HRF estimation at the whole-brain voxel level, without assuming a response profile at the individual level. We then employ a roughness penalty at the population level to estimate the response curve, aiming to enhance predictive accuracy, inferential efficiency, and cross-study reproducibility. By examining a fast event-related FMRI dataset, we demonstrate the shortcomings and information loss associated with adopting the canonical approach. Furthermore, we address the following key questions: 1) To what extent does the HRF shape vary across different regions, conditions, and participant groups? 2) Does the data-driven approach improve detection sensitivity compared to the canonical approach? 3) Can analyzing the HRF shape help validate the presence of an effect in conjunction with statistical evidence? 4) Does analyzing the HRF shape offer evidence for whole-brain response during a simple task?

Highlight results, don't hide them: Enhance interpretation, reduce biases and improve reproducibility

Most neuroimaging studies display results that represent only a tiny fraction of the collected data. While it is conventional to present "only the significant results" to the reader, here we suggest that this practice has several negative consequences for both reproducibility and understanding. This practice hides away most of the results of the dataset and leads to problems of selection bias and irreproducibility, both of which have been recognized as major issues in neuroimaging studies recently. Opaque, all-or-nothing thresholding, even if well-intentioned, places undue influence on arbitrary filter values, hinders clear communication of scientific results, wastes data, is antithetical to good scientific practice, and leads to conceptual inconsistencies. It is also inconsistent with the properties of the acquired data and the underlying biology being studied. Instead of presenting only a few statistically significant locations and hiding away the remaining results, studies should "highlight" the former while also showing as much as possible of the rest. This is distinct from but complementary to utilizing data sharing repositories: the initial presentation of results has an enormous impact on the interpretation of a study. We present practical examples and extensions of this approach for voxelwise, regionwise and cross-study analyses using publicly available data that was analyzed previously by 70 teams (NARPS; Botvinik-Nezer, et al., 2020), showing that it is possible to balance the goals of displaying a full set of results with providing the reader reasonably concise and "digestible" findings. In particular, the highlighting approach sheds useful light on the kind of variability present among the NARPS teams' results, which is primarily a varied strength of agreement rather than disagreement. Using a meta-analysis built on the informative "highlighting" approach shows this relative agreement, while one using the standard "hiding" approach does not. We describe how this simple but powerful change in practice-focusing on highlighting results, rather than hiding all but the strongest ones-can help address many large concerns within the field, or at least to provide more complete information about them. We include a list of practical suggestions for results reporting to improve reproducibility, cross-study comparisons and meta-analyses.

The art and science of using quality control to understand and improve fMRI data

Designing and executing a good quality control (QC) process is vital to robust and reproducible science and is often taught through hands on training. As FMRI research trends toward studies with larger sample sizes and highly automated processing pipelines, the people who analyze data are often distinct from those who collect and preprocess the data. While there are good reasons for this trend, it also means that important information about how data were acquired, and their quality, may be missed by those working at later stages of these workflows. Similarly, an abundance of publicly available datasets, where people (not always correctly) assume others already validated data quality, makes it easier for trainees to advance in the field without learning how to identify problematic data. This manuscript is designed as an introduction for researchers who are already familiar with fMRI, but who did not get hands on QC training or who want to think more deeply about QC. This could be someone who has analyzed fMRI data but is planning to personally acquire data for the first time, or someone who regularly uses openly shared data and wants to learn how to better assess data quality. We describe why good QC processes are important, explain key priorities and steps for fMRI QC, and as part of the FMRI Open QC Project, we demonstrate some of these steps by using AFNI software and AFNI’s QC reports on an openly shared dataset. A good QC process is context dependent and should address whether data have the potential to answer a scientific question, whether any variation in the data has the potential to skew or hide key results, and whether any problems can potentially be addressed through changes in acquisition or data processing. Automated metrics are essential and can often highlight a possible problem, but human interpretation at every stage of a study is vital for understanding causes and potential solutions.

Quality control practices in FMRI analysis: Philosophy, methods and examples using AFNI

Quality control (QC) is a necessary, but often an under-appreciated, part of FMRI processing. Here we describe procedures for performing QC on acquired or publicly available FMRI datasets using the widely used AFNI software package. This work is part of the Research Topic, "Demonstrating Quality Control (QC) Procedures in fMRI." We used a sequential, hierarchical approach that contained the following major stages: (1) GTKYD (getting to know your data, esp. its basic acquisition properties), (2) APQUANT (examining quantifiable measures, with thresholds), (3) APQUAL (viewing qualitative images, graphs, and other information in systematic HTML reports) and (4) GUI (checking features interactively with a graphical user interface); and for task data, and (5) STIM (checking stimulus event timing statistics). We describe how these are complementary and reinforce each other to help researchers stay close to their data. We processed and evaluated the provided, publicly available resting state data collections (7 groups, 139 total subjects) and task-based data collection (1 group, 30 subjects). As specified within the Topic guidelines, each subject's dataset was placed into one of three categories: Include, exclude or uncertain. The main focus of this paper, however, is the detailed description of QC procedures: How to understand the contents of an FMRI dataset, to check its contents for appropriateness, to verify processing steps, and to examine potential quality issues. Scripts for the processing and analysis are freely available.

Quality control in functional MRI studies with MRIQC and fMRIPrep

The implementation of adequate quality assessment (QA) and quality control (QC) protocols within the magnetic resonance imaging (MRI) research workflow is resource- and time-consuming and even more so is their execution. As a result, QA/QC practices highly vary across laboratories and "MRI schools", ranging from highly specialized knowledge spots to environments where QA/QC is considered overly onerous and costly despite evidence showing that below-standard data increase the false positive and false negative rates of the final results. Here, we demonstrate a protocol based on the visual assessment of images one-by-one with reports generated by MRIQC and fMRIPrep, for the QC of data in functional (blood-oxygen dependent-level; BOLD) MRI analyses. We particularize the proposed, open-ended scope of application to whole-brain voxel-wise analyses of BOLD to correspondingly enumerate and define the exclusion criteria applied at the QC checkpoints. We apply our protocol on a composite dataset (n = 181 subjects) drawn from open fMRI studies, resulting in the exclusion of 97% of the data (176 subjects). This high exclusion rate was expected because subjects were selected to showcase artifacts. We describe the artifacts and defects more commonly found in the dataset that justified exclusion. We moreover release all the materials we generated in this assessment and document all the QC decisions with the expectation of contributing to the standardization of these procedures and engaging in the discussion of QA/QC by the community.

Demonstrating quality control procedures for fMRI in DPABI

Quality control (QC) is an important stage for functional magnetic resonance imaging (fMRI) studies. The methods for fMRI QC vary in different fMRI preprocessing pipelines. The inflating sample size and number of scanning sites for fMRI studies further add to the difficulty and working load of the QC procedure. Therefore, as a constituent part of the Demonstrating Quality Control Procedures in fMRI research topic in Frontiers, we preprocessed a well-organized open-available dataset using DPABI pipelines to illustrate the QC procedure in DPABI. Six categories of DPABI-derived reports were used to eliminate images without adequate quality. After the QC procedure, twelve participants (8.6%) were categorized as excluded and eight participants (5.8%) were categorized as uncertain. More automatic QC tools were needed in the big-data era while visually inspecting images was still indispensable now.

A comprehensive macaque fMRI pipeline and hierarchical atlas

Functional neuroimaging research in the non-human primate (NHP) has been advancing at a remarkable rate. The increase in available data establishes a need for robust analysis pipelines designed for NHP neuroimaging and accompanying template spaces to standardize the localization of neuroimaging results. Our group recently developed the NIMH Macaque Template (NMT), a high-resolution population average anatomical template and associated neuroimaging resources, providing researchers with a standard space for macaque neuroimaging . Here, we release NMT v2, which includes both symmetric and asymmetric templates in stereotaxic orientation, with improvements in spatial contrast, processing efficiency, and segmentation. We also introduce the Cortical Hierarchy Atlas of the Rhesus Macaque (CHARM), a hierarchical parcellation of the macaque cerebral cortex with varying degrees of detail. These tools have been integrated into the neuroimaging analysis software AFNI to provide a comprehensive and robust pipeline for fMRI processing, visualization and analysis of NHP data. AFNI's new @animal_warper program can be used to efficiently align anatomical scans to the NMT v2 space, and afni_proc.py integrates these results with full fMRI processing using macaque-specific parameters: from motion correction through regression modeling. Taken together, the NMT v2 and AFNI represent an all-in-one package for macaque functional neuroimaging analysis, as demonstrated with available demos for both task and resting state fMRI.

Different activation signatures in the primary sensorimotor and higher-level regions for haptic three-dimensional curved surface exploration

Haptic object perception begins with continuous exploratory contact, and the human brain needs to accumulate sensory information continuously over time. However, it is still unclear how the primary sensorimotor cortex (PSC) interacts with these higher-level regions during haptic exploration over time. This functional magnetic resonance imaging (fMRI) study investigates time-dependent haptic object processing by examining brain activity during haptic 3D curve and roughness estimations. For this experiment, we designed sixteen haptic stimuli (4 kinds of curves × 4 varieties of roughness) for the haptic curve and roughness estimation tasks. Twenty participants were asked to move their right index and middle fingers along the surface twice and to estimate one of the two features-roughness or curvature-depending on the task instruction. We found that the brain activity in several higher-level regions (e.g., the bilateral posterior parietal cortex) linearly increased as the number of curves increased during the haptic exploration phase. Surprisingly, we found that the contralateral PSC was parametrically modulated by the number of curves only during the late exploration phase but not during the early exploration phase. In contrast, we found no similar parametric modulation activity patterns during the haptic roughness estimation task in either the contralateral PSC or in higher-level regions. Thus, our findings suggest that haptic 3D object perception is processed across the cortical hierarchy, whereas the contralateral PSC interacts with other higher-level regions across time in a manner that is dependent upon the features of the object.