Statistical methods of estimation and inference for functional MR image analysis

Interpretability of Multivariate Brain Maps in Linear Brain Decoding: Definition, and Heuristic Quantification in Multivariate Analysis of MEG Time-Locked Effects

Brain decoding is a popular multivariate approach for hypothesis testing in neuroimaging. Linear classifiers are widely employed in the brain decoding paradigm to discriminate among experimental conditions. Then, the derived linear weights are visualized in the form of multivariate brain maps to further study spatio-temporal patterns of underlying neural activities. It is well known that the brain maps derived from weights of linear classifiers are hard to interpret because of high correlations between predictors, low signal to noise ratios, and the high dimensionality of neuroimaging data. Therefore, improving the interpretability of brain decoding approaches is of primary interest in many neuroimaging studies. Despite extensive studies of this type, at present, there is no formal definition for interpretability of multivariate brain maps. As a consequence, there is no quantitative measure for evaluating the interpretability of different brain decoding methods. In this paper, first, we present a theoretical definition of interpretability in brain decoding; we show that the interpretability of multivariate brain maps can be decomposed into their reproducibility and representativeness. Second, as an application of the proposed definition, we exemplify a heuristic for approximating the interpretability in multivariate analysis of evoked magnetoencephalography (MEG) responses. Third, we propose to combine the approximated interpretability and the generalization performance of the brain decoding into a new multi-objective criterion for model selection. Our results, for the simulated and real MEG data, show that optimizing the hyper-parameters of the regularized linear classifier based on the proposed criterion results in more informative multivariate brain maps. More importantly, the presented definition provides the theoretical background for quantitative evaluation of interpretability, and hence, facilitates the development of more effective brain decoding algorithms in the future.

INVESTIGATING DIFFERENCES IN BRAIN FUNCTIONAL NETWORKS USING HIERARCHICAL COVARIATE-ADJUSTED INDEPENDENT COMPONENT ANALYSIS

Human brains perform tasks via complex functional networks consisting of separated brain regions. A popular approach to characterize brain functional networks in fMRI studies is independent component analysis (ICA), which is a powerful method to reconstruct latent source signals from their linear mixtures. In many fMRI studies, an important goal is to investigate how brain functional networks change according to specific clinical and demographic variabilities. Existing ICA methods, however, cannot directly incorporate covariate effects in ICA decomposition. Heuristic post-ICA analysis to address this need can be inaccurate and inefficient. In this paper, we propose a hierarchical covariate-adjusted ICA (hc-ICA) model that provides a formal statistical framework for estimating covariate effects and testing differences between brain functional networks. Our method provides a more reliable and powerful statistical tool for evaluating group differences in brain functional networks while appropriately controlling for potential confounding factors. We present an analytically tractable EM algorithm to obtain maximum likelihood estimates of our model. We also develop a subspace-based approximate EM that runs significantly faster while retaining high accuracy. To test the differences in functional networks, we introduce a voxel-wise approximate inference procedure which eliminates the need of computationally expensive covariance matrix estimation and inversion. We demonstrate the advantages of our methods over the existing method via simulation studies. We apply our method to an fMRI study to investigate differences in brain functional networks associated with post-traumatic stress disorder (PTSD).

Potential pitfalls when denoising resting state fMRI data using nuisance regression

In resting state fMRI, it is necessary to remove signal variance associated with noise sources, leaving cleaned fMRI time-series that more accurately reflect the underlying intrinsic brain fluctuations of interest. This is commonly achieved through nuisance regression, in which the fit is calculated of a noise model of head motion and physiological processes to the fMRI data in a General Linear Model, and the "cleaned" residuals of this fit are used in further analysis. We examine the statistical assumptions and requirements of the General Linear Model, and whether these are met during nuisance regression of resting state fMRI data. Using toy examples and real data we show how pre-whitening, temporal filtering and temporal shifting of regressors impact model fit. Based on our own observations, existing literature, and statistical theory, we make the following recommendations when employing nuisance regression: pre-whitening should be applied to achieve valid statistical inference of the noise model fit parameters; temporal filtering should be incorporated into the noise model to best account for changes in degrees of freedom; temporal shifting of regressors, although merited, should be achieved via optimisation and validation of a single temporal shift. We encourage all readers to make simple, practical changes to their fMRI denoising pipeline, and to regularly assess the appropriateness of the noise model used. By negotiating the potential pitfalls described in this paper, and by clearly reporting the details of nuisance regression in future manuscripts, we hope that the field will achieve more accurate and precise noise models for cleaning the resting state fMRI time-series.

Comparison of fMRI analysis methods for heterogeneous BOLD responses in block design studies

A large number of fMRI studies have shown that the temporal dynamics of evoked BOLD responses can be highly heterogeneous. Failing to model heterogeneous responses in statistical analysis can lead to significant errors in signal detection and characterization and alter the neurobiological interpretation. However, to date it is not clear that, out of a large number of options, which methods are robust against variability in the temporal dynamics of BOLD responses in block-design studies. Here, we used rodent optogenetic fMRI data with heterogeneous BOLD responses and simulations guided by experimental data as a means to investigate different analysis methods' performance against heterogeneous BOLD responses. Evaluations are carried out within the general linear model (GLM) framework and consist of standard basis sets as well as independent component analysis (ICA). Analyses show that, in the presence of heterogeneous BOLD responses, conventionally used GLM with a canonical basis set leads to considerable errors in the detection and characterization of BOLD responses. Our results suggest that the 3rd and 4th order gamma basis sets, the 7th to 9th order finite impulse response (FIR) basis sets, the 5th to 9th order B-spline basis sets, and the 2nd to 5th order Fourier basis sets are optimal for good balance between detection and characterization, while the 1st order Fourier basis set (coherence analysis) used in our earlier studies show good detection capability. ICA has mostly good detection and characterization capabilities, but detects a large volume of spurious activation with the control fMRI data.

Characterizing and differentiating task-based and resting state fMRI signals via two-stage sparse representations

A relatively underexplored question in fMRI is whether there are intrinsic differences in terms of signal composition patterns that can effectively characterize and differentiate task-based or resting state fMRI (tfMRI or rsfMRI) signals. In this paper, we propose a novel two-stage sparse representation framework to examine the fundamental difference between tfMRI and rsfMRI signals. Specifically, in the first stage, the whole-brain tfMRI or rsfMRI signals of each subject were composed into a big data matrix, which was then factorized into a subject-specific dictionary matrix and a weight coefficient matrix for sparse representation. In the second stage, all of the dictionary matrices from both tfMRI/rsfMRI data across multiple subjects were composed into another big data-matrix, which was further sparsely represented by a cross-subjects common dictionary and a weight matrix. This framework has been applied on the recently publicly released Human Connectome Project (HCP) fMRI data and experimental results revealed that there are distinctive and descriptive atoms in the cross-subjects common dictionary that can effectively characterize and differentiate tfMRI and rsfMRI signals, achieving 100% classification accuracy. Moreover, our methods and results can be meaningfully interpreted, e.g., the well-known default mode network (DMN) activities can be recovered from the very noisy and heterogeneous aggregated big-data of tfMRI and rsfMRI signals across all subjects in HCP Q1 release.

fMRI single trial discovery of spatio-temporal brain activity patterns

There is growing interest in the description of short-lived patterns in the spatiotemporal cortical activity monitored via neuroimaging. Most traditional analysis methods, designed to estimate relatively long-term brain dynamics, are not always appropriate to capture these patterns. Here we introduce a novel data-driven approach for detecting short-lived fMRI brain activity patterns. Exploiting Density Peak Clustering (Rodriguez and Laio [2014]), our approach reveals well localized clusters by identifying and grouping together voxels whose time-series are similar, irrespective of their brain location, even when very short time windows (∼10 volumes) are used. The method, which we call Coherence Density Peak Clustering (CDPC), is first tested on simulated data and compared with a standard unsupervised approach for fMRI analysis, independent component analysis (ICA). CDPC identifies activated voxels with essentially no false-positives and proves more reliable than ICA, which is troubled by a number of false positives comparable to that of true positives. The reliability of the method is demonstrated on real fMRI data from a simple motor task, containing brief iterations of the same movement. The clusters identified are found in regions expected to be involved in the task, and repeat synchronously with the paradigm. The methodology proposed is especially suitable for the study of short-time brain dynamics and single trial experiments, where the event or task of interest cannot be repeated for the same subject, as happens, for instance, in problem-solving, learning and decision-making. A GUI implementation of our method is available for download at https://github.com/micheleallegra/CDPC. Hum Brain Mapp 38:1421-1437, 2017. © 2016 Wiley Periodicals, Inc.

Detecting Spatio-Temporal Modes in Multivariate Data by Entropy Field Decomposition

A new data analysis method that addresses a general problem of detecting spatio-temporal variations in multivariate data is presented. The method utilizes two recent and complimentary general approaches to data analysis, information field theory (IFT) and entropy spectrum pathways (ESP). Both methods reformulate and incorporate Bayesian theory, thus use prior information to uncover underlying structure of the unknown signal. Unification of ESP and IFT creates an approach that is non-Gaussian and non-linear by construction and is found to produce unique spatio-temporal modes of signal behavior that can be ranked according to their significance, from which space-time trajectories of parameter variations can be constructed and quantified. Two brief examples of real world applications of the theory to the analysis of data bearing completely different, unrelated nature, lacking any underlying similarity, are also presented. The first example provides an analysis of resting state functional magnetic resonance imaging (rsFMRI) data that allowed us to create an efficient and accurate computational method for assessing and categorizing brain activity. The second example demonstrates the potential of the method in the application to the analysis of a strong atmospheric storm circulation system during the complicated stage of tornado development and formation using data recorded by a mobile Doppler radar. Reference implementation of the method will be made available as a part of the QUEST toolkit that is currently under development at the Center for Scientific Computation in Imaging.

Semiparametric Bayesian models for human brain mapping

Functional magnetic resonance imaging (fMRI) has led to enormous progress in human brain mapping. Adequate analysis of the massive spatiotemporal data sets generated by this imaging technique, combining parametric and non-parametric components, imposes challenging problems in statistical modelling. Complex hierarchical Bayesian models in combination with computer-intensive Markov chain Monte Carlo inference are promising tools.

The purpose of this paper is twofold. First, it provides a review of general semiparametric Bayesian models for the analysis of fMRI data. Most approaches focus on important but separate temporal or spatial aspects of the overall problem, or they proceed by stepwise procedures. Therefore, as a second aim, we suggest a complete spatiotemporal model for analysing fMRI data within a unified semiparametric Bayesian framework. An application to data from a visual stimulation experiment illustrates our approach and demonstrates its computational feasibility.

Wavelets and statistical analysis of functional magnetic resonance images of the human brain

Wavelets provide an orthonormal basis for multiresolution analysis and decorrelation or ‘whitening’ of nonstationary time series and spatial processes. Wavelets are particularly well suited to analysis of biological signals and images, such as human brain imaging data, which often have fractal or scale-invariant properties. We briefly define some key properties of the discrete wavelet transform (DWT) and review its applications to statistical analysis of functional magnetic resonance imaging (fMRI) data. We focus on time series resampling by ‘wavestrapping’ of wavelet coefficients, methods for efficient linear model estimation in the wavelet domain, and wavelet-based methods for multiple hypothesis testing, all of which are somewhat simplified by the decorrelating property of the DWT.

Testing for signals with unknown location and scale in a χ2random field, with an application to fMRI

Siegmund and Worsley (1995) considered the problem of testing for signals with unknown location and scale in a Gaussian random field defined on ℝN. The test statistic was the maximum of a Gaussian random field in anN+1 dimensional ‘scale space’,Ndimensions for location and 1 dimension for the scale of a smoothing filter. Scale space is identical to a continuous wavelet transform with a kernel smoother as the wavelet, though the emphasis here is on signal detection rather than image compression or enhancement. Two methods were used to derive an approximate null distribution forN=2 andN=3: one based on the method of volumes of tubes, the other based on the expected Euler characteristic of the excursion set. The purpose of this paper is two-fold: to show how the latter method can be extended to higher dimensions, and to apply this more general result to χ2fields. The result of Siegmund and Worsley (1995) then follows as a special case. In this paper the results are applied to the problem of searching for activation in brain images obtained by functional magnetic resonance imaging (fMRI).

Testing for signals with unknown location and scale in a χ2 random field, with an application to fMRI

Siegmund and Worsley (1995) considered the problem of testing for signals with unknown location and scale in a Gaussian random field defined on ℝN. The test statistic was the maximum of a Gaussian random field in anN+1 dimensional ‘scale space’,Ndimensions for location and 1 dimension for the scale of a smoothing filter. Scale space is identical to a continuous wavelet transform with a kernel smoother as the wavelet, though the emphasis here is on signal detection rather than image compression or enhancement. Two methods were used to derive an approximate null distribution forN=2 andN=3: one based on the method of volumes of tubes, the other based on the expected Euler characteristic of the excursion set. The purpose of this paper is two-fold: to show how the latter method can be extended to higher dimensions, and to apply this more general result to χ2fields. The result of Siegmund and Worsley (1995) then follows as a special case. In this paper the results are applied to the problem of searching for activation in brain images obtained by functional magnetic resonance imaging (fMRI).

Progressive Bidirectional Age-Related Changes in Default Mode Network Effective Connectivity across Six Decades

The default mode network (DMN) is a set of regions that is tonically engaged during the resting state and exhibits task-related deactivation that is readily reproducible across a wide range of paradigms and modalities. The DMN has been implicated in numerous disorders of cognition and, in particular, in disorders exhibiting age-related cognitive decline. Despite these observations, investigations of the DMN in normal aging are scant. Here, we used blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) acquired during rest to investigate age-related changes in functional connectivity of the DMN in 120 healthy normal volunteers comprising six, 20-subject, decade cohorts (from 20–29 to 70–79). Structural equation modeling (SEM) was used to assess age-related changes in inter-regional connectivity within the DMN. SEM was applied both using a previously published, meta-analytically derived, node-and-edge model, and using exploratory modeling searching for connections that optimized model fit improvement. Although the two models were highly similar (only 3 of 13 paths differed), the sample demonstrated significantly better fit with the exploratory model. For this reason, the exploratory model was used to assess age-related changes across the decade cohorts. Progressive, highly significant changes in path weights were found in 8 (of 13) paths: four rising, and four falling (most changes were significant by the third or fourth decade). In all cases, rising paths and falling paths projected in pairs onto the same nodes, suggesting compensatory increases associated with age-related decreases. This study demonstrates that age-related changes in DMN physiology (inter-regional connectivity) are bidirectional, progressive, of early onset and part of normal aging.

Mixed Spectrum Analysis on fMRI Time-Series

Temporal autocorrelation present in functional magnetic resonance image (fMRI) data poses challenges to its analysis. The existing approaches handling autocorrelation in fMRI time-series often presume a specific model of autocorrelation such as an auto-regressive model. The main limitation here is that the correlation structure of voxels is generally unknown and varies in different brain regions because of different levels of neurogenic noises and pulsatile effects. Enforcing a universal model on all brain regions leads to bias and loss of efficiency in the analysis. In this paper, we propose the mixed spectrum analysis of the voxel time-series to separate the discrete component corresponding to input stimuli and the continuous component carrying temporal autocorrelation. A mixed spectral analysis technique based on M-spectral estimator is proposed, which effectively removes autocorrelation effects from voxel time-series and identify significant peaks of the spectrum. As the proposed method does not assume any prior model for the autocorrelation effect in voxel time-series, varying correlation structure among the brain regions does not affect its performance. We have modified the standard M-spectral method for an application on a spatial set of time-series by incorporating the contextual information related to the continuous spectrum of neighborhood voxels, thus reducing considerably the computation cost. Likelihood of the activation is predicted by comparing the amplitude of discrete component at stimulus frequency of voxels across the brain by using normal distribution and modeling spatial correlations among the likelihood with a conditional random field. We also demonstrate the application of the proposed method in detecting other desired frequencies.

Neural circuits underlying mother's voice perception predict social communication abilities in children

The human voice is a critical social cue, and listeners are extremely sensitive to the voices in their environment. One of the most salient voices in a child's life is mother's voice: Infants discriminate their mother's voice from the first days of life, and this stimulus is associated with guiding emotional and social function during development. Little is known regarding the functional circuits that are selectively engaged in children by biologically salient voices such as mother's voice or whether this brain activity is related to children's social communication abilities. We used functional MRI to measure brain activity in 24 healthy children (mean age, 10.2 y) while they attended to brief (<1 s) nonsense words produced by their biological mother and two female control voices and explored relationships between speech-evoked neural activity and social function. Compared to female control voices, mother's voice elicited greater activity in primary auditory regions in the midbrain and cortex; voice-selective superior temporal sulcus (STS); the amygdala, which is crucial for processing of affect; nucleus accumbens and orbitofrontal cortex of the reward circuit; anterior insula and cingulate of the salience network; and a subregion of fusiform gyrus associated with face perception. The strength of brain connectivity between voice-selective STS and reward, affective, salience, memory, and face-processing regions during mother's voice perception predicted social communication skills. Our findings provide a novel neurobiological template for investigation of typical social development as well as clinical disorders, such as autism, in which perception of biologically and socially salient voices may be impaired.

Spatial regression analysis of serial DTI for subject-specific longitudinal changes of neurodegenerative disease

Quantitative measurement of localized longitudinal changes in brain abnormalities at an individual level may offer critical information for disease diagnosis and treatment. The voxel-wise permutation-based method SPREAD/iSPREAD, which combines resampling and spatial regression of neighboring voxels, provides an effective and robust method for detecting subject-specific longitudinal changes within the whole brain, especially for longitudinal studies with a limited number of scans. As an extension of SPREAD/iSPREAD, we present a general method that facilitates analysis of serial Diffusion Tensor Imaging (DTI) measurements (with more than two time points) for testing localized changes in longitudinal studies. Two types of voxel-level test statistics (model-free test statistics, which measure intra-subject variability across time, and test statistics based on general linear model that incorporate specific lesion evolution models) were estimated and tested against the null hypothesis among groups of DTI data across time. The implementation and utility of the proposed statistical method were demonstrated by both Monte Carlo simulations and applications on clinical DTI data from human brain in vivo. By a design of test statistics based on the disease progression model, it was possible to apportion the true significant voxels attributed to the disease progression and those caused by underlying anatomical differences that cannot be explained by the model, which led to improvement in false positive (FP) control in the results. Extension of the proposed method to include other diseases or drug effect models, as well as the feasibility of global statistics, was discussed. The proposed statistical method can be extended to a broad spectrum of longitudinal studies with carefully designed test statistics, which helps to detect localized changes at the individual level.

•

A nonparametric method for detecting subject-specific localized changes in longitudinal DTI

•

Obtain sufficient statistical power even with limited scans available

•

Various voxel-level test statistics for hypothesis tests among groups across time

•

Excellent at controlling false positive ratio with the model-based test statistics

Interneurons contribute to the hemodynamic/metabolic response to epileptiform discharges

Interpretation of hemodynamic responses in epilepsy is hampered by an incomplete understanding of the underlying neurovascular coupling, especially the contributions of excitation and inhibition. We made simultaneous multimodal recordings of local field potentials (LFPs), firing of individual neurons, blood flow, and oxygen level in the somatosensory cortex of anesthetized rats. Epileptiform discharges induced by bicuculline injections were used to trigger large local events. LFP and blood flow were robustly coupled, as were LFP and tissue oxygen. In a parametric linear model, LFP and the baseline activities of cerebral blood flow and tissue partial oxygen tension contributed significantly to blood flow and oxygen responses. In an analysis of recordings from 402 neurons, blood flow/tissue oxygen correlated with the discharge of putative interneurons but not of principal cells. Our results show that interneuron activity is important in the vascular and metabolic responses during epileptiform discharges.

Social hierarchy modulates neural responses of empathy for pain

Recent evidence indicates that empathic responses to others' pain are modulated by various situational and individual factors. However, few studies have examined how empathy and underlying brain functions are modulated by social hierarchies, which permeate human society with an enormous impact on social behavior and cognition. In this study, social hierarchies were established based on incidental skill in a perceptual task in which all participants were mediumly ranked. Afterwards, participants were scanned with functional magnetic resonance imaging while watching inferior-status or superior-status targets receiving painful or non-painful stimulation. The results revealed that painful stimulation applied to inferior-status targets induced higher activations in the anterior insula (AI) and anterior medial cingulate cortex (aMCC), whereas these empathic brain activations were significantly attenuated in response to superior-status targets' pain. Further, this neural empathic bias to inferior-status targets was accompanied by stronger functional couplings of AI with brain regions important in emotional processing (i.e. thalamus) and cognitive control (i.e. middle frontal gyrus). Our findings indicate that emotional sharing with others' pain is shaped by relative positions in a social hierarchy such that underlying empathic neural responses are biased toward inferior-status compared with superior-status individuals.

A wavelet-based estimator of the degrees of freedom in denoised fMRI time series for probabilistic testing of functional connectivity and brain graphs

Connectome mapping using techniques such as functional magnetic resonance imaging (fMRI) has become a focus of systems neuroscience. There remain many statistical challenges in analysis of functional connectivity and network architecture from BOLD fMRI multivariate time series. One key statistic for any time series is its (effective) degrees of freedom, df, which will generally be less than the number of time points (or nominal degrees of freedom, N). If we know the df, then probabilistic inference on other fMRI statistics, such as the correlation between two voxel or regional time series, is feasible. However, we currently lack good estimators of df in fMRI time series, especially after the degrees of freedom of the “raw” data have been modified substantially by denoising algorithms for head movement. Here, we used a wavelet-based method both to denoise fMRI data and to estimate the (effective) df of the denoised process. We show that seed voxel correlations corrected for locally variable df could be tested for false positive connectivity with better control over Type I error and greater specificity of anatomical mapping than probabilistic connectivity maps using the nominal degrees of freedom. We also show that wavelet despiked statistics can be used to estimate all pairwise correlations between a set of regional nodes, assign a P value to each edge, and then iteratively add edges to the graph in order of increasing P. These probabilistically thresholded graphs are likely more robust to regional variation in head movement effects than comparable graphs constructed by thresholding correlations. Finally, we show that time-windowed estimates of df can be used for probabilistic connectivity testing or dynamic network analysis so that apparent changes in the functional connectome are appropriately corrected for the effects of transient noise bursts. Wavelet despiking is both an algorithm for fMRI time series denoising and an estimator of the (effective) df of denoised fMRI time series. Accurate estimation of df offers many potential advantages for probabilistically thresholding functional connectivity and network statistics tested in the context of spatially variant and non-stationary noise. Code for wavelet despiking, seed correlational testing and probabilistic graph construction is freely available to download as part of the BrainWavelet Toolbox at www.brainwavelet.org.

•

We lack good estimators of degrees of freedom (df) in denoised fMRI time series.

•

Wavelet despiking can simultaneously denoise and estimate effective df post-denoising.

•

We show Type I error controlled single-subject probabilistic inference for seed connectivity.

•

We describe new methods for building and probabilistically thresholding brain graphs.

•

We extend these methods to probabilistic inference for time-varying connectomics.

Response inhibition failure to visual stimuli paired with a "single-type" stressor in PTSD patients: an fMRI pilot study

Patients with post-traumatic stress disorder (PTSD) tend to misinterpret innocuous stimuli as potential threats, possibly due to a conditioning provoked by traumatic episodes. Previous neuroimaging evidence has shown an abnormal activation of the amygdala and prefrontal cortex in PTSD patients during fear conditioning and extinction. Nevertheless, the effects of a single-type adverse stressor on that circuit remain poorly explored. We tested the hypothesis that a single-type adverse episode is able to affect the prefrontal cortex and amygdala response to conditioned stimuli. To test this hypothesis, fMRI recordings were performed in PTSD patients and trauma-exposed controls during the observation of neutral and negative paired or non-paired pictures with an adverse stimulus by means of a single association. Results showed that left amygdala activation during negative reinforced stimuli was correlated with the score of PTSD clinical scale across all subjects. Furthermore, in the traumatized non-PTSD group, the activation of the dorso-medial prefrontal cortex and bilateral amygdala was lower during the observation of the reinforced (CS(+)) versus non-reinforced pictures (CS(-)) in response to emotionally negative stimuli. This was not the case in the PTSD patients. These results suggest that in PTSD patients, a single-episode conditioning unveils the failure of an inhibitory mechanism moderating the activity of the prefrontal cortex and amygdala in response to adverse and neutral stimuli.

Intersubject information mapping: revealing canonical representations of complex natural stimuli

Abstract Real-world time-continuous stimuli such as video promise greater naturalism for studies of brain function. However, modeling the stimulus variation is challenging and introduces a bias in favor of particular descriptive dimensions. Alternatively, we can look for brain regions whose signal is correlated between subjects, essentially using one subject to model another. Intersubject correlation mapping (ICM) allows us to find brain regions driven in a canonical manner across subjects by a complex natural stimulus. However, it requires a direct voxel-to-voxel match between the spatiotemporal activity patterns and is thus only sensitive to common activations sufficiently extended to match up in Talairach space (or in an alternative, e.g. cortical-surface-based, common brain space). Here we introduce the more general approach of intersubject information mapping (IIM). For each brain region, IIM determines how much information is shared between the subjects' local spatiotemporal activity patterns. We estimate the intersubject mutual information using canonical correlation analysis applied to voxels within a spherical searchlight centered on each voxel in turn. The intersubject information estimate is invariant to linear transforms including spatial rearrangement of the voxels within the searchlight. This invariance to local encoding will be crucial in exploring fine-grained brain representations, which cannot be matched up in a common space and, more fundamentally, might be unique to each individual – like fingerprints. IIM yields a continuous brain map, which reflects intersubject information in fine-grained patterns. Performed on data from functional magnetic resonance imaging (fMRI) of subjects viewing the same television show, IIM and ICM both highlighted sensory representations, including primary visual and auditory cortices. However, IIM revealed additional regions in higher association cortices, namely temporal pole and orbitofrontal cortex. These regions appear to encode the same information across subjects in their fine-grained patterns, although their spatial-average activation was not significantly correlated between subjects.

Beyond BOLD: optimizing functional imaging in stroke populations

Blood oxygenation level-dependent (BOLD) signal changes are often assumed to directly reflect neural activity changes. Yet the real relationship is indirect, reliant on numerous assumptions, and subject to several sources of noise. Deviations from the core assumptions of BOLD contrast functional magnetic resonance imaging (fMRI), and their implications, have been well characterized in healthy populations, but are frequently neglected in stroke populations. In addition to conspicuous local structural and vascular changes after stroke, there are many less obvious challenges in the imaging of stroke populations. Perilesional ischemic changes, remodeling in regions distant to lesion sites, and diffuse perfusion changes all complicate interpretation of BOLD signal changes in standard fMRI protocols. Most stroke patients are also older than the young populations on which assumptions of neurovascular coupling and the typical analysis pipelines are based. We present a review of the evidence to show that the basic assumption of neurovascular coupling on which BOLD-fMRI relies does not capture the complex changes arising from stroke, both pathological and recovery related. As a result, estimating neural activity using the canonical hemodynamic response function is inappropriate in a number of contexts. We review methods designed to better estimate neural activity in stroke populations. One promising alternative to event-related fMRI is a resting-state-derived functional connectivity approach. Resting-state fMRI is well suited to stroke populations because it makes no performance demands on patients and is capable of revealing network-based pathology beyond the lesion site.

Functional neuroanatomy of executive function after neonatal brain injury in adults who were born very preterm

Individuals who were born very preterm (VPT; <33 gestational weeks) are at risk of experiencing deficits in tasks involving executive function in childhood and beyond. In addition, the type and severity of neonatal brain injury associated with very preterm birth may exert differential effects on executive functioning by altering its neuroanatomical substrates. Here we addressed this question by investigating with functional magnetic resonance imaging (fMRI) the haemodynamic response during executive-type processing using a phonological verbal fluency and a working memory task in VPT-born young adults who had experienced differing degrees of neonatal brain injury. 12 VPT individuals with a history of periventricular haemorrhage and ventricular dilatation (PVH+VD), 17 VPT individuals with a history of uncomplicated periventricular haemorrhage (UPVH), 13 VPT individuals with no history of neonatal brain injury and 17 controls received an MRI scan whilst completing a verbal fluency task with two cognitive loads (‘easy’ and ‘hard’ letters). Two groups of VPT individuals (PVH+VD; n = 10, UPVH; n = 8) performed an n-back task with three cognitive loads (1-, 2-, 3-back). Results demonstrated that VPT individuals displayed hyperactivation in frontal, temporal, and parietal cortices and in caudate nucleus, insula and thalamus compared to controls, as demands of the verbal fluency task increased, regardless of type of neonatal brain injury. On the other hand, during the n-back task and as working memory load increased, the PVH+VD group showed less engagement of the frontal cortex than the UPVH group. In conclusion, this study suggests that the functional neuroanatomy of different executive-type processes is altered following VPT birth and that neural activation associated with specific aspects of executive function (i.e., working memory) may be particularly sensitive to the extent of neonatal brain injury.

Mapping tonotopic organization in human temporal cortex: representational similarity analysis in EMEG source space

A wide variety of evidence, from neurophysiology, neuroanatomy, and imaging studies in humans and animals, suggests that human auditory cortex is in part tonotopically organized. Here we present a new means of resolving this spatial organization using a combination of non-invasive observables (EEG, MEG, and MRI), model-based estimates of spectrotemporal patterns of neural activation, and multivariate pattern analysis. The method exploits both the fine-grained temporal patterning of auditory cortical responses and the millisecond scale temporal resolution of EEG and MEG. Participants listened to 400 English words while MEG and scalp EEG were measured simultaneously. We estimated the location of cortical sources using the MRI anatomically constrained minimum norm estimate (MNE) procedure. We then combined a form of multivariate pattern analysis (representational similarity analysis) with a spatiotemporal searchlight approach to successfully decode information about patterns of neuronal frequency preference and selectivity in bilateral superior temporal cortex. Observed frequency preferences in and around Heschl's gyrus matched current proposals for the organization of tonotopic gradients in primary acoustic cortex, while the distribution of narrow frequency selectivity similarly matched results from the fMRI literature. The spatial maps generated by this novel combination of techniques seem comparable to those that have emerged from fMRI or ECOG studies, and a considerable advance over earlier MEG results.

Healthy individuals treated with clomipramine: an fMRI study of brain activity during autobiographical recall of emotions

Various functional magnetic resonance imaging studies addressed the effects of antidepressant drugs on brain functioning in healthy subjects; however, none specifically investigated positive mood changes to antidepressant drug. Sixteen subjects with no personal or family history of psychiatric disorders were selected from an ongoing 4-week open trial of small doses of clomipramine. Follow-up interviews documented clear positive treatment effects in six subjects, with reduced irritability and tension in social interactions, improved decision making, higher self-confidence and brighter mood. These subjects were then included in a placebo-controlled confirmatory trial and were scanned immediately after 4 weeks of clomipramine use and again 4 weeks after the last dose of clomipramine. The functional magnetic resonance imaging (fMRI) scans were run during emotion-eliciting stimuli. Repeated-measures analysis of variance of brain activity patterns showed significant interactions between group and treatment status during induced irritability (P<0.005 cluster-based) but not during happiness. Individuals displaying a positive subjective response do clomipramine had higher frontoparietal cortex activity during irritability than during happiness and neutral emotion, and higher temporo-parieto-occipital cortex activity during irritability than during happiness. We conclude that antidepressants not only induce positive mood responses but also act upon autobiographical recall of negative emotions.

Permutation inference for the general linear model

Permutation methods can provide exact control of false positives and allow the use of non-standard statistics, making only weak assumptions about the data. With the availability of fast and inexpensive computing, their main limitation would be some lack of flexibility to work with arbitrary experimental designs. In this paper we report on results on approximate permutation methods that are more flexible with respect to the experimental design and nuisance variables, and conduct detailed simulations to identify the best method for settings that are typical for imaging research scenarios. We present a generic framework for permutation inference for complex general linear models (GLMS) when the errors are exchangeable and/or have a symmetric distribution, and show that, even in the presence of nuisance effects, these permutation inferences are powerful while providing excellent control of false positives in a wide range of common and relevant imaging research scenarios. We also demonstrate how the inference on GLM parameters, originally intended for independent data, can be used in certain special but useful cases in which independence is violated. Detailed examples of common neuroimaging applications are provided, as well as a complete algorithm - the "randomise" algorithm - for permutation inference with the GLM.

Group-wise FMRI activation detection on DICCCOL landmarks

Group-wise activation detection in task-based fMRI has been widely used because of its robustness to noises and its capacity to deal with variability of individual brains. However, current group-wise fMRI activation detection methods typically rely on the co-registration of individual brains' fMRI images, which has difficulty in dealing with the remarkable anatomic variation of different brains. As a consequence, the resulted misalignments could significantly degrade the required inter-subject correspondences, thus substantially reducing the sensitivity and specificity of group-wise fMRI activation detection. To deal with these challenges, this paper presents a novel approach to detecting group-wise fMRI activation on our recently developed and validated Dense Individualized and Common Connectivity-based Cortical Landmarks (DICCCOL). The basic idea here is that the first-level general linear model (GLM) analysis is first performed on the fMRI signal of each corresponding DICCCOL landmark in individual brain's own space, and then the estimated effect sizes of the same landmark from a group of subjects are statistically assessed with the mixed-effect model at the group level. Finally, the consistently activated DICCCOL landmarks are determined and declared in a group-wise fashion in response to external block-based stimuli. Our experimental results have demonstrated that the proposed approach can detect meaningful activations.

A wavelet method for modeling and despiking motion artifacts from resting-state fMRI time series

The impact of in-scanner head movement on functional magnetic resonance imaging (fMRI) signals has long been established as undesirable. These effects have been traditionally corrected by methods such as linear regression of head movement parameters. However, a number of recent independent studies have demonstrated that these techniques are insufficient to remove motion confounds, and that even small movements can spuriously bias estimates of functional connectivity. Here we propose a new data-driven, spatially-adaptive, wavelet-based method for identifying, modeling, and removing non-stationary events in fMRI time series, caused by head movement, without the need for data scrubbing. This method involves the addition of just one extra step, the Wavelet Despike, in standard pre-processing pipelines. With this method, we demonstrate robust removal of a range of different motion artifacts and motion-related biases including distance-dependent connectivity artifacts, at a group and single-subject level, using a range of previously published and new diagnostic measures. The Wavelet Despike is able to accommodate the substantial spatial and temporal heterogeneity of motion artifacts and can consequently remove a range of high and low frequency artifacts from fMRI time series, that may be linearly or non-linearly related to physical movements. Our methods are demonstrated by the analysis of three cohorts of resting-state fMRI data, including two high-motion datasets: a previously published dataset on children (N=22) and a new dataset on adults with stimulant drug dependence (N=40). We conclude that there is a real risk of motion-related bias in connectivity analysis of fMRI data, but that this risk is generally manageable, by effective time series denoising strategies designed to attenuate synchronized signal transients induced by abrupt head movements. The Wavelet Despiking software described in this article is freely available for download at www.brainwavelet.org.

Resting-state fMRI using passband balanced steady-state free precession

OBJECTIVE

Resting-state functional MRI (rsfMRI) has been increasingly used for understanding brain functional architecture. To date, most rsfMRI studies have exploited blood oxygenation level-dependent (BOLD) contrast using gradient-echo (GE) echo planar imaging (EPI), which can suffer from image distortion and signal dropout due to magnetic susceptibility and inherent long echo time. In this study, the feasibility of passband balanced steady-state free precession (bSSFP) imaging for distortion-free and high-resolution rsfMRI was investigated.

METHODS

rsfMRI was performed in humans at 3 T and in rats at 7 T using bSSFP with short repetition time (TR = 4/2.5 ms respectively) in comparison with conventional GE-EPI. Resting-state networks (RSNs) were detected using independent component analysis.

RESULTS AND SIGNIFICANCE

RSNs derived from bSSFP images were shown to be spatially and spectrally comparable to those derived from GE-EPI images with considerable intra- and inter-subject reproducibility. High-resolution bSSFP images corresponded well to the anatomical images, with RSNs exquisitely co-localized to the gray matter. Furthermore, RSNs at areas of severe susceptibility such as human anterior prefrontal cortex and rat piriform cortex were proved accessible. These findings demonstrated for the first time that passband bSSFP approach can be a promising alternative to GE-EPI for rsfMRI. It offers distortion-free and high-resolution RSNs and is potentially suited for high field studies.

Application of mean-shift clustering to blood oxygen level dependent functional MRI activation detection

Background

Functional magnetic resonance imaging (fMRI) analysis is commonly done with cross-correlation analysis (CCA) and the General Linear Model (GLM). Both CCA and GLM techniques, however, typically perform calculations on a per-voxel basis and do not consider relationships neighboring voxels may have. Clustered voxel analyses have then been developed to improve fMRI signal detections by taking advantages of relationships of neighboring voxels. Mean-shift clustering (MSC) is another technique which takes into account properties of neighboring voxels and can be considered for enhancing fMRI activation detection.

Methods

This study examines the adoption of MSC to fMRI analysis. MSC was applied to a Statistical Parameter Image generated with the CCA technique on both simulated and real fMRI data. The MSC technique was then compared with CCA and CCA plus cluster analysis. A range of kernel sizes were used to examine how the technique behaves.

Results

Receiver Operating Characteristic curves shows an improvement over CCA and Cluster analysis. False positive rates are lower with the proposed technique. MSC allows the use of a low intensity threshold and also does not require the use of a cluster size threshold, which improves detection of weak activations and highly focused activations.

Conclusion

The proposed technique shows improved activation detection for both simulated and real Blood Oxygen Level Dependent fMRI data. More detailed studies are required to further develop the proposed technique.

Direct cortical hemodynamic mapping of somatotopy of pig nostril sensation by functional near-infrared cortical imaging (fNCI)

Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique for the noninvasive monitoring of human brain activation states utilizing the coupling between neural activity and regional cerebral hemodynamics. Illuminators and detectors, together constituting optodes, are placed on the scalp, but due to the presence of head tissues, an inter-optode distance of more than 2.5cm is necessary to detect cortical signals. Although direct cortical monitoring with fNIRS has been pursued, a high-resolution visualization of hemodynamic changes associated with sensory, motor and cognitive neural responses directly from the cortical surface has yet to be realized. To acquire robust information on the hemodynamics of the cortex, devoid of signal complications in transcranial measurement, we devised a functional near-infrared cortical imaging (fNCI) technique. Here we demonstrate the first direct functional measurement of temporal and spatial patterns of cortical hemodynamics using the fNCI technique. For fNCI, inter-optode distance was set at 5mm, and light leakage from illuminators was prevented by a special optode holder made of a light-shielding rubber sheet. fNCI successfully detected the somatotopy of pig nostril sensation, as assessed in comparison with concurrent and sequential somatosensory-evoked potential (SEP) measurements on the same stimulation sites. Accordingly, the fNCI system realized a direct cortical hemodynamic measurement with a spatial resolution comparable to that of SEP mapping on the rostral region of the pig brain. This study provides an important initial step toward realizing functional cortical hemodynamic monitoring during neurosurgery of human brains.

Dopamine precursor depletion impairs structure and efficiency of resting state brain functional networks

Spatial patterns of functional connectivity derived from resting brain activity may be used to elucidate the topological properties of brain networks. Such networks are amenable to study using graph theory, which shows that they possess small world properties and can be used to differentiate healthy subjects and patient populations. Of particular interest is the possibility that some of these differences are related to alterations in the dopamine system. To investigate the role of dopamine in the topological organization of brain networks at rest, we tested the effects of reducing dopamine synthesis in 13 healthy subjects undergoing functional magnetic resonance imaging. All subjects were scanned twice, in a resting state, following ingestion of one of two amino acid drinks in a randomized, double-blind manner. One drink was a nutritionally balanced amino acid mixture, and the other was tyrosine and phenylalanine deficient. Functional connectivity between 90 cortical and subcortical regions was estimated for each individual subject under each dopaminergic condition. The lowered dopamine state caused the following network changes: reduced global and local efficiency of the whole brain network, reduced regional efficiency in limbic areas, reduced modularity of brain networks, and greater connection between the normally anti-correlated task-positive and default-mode networks. We conclude that dopamine plays a role in maintaining the efficient small-world properties and high modularity of functional brain networks, and in segregating the task-positive and default-mode networks. This article is part of the Special Issue Section entitled 'Neuroimaging in Neuropharmacology'.

Brain Imaging Analysis

The increasing availability of brain imaging technologies has led to intense neuroscientific inquiry into the human brain. Studies often investigate brain function related to emotion, cognition, language, memory, and numerous other externally induced stimuli as well as resting-state brain function. Studies also use brain imaging in an attempt to determine the functional or structural basis for psychiatric or neurological disorders and, with respect to brain function, to further examine the responses of these disorders to treatment. Neuroimaging is a highly interdisciplinary field, and statistics plays a critical role in establishing rigorous methods to extract information and to quantify evidence for formal inferences. Neuroimaging data present numerous challenges for statistical analysis, including the vast amounts of data collected from each individual and the complex temporal and spatial dependence present. We briefly provide background on various types of neuroimaging data and analysis objectives that are commonly targeted in the field. We present a survey of existing methods targeting these objectives and identify particular areas offering opportunities for future statistical contribution.

A hierarchical model for probabilistic independent component analysis of multi-subject fMRI studies

An important goal in fMRI studies is to decompose the observed series of brain images to identify and characterize underlying brain functional networks. Independent component analysis (ICA) has been shown to be a powerful computational tool for this purpose. Classic ICA has been successfully applied to single-subject fMRI data. The extension of ICA to group inferences in neuroimaging studies, however, is challenging due to the unavailability of a pre-specified group design matrix. Existing group ICA methods generally concatenate observed fMRI data across subjects on the temporal domain and then decompose multi-subject data in a similar manner to single-subject ICA. The major limitation of existing methods is that they ignore between-subject variability in spatial distributions of brain functional networks in group ICA. In this article, we propose a new hierarchical probabilistic group ICA method to formally model subject-specific effects in both temporal and spatial domains when decomposing multi-subject fMRI data. The proposed method provides model-based estimation of brain functional networks at both the population and subject level. An important advantage of the hierarchical model is that it provides a formal statistical framework to investigate similarities and differences in brain functional networks across subjects, for example, subjects with mental disorders or neurodegenerative diseases such as Parkinson's as compared to normal subjects. We develop an EM algorithm for model estimation where both the E-step and M-step have explicit forms. We compare the performance of the proposed hierarchical model with that of two popular group ICA methods via simulation studies. We illustrate our method with application to an fMRI study of Zen meditation.

Post-hoc power estimation for topological inference in fMRI

When analyzing functional MRI data, several thresholding procedures are available to account for the huge number of volume units or features that are tested simultaneously. The main focus of these methods is to prevent an inflation of false positives. However, this comes with a serious decrease in power and leads to a problematic imbalance between type I and type II errors. In this paper, we show how estimating the number of activated peaks or clusters enables one to estimate post-hoc how powerful the selection procedure performs. This procedure can be used in real studies as a diagnostics tool, and raises awareness on how much activation is potentially missed. The method is evaluated and illustrated using simulations and a real data example. Our real data example illustrates the lack of power in current fMRI research.

Language-specific cortical activation patterns for verbal fluency tasks in Japanese as assessed by multichannel functional near-infrared spectroscopy

In Japan, verbal fluency tasks are commonly utilized as a standard paradigm for neuropsychological testing of cognitive and linguistic abilities. The Japanese "letter fluency task" is a mora/letter fluency task based on the phonological and orthographical characteristics of the Japanese language. Whether there are similar activation patterns across languages or a Japanese-specific mora/letter fluency pattern is not certain. We investigated the neural correlates of overt mora/letter and category fluency tasks in healthy Japanese. The category fluency task activated the bilateral fronto-temporal language-related regions with left-superior lateralization, while the mora/letter fluency task led to wider activation including the inferior parietal regions (left and right supramarginal gyrus). Specific bilateral supramarginal activation during the mora/letter fluency task in Japanese was distinct from that of similar letter fluency tasks in syllable-alphabet-based languages: this might be due to the requirement of additional phonological processing and working memory, or due to increased cognitive load in general.

Correcting brain-wide correlation differences in resting-state FMRI

Brain function in "resting" state has been extensively studied with functional magnetic resonance imaging (FMRI). However, drawing valid inferences, particularly for group comparisons, is fraught with pitfalls. Differing levels of brain-wide correlations can confound group comparisons. Global signal regression (GSReg) attempts to reduce this confound and is commonly used, even though it differentially biases correlations over brain regions, potentially leading to false group differences. We propose to use average brain-wide correlations as a measure of global correlation (GCOR), and examine the circumstances under which it can be used to identify or correct for differences in global fluctuations. In the process, we show the bias induced by GSReg to be a function only of the data's covariance matrix, and use simulations to compare corrections with GCOR as covariate to GSReg under various scenarios. We find that unlike GSReg, GCOR is a conservative approach that can reduce global variations, while avoiding the introduction of false significant differences, as GSReg can. However, as with GSReg, one cannot escape the interaction effect between the grouping variable and GCOR covariate on effect size. While GCOR is a complementary measure for resting state-FMRI applicable to legacy data, it is a lesser substitute for proper level-I denoising. We also assess the applicability of GCOR to empirical data with motion-based subject grouping and compare group differences to those using GSReg. We find that, while GCOR reduced correlation differences between high and low movers, it is doubtful that motion was the sole driver behind the differences in the first place.

Human left ventral premotor cortex mediates matching of hand posture to object use

Visuomotor transformations for grasping have been associated with a fronto-parietal network in the monkey brain. The human homologue of the parietal monkey region (AIP) has been identified as the anterior part of the intraparietal sulcus (aIPS), whereas the putative human equivalent of the monkey frontal region (F5) is located in the ventral part of the premotor cortex (vPMC). Results from animal studies suggest that monkey F5 is involved in the selection of appropriate hand postures relative to the constraints of the task. In humans, the functional roles of aIPS and vPMC appear to be more complex and the relative contribution of each region to grasp selection remains uncertain. The present study aimed to identify modulation in brain areas sensitive to the difficulty level of tool object - hand posture matching. Seventeen healthy right handed participants underwent fMRI while observing pictures of familiar tool objects followed by pictures of hand postures. The task was to decide whether the hand posture matched the functional use of the previously shown object. Conditions were manipulated for level of difficulty. Compared to a picture matching control task, the tool object – hand posture matching conditions conjointly showed increased modulation in several left hemispheric regions of the superior and inferior parietal lobules (including aIPS), the middle occipital gyrus, and the inferior temporal gyrus. Comparison of hard versus easy conditions selectively modulated the left inferior frontal gyrus with peak activity located in its opercular part (Brodmann area (BA) 44). We suggest that in the human brain, vPMC/BA44 is involved in the matching of hand posture configurations in accordance with visual and functional demands.

Depth-compensated diffuse optical tomography enhanced by general linear model analysis and an anatomical atlas of human head

One of the main challenges in functional diffuse optical tomography (DOT) is to accurately recover the depth of brain activation, which is even more essential when differentiating true brain signals from task-evoked artifacts in the scalp. Recently, we developed a depth-compensated algorithm (DCA) to minimize the depth localization error in DOT. However, the semi-infinite model that was used in DCA deviated significantly from the realistic human head anatomy. In the present work, we incorporated depth-compensated DOT (DC-DOT) with a standard anatomical atlas of human head. Computer simulations and human measurements of sensorimotor activation were conducted to examine and prove the depth specificity and quantification accuracy of brain atlas-based DC-DOT. In addition, node-wise statistical analysis based on the general linear model (GLM) was also implemented and performed in this study, showing the robustness of DC-DOT that can accurately identify brain activation at the correct depth for functional brain imaging, even when co-existing with superficial artifacts.

Spatiotemporal Segregation of Neural Response to Auditory Stimulation: An fMRI Study Using Independent Component Analysis and Frequency-Domain Analysis

Although auditory processing has been widely studied with conventional parametric methods, there have been a limited number of independent component analysis (ICA) applications in this area. The purpose of this study was to examine spatiotemporal behavior of brain networks in response to passive auditory stimulation using ICA. Continuous broadband noise was presented binaurally to 19 subjects with normal hearing. ICA was performed to segregate spatial networks, which were subsequently classified according to their temporal relation to the stimulus using power spectrum analysis. Classification of separated networks resulted in 3 stimulus-activated, 9 stimulus-deactivated, 2 stimulus-neutral (stimulus-dependent but not correlated with the stimulation timing), and 2 stimulus-unrelated (fluctuations that did not follow the stimulus cycles) components. As a result of such classification, spatiotemporal subdivisions were observed in a number of cortical structures, namely auditory, cingulate, and sensorimotor cortices, where parts of the same cortical network responded to the stimulus with different temporal patterns. The majority of the classified networks seemed to comprise subparts of the known resting-state networks (RSNs); however, they displayed different temporal behavior in response to the auditory stimulus, indicating stimulus-dependent temporal segregation of RSNs. Only one of nine deactivated networks coincided with the “classic” default-mode network, suggesting the existence of a stimulus-dependent default-mode network, different from that commonly accepted.

Functional activation during reading in relation to the sulci of the angular gyrus region

Neurological studies suggest that the angular gyrus region of the inferior parietal lobule may be critical for reading. However, unambiguous demonstration of angular gyrus involvement from lesion and functional neuroimaging studies is lacking, partly because of the absence of detailed morphological descriptions of this region. On the basis of our recent anatomical examination of this region and a tightly controlled functional magnetic resonance imaging paradigm, the present investigation demonstrated reading-related activity in the region of the angular gyrus that lies between the central and posterior branches of the caudal superior temporal sulcus, namely cytoarchitectonic area PG. Analysis of functional connectivity showed increased functional coupling during reading of area PG with the language areas of Broca and Wernicke, and a region previously identified as the visual word form area. Thus, the parietal reading area has been precisely localized, and its interactions with other cortical areas during reading have been demonstrated.

Statistical analysis of high density diffuse optical tomography

High density diffuse optical tomography (HD-DOT) is a noninvasive neuroimaging modality with moderate spatial resolution and localization accuracy. Due to portability and wear-ability advantages, HD-DOT has the potential to be used in populations that are not amenable to functional magnetic resonance imaging (fMRI), such as hospitalized patients and young children. However, whereas the use of event-related stimuli designs, general linear model (GLM) analysis, and imaging statistics are standardized and routine with fMRI, such tools are not yet common practice in HD-DOT. In this paper we adapt and optimize fundamental elements of fMRI analysis for application to HD-DOT. We show the use of event-related protocols and GLM de-convolution analysis in un-mixing multi-stimuli event-related HD-DOT data. Statistical parametric mapping (SPM) in the framework of a general linear model is developed considering the temporal and spatial characteristics of HD-DOT data. The statistical analysis utilizes a random field noise model that incorporates estimates of the local temporal and spatial correlations of the GLM residuals. The multiple-comparison problem is addressed using a cluster analysis based on non-stationary Gaussian random field theory. These analysis tools provide access to a wide range of experimental designs necessary for the study of the complex brain functions. In addition, they provide a foundation for understanding and interpreting HD-DOT results with quantitative estimates for the statistical significance of detected activation foci.