Changes in the subject’s breathing rate or depth, such as a breath-hold challenge, can cause significant MRI signal changes. However, the response function that best models breath-holding-induced ...signal changes, as well as those resulting from a wider range of breathing variations including those occurring during rest, has not yet been determined. Respiration related signal changes appear to be slower than neuronally induced BOLD signal changes and are not modeled accurately using the typical hemodynamic response functions used in fMRI. In this study, we derive a new response function to model the average MRI signal changes induced by variations in the respiration volume (breath-to-breath changes in the respiration depth and rate). This was done by averaging the response to a series of single deep breaths performed once every 40 s amongst otherwise constant breathing. The new “respiration response function” consists of an early overshoot followed by a later undershoot (peaking at approximately 16 s), and accurately models the MRI signal changes resulting from breath-holding as well as cued depth and rate changes.
An artificial neural network with multiple hidden layers (known as a deep neural network, or DNN) was employed as a predictive model (DNNp) for the first time to predict emotional responses using ...whole-brain functional magnetic resonance imaging (fMRI) data from individual subjects. During fMRI data acquisition, 10 healthy participants listened to 80 International Affective Digital Sound stimuli and rated their own emotions generated by each sound stimulus in terms of the arousal, dominance, and valence dimensions. The whole-brain spatial patterns from a general linear model (i.e., beta-valued maps) for each sound stimulus and the emotional response ratings were used as the input and output for the DNNP, respectively. Based on a nested five-fold cross-validation scheme, the paired input and output data were divided into training (three-fold), validation (one-fold), and test (one-fold) data. The DNNP was trained and optimized using the training and validation data and was tested using the test data. The Pearson's correlation coefficients between the rated and predicted emotional responses from our DNNP model with weight sparsity optimization (mean ± standard error 0.52 ± 0.02 for arousal, 0.51 ± 0.03 for dominance, and 0.51 ± 0.03 for valence, with an input denoising level of 0.3 and a mini-batch size of 1) were significantly greater than those of DNN models with conventional regularization schemes including elastic net regularization (0.15 ± 0.05, 0.15 ± 0.06, and 0.21 ± 0.04 for arousal, dominance, and valence, respectively), those of shallow models including logistic regression (0.11 ± 0.04, 0.10 ± 0.05, and 0.17 ± 0.04 for arousal, dominance, and valence, respectively; average of logistic regression and sparse logistic regression), and those of support vector machine-based predictive models (SVMps; 0.12 ± 0.06, 0.06 ± 0.06, and 0.10 ± 0.06 for arousal, dominance, and valence, respectively; average of linear and non-linear SVMps). This difference was confirmed to be significant with a Bonferroni-corrected p-value of less than 0.001 from a one-way analysis of variance (ANOVA) and subsequent paired t-test. The weights of the trained DNNPs were interpreted and input patterns that maximized or minimized the output of the DNNPs (i.e., the emotional responses) were estimated. Based on a binary classification of each emotion category (e.g., high arousal vs. low arousal), the error rates for the DNNP (31.2% ± 1.3% for arousal, 29.0% ± 1.7% for dominance, and 28.6% ± 3.0% for valence) were significantly lower than those for the linear SVMP (44.7% ± 2.0%, 50.7% ± 1.7%, and 47.4% ± 1.9% for arousal, dominance, and valence, respectively) and the non-linear SVMP (48.8% ± 2.3%, 52.2% ± 1.9%, and 46.4% ± 1.3% for arousal, dominance, and valence, respectively), as confirmed by the Bonferroni-corrected p < 0.001 from the one-way ANOVA. Our study demonstrates that the DNNp model is able to reveal neuronal circuitry associated with human emotional processing – including structures in the limbic and paralimbic areas, which include the amygdala, prefrontal areas, anterior cingulate cortex, insula, and caudate. Our DNNp model was also able to use activation patterns in these structures to predict and classify emotional responses to stimuli.
•Deep neural network (DNN) was trained to predict the human emotion measured from fMRI.•The prediction performance of the DNN was superior to that of the support vector machine.•Weight representation and input pattern estimation were introduced to interpret the trained DNN.•Brain regions related to emotional processing were identified from the DNN interpretation.•Representations of neuronal activations were readily separable at the higher hidden layer of the DNN.
The human brain coordinates a wide variety of motor activities. On a large scale, the cortical motor system is topographically organized such that neighboring body parts are represented by ...neighboring brain areas. This homunculus-like somatotopic organization along the central sulcus has been observed using neuroimaging for large body parts such as the face, hands and feet. However, on a finer scale, invasive electrical stimulation studies show deviations from this somatotopic organization that suggest an organizing principle based on motor actions rather than body part moved. It has not been clear how the action-map organization principle of the motor cortex in the mesoscopic (sub-millimeter) regime integrates into a body map organization principle on a macroscopic scale (cm). Here we developed and applied advanced mesoscopic (sub-millimeter) fMRI and analysis methodology to non-invasively investigate the functional organization topography across columnar and laminar structures in humans. Compared to previous methods, in this study, we could capture locally specific blood volume changes across entire brain regions along the cortical curvature. We find that individual fingers have multiple mirrored representations in the primary motor cortex depending on the movements they are involved in. We find that individual digits have cortical representations up to 3 mm apart from each other arranged in a column-like fashion. These representations are differentially engaged depending on whether the digits’ muscles are used for different motor actions such as flexion movements, like grasping a ball or retraction movements like releasing a ball. This research provides a starting point for non-invasive investigation of mesoscale topography across layers and columns of the human cortex and bridges the gap between invasive electrophysiological investigations and large coverage non-invasive neuroimaging.
Display omitted
•A sub-millimeter fMRI method is developed to image neural microcircuitry in humans.•The method can capture large FOVs with thin slices for ‛columnar’ and ‛laminar’ mapping.•An analysis pipeline is developed to investigate topographical representations that have only been visible in animals so far.•Novel findings include a mirrored finger representation in the human motor cortex.
Aerobic activity is a powerful stimulus for improving mental health and for generating structural changes in the brain. We review the literature documenting these structural changes and explore ...exactly where in the brain these changes occur as well as the underlying substrates of the changes including neural, glial, and vasculature components. Aerobic activity has been shown to produce different types of changes in the brain. The presence of novel experiences or learning is an especially important component in how these changes are manifest. We also discuss the distinct time courses of structural brain changes with both aerobic activity and learning as well as how these effects might differ in diseased and elderly groups.
The brain is the body's largest energy consumer, even in the absence of demanding tasks. Electrophysiologists report on-going neuronal firing during stimulation or task in regions beyond those of ...primary relationship to the perturbation. Although the biological origin of consciousness remains elusive, it is argued that it emerges from complex, continuous whole-brain neuronal collaboration. Despite converging evidence suggesting the whole brain is continuously working and adapting to anticipate and actuate in response to the environment, over the last 20 y, task-based functional MRI (fMRI) have emphasized a localizationist view of brain function, with fMRI showing only a handful of activated regions in response to task/stimulation. Here, we challenge that view with evidence that under optimal noise conditions, fMRI activations extend well beyond areas of primary relationship to the task; and blood-oxygen level-dependent signal changes correlated with task-timing appear in over 95% of the brain for a simple visual stimulation plus attention control task. Moreover, we show that response shape varies substantially across regions, and that whole-brain parcellations based on those differences produce distributed clusters that are anatomically and functionally meaningful, symmetrical across hemispheres, and reproducible across subjects. These findings highlight the exquisite detail lying in fMRI signals beyond what is normally examined, and emphasize both the pervasiveness of false negatives, and how the sparseness of fMRI maps is not a result of localized brain function, but a consequence of high noise and overly strict predictive response models.
Earlier research in cats has shown that both cerebral blood volume (CBV) and cerebral blood flow (CBF) can be used to identify layer-dependent fMRI activation with spatial specificity superior to ...gradient-echo blood-oxygen-level-dependent (BOLD) contrast (Jin and Kim, 2008a). CBF contrast of perfusion fMRI at ultra-high field has not been widely applied in humans to measure laminar activity due to its low sensitivity, while CBV contrast for fMRI using vascular space occupancy (VASO) has been successfully used. However, VASO can be compromised by interference of blood in-flow effects and a temporally limited acquisition window around the blood-nulling time point. Here, we proposed to use DANTE (Delay Alternating with Nutation for Tailored Excitation) pulse trains combined with 3D-EPI to acquire an integrated VASO and perfusion (VAPER) contrast. The signal origin of the VAPER contrast was theoretically evaluated with respect to its CBV and CBF contributions using a four-compartment simulation model. The feasibility of VAPER to measure layer-dependent activity was empirically investigated in human primary motor cortex at 7 T. We demonstrated this new tool, with its highly specified functional layer profile, robust reproducibility, and improved sensitivity, to allow investigation of layer-specific cortical functions.
Wakefulness levels modulate estimates of functional connectivity (FC), and, if unaccounted for, can become a substantial confound in resting-state fMRI. Unfortunately, wakefulness is rarely monitored ...due to the need for additional concurrent recordings (e.g., eye tracking, EEG). Recent work has shown that strong fluctuations around 0.05Hz, hypothesized to be CSF inflow, appear in the fourth ventricle (FV) when subjects fall asleep, and that they correlate significantly with the global signal. The analysis of these fluctuations could provide an easy way to evaluate wakefulness in fMRI-only data and improve our understanding of FC during sleep. Here we evaluate this possibility using the 7T resting-state sample from the Human Connectome Project (HCP). Our results replicate the observation that fourth ventricle ultra-slow fluctuations (∼0.05Hz) with inflow-like characteristics (decreasing in intensity for successive slices) are present in scans during which subjects did not comply with instructions to keep their eyes open (i.e., drowsy scans). This is true despite the HCP data not being optimized for the detection of inflow-like effects. In addition, time-locked BOLD fluctuations of the same frequency could be detected in large portions of grey matter with a wide range of temporal delays and contribute in significant ways to our understanding of how FC changes during sleep. First, these ultra-slow fluctuations explain half of the increase in global signal that occurs during descent into sleep. Similarly, global shifts in FC between awake and sleep states are driven by changes in this slow frequency band. Second, they can influence estimates of inter-regional FC. For example, disconnection between frontal and posterior components of the Defulat Mode Network (DMN) typically reported during sleep were only detectable after regression of these ultra-slow fluctuations. Finally, we report that the temporal evolution of the power spectrum of these ultra-slow FV fluctuations can help us reproduce sample-level sleep patterns (e.g., a substantial number of subjects descending into sleep 3 minutes following scanning onset), partially rank scans according to overall drowsiness levels, and predict individual segments of elevated drowsiness (at 60 seconds resolution) with 71% accuracy.
•Most studies do not present all results of their analysis, hiding subthreshold ones.•Hiding results negatively affects the interpretation and understanding of the study.•Neuroimagers should present ...all results of their study, highlighting key ones.•Using the public NARPS data, we show several benefits of the "highlighting" approach.•The highlighting approach improves individual studies and meta-analyses.
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.
LayNii: A software suite for layer-fMRI Huber, Laurentius (Renzo); Poser, Benedikt A.; Bandettini, Peter A. ...
NeuroImage (Orlando, Fla.),
08/2021, Volume:
237
Journal Article
Peer reviewed
Open access
•A new software toolbox is introduced for layer-specific functional MRI: LayNii.•LayNii is a suite of command-line executable C++ programs for Linux, Windows, and macOS.•LayNii is designed for ...layer-fMRI data that suffer from SNR and coverage constraints.•LayNii performs layerification in the native voxel space of functional data.•LayNii performs layer-smoothing, GE-BOLD deveining, QA, and VASO analysis.
High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain ‘layerification’ and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.
Display omitted