Quality and denoising in real-time functional magnetic resonance imaging neurofeedback: A methods review

Neurofeedback training using real-time functional magnetic resonance imaging (rtfMRI-NF) allows subjects voluntary control of localised and distributed brain activity. It has sparked increased interest as a promising non-invasive treatment option in neuropsychiatric and neurocognitive disorders, although its efficacy and clinical significance are yet to be determined. In this work, we present the first extensive review of acquisition, processing and quality control methods available to improve the quality of the neurofeedback signal. Furthermore, we investigate the state of denoising and quality control practices in 128 recently published rtfMRI-NF studies. We found: (a) that less than a third of the studies reported implementing standard real-time fMRI denoising steps, (b) significant room for improvement with regards to methods reporting and (c) the need for methodological studies quantifying and comparing the contribution of denoising steps to the neurofeedback signal quality. Advances in rtfMRI-NF research depend on reproducibility of methods and results. Notably, a systematic effort is needed to build up evidence that disentangles the various mechanisms influencing neurofeedback effects. To this end, we recommend that future rtfMRI-NF studies: (a) report implementation of a set of standard real-time fMRI denoising steps according to a proposed COBIDAS-style checklist (https://osf.io/kjwhf/), (b) ensure the quality of the neurofeedback signal by calculating and reporting community-informed quality metrics and applying offline control checks and (c) strive to adopt transparent principles in the form of methods and data sharing and support of open-source rtfMRI-NF software. Code and data for reproducibility, as well as an interactive environment to explore the study data, can be accessed at https://github.com/jsheunis/quality-and-denoising-in-rtfmri-nf.

Multi-dimensional connectivity: a conceptual and mathematical review

The estimation of functional connectivity between regions of the brain, for example based on statistical dependencies between the time series of activity in each region, has become increasingly important in neuroimaging. Typically, multiple time series (e.g. from each voxel in fMRI data) are first reduced to a single time series that summarises the activity in a region of interest, e.g. by averaging across voxels or by taking the first principal component; an approach we call one-dimensional connectivity. However, this summary approach ignores potential multi-dimensional connectivity between two regions, and a number of recent methods have been proposed to capture such complex dependencies. Here we review the most common multi-dimensional connectivity methods, from an intuitive perspective, from a formal (mathematical) point of view, and through a number of simulated and real (fMRI and MEG) data examples that illustrate the strengths and weaknesses of each method. The paper is accompanied with both functions and scripts, which implement each method and reproduce all the examples.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

Hemodynamic responses are abnormal in isolated cervical dystonia

Neuroimaging studies using functional magnetic resonance imaging (fMRI), which measures brain activity by detecting the changes in blood oxygenation levels, are advancing our understanding of the pathophysiology of dystonia. Neurobiological disturbances in dystonia, however, may affect neurovascular coupling and impact the interpretability of fMRI studies. We evaluated here whether the hemodynamic response patterns during a behaviorally matched motor task are altered in isolated cervical dystonia (CD). Twenty-five CD patients and 25 healthy controls (HCs) underwent fMRI scanning during a paced finger tapping task (nondystonic task in patients). Imaging data were analyzed using a constrained principal component analysis-a statistical method that combines regression analysis and principal component analysis and enables the extraction of task-related functional networks and determination of the spatial and temporal hemodynamic response patterns associated with the task performance. Data from three patients and two controls were removed due to excessive movement. No significant differences in demographics or motor performance were observed. Three task-associated functional brain networks were identified. During task performance, reduced hemodynamic responses were seen in a sensorimotor network and in a network that included key nodes of the default mode, executive control and visual networks. During rest, reductions in hemodynamic responses were seen in the cognitive/visual network. Lower hemodynamic responses within the primary sensorimotor network in patients were correlated with the increased dystonia severity. Pathophysiological disturbances in isolated CD, such as alterations in inhibitory signaling and dopaminergic neurotransmission, may impact neurovascular coupling. Not accounting for hemodynamic response differences in fMRI studies of dystonia could lead to inaccurate results and interpretations.

MALINI (Machine Learning in NeuroImaging): A MATLAB toolbox for aiding clinical diagnostics using resting-state fMRI data

Resting-state functional Magnetic Resonance Imaging (rs-fMRI) has been extensively used for diagnostic classification because it does not require task compliance and is easier to pool data from multiple imaging sites, thereby increasing the sample size. A MATLAB-based toolbox called Machine Learning in NeuroImaging (MALINI) for feature extraction and disease classification is presented. The MALINI toolbox extracts functional and effective connectivity features from preprocessed rs-fMRI data and performs classification between healthy and disease groups using any of 18 popular and widely used machine learning algorithms that are based on diverse principles. A consensus classifier combining the power of multiple classifiers is also presented. The utility of the toolbox is illustrated by accompanying data consisting of resting-state functional connectivity features from healthy controls and subjects with various brain-based disorders: autism spectrum disorder from autism brain imaging data exchange (ABIDE), Alzheimer's disease and mild cognitive impairment from Alzheimer's disease neuroimaging initiative (ADNI), attention deficit hyperactivity disorder from ADHD-200, and post-traumatic stress disorder and post-concussion syndrome acquired in-house. Results of classification performed on the above datasets can be obtained from the main article titled “Supervised machine learning for diagnostic classification from large-scale neuroimaging datasets” [1]. The data was divided into homogeneous and heterogeneous splits, such that 80% could be used for training, model building and cross-validation, while the remaining 20% of the data could be used as a hold-out independent test data for replication of the classification performance, to ensure the robustness of the classifiers to population variance in image acquisition site and age of the sample.

Individual structural features constrain the mouse functional connectome

Whole brain dynamics intuitively depend upon the internal wiring of the brain; but to which extent the individual structural connectome constrains the corresponding functional connectome is unknown, even though its importance is uncontested. After acquiring structural data from individual mice, we virtualized their brain networks and simulated in silico functional MRI data. Theoretical results were validated against empirical awake functional MRI data obtained from the same mice. We demonstrate that individual structural connectomes predict the functional organization of individual brains. Using a virtual mouse brain derived from the Allen Mouse Brain Connectivity Atlas, we further show that the dominant predictors of individual structure-function relations are the asymmetry and the weights of the structural links. Model predictions were validated experimentally using tracer injections, identifying which missing connections (not measurable with diffusion MRI) are important for whole brain dynamics in the mouse. Individual variations thus define a specific structural fingerprint with direct impact upon the functional organization of individual brains, a key feature for personalized medicine.

Deterioration from healthy to mild cognitive impairment and Alzheimer's disease mirrored in corresponding loss of centrality in directed brain networks

OBJECTIVE

It is important to identify brain-based biomarkers that progressively deteriorate from healthy to mild cognitive impairment (MCI) to Alzheimer's disease (AD). Cortical thickness, amyloid-ß deposition, and graph measures derived from functional connectivity (FC) networks obtained using functional MRI (fMRI) have been previously identified as potential biomarkers. Specifically, in the latter case, betweenness centrality (BC), a nodal graph measure quantifying information flow, is reduced in both AD and MCI. However, all such reports have utilized BC calculated from undirected networks that characterize synchronization rather than information flow, which is better characterized using directed networks.

METHODS

Therefore, we estimated BC from directed networks using Granger causality (GC) on resting-state fMRI data (N = 132) to compare the following populations (p < 0.05, FDR corrected for multiple comparisons): normal control (NC), early MCI (EMCI), late MCI (LMCI) and AD. We used an additional metric called middleman power (MP), which not only characterizes nodal information flow as in BC, but also measures nodal power critical for information flow in the entire network.

RESULTS

MP detected more brain regions than BC that progressively deteriorated from NC to EMCI to LMCI to AD, as well as exhibited significant associations with behavioral measures. Additionally, graph measures obtained from conventional FC networks could not identify a single node, underscoring the relevance of GC.

CONCLUSION

Our findings demonstrate the superiority of MP over BC as well as GC over FC in our case. MP obtained from GC networks could serve as a potential biomarker for progressive deterioration of MCI and AD.

Advancing functional connectivity research from association to causation

Cognition and behavior emerge from brain network interactions, such that investigating causal interactions should be central to the study of brain function. Approaches that characterize statistical associations among neural time series-functional connectivity (FC) methods-are likely a good starting point for estimating brain network interactions. Yet only a subset of FC methods ('effective connectivity') is explicitly designed to infer causal interactions from statistical associations. Here we incorporate best practices from diverse areas of FC research to illustrate how FC methods can be refined to improve inferences about neural mechanisms, with properties of causal neural interactions as a common ontology to facilitate cumulative progress across FC approaches. We further demonstrate how the most common FC measures (correlation and coherence) reduce the set of likely causal models, facilitating causal inferences despite major limitations. Alternative FC measures are suggested to immediately start improving causal inferences beyond these common FC measures.

Transitions in information processing dynamics at the whole-brain network level are driven by alterations in neural gain

A key component of the flexibility and complexity of the brain is its ability to dynamically adapt its functional network structure between integrated and segregated brain states depending on the demands of different cognitive tasks. Integrated states are prevalent when performing tasks of high complexity, such as maintaining items in working memory, consistent with models of a global workspace architecture. Recent work has suggested that the balance between integration and segregation is under the control of ascending neuromodulatory systems, such as the noradrenergic system, via changes in neural gain (in terms of the amplification and non-linearity in stimulus-response transfer function of brain regions). In a previous large-scale nonlinear oscillator model of neuronal network dynamics, we showed that manipulating neural gain parameters led to a 'critical' transition in phase synchrony that was associated with a shift from segregated to integrated topology, thus confirming our original prediction. In this study, we advance these results by demonstrating that the gain-mediated phase transition is characterized by a shift in the underlying dynamics of neural information processing. Specifically, the dynamics of the subcritical (segregated) regime are dominated by information storage, whereas the supercritical (integrated) regime is associated with increased information transfer (measured via transfer entropy). Operating near to the critical regime with respect to modulating neural gain parameters would thus appear to provide computational advantages, offering flexibility in the information processing that can be performed with only subtle changes in gain control. Our results thus link studies of whole-brain network topology and the ascending arousal system with information processing dynamics, and suggest that the constraints imposed by the ascending arousal system constrain low-dimensional modes of information processing within the brain.

Increasing robustness of pairwise methods for effective connectivity in magnetic resonance imaging by using fractional moment series of BOLD signal distributions

Estimating causal interactions in the brain from functional magnetic resonance imaging (fMRI) data remains a challenging task. Multiple studies have demonstrated that all current approaches to determine direction of connectivity perform poorly when applied to synthetic fMRI datasets. Recent advances in this field include methods for pairwise inference, which involve creating a sparse connectome in the first step, and then using a classifier in order to determine the directionality of connection between every pair of nodes in the second step. In this work, we introduce an advance to the second step of this procedure, by building a classifier based on fractional moments of the BOLD distribution combined into cumulants. The classifier is trained on datasets generated under the dynamic causal modeling (DCM) generative model. The directionality is inferred based on statistical dependencies between the two-node time series, for example, by assigning a causal link from time series of low variance to time series of high variance. Our approach outperforms or performs as well as other methods for effective connectivity when applied to the benchmark datasets. Crucially, it is also more resilient to confounding effects such as differential noise level across different areas of the connectome.

Dynamics of Segregation and Integration in Directional Brain Networks: Illustration in Soldiers With PTSD and Neurotrauma

Brain functioning relies on various segregated/specialized neural regions functioning as an integrated-interconnected network (i.e., metastability). Various psychiatric and neurologic disorders are associated with aberrant functioning of these brain networks. In this study, we present a novel framework integrating the strength and temporal variability of metastability in brain networks. We demonstrate that this approach provides novel mechanistic insights which enables better imaging-based predictions. Using whole-brain resting-state fMRI and a graph-theoretic framework, we integrated strength and temporal-variability of complex-network properties derived from effective connectivity networks, obtained from 87 U.S. Army soldiers consisting of healthy combat controls (n = 28), posttraumatic stress disorder (PTSD; n = 17), and PTSD with comorbid mild-traumatic brain injury (mTBI; n = 42). We identified prefrontal dysregulation of key subcortical and visual regions in PTSD/mTBI, with all network properties exhibiting lower variability over time, indicative of poorer flexibility. Larger impairment in the prefrontal-subcortical pathway but not prefrontal-visual pathway differentiated comorbid PTSD/mTBI from the PTSD group. Network properties of the prefrontal-subcortical pathway also had significant association (R ² = 0.56) with symptom severity and neurocognitive performance; and were also found to possess high predictive ability (81.4% accuracy in classifying the disorders, explaining 66-72% variance in symptoms), identified through machine learning. Our framework explained 13% more variance in behaviors compared to the conventional framework. These novel insights and better predictions were made possible by our novel framework using static and time-varying network properties in our three-group scenario, advancing the mechanistic understanding of PTSD and comorbid mTBI. Our contribution has wide-ranging applications for network-level characterization of healthy brains as well as mental disorders.

Investigating the impact of autocorrelation on time-varying connectivity

In recent years, a number of studies have reported on the existence of time-varying functional connectivity (TVC) in resting-state functional magnetic resonance imaging (rs-fMRI) data. The sliding-window technique is currently one of the most commonly used methods to estimate TVC. Although previous studies have shown that autocorrelation can negatively impact estimates of static functional connectivity, its impact on TVC estimates is not well known at this time. In this paper, we show both theoretically and empirically that the existence of autocorrelation within a time series can inflate the sampling variability of TVC estimated using the sliding-window technique. This can in turn increase the risk of misinterpreting noise as true TVC and negatively impact subsequent estimation of whole-brain time-varying FC profiles, or "brain states". The latter holds as more variable input measures lead to more variable output measures in the state estimation procedure. Finally, we demonstrate that prewhitening the data prior to analysis can lower the variance of the estimated TVC and improve brain state estimation. These results suggest that careful consideration is required when making inferences on TVC.

A deconvolution algorithm for multi-echo functional MRI: Multi-echo Sparse Paradigm Free Mapping

This work introduces a novel algorithm for deconvolution of the BOLD signal in multi-echo fMRI data: Multi-echo Sparse Paradigm Free Mapping (ME-SPFM). Assuming a linear dependence of the BOLD percent signal change on the echo time (TE) and using sparsity-promoting regularized least squares estimation, ME-SPFM yields voxelwise time-varying estimates of the changes in the apparent transverse relaxation (ΔR₂^⁎) without prior knowledge of the timings of individual BOLD events. Our results in multi-echo fMRI data collected during a multi-task event-related paradigm at 3 Tesla demonstrate that the maps of R₂^⁎ changes obtained with ME-SPFM at the times of the stimulus trials show high spatial and temporal concordance with the activation maps and BOLD signals obtained with standard model-based analysis. This method yields estimates of ΔR₂^⁎ having physiologically plausible values. Owing to its ability to blindly detect events, ME-SPFM also enables us to map ΔR₂^⁎ associated with spontaneous, transient BOLD responses occurring between trials. This framework is a step towards deciphering the dynamic nature of brain activity in naturalistic paradigms, resting-state or experimental paradigms with unknown timing of the BOLD events.

Understanding white matter structural connectivity differences between cognitively impaired and nonimpaired active professional fighters

Long-term traumatic brain injury due to repeated head impacts (RHI) has been shown to be a risk factor for neurodegenerative disorders, characterized by a loss in cognitive performance. Establishing the correlation between changes in the white matter (WM) structural connectivity measures and neuropsychological test scores might help to identify the neural correlates of the scores that are used in daily clinical setting to investigate deficits due to repeated head blows. Hence, in this study, we utilized high angular diffusion MRI (dMRI) of 69 cognitively impaired and 70 nonimpaired active professional fighters from the Professional Fighters Brain Health Study, and constructed structural connectomes to understand: (a) whether there is a difference in the topological WM organization between cognitively impaired and nonimpaired active professional fighters, and (b) whether graph-theoretical measures exhibit correlations with neuropsychological scores in these groups. A dMRI derived structural connectome was constructed for every participant using brain regions defined in AAL atlas as nodes, and the product of fiber number and average fractional anisotropy of the tracts connecting the nodes as edges. Our study identified a topological WM reorganization due to RHI in fighters prone to cognitive decline that was correlated with neuropsychological scores. Furthermore, graph-theoretical measures were correlated differentially with neuropsychological scores between groups. We also found differentiated WM connectivity involving regions of hippocampus, precuneus, and insula within our cohort of cognitively impaired fighters suggesting that there is a discernible WM topological reorganization in fighters prone to cognitive decline.

Differential coupling between subcortical calcium and BOLD signals during evoked and resting state through simultaneous calcium fiber photometry and fMRI

Task based and resting state fMRI has been widely utilized to study brain functions. As the foundation of fMRI, the underlying neural basis of the BOLD signal has been extensively studied, but the detailed mechanism remains elusive, particularly during the resting state. To examine the neurovascular coupling, it is important to simultaneously record neural and vascular signals. Here we developed a novel setup of camera based, scalable simultaneous calcium fiber photometry and fMRI in rats. Using this setup, we recorded calcium signals of superior colliculus (SC) and lateral geniculate nucleus (LGN) and fMRI simultaneously during visual stimulation and the resting state. Our results revealed robust, region-specific coupling between calcium and BOLD signals in the task state and weaker, whole brain correlation in the resting state. Interestingly, the spatial specificity of such correlation in the resting state was improved upon regression of white matter, ventricle signals and global signals in fMRI data. Overall, our results suggest differential coupling of calcium and BOLD signals for subcortical regions between evoked and resting states, and the coupling relationship in the resting state was related with resting state BOLD preprocessing strategies.

Modeling of the BOLD signal using event-related simultaneous EEG-fMRI and convolutional sparse coding analysis

In this work, we employ simultaneous EEG-fMRI data acquired during a visually-guided attention task along with convolutional sparse coding (CSC) analysis to extract transient events from the EEG. Subsequently, we use these events in a standard voxel-wise fMRI analysis and compare the resultant activation maps with maps obtained using the subjects' response time (RT) in detection of visual target stimuli. We also employ FIR models to obtain HRF estimates using the detected CSC events. Our results show concordance between the resultant activation maps and consistent HRF shapes for most of the subjects, suggesting that CSC can be used as a tool for the detection of reliable events in the EEG.

A Functional Data Method for Causal Dynamic Network Modeling of Task-Related fMRI

Functional MRI (fMRI) is a popular approach to investigate brain connections and activations when human subjects perform tasks. Because fMRI measures the indirect and convoluted signals of brain activities at a lower temporal resolution, complex differential equation modeling methods (e.g., Dynamic Causal Modeling) are usually employed to infer the neuronal processes and to fit the resulting fMRI signals. However, this modeling strategy is computationally expensive and remains to be mostly a confirmatory or hypothesis-driven approach. One major statistical challenge here is to infer, in a data-driven fashion, the underlying differential equation models from fMRI data. In this paper, we propose a causal dynamic network (CDN) method to estimate brain activations and connections simultaneously. Our method links the observed fMRI data with the latent neuronal states modeled by an ordinary differential equation (ODE) model. Using the basis function expansion approach in functional data analysis, we develop an optimization-based criterion that combines data-fitting errors and ODE fitting errors. We also develop and implement a block coordinate-descent algorithm to compute the ODE parameters efficiently. We illustrate the numerical advantages of our approach using data from realistic simulations and two task-related fMRI experiments. Compared with various effective connectivity methods, our method achieves higher estimation accuracy while improving the computational speed by from tens to thousands of times. Though our method is developed for task-related fMRI, we also demonstrate the potential applicability of our method (with a simple modification) to resting-state fMRI, by analyzing both simulated and real data from medium-sized networks.

Robust Recovery of Temporal Overlap Between Network Activity Using Transient-Informed Spatio-Temporal Regression

Functional magnetic resonance imaging is a non-invasive tomographic imaging modality that has provided insights into system-level brain function. New analysis methods are emerging to study the dynamic behavior of brain activity. The innovation-driven co-activation pattern (iCAP) approach is one such approach that relies on the detection of timepoints with a significant transient activity to subsequently retrieve spatially and temporally overlapping large-scale brain networks. To recover temporal profiles of the iCAPs for further time-resolved analysis, spatial patterns are fitted back to the activity-inducing signals. In this crucial step, spatial dependences can hinder the recovery of temporal overlapping activity. To overcome this effect, we propose a novel back-projection method that optimally fits activity-inducing signals given a set of transient timepoints and spatial maps of iCAPs, thus taking into account both spatial and temporal constraints. Validation on simulated data shows that transient-based constraints improve the quality of fitted time courses. Further evaluation on experimental data demonstrates that overfitting and underfitting are prevented by the use of optimized spatio-temporal constraints. Spatial and temporal properties of resulting iCAPs support that brain activity is characterized by the recurrent co-activation and co-deactivation of spatially overlapping large-scale brain networks. This new approach opens new avenues to explore the brain's dynamic core.

Investigating the Correspondence of Clinical Diagnostic Grouping With Underlying Neurobiological and Phenotypic Clusters Using Unsupervised Machine Learning

Many brain-based disorders are traditionally diagnosed based on clinical interviews and behavioral assessments, which are recognized to be largely imperfect. Therefore, it is necessary to establish neuroimaging-based biomarkers to improve diagnostic precision. Resting-state functional magnetic resonance imaging (rs-fMRI) is a promising technique for the characterization and classification of varying disorders. However, most of these classification methods are supervised, i.e., they require a priori clinical labels to guide classification. In this study, we adopted various unsupervised clustering methods using static and dynamic rs-fMRI connectivity measures to investigate whether the clinical diagnostic grouping of different disorders is grounded in underlying neurobiological and phenotypic clusters. In order to do so, we derived a general analysis pipeline for identifying different brain-based disorders using genetic algorithm-based feature selection, and unsupervised clustering methods on four different datasets; three of them-ADNI, ADHD-200, and ABIDE-which are publicly available, and a fourth one-PTSD and PCS-which was acquired in-house. Using these datasets, the effectiveness of the proposed pipeline was verified on different disorders: Attention Deficit Hyperactivity Disorder (ADHD), Alzheimer's Disease (AD), Autism Spectrum Disorder (ASD), Post-Traumatic Stress Disorder (PTSD), and Post-Concussion Syndrome (PCS). For ADHD and AD, highest similarity was achieved between connectivity and phenotypic clusters, whereas for ASD and PTSD/PCS, highest similarity was achieved between connectivity and clinical diagnostic clusters. For multi-site data (ABIDE and ADHD-200), we report site-specific results. We also reported the effect of elimination of outlier subjects for all four datasets. Overall, our results suggest that neurobiological and phenotypic biomarkers could potentially be used as an aid by the clinician, in additional to currently available clinical diagnostic standards, to improve diagnostic precision. Data and source code used in this work is publicly available at https://github.com/xinyuzhao/identification-of-brain-based-disorders.git.

Zinc Nanoparticles Enhance Brain Connectivity in the Canine Olfactory Network: Evidence From an fMRI Study in Unrestrained Awake Dogs

Prior functional Magnetic Resonance Imaging (fMRI) studies have indicated increased neural activation when zinc nanoparticles are added to odorants in canines. Here we demonstrate that zinc nanoparticles up-regulate directional brain connectivity in parts of the canine olfactory network. This provides an explanation for previously reported enhancement in the odor detection capability of the dogs in the presence of zinc nanoparticles. In this study, we obtained fMRI data from awake and unrestrained dogs while they were being exposed to odorants with and without zinc nanoparticles, zinc nanoparticles suspended in water vapor, as well as just water vapor alone. We obtained directional connectivity between the brain regions of the olfactory network that were significantly stronger for the condition of odorant + zinc nanoparticles compared to just odorants, water vapor + zinc nanoparticles and water vapor alone. We observed significant strengthening of the paths of the canine olfactory network in the presence of zinc nanoparticles. This result indicates that zinc nanoparticles could potentially be used to increase canine detection capabilities in the environments of very low concentrations of the odorants, which would have otherwise been undetected.

Estimated hemodynamic response function parameters obtained from resting state BOLD fMRI signals in subjects with autism spectrum disorder and matched healthy subjects

In Functional magnetic resonance imaging (fMRI), the blood oxygen level dependent (BOLD) signal is modeled as a convolution of the hemodynamic response function (HRF) and the unmeasured latent neural signal. Although most cortical and subcortical brain regions share the canonical shape of the HRF, the temporal structure of HRFs are variable across brain regions and subjects. This variability is induced by both neural and non-neural factors. The variability between subjects can be examined by three parameters that characterize the HRF: response height (RH), time-to-peak (TTP) and full-width at half-max (FWHM). This data provides three HRF parameters at every voxel, obtained from Autism Spectrum Disorder (ASD) patients (N = 531), and matched healthy controls (N = 571). Since ongoing studies suggest that non-standard populations have important differences in their HRFs when compared with healthy control, this data set is valuable in studying variability of HRF in ASD group and inferring the underlying pathology that also affects the HRF. It also has implications for fMRI analyses like resting-sate connectivity analysis.

Aberrant hemodynamic responses in autism: Implications for resting state fMRI functional connectivity studies

Functional MRI (fMRI) is modeled as a convolution of the hemodynamic response function (HRF) and an unmeasured latent neural signal. However, HRF itself is variable across brain regions and subjects. This variability is induced by both neural and non-neural factors. Aberrations in underlying neurochemical mechanisms, which control HRF shape, have been reported in autism spectrum disorders (ASD). Therefore, we hypothesized that this will lead to voxel-specific, yet systematic differences in HRF shape between ASD and healthy controls. As a corollary, we also hypothesized that such alterations will lead to differences in estimated functional connectivity in fMRI space compared to latent neural space. To test these hypotheses, we performed blind deconvolution of resting-state fMRI time series acquired from large number of ASD and control subjects obtained from the Autism Brain Imaging Data Exchange (ABIDE) database (N = 1102). Many brain regions previously implicated in autism showed systematic differences in HRF shape in ASD. Specifically, we found that precuneus had aberrations in all HRF parameters. Consequently, we obtained precuneus-seed-based functional connectivity differences between ASD and controls using fMRI as well as using latent neural signals. We found that non-deconvolved fMRI data failed to detect group differences in connectivity between precuneus and certain brain regions that were instead observed in deconvolved data. Our results are relevant for the understanding of hemodynamic and neurochemical aberrations in ASD, as well as have methodological implications for resting-state functional connectivity studies in Autism, and more generally in disorders that are accompanied by neurochemical alterations that may impact HRF shape.

Parameterized hemodynamic response function data of healthy individuals obtained from resting-state functional MRI in a 7T MRI scanner

Functional magnetic resonance imaging (fMRI), being an indirect measure of brain activity, is mathematically defined as a convolution of the unmeasured latent neural signal and the hemodynamic response function (HRF). The HRF is known to vary across the brain and across individuals, and it is modulated by neural as well as non-neural factors. Three parameters characterize the shape of the HRF, which is obtained by performing deconvolution on resting-state fMRI data: response height, time-to-peak and full-width at half-max. The data provided here, obtained from 47 healthy adults, contains these three HRF parameters at every voxel in the brain, as well as HRF parameters from the default-mode network (DMN). In addition, we have provided functional connectivity (FC) data from the same DMN regions, obtained for two cases: data with deconvolution (HRF variability minimized) and data with no deconvolution (HRF variability corrupted). This would enable researchers to compare regional changes in HRF with corresponding FC differences, to assess the impact of HRF variability on FC. Importantly, the data was obtained in a 7T MRI scanner. While most fMRI studies are conducted at lower field strengths, like 3T, ours is the first study to report HRF data obtained at 7T. FMRI data at ultra-high fields contains larger contributions from small vessels, consequently HRF variability is lower for small vessels at higher field strengths. This implies that findings made from this data would be more conservative than from data acquired at lower fields, such as 3T. Results obtained with this data and further interpretations are available in our recent research study (Rangaprakash et al., in press) [1]. This is a valuable dataset for studying HRF variability in conjunction with FC, and for developing the HRF profile in healthy individuals, which would have direct implications for fMRI data analysis, especially resting-state connectivity modeling. This is the first public HRF data at 7T.