Search

Diffusion MRI is a powerful tool for imaging brain structure, but it is challenging to discern the biological underpinnings of plasticity inferred from these and other non-invasive MR measurements. Biophysical modeling of the diffusion signal aims to render a more biologically rich image of tissue microstructure, but the application of these models comes with important caveats. A separate approach for gaining biological specificity has been to seek converging evidence from multi-modal datasets. Here we use metrics derived from diffusion kurtosis imaging (DKI) and the white matter tract integrity (WMTI) model along with quantitative MRI measurements of T1 relaxation to characterize changes throughout the white matter during an 8-week, intensive reading intervention (160 total hours of instruction). Behavioral measures, multi-shell diffusion MRI data, and quantitative T1 data were collected at regular intervals during the intervention in a group of 33 children with reading difficulties (7-12 years old), and over the same period in an age-matched non-intervention control group. Throughout the white matter, mean ‘extra-axonal’ diffusivity was inversely related to intervention time. In contrast, model estimated axonal water fraction (AWF), overall diffusion kurtosis, and T1 relaxation time showed no significant change over the intervention period. Both diffusion and quantitative T1 based metrics were correlated with pre-intervention reading performance, albeit with distinct anatomical distributions. These results are consistent with the view that rapid changes in diffusion properties reflect phenomena other than widespread changes in myelin density. We discuss this result in light of recent work highlighting non-axonal factors in experience-dependent plasticity and learning.HighlightsDiffusion MRI measurements in white matter show changes linked to an educational intervention.Tissue modeling results point to changes within the extra-axonal space.Complementary MRI measurements fail to suggest a widespread change in white matter in myelination over the intervention period.Both diffusion and quantitative T1 measures correlate with pre-intervention reading skill.

The ‘sensitive period’ for phonetic learning (∼6-12 months) is one of the earliest milestones in language acquisition where infants start to become specialized in processing speech sounds in their native language. In the last decade, advancements in neuroimaging technologies for infants are starting to shed light on the underlying neural mechanisms supporting this important learning period. The current study reports on the largest longitudinal dataset to date with the aim to replicate and extend on two important questions: 1) what are the developmental changes during the ‘sensitive period’ for native and nonnative speech processing? 2) how does native and nonnative speech processing in infants predict later language outcomes? Fifty-four infants were recruited at 7 months of age and their neural processing of speech was measured using Magnetoencephalography (MEG). Specifically, the neural sensitivity to a native and a nonnative speech contrast was indexed by the mismatch response (MMR). They repeated the measurement again at 11 months of age and their language development was further tracked from 12 months to 30 months of age using the MacArthur-Bates Communicative Development Inventory (CDI). Using an a prior Region-of-Interest (ROI) approach, we observed significant increases for the Native MMR in the left inferior frontal region (IF) and superior temporal region (ST) from 7 to 11 months, but not for the Nonnative MMR. Complementary whole brain comparison revealed more widespread developmental changes for both contrasts. However, only individual differences in the left IF and ST for the Nonnative MMR at 11 months of age were significant predictors of individual vocabulary growth up to 30 months of age. An exploratory machine-learning based analysis further revealed that whole brain MMR for both Native and Nonnative contrasts can robustly predict later outcomes, but with very different underlying spatial-temporal patterns. The current study extends our current knowledge and suggests that native and nonnative speech processing may follow different developmental trajectories and utilize different mechanisms that are relevant for later language skills.

Gestures are an integral part of in-person conversations and complement the meaning of the speech they accompany. The neural processing of co-speech gestures is supported by a mostly left-lateralized network of fronto-temporal regions. However, in contrast to iconic gestures, metaphoric as well as unrelated gestures have been found to more strongly engage the left and right inferior frontal gyrus (IFG), respectively. With this study, we conducted the first systematic comparison of all three types of gestures and resulting potential laterality effects. During collection of functional imaging data, 74 subjects were presented with 5 s videos of abstract speech with related metaphoric gestures, concrete speech with related iconic gestures and concrete speech with unrelated gestures. They were asked to judge whether the content of the speech and gesture matched or not. Differential contrasts revealed that both abstract related and concrete unrelated compared to concrete related stimuli elicited stronger activation of the bilateral IFG. Analyses of lateralization indices for IFG activation further showed a left hemispheric dominance for metaphoric gestures and a right hemispheric dominance for unrelated gestures. Our results give support to the hypothesis that the bilateral IFG is activated specifically when processing load for speech-gesture combinations is high. In addition, laterality effects indicate a stronger involvement of the right IFG in mismatch detection and conflict processing, whereas the left IFG performs the actual integration of information from speech and gesture.

In the last few years, a significant amount of work has aimed to characterize maturational trajectories of cortical development. The role of pericortical microstructure putatively characterized as the gray-white matter contrast (GWC) at the pericortical gray-white matter boundary and its relationship to more traditional morphological measures of cortical morphometry has emerged as a means to examine finer grained neuroanatomical underpinnings of cortical changes. In this work, we characterize the GWC developmental trajectories in a representative sample (n = 394) of children and adolescents (~4 to ~22 years of age), with repeated scans (1–3 scans per subject, total scans n = 819). We tested whether linear, quadratic, or cubic trajectories of contrast development best described changes in GWC. A best-fit model was identified vertex-wise across the whole cortex via the Akaike Information Criterion (AIC). GWC across nearly the whole brain was found to significantly change with age. Cubic trajectories were likeliest for 63% of vertices, quadratic trajectories were likeliest for 20% of vertices, and linear trajectories were likeliest for 16% of vertices. A main effect of sex was observed in some regions, where males had a higher GWC than females. However, no sex by age interactions were found on GWC. In summary, our results suggest a progressive decrease in GWC at the pericortical boundary throughout childhood and adolescence. This work contributes to efforts seeking to characterize typical, healthy brain development and, by extension, can help elucidate aberrant developmental trajectories.

Technological advancements in fMRI have afforded the opportunity to conduct submillimeter investigations into human brain function. The ability to do cortical depth dependent (or layer-specific) fMRI could allow probing intrinsic neuronal organizations and inter-connections, including the directionality of interregional information flow, while ultimately shedding light on uniquely human behaviors. The methodological development and applications of cortical depth dependent fMRI has been ongoing for nearly a decade, yet a consensus on protocols, analysis pipelines and interpretations of data has yet to be reached. In this article we discuss the current status and understandings of high resolution fMRI data and methods. In particular, the spatial signal characteristics of the BOLD contrast are examined in relation to the intra-cortical vasculature and associated underlying contrast mechanisms for three pulse sequences used in laminar fMRI: gradient echo (GE), spin echo (SE) and 3D-GRASE. We conclude with a short overview of factors to be taken into account when designing a specific study, including the choice of pulse sequence, experimental design, and data analysis strategy.

Previous research suggests that default-mode network (DMN) and dorsal attention network (DAN) are involved in internally- and externally-directed attention, respectively, through interactions with salience network (SN) and frontoparietal network (FPCN). Performing a task requiring external attention is often accompanied by a down-regulation of attention to internal thoughts, and vice-versa. In contrast, we often divide our attention between the external environment and internal thoughts to pursue distinct goals, yet virtually no prior research has examined how brain networks support this functionally critical neurocognitive process. In the current study, participants planned their responses for an upcoming alternate uses divergent thinking task (AUT-Condition), indicated whether arrows were pointing left or right (Arrows-Condition) or performed both tasks simultaneously (Dual-Task condition). Behaviorally, the Dual-Task condition was associated with equivalent generation of alternate uses but increased RT variability compared to the single-task conditions. Static connectivity analyses indicated that FPCN and SN increased their connectivity to DMN and reduced their connectivity to DAN during the Dual-Task condition and the AUT-Condition compared to the Arrows-Condition. Furthermore, DAN-SN connectivity was highest during the Arrows-Condition, intermediate during the Dual-Task condition and lowest during the AUT-Condition. Finally, time-varying connectivity analyses indicated that individuals who reported spending less time thinking of alternate uses during the Dual-Task condition spent more time in a state associated with performing the Arrows-Condition. Overall, our results suggest that interactions between DMN, FPCN, SN and DAN allow internal-external dual-tasking, and that time-varying functional connectivity between these networks is sensitive to attentional fluctuations between tasks during dual-tasking.

Functional magnetic resonance imaging (fMRI) using the blood oxygenation level-dependent (BOLD) contrast indirectly probes neuronal activity changes via evoked cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen (CMRO2) changes. The gradient-echo BOLD signal is mostly sensitive to ascending veins in the tissue and to pial veins. Thereby, the achievable spatial specificity to neuronal activation is limited. Furthermore, the non-linear interaction of CBF, CBV and CMRO2 can hamper quantitative interpretations of the BOLD signal across cortical depths with different baseline physiology. Measuring CBF, CBV or CMRO2 directly on a depth-dependent level has the potential to overcome these limitations. Here, we review these candidates of physiologically well-defined contrasts with the particular focus on arterial spin labeling (ASL), vascular space occupancy (VASO) and calibrated fMRI. These methods are reviewed with respect to their fMRI sequence parameter space and the applicability for neuroscientific studies in humans. We show representative results of depth-dependent 'non-BOLD-fMRI' in humans and their spatiotemporal characteristics. We conclude that non-BOLD methods are promising alternatives compared to conventional fMRI as they can provide improved spatial specificity, quantifiability and, hence, physiological interpretability as a function of cortical depth. At submillimeter resolution with inherently low signal-to-noise ratio (SNR), however, their use is still challenging. Nevertheless, we believe that 'non-BOLD-fMRI' is a useful alternative for depth-dependent investigations, by providing valuable insights into neurovascular coupling models that facilitate the interpretability of fMRI for neuroscientific applications.

BackgroundThe use of functional neuroimaging has been an extremely fruitful avenue for investigating the neural basis of human reward function. This approach has included identification of potential neurobiological mechanisms of psychiatric disease and examination of environmental, experiential, and biological factors that may contribute to disease risk via effects on the reward system. However, a central and largely unexamined assumption of much of this research is that neural reward function is an individual difference characteristic that is relatively stable over time.MethodsIn two independent samples of adolescents and young adults studied longitudinally (Ns= 145 & 153, 100% female & 100% male, ages 15-21 & 20-22, 2-4 scans & 2 scans respectively), we tested within-person stability of reward-task BOLD activation, with a median of 1 and 2 years between scans. We examined multiple commonly used contrasts of active states and baseline in both the anticipation and feedback phases of a card-guessing reward task. We examined the effects of cortical parcellation resolution, contrast, network (reward regions and resting-state networks), region-size, and activation strength and variability on the stability of reward-related activation.ResultsOverall, stability (ICC; intra-class correlation) across 1-2 years was modest. In both samples, contrasts of an active state relative to a baseline were more stable (e.g., Win>Baseline; mean ICC = 0.13 – 0.33) than contrasts of two active states (e.g., Win>Loss; mean ICC = 0.048 – 0.05). Additionally, activation in reward regions was less stable than in many non-task networks (e.g., dorsal attention), and activation in regions with greater between-subject variability showed higher stability in both samples.ConclusionsThese results show that functional neuroimaging activation to reward has modest stability over 1-2 years. Notably, results suggest that contrasts intended to map cognitive function and show robust group-level effects (i.e. Win > Loss) may be less effective in studies of individual differences and disease risk. The robustness of group-level activation should be weighed against other factors when selecting regions of interest in individual difference fMRI studies.

We address the problem of estimating how different parts of the brain develop and change throughout the lifespan, and how these trajectories are affected by genetic and environmental factors. Estimation of these lifespan trajectories is statistically challenging, since their shapes are typically highly nonlinear, and although true change can only be quantified by longitudinal examinations, as follow-up intervals in neuroimaging studies typically cover less than 10 \% of the lifespan, use of cross-sectional information is necessary. Linear mixed models (LMMs) and structural equation models (SEMs) commonly used in longitudinal analysis rely on assumptions which are typically not met with lifespan data, in particular when the data consist of observations combined from multiple studies. While LMMs require a priori specification of a polynomial functional form, SEMs do not easily handle data with unstructured time intervals between measurements. Generalized additive mixed models (GAMMs) offer an attractive alternative, and in this paper we propose various ways of formulating GAMMs for estimation of lifespan trajectories of 12 brain regions, using a large longitudinal dataset and realistic simulation experiments. We show that GAMMs are able to more accurately fit lifespan trajectories, distinguish longitudinal and cross-sectional effects, and estimate effects of genetic and environmental exposures. Finally, we discuss and contrast questions related to lifespan research which strictly require repeated measures data and questions which can be answered with a single measurement per participant, and in the latter case, which simplifying assumptions that need to be made. The examples are accompanied with R code, providing a tutorial for researchers interested in using GAMMs.

Brain networks are flexible and reconfigure over time to support ongoing cognitive processes. However, tracking statistically meaningful reconfigurations across time has proven difficult. This has to do largely with issues related to sampling variability, making instantaneous estimation of network organization difficult, along with increased reliance on task-free (cognitively unconstrained) experimental paradigms, limiting the ability to interpret the origin of changes in network structure over time. Here, we address these challenges using time-varying network analysis in conjunction with a naturalistic viewing paradigm. Specifically, we developed a measure of inter-subject network similarity and used this measure as a coincidence filter to identify synchronous fluctuations in network organization across individuals. Applied to movie-watching data, we found that periods of high inter-subject similarity coincided with reductions in network modularity and increased connectivity between cognitive systems. In contrast, low inter-subject similarity was associated with increased system segregation and more rest-like architectures. We then used a data-driven approach to uncover clusters of functional connections that follow similar trajectories over time and are more strongly correlated during movie-watching than at rest. Finally, we show that synchronous fluctuations in network architecture over time can be linked to a subset of features in the movie. Our findings link dynamic fluctuations in network integration and segregation to patterns of inter-subject similarity, and suggest that moment-to-moment fluctuations in functional connectivity reflect shared cognitive processing across individuals.

Rapidly and accurately processing information from faces is a critical human function that is known to improve with developmental age. Understanding the underlying drivers of this improvement remains a contentious question, with debate continuing as to the presence of early vs. late maturation of face-processing mechanisms. Recent behavioural evidence suggests an important ‘hallmark’ of expert face processing – the face inversion effect – is present in very young children, yet neural support for this remains unclear. To address this, we conducted a detailed investigation of the neural dynamics of face-selective processing in children spanning a range of ages (6 – 11 years) and adults. Uniquely, we applied multivariate pattern analysis (MVPA) to the electroencephalogram signal (EEG) to test for the presence of a distinct neural profile associated with canonical upright faces when compared both to other objects (houses) and to inverted faces. Results revealed robust discrimination profiles, at the individual level, of differentiated neural activity associated with broad face categorization and further with its expert processing, as indexed by the face inversion effect, from the youngest ages tested. This result is consistent with an early functional maturation of broad face processing mechanisms. Yet, clear quantitative differences between the response profile of children and adults is suggestive of age-related refinement of this system with developing face and general expertise. Standard ERP analysis also provides some support for qualitative differences in the neural response to inverted faces in children in contrast to adults. This neural profile is in line with recent behavioural studies that have reported impressively expert early face abilities during childhood, while also providing novel evidence of the ongoing neural specialisation between child and adulthood.

Reward Prediction Errors (RPEs), defined as the difference between the expected and received outcomes, are integral to reinforcement learning models and play an important role in development and psychopathology. In humans, RPE encoding can be estimated using fMRI recordings, however, a basic measurement property of RPE signals, their test-retest reliability across different time scales, remains an open question. In this paper, we examine the 3-month and 3-year reliability of RPE encoding in youth (mean age at baseline = 10.6 ± 0.3 years), a period of developmental transitions in reward processing. We show that RPE encoding is differentially distributed between the positive values being encoded predominantly in the striatum and negative RPEs primarily encoded in the insula. The encoding of negative RPE values is highly reliable in the right insula, across both the long and the short time intervals. Insula reliability for RPE encoding is the most robust finding, while other regions, such as the striatum, are less consistent. Striatal reliability appeared significant as well once covarying for factors, which were possibly confounding the signal to noise ratio. By contrast, task activation during feedback in the striatum is highly reliable across both time intervals. These results demonstrate the valence-dependent differential encoding of RPE signals between the insula and striatum, and the consistency of RPE signals or lack thereof, during childhood and into adolescence. Characterizing the regions where the RPE signal in BOLD fMRI is a reliable marker is key for estimating reward-processing alterations in longitudinal designs, such as developmental or treatment studies.

Multiple studies have shown that data quality is a critical confound in the construction of brain networks derived from functional MRI. This problem is particularly relevant for studies of human brain development where important variables (such as participant age) are correlated with data quality. Nevertheless, the impact of head motion on estimates of structural connectivity derived from diffusion tractography methods remains poorly characterized. Here, we evaluated the impact of in-scanner head motion on structural connectivity using a sample of 949 participants (ages 8-23 years old) who passed a rigorous quality assessment protocol for diffusion magnetic resonance imaging (dMRI) acquired as part of the Philadelphia Neurodevelopmental Cohort. Structural brain networks were constructed for each participant using both deterministic and probabilistic tractography. We hypothesized that subtle variation in head motion would systematically bias estimates of structural connectivity and confound developmental inference, as observed in previous studies of functional connectivity. Even following quality assurance and retrospective correction for head motion, eddy currents, and field distortions, in-scanner head motion significantly impacted the strength of structural connectivity in a consistency- and length-dependent manner. Specifically, increased head motion was associated with reduced estimates of structural connectivity for network edges with high inter-subject consistency, which included both short- and long-range connections. In contrast, motion inflated estimates of structural connectivity for low-consistency network edges that were primarily shorter-range. Finally, we demonstrate that age-related differences in head motion can both inflate and obscure developmental inferences on structural connectivity. Taken together, these data delineate the systematic impact of head motion on structural connectivity, and provide a critical context for identifying motion-related confounds in studies of structural brain network development.

Reports linking a 'jumping-to-conclusions' bias to delusions have led to growing interest in the neurobiological correlates of probabilistic reasoning. Several brain areas have been implicated in probabilistic reasoning; however, findings are difficult to integrate into a coherent account. The present study aimed to provide additional evidence by investigating, for the first time, effective connectivity among brain areas involved in different stages of evidence gathering. We investigated evidence gathering in 25 healthy individuals using fMRI and a new paradigm (Box Task) designed such as to minimize the effects of cognitive effort and reward processing. Decisions to collect more evidence ('draws') were contrasted to decisions to reach a final choice ('conclusions') with respect to BOLD activity. Psychophysiological interaction analysis was used to investigate effective connectivity. Conclusion events were associated with extensive brain activations in widely distributed brain areas associated with the task-positive network. In contrast, draw events were characterized by higher activation in areas assumed to be part of the task-negative network. Effective connectivity between the two networks decreased during draws and increased during conclusion events. Our findings indicate that probabilistic reasoning may depend on the balance between the task-positive and task-negative network, and that shifts in connectivity between the two may be crucial for evidence gathering. Thus, abnormal connectivity between the two systems may significantly contribute to the jumping-to-conclusions bias.

Large, longitudinal, multi-center MR neuroimaging studies require comprehensive quality assurance (QA) protocols for assessing the general quality of the compiled data, indicating potential malfunctions in the scanning equipment, and evaluating inter-site differences that need to be accounted for in subsequent analyses. We describe the implementation of a QA protocol for functional magnet resonance imaging (fMRI) data based on the regular measurement of an MRI phantom and an extensive variety of currently published QA statistics. The protocol is implemented in the MACS (Marburg-Münster Affective Disorders Cohort Study, http://for2107.de/), a two-center research consortium studying the neurobiological foundations of affective disorders. Between February 2015 and October 2016, 1214 phantom measurements have been acquired using a standard fMRI protocol. Using 444 healthy control subjects which have been measured between 2014 and 2016 in the cohort, we investigate the extent of between-site differences in contrast to the dependence on subject-specific covariates (age and sex) for structural MRI, fMRI, and diffusion tensor imaging (DTI) data. We show that most of the presented QA statistics differ severely not only between the two scanners used for the cohort but also between experimental settings (e.g. hardware and software changes), demonstrate that some of these statistics depend on external variables (e.g. time of day, temperature), highlight their strong dependence on proper handling of the MRI phantom, and show how the use of a phantom holder may balance this dependence. Site effects, however, do not only exist for the phantom data, but also for human MRI data. Using T1-weighted structural images, we show that total intracranial (TIV), grey matter (GMV), and white matter (WMV) volumes significantly differ between the MR scanners, showing large effect sizes. Voxel-based morphometry (VBM) analyses show that these structural differences observed between scanners are most pronounced in the bilateral basal ganglia, thalamus, and posterior regions. Using DTI data, we also show that fractional anisotropy (FA) differs between sites in almost all regions assessed. When pooling data from multiple centers, our data show that it is a necessity to account not only for inter-site differences but also for hardware and software changes of the scanning equipment. Also, the strong dependence of the QA statistics on the reliable placement of the MRI phantom shows that the use of a phantom holder is recommended to reduce the variance of the QA statistics and thus to increase the probability of detecting potential scanner malfunctions.

Visual expertise induces changes in neural processing for many different domains of expertise. However, it is unclear how expertise effects for different domains of expertise are related. In the present fMRI study, we combine large-scale univariate and multi-voxel analyses to contrast the expertise-related neural changes associated with two different domains of expertise, bird expertise (ornithology) and mineral expertise (mineralogy). Results indicated distributed expertise-related neural changes, with effects for both domains of expertise in high-level visual cortex and effects for bird expertise even extending to low-level visual regions and the frontal lobe. Importantly, a multivariate generalization analysis showed that effects in high-level visual cortex were specific to the domain of expertise. In contrast, the neural changes in the frontal lobe relating to expertise showed significant generalization, signaling the presence of domain-independent expertise effects. In conclusion, expertise is related to a combination of domain-specific and domain-general changes in neural processing. ispartof: NeuroImage vol:169 pages:80-93 ispartof: location:United States status: published

Following rapid technological advances, ultra-high field functional MRI (fMRI) enables exploring correlates of neuronal population activity at an increasing spatial resolution. However, as the fMRI blood-oxygenation-level-dependent (BOLD) contrast is a vascular signal, the spatial specificity of fMRI data is ultimately determined by the characteristics of the underlying vasculature. At 7T, fMRI measurement parameters determine the relative contribution of the macro- and microvasculature to the acquired signal. Here we investigate how these parameters affect relevant high-end fMRI analyses such as encoding, decoding, and submillimeter mapping of voxel preferences in the human auditory cortex. Specifically, we compare a T2* weighted fMRI dataset, obtained with 2D gradient echo (GE) EPI, to a predominantly T2 weighted dataset obtained with 3D GRASE. We first investigated the decoding accuracy based on two encoding models that represented different hypotheses about auditory cortical processing. This encoding/decoding analysis profited from the large spatial coverage and sensitivity of the T2* weighted acquisitions, as evidenced by a significantly higher prediction accuracy in the GE-EPI dataset compared to the 3D GRASE dataset for both encoding models. The main disadvantage of the T2* weighted GE-EPI dataset for encoding/decoding analyses was that the prediction accuracy exhibited cortical depth dependent vascular biases. However, we propose that the comparison of prediction accuracy across the different encoding models may be used as a post processing technique to salvage the spatial interpretability of the GE-EPI cortical depth-dependent prediction accuracy. Second, we explored the mapping of voxel preferences. Large-scale maps of frequency preference (i.e., tonotopy) were similar across datasets, yet the GE-EPI dataset was preferable due to its larger spatial coverage and sensitivity. However, submillimeter tonotopy maps revealed biases in assigned frequency preference and selectivity for the GE-EPI dataset, but not for the 3D GRASE dataset. Thus, a T2 weighted acquisition is recommended if high specificity in tonotopic maps is required. In conclusion, different fMRI acquisitions were better suited for different analyses. It is therefore critical that any sequence parameter optimization considers the eventual intended fMRI analyses and the nature of the neuroscience questions being asked.

Despite the multidimensional and temporally fleeting nature of auditory signals we quickly learn to assign novel sounds to behaviorally relevant categories. The neural systems underlying the learning and representation of novel auditory categories are far from understood. Current models argue for a rigid specialization of hierarchically organized core regions that are fine-tuned to extracting and mapping relevant auditory dimensions to meaningful categories. Scaffolded within a dual-learning systems approach, we test a competing hypothesis: the spatial and temporal dynamics of emerging auditory-category representations are not driven by the underlying dimensions but are constrained by category structure and learning strategies. To test these competing models, we used functional Magnetic Resonance Imaging (fMRI) to assess representational dynamics during the feedback-based acquisition of novel non-speech auditory categories with identical dimensions but differing category structures: rule-based (RB) categories, hypothesized to involve an explicit sound-to-rule mapping network, and information integration (II) based categories, involving pre-decisional integration of dimensions via a procedural-based sound-to-reward mapping network. Adults were assigned to either the RB (n = 30, 19 females) or II (n = 30, 22 females) learning tasks. Despite similar behavioral learning accuracies, learning strategies derived from computational modeling and involvements of corticostriatal systems during feedback processing differed across tasks. Spatiotemporal multivariate representational similarity analysis revealed an emerging representation within an auditory sensory-motor pathway exclusively for the II learning task, prominently involving the superior temporal gyrus (STG), inferior frontal gyrus (IFG), and posterior precentral gyrus. In contrast, the RB learning task yielded distributed neural representations within regions involved in cognitive-control and attentional processes that emerged at different time points of learning. Our results unequivocally demonstrate that auditory learners’ neural systems are highly flexible and show distinct spatial and temporal patterns that are not dimension-specific but reflect underlying category structures and learning strategies.

In magnetic resonance (MR) imaging, a lack of standardization in acquisition often causes pulse sequence-based contrast variations in MR images from site to site, which impedes consistent measurements in automatic analyses. In this paper, we propose an unsupervised MR image harmonization approach, CALAMITI (Contrast Anatomy Learning and Analysis for MR Intensity Translation and Integration), which aims to alleviate contrast variations in multi-site MR imaging. Designed using information bottleneck theory, CALAMITI learns a globally disentangled latent space containing both anatomical and contrast information, which permits harmonization. In contrast to supervised harmonization methods, our approach does not need a sample population to be imaged across sites. Unlike traditional unsupervised harmonization approaches which often suffer from geometry shifts, CALAMITI better preserves anatomy by design. The proposed method is also able to adapt to a new testing site with a straightforward fine-tuning process. Experiments on MR images acquired from ten sites show that CALAMITI achieves superior performance compared with other harmonization approaches.

Decades of research have led to several competing theories regarding the neural contributors to impaired reading. But how can we know which theory (or theories) identifies the types of markers that indeed differentiate between individuals with reading disabilities (RD) and their typically developing (TD) peers? To answer this question, we propose a new analytical tool for theory evaluation and comparison, grounded in the Bayesian latent-mixture modeling framework. We start by constructing a series of latent-mixture classification models, each reflecting one existing theoretical claim regarding the neurofunctional markers of RD (highlighting network-level differences in either mean activation, inter-subject heterogeneity, inter-region variability, or connectivity). Then, we run each model on fMRI data alone (i.e., while models are blind to participants' behavioral status), which enables us to interpret the fit between a model's classification of participants and their behavioral (known) RD/TD status as an estimate of its explanatory power. Results from n=127 adolescents and young adults (RD: n=59; TD: n=68) show that models based on network-level differences in mean activation and heterogeneity failed to differentiate between TD and RD individuals. In contrast, classifications based on variability and connectivity were significantly associated with participants' behavioral status. These findings suggest that differences in inter-region variability and connectivity may be better network-level markers of RD than mean activation or heterogeneity (at least in some populations and tasks). More broadly, the results demonstrate the promise of latent-mixture modeling as a theory-driven tool for evaluating different theoretical claims regarding neural contributors to language disorders and other cognitive traits.

The functional architecture of spontaneous BOLD fluctuations has been characterized in detail by numerous studies, demonstrating its potential relevance as a biomarker. However, the systematic investigation of its consistency is still in its infancy. Here, we analyze both the within- and between-subject variability as well as the test-retest reliability of resting-state functional connectivity (FC) estimates in a unique data set comprising multiple fMRI scans (42) from 5 subjects, and 50 single scans from 50 subjects. To this aim we adopted a statistical framework enabling us to disentangle the contribution of different sources of variability and their dependence on scan duration, and showed that the low reliability of single links can be largely improved using multiple scans per subject. Moreover, we show that practically all observed inter-region variability (at the link-level) is not significant and due to the statistical uncertainty of the estimator itself rather than to genuine variability among areas. Finally, we use the proposed statistical framework to demonstrate that, despite the poor consistency of single links, the information carried by the whole-brain spontaneous correlation structure is indeed robust, and can in fact be used as a functional fingerprint.

Longitudinal brain morphometry probes time-related brain morphometric patterns. We propose a method called dynamic network modeling with continuous valued nodes to generate a dynamic brain network from continuous valued longitudinal morphometric data. The mathematical framework of this method is based on state-space modeling. We use a bootstrap-enhanced least absolute shrinkage operator to solve the network-structure generation problem. In contrast to discrete dynamic Bayesian network modeling, the proposed method enables network generation directly from continuous valued high-dimensional short sequence data, being free from any discretization process. We applied the proposed method to a study of normal brain development.

While there is substantial overlap in the neural systems underlying empathy for people we know as opposed to strangers, social distance has been shown to significantly moderate empathic neural responses towards the negative experiences of others. Intriguingly however, variance in empathic neural responses towards known and unknown targets has not been reflected by behavioral differences as indexed by self-reported empathic ratings. One explanation for this disconnect is that empathic evaluations of known and unknown individuals draw on different bases (e.g. target identity/reactions) within the empathic process. To test this hypothesis, we utilized high density EEG to assess how individuating targets with personal names moderated the link between behavioral pain ratings and attentional processing oriented towards (a) initial target processing and (b) subsequent expressions target discomfort. Consistent with prior findings, no differences in pain ratings between individuated and unindividuated targets was observed. However, individual mean pain rating differences for individuated targets was strongly positively related to attentional processing levels, indexed by the P300, during the initial presentation of those targets, a relationship absent for unindividuated targets. In contrast, pain ratings for unindividuated targets was positively related to levels of attentional processing, indexed by the Late Positive Potential (LPP), during the subsequent discomfort expression stage. Furthermore, the LPP response to individuated target discomfort was positively linked to behavioral measures of emotional expressivity whereas the LPP response to unindividuated target discomfort was positively associated with cognitive appraisal. These findings suggest that individuation can significantly shift the bases of empathic responding.

Effort discounting theory suggests that the value of a reward should be lower if it was effortful to obtain, whereas contrast theory suggests that the contrast between the costly effort and the reward makes the reward seem more valuable. To test these alternative hypotheses, we used functional magnetic resonance imaging (fMRI) as participants engaged in feedback-based learning that required low or high cognitive effort to obtain positive feedback, while the objective amount of information provided by feedback remained constant. In the low effort condition, a single image was presented with four response options. In the high effort condition, two images were presented, each with two response options, and correct feedback was presented only when participants responded correctly to both of the images. Accuracy was significantly lower for the high effort condition, and all participants reported that the high effort condition was more difficult. A region of the ventral striatum selected for sensitivity to feedback value also showed increased activation to feedback presentation associated with the high effort condition relative to the low effort condition, when controlling for activation from corresponding control conditions where feedback was random. These results suggest that increased cognitive effort produces corresponding increases in positive feedback-related ventral striatum activity, in line with the predictions made by contrast theory. The accomplishment of obtaining a hard-earned intrinsic reward, such as positive feedback, may be particularly likely to promote reward-related brain activity.

Research on whether perception or other processes depend on the phase of neural oscillations is rapidly gaining popularity. However, it is unknown which methods are optimally suited to evaluate the hypothesized phase effect. Using a simulation approach, we here test the ability of different methods to detect such an effect on dichotomous (e.g., “hit” vs “miss”) and continuous (e.g., scalp potentials) response variables. We manipulated parameters that characterise the phase effect or define the experimental approach to test for this effect. For each parameter combination and response variable, we identified an optimal method. We found that methods regressing single-trial responses on circular (sine and cosine) predictors perform best for all of the simulated parameters, regardless of the nature of the response variable (dichotomous or continuous). In sum, our study lays a foundation for optimized experimental designs and analyses in future studies investigating the role of phase for neural and behavioural responses. We provide MATLAB code for the statistical methods tested., Highlights • The phase of neural oscillations modulates neural and behavioural responses. • We test the ability of different methods to detect such a phase effect. • We test how this ability changes with nature of effect and experimental parameters. • Results represent groundwork for future investigations of phase effects.

Encoding models based on deep convolutional neural networks (DCNN) predict BOLD responses to natural scenes in the human visual system more accurately than many other currently available models. However, DCNN-based encoding models fail to predict a significant amount of variance in the activity of most voxels in all visual areas. This failure could reflect limitations in the data (e.g., a noise ceiling), or could reflect limitations of the DCNN as a model of computation in the brain. Understanding the source and structure of the unexplained variance could therefore provide helpful clues for improving models of brain computation. Here, we characterize the structure of the variance that DCNN-based encoding models cannot explain. Using a publicly available dataset of BOLD responses to natural scenes, we determined if the source of unexplained variance was shared across voxels, individual brains, retinotopic locations, and hierarchically distant visual brain areas. We answered these questions using voxel-to-voxel (vox2vox) models that predict activity in a target voxel given activity in a population of source voxels. We found that simple linear vox2vox models increased within-subject prediction accuracy over DCNN-based models for any pair of source/target visual areas, clearly demonstrating that the source of unexplained variance is widely shared within and across visual brain areas. However, vox2vox models were not more accurate than DCNN-based encoding models when source and target voxels came from different brains, demonstrating that the source of unexplained variance was not shared across brains. Importantly, control analyses demonstrated that the source of unexplained variance was not encoded in the mean activity of source voxels, or the activity of voxels in white matter. Interestingly, the weights of vox2vox models revealed preferential connection of target voxel activity to source voxels with adjacent receptive fields, even when source and target voxels were in different functional brain areas. Finally, we found that the prediction accuracy of the vox2vox models decayed with hierarchical distance between the source and target voxels but showed detailed patterns of dependence on hierarchical relationships that we did not observe in DCNNs. Given these results, we argue that the structured variance unexplained by DCNN-based encoding models is unlikely to be entirely caused by non-neural artifacts (e.g., spatially correlated measurement noise) or a failure of DCNNs to approximate the features encoded in brain activity; rather, our results point to a need for brain models that provide both mechanistic and computational explanations for structured ongoing activity in the brain. Keywords: fMRI, encoding models, deep neural networks, functional connectivity.

The deep grey matter (DGM) nuclei of the brain play a crucial role in learning, behaviour, cognition, movement and memory. Although automated segmentation strategies can provide insight into the impact of multiple neurological conditions affecting these structures, such as Multiple Sclerosis (MS), Huntington's disease (HD), Alzheimer's disease (AD), Parkinson's disease (PD) and Cerebral Palsy (CP), there are a number of technical challenges limiting an accurate automated segmentation of the DGM. Namely, the insufficient contrast of T1 sequences to completely identify the boundaries of these structures, as well as the presence of iso-intense white matter lesions or extensive tissue loss caused by brain injury. Therefore in this systematic review, 269 eligible studies were analysed and compared to determine the optimal approaches for addressing these technical challenges. The automated approaches used among the reviewed studies fall into three broad categories, atlas-based approaches focusing on the accurate alignment of atlas priors, algorithmic approaches which utilise intensity information to a greater extent, and learning-based approaches that require an annotated training set. Studies that utilise freely available software packages such as FIRST, FreeSurfer and LesionTOADS were also eligible, and their performance compared. Overall, deep learning approaches achieved the best overall performance, however these strategies are currently hampered by the lack of large-scale annotated data. Improving model generalisability to new datasets could be achieved in future studies with data augmentation and transfer learning. Multi-atlas approaches provided the second-best performance overall, and may be utilised to construct a "silver standard" annotated training set for deep learning. To address the technical challenges, providing robustness to injury can be improved by using multiple channels, highly elastic diffeomorphic transformations such as LDDMM, and by following atlas-based approaches with an intensity driven refinement of the segmentation, which has been done with the Expectation Maximisation (EM) and level sets methods. Accounting for potential lesions should be achieved with a separate lesion segmentation approach, as in LesionTOADS. Finally, to address the issue of limited contrast, R2*, T2* and QSM sequences could be used to better highlight the DGM due to its higher iron content. Future studies could look to additionally acquire these sequences by retaining the phase information from standard structural scans, or alternatively acquiring these sequences for only a training set, allowing models to learn the "improved" segmentation from T1-sequences alone.

Neuroimaging techniques such as MRI have been widely used to explore the associations between brain areas. Structural connectivity (SC) captures the anatomical pathways across the brain and functional connectivity (FC) measures the correlation between the activity of brain regions. These connectivity measures have been much studied using network theory in order to uncover the distributed organization of brain structures, in particular FC for task-specific brain communication. However, the application of network theory to study FC matrices is often “static” despite the dynamic nature of time series obtained from fMRI. The present study aims to overcome this limitation by introducing a network-oriented analysis applied to whole-brain effective connectivity (EC) useful to interpret the brain dynamics. Technically, we tune a multivariate Ornstein-Uhlenbeck (MOU) process to reproduce the statistics of the whole-brain resting-state fMRI signals, which provides estimates for MOU-EC as well as input properties (similar to local excitabilities). The network analysis is then based on the Green function (or network impulse response) that describes the interactions between nodes across time for the estimated dynamics. This model-based approach provides time-dependent graph-like descriptor, named communicability, that characterize the roles that either nodes or connections play in the propagation of activity within the network. They can be used at both global and local levels, and also enables the comparison of estimates from real data with surrogates (e.g. random network or ring lattice). In contrast to classical graph approaches to study SC or FC, our framework stresses the importance of taking the temporal aspect of fMRI signals into account. Our results show a merging of functional communities over time (in which input properties play a role), moving from segregated to global integration of the network activity. Our formalism sets a solid ground for the analysis and interpretation of fMRI data, including task-evoked activity.

Connectomics is essential for understanding large-scale brain networks but requires that individual connection estimates are neurobiologically interpretable. In particular, a principle of brain organization is that reciprocal connections between cortical areas are functionally asymmetric. This is a challenge for fMRI-based connectomics in humans where only undirected functional connectivity estimates are routinely available. By contrast, whole-brain estimates of effective (directed) connectivity are computationally challenging, and emerging methods require empirical validation. Here, using a motor task at 7T, we demonstrate that a novel generative model can infer known connectivity features in a whole-brain network (>200 regions, >40,000 connections) highly efficiently. Furthermore, graph-theoretical analyses of directed connectivity estimates identify functional roles of motor areas more accurately than undirected functional connectivity estimates. These results, which can be achieved in an entirely unsupervised manner, demonstrate the feasibility of inferring directed connections in whole-brain networks and open new avenues for human connectomics., NeuroImage, 225, ISSN:1053-8119, ISSN:1095-9572

The present study examined the neural and behavioral correlates of selective stopping, a form of inhibition that has scarcely been investigated. The selectivity of the inhibitory process is needed when individuals have to deal with an environment filled with multiple stimuli, some of which require inhibition and some of which do not. The stimulus-selective stop-signal task has been used to explore this issue assuming that all participants interrupt their ongoing responses selectively to stop but not to ignore signals. However, recent behavioral evidence suggests that some individuals do not carry out the task as experimenters expect, since they seemed to interrupt their response non-selectively to both signals. In the present study, we detected and controlled the cognitive strategy adopted by participants (n=57) when they performed a stimulus-selective stop-signal task before comparing brain activation between conditions. In order to determine both the onset and the end of the response cancellation process underlying each strategy and to fully take advantage of the precise temporal resolution of event-related potentials, we used a mass univariate approach. Source localization techniques were also employed to estimate the neural underpinnings of the effects observed at the scalp level. Our results from scalp and source level analysis support the behavioral-based strategy classification. Specific effects were observed depending on the strategy adopted by participants. Thus, when contrasting successful stop versus ignore conditions, increased activation was only evident for subjects who were classified as using a strategy whereby the response interruption process was selective to stop trials. This increased activity was observed during the P3 time window in several left-lateralized brain regions, including middle and inferior frontal gyri, as well as parietal and insular cortices. By contrast, in those participants who used a strategy characterized by stopping non-selectively, no activation differences between successful stop and ignore conditions were observed at the estimated time at which response interruption process occurs. Overall, results from the current study highlight the importance of controlling for the different strategies adopted by participants to perform selective stopping tasks before analyzing brain activation patterns.

There has been a long debate about the neural mechanism of event-related potentials (ERPs). Previously, no evidence or method was apparent to validate the two competing models, the evoked model and the oscillation model. One argument is whether the pre-stimulus brain oscillation could influence the following ERP. This study carried out an innovative visual oddball task experiment to investigate the dynamic process of visual evoked potentials. A period of stable oscillations of specified dominant frequencies and initial phases, i.e. the steady-state baseline, would be induced before responses to transient stimuli of different contrasts, which could overcome the artifact problem caused by the 'sorting' method. The result first revealed a 'three-period-transition' for the generation of visual evoked potentials by an objective decomposition. The ERP almost retained the preceding oscillation during the first period, provided an unstable negative potential in the second period, and generated the N1 component in the third period. The cross term analysis showed that the evoked model couldn't be the whole explanation for the ERP generation. Furthermore, the component analysis revealed that the N1 latency was sensitive to the initial phase under the low stimulus contrast (supporting the oscillation model) but not under the high stimulus contrast (supporting the evoked model). It demonstrated that the external stimulus contrast is a significant factor deciding the explicit model for ERPs. Our method and preliminary results may help reconcile the previous, seemly contradictory findings on the ERP mechanism.

We investigate the relationship of resting-state fMRI functional connectivity estimated over long periods of time with time-varying functional connectivity estimated over shorter time intervals. We show that using Pearson's correlation to estimate functional connectivity implies that the range of fluctuations of functional connections over short time scales is subject to statistical constraints imposed by their connectivity strength over longer scales. We present a method for estimating time-varying functional connectivity that is designed to mitigate this issue and allows us to identify episodes where functional connections are unexpectedly strong or weak. We apply this method to data recorded from $N=80$ participants, and show that the number of unexpectedly strong/weak connections fluctuates over time, and that these variations coincide with intermittent periods of high and low modularity in time-varying functional connectivity. We also find that during periods of relative quiescence regions associated with default mode network tend to join communities with attentional, control, and primary sensory systems. In contrast, during periods where many connections are unexpectedly strong/weak, default mode regions dissociate and form distinct modules. Finally, we go on to show that, while all functional connections can at times manifest stronger (more positively correlated) or weaker (more negatively correlated) than expected, a small number of connections, mostly within the visual and somatomotor networks, do so a disproportional number of times. Our statistical approach allows the detection of functional connections that fluctuate more or less than expected based on their long-time averages and may be of use in future studies characterizing the spatio-temporal patterns of time-varying functional connectivity, Comment: 47 Pages, 8 Figures, 4 Supplementary Figures

Speech recognition is rapid, automatic and amazingly robust. How the brain is able to decode speech from noisy acoustic inputs is unknown. We show that the brain recognizes speech by integrating bottom-up acoustic signals with top-down predictions. Subjects listened to intelligible normal and unintelligible fine structure speech that lacked the predictability of the temporal envelope and did not enable access to higher linguistic representations. Their top-down predictions were manipulated using priming. Activation for unintelligible fine structure speech was confined to primary auditory cortices, but propagated into posterior middle temporal areas when fine structure speech was made intelligible by top-down predictions. By contrast, normal speech engaged posterior middle temporal areas irrespective of subjects’ predictions. Critically, when speech violated subjects’ expectations, activation increases in anterior temporal gyri/sulci signalled a prediction error and the need for new semantic integration. In line with predictive coding, our findings compellingly demonstrate that top-down predictions determine whether and how the brain translates bottom-up acoustic inputs into intelligible speech.

Functional brain connectivity is increasingly being seen as critical for cognition, perception and motor control.Magnetoencephalography and electroencephalography are modalities that offer noninvasive mapping of electrophysiological interactions among brain regions, yet suffer from signal leakage and signal cancellation when estimating brain activity. This leads to biased connectivity values which complicate interpretation. In this study, we test the hypothesis that a Multiple Constrained Minimum Variance beamformer (MCMV) outperforms the more traditional Linearly Constrained Minimum Variance beamformer (LCMV) for estimation of electrophysiological connectivity. To this end, MCMV and LCMV performance is compared in task related analyses with both simulated data and human MEG recordings of visual steady state signals, and in resting state analyses with simulated data and human MEG data of 89 subjects. In task related scenarios connectivity was estimated using coherence and phase locking values, whereas envelope correlations were used for the resting state data. We also introduce a novel Augmented Pairwise MCMV (APW-MCMV) approach for signal leakage suppression in resting state analyses and assess its performance against LCMV and more conventional MCMV approaches. We demonstrate that with MCMV effects of signal mixing and coherent source cancellation are greatly reduced in both task related and resting state conditions, while in contrast to other approaches 0-and short time lag interactions are preserved. In addition, we demonstrate that in resting state analyses, APW-MCMV strongly reduces spurious connections while better controlling for false negatives compared to more conservative measures such as symmetrical orthogonalization.

Recently, resting functional MRI has provided invaluable insight into the brain developmental processes of early infancy and childhood. A common feature of previous functional development studies is the use of age to separate subjects into different cohorts for group comparisons. However, functional maturation paces vary tremendously from subject to subject. Since this is particularly true for the first years of life, an alternative to physical age alone is needed for cluster analysis. Here, a data-driven approach based on individual brain functional connectivity was employed to cluster typically developing children who were longitudinally imaged using MRI without sedation for the first two years of life. Specifically, three time periods were determined based on the distinction of brain functional connectivity patterns, including 0-1 month (group 1), 2-7 months (group 2), and 8-24 (group 3) of age, respectively. From groups 1 to 2, connection density increased by almost two-fold, local efficacy (LE) is significantly improved, and there was no change in global efficiency (GE). From groups 2 to 3, connection density increased slightly, LE showed no change, and a significant increase in GE were observed. Furthermore, 27 core brain regions were identified which yielded clustering results that resemble those obtained using all brain regions. These core regions were largely associated with the motor, visual and language functional domains as well as regions associated with higher order cognitive functional domains. Both visual and language functional domains exhibited a persistent and significant increase within domain connection from groups 1 to 3, while no changes were observed for the motor domain. In contrast, while a reduction of inter-domain connection was the general developmental pattern, the motor domain exhibited an interesting "V" shape pattern in its relationship to visual and language associated areas, showing a decrease from groups 1 to 2, followed by an increase from groups 2 to 3. In summary, our results offer new insights into functional brain development and identify 27 core brain regions critically important for early brain development.

How the brain encodes abstract numerical symbols is a fundamental question in philosophy and cognitive neuroscience alike. Here we probe the nature of symbolic number representation in the brain by characterizing the neural similarity space for symbolic quantities in regions sensitive to their semantic content. In parietal and occipital regions, the similarity space of number symbols was positively predicted by the lexical frequency of numerals in parietal and occipital areas, and was unrelated to numerical ratio. These results are more consistent with a categorical, frequency-based account of symbolic quantity encoding. In contrast, the similarity space of analog quantities was positively predicted by ratio in prefrontal, parietal and occipital regions. We thus provide an explanation for why previous work has indicated that symbolic and analog quantities are distinct: number symbols operate primarily like discrete categories sensitive to input frequency, while analog quantities operate more like approximate perceptual magnitudes. In addition, we find substantial evidence for related patterns of activity across formats in prefrontal, parietal and occipital regions. Crucially however, between-format relations were not specific to individual quantities, indicating common processing as opposed to common representation. Moreover, evidence for between-format processing was strongest for quantities that could be represented as exact, discrete values in both systems (quantities in the 'subitizing' range: 1-4). In sum, converging evidence presented here indicates that symbolic quantities are coded in the brain as discrete categories sensitive to input frequency and largely independent of approximate, analog quantities.

Purpose The high chemical shift separation at 9.4 T allows for selective saturation of proton pools in exchange with water protons. For the first time, highly selective and comprehensive chemical exchange saturation transfer (CEST) experiments were performed in the human brain at 9.4 T. This work provides insight into CEST signals in the human brain in comparison with existing animal studies, as well as with CEST effects in vivo at lower field strengths. Methods A novel snapshot-CEST method for human brain scans at 9.4 T was optimized and employed for highly-spectrally-resolved (95 offsets) CEST measurements in healthy subjects and one brain tumor patient. Reproducibility and stability between scans was verified in grey and white matter after B0, B1, and motion correction of the acquired 3D CEST volumes. Two-step Lorentzian fitting was used to further improve separation of spectrally discernible signals to create known and novel CEST contrast maps at 9.4 T. Results At a saturation power of B1 = 0.5 μT most selective CEST effects could be obtained in the human brain with high inter-scan reproducibility. While contrast behavior of previously measured signals at lower field, namely amide-, guanidyl- and NOE-CEST effects, could be reproduced, novel signals at 2.7 ppm, and −1.6 ppm could be verified in healthy subjects and in a brain tumor patient for the first time. Conclusion High spectral resolution chemical exchange saturation transfer at 9.4 T allows deeper insights into the Z-spectrum structure of the human brain, and provides many different contrasts showing different correlations in healthy tissue and in tumor-affected areas of the brain, generating hypotheses for future investigations of in-vivo-CEST at UHF.

Functional connectivity provides an informative and powerful framework for exploring brain organization. Despite this, few statistical methods are available for the accurate estimation of dynamic changes in functional network architecture. To date, the majority of existing statistical techniques have assumed that connectivity structure is stationary, which is in direct contrast to emerging data that suggests that the strength of connectivity between regions is variable over time. Therefore, the development of statistical methods that enable exploration of dynamic changes in functional connectivity is currently of great importance to the neuroscience community. In this paper, we introduce the 'Multiplication of Temporal Derivatives' (MTD) and then demonstrate the utility of this metric to: (i) detect dynamic changes in connectivity using data from a novel state-switching simulation; (ii) accurately estimate graph structure in a previously-described 'ground-truth' simulated dataset; and (iii) identify task-driven alterations in functional connectivity. We show that the MTD is more sensitive than existing sliding-window methods in detecting dynamic alterations in connectivity structure across a range of correlation strengths and window lengths in simulated data. In addition to the temporal precision offered by MTD, we demonstrate that the metric is also able to accurately estimate stationary network structure in both simulated and real task-based data, suggesting that the method may be used to identify dynamic changes in network structure as they evolve through time.

Independent component analysis (ICA) and seed-based analyses are widely used techniques for studying intrinsic neuronal activity in task-based or resting scans. In this work, we show there is a direct link between the two, and show that there are some important differences between the two approaches in terms of what information they capture. We developed an enhanced connectivity-matrix independent component analysis (cmICA) for calculating whole brain voxel maps of functional connectivity, which reduces the computational complexity of voxel-based connectivity analysis on performing many temporal correlations. We also show there is a mathematical equivalency between parcellations on voxel-to-voxel functional connectivity and simplified cmICA. Next, we used this cost-efficient data-driven method to examine the resting state fMRI connectivity in schizophrenia patients (SZ) and healthy controls (HC) on a whole brain scale and further quantified the relationship between brain functional connectivity and cognitive performances measured by the Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) battery. Current results suggest that SZ exhibit a wide-range abnormality, primarily a decrease, in functional connectivity both between networks and within different network hubs. Specific functional connectivity decreases were associated with MATRICS performance deficits. In addition, we found that resting state functional connectivity decreases was extensively associated with aging regardless of groups. In contrast, there was no relationship between positive and negative symptoms in the patients and functional connectivity. In sum, we have developed a novel mathematical relationship between ICA and seed-based connectivity that reduces computational complexity, which has broad applicability, and showed a specific application of this approach to characterize connectivity changes associated with cognitive scores in SZ.

Recent studies have suggested that individuals can form multiple motor memories when simultaneously adapting to multiple, but oppositely-oriented perturbations. These findings predict that individuals detect the change in learning context, allowing the selective initialization and update of motor memories. However, previous electrophysiological studies of sensorimotor adaptation have not identified a neural mechanism supporting the detection of a context switch and adaptation to separate contexts. Here, we tested the hypothesis that such a mechanism is identifiable through neural oscillations measured through EEG. Human participants learned to manipulate an object in two opposite contexts (mass distribution). This task was designed based on previous work showing that people can adapt to both contexts. We found that sensorimotor α and β, and medial frontal θ frequency bands all exhibited different response patterns with respect to the error in each context. To determine whether any frequency's responses to error were distinctly related to a switch in context, we predicted single-trial EEG data from a computational learning model that can adapt to multiple contexts simultaneously based on a switching mechanism. This analysis revealed that only medial frontal θ was predicted by a component of the model state that adapts to errors based on a context switch. In contrast, α and β were predicted by a model state that was updated from performance errors independent of the context. These findings provide novel evidence showing that sensorimotor and medial frontal oscillations are predicted by different adaptation processes, and that changes in medial frontal activity may indicate the formation of motor memories by responding to changes in learning context.

Resting-state functional connectivity is a powerful tool for studying human functional brain networks. Temporal fluctuations in functional connectivity, i.e.,dynamicfunctional connectivity (dFC), are thought to reflect dynamic changes in brain organization andnon-stationaryswitching of discrete brain states. However, recent studies have suggested that dFC might be attributed to sampling variability of static FC. Despite this controversy, a detailed exposition of stationarity and statistical testing of dFC is lacking in the literature. This article seeks an in-depth exploration of these statistical issues at a level appealing to both neuroscientists and statisticians.We first review the statistical notion of stationarity, emphasizing its reliance on ensemble statistics. In contrast, all FC measures depend on sample statistics. An important consequence is that the space of stationary signals is much broader than expected, e.g., encompassing hidden markov models (HMM) widely used to extract discrete brain states. In other words, stationarity does not imply the absence of brain states. We then expound the assumptions underlying the statistical testing of dFC. It turns out that the two popular frameworks - phase randomization (PR) and autoregressive randomization (ARR) - generate stationary, linear, Gaussian null data. Therefore, statistical rejection can be due to non-stationarity, nonlinearity and/or non-Gaussianity. For example, the null hypothesis can be rejected for the stationary HMM due to nonlinearity and non-Gaussianity. Finally, we show that a common form of ARR (bivariate ARR) is susceptible to false positives compared with PR and an adapted version of ARR (multivariate ARR).Application of PR and multivariate ARR to Human Connectome Project data suggests that the stationary, linear, Gaussian null hypothesis cannot be rejected for most participants. However, failure to reject the null hypothesis does not imply that static FC can fully explain dFC. We find that first order AR models explain temporal FC fluctuations significantly better than static FC models. Since first order AR models encode both static FC and one-lag FC, this suggests the presence of dynamical information beyond static FC. Furthermore, even in subjects where the null hypothesis was rejected, AR models explain temporal FC fluctuations significantly better than a popular HMM, suggesting the lack of discrete states (as measured by resting-state fMRI). Overall, our results suggest that AR models are not only useful as a means for generating null data, but may be a powerful tool for exploring the dynamical properties of resting-state fMRI. Finally, we discuss how apparent contradictions in the growing dFC literature might be reconciled.

Scientific and mathematical thinking relies on the ability to evaluate whether conclusions drawn from conditional (if-then) arguments are logically valid. Yet, the neural development of this ability -- termed deductive reasoning -- is largely unknown. Here we aimed to identify the neural mechanisms that underlie the emergence of deductive reasoning with conditional rules in children. We further tested whether these mechanisms have their roots in the neural mechanisms involved in judging the likelihood of conclusions. In a functional Magnetic Resonance Imaging (fMRI) scanner, 8- to 13-year-olds were presented with causal conditional problems such as “If a baby is hungry then he will start crying; The baby is crying; Is the baby hungry?”. In Validity trials, children were asked to indicate whether the conclusion followed out of necessity from the premises. In Likelihood trials, they indicated the degree of likelihood of the conclusion. We found that children who made accurate judgments of logical validity (as compared to those who did not) exhibited enhanced activity in left and medial frontal regions. In contrast, differences in likelihood ratings between children were related to differences of activity in right frontal and bilateral parietal regions. There was no overlap between the brain regions underlying validity and likelihood judgments. Therefore, our results suggest that the ability to evaluate the logical validity of conditional arguments emerges from brain mechanisms that qualitatively differ from those involved in evaluating the likelihood of these arguments in children.

Although analogical reasoning (AR) plays a central role in higher-level cognition and constitutes a key source of individual differences in intellectual ability, the neural mechanisms that account for individual differences in AR remain to be well characterized. Here we investigated individual differences in AR within a large sample (n = 229), using multivariate fMRI analysis (a simple multiple kernel learning machine). The individual AR capability was positively correlated with activation level in a prefrontal executive network and a visuospatial network. Notably, the best predictors of individual differences in AR within these networks were activation in the dorsomedial prefrontal cortex (response selection) and the lingual gyrus (visual feature mapping). In contrast, AR capability was negatively correlated with activation in the default mode network. The implications of the reported findings are twofold: (i) Individual differences in AR depend on multiple executive and visuospatial brain regions, where their respective contributions are contingent upon the individuals' cognitive skills; (ii) Brain regions associated with individual differences in AR only partially overlap with brain regions sensitive to the associated task demands (i.e., brain regions sensitive to the analogy relational complexity, at the group-level). We discuss implications of such brain organization supporting AR as an example for brain architecture underlying higher-level cognitive processes.

Structural and functional connectomes are emerging as important instruments in the study of normal brain function and in the development of new biomarkers for a variety of brain disorders. In contrast to single-network studies that presently dominate the (non-connectome) network literature, connectome analyses typically examine groups of empirical networks and then compare these against standard (stochastic) network models. The current practice in connectome studies is to employ stochastic network models derived from social science and engineering contexts as the basis for the comparison. However, these are not necessarily best suited for the analysis of connectomes, which often contain groups of very closely related networks, such as occurs with a set of controls or a set of patients with a specific disorder. This paper studies important extensions of standard stochastic models that make them better adapted for analysis of connectomes, and develops new statistical fitting methodologies that account for inter-subject variations. The extensions explicitly incorporate geometric information about a network based on distances and inter/intra hemispherical asymmetries (to supplement ordinary degree-distribution information), and utilize a stochastic choice of network density levels (for fixed threshold networks) to better capture the variance in average connectivity among subjects. The new statistical tools introduced here allow one to compare groups of networks by matching both their average characteristics and the variations among them. A notable finding is that connectomes have high "smallworldness" beyond that arising from geometric and degree considerations alone.

Current findings suggest that confidence emerges only after decision making. However, the temporal and neural dynamics of the emergence of post-decision confidence--a metacognitive judgement--are not fully explored. To gain insight into the dynamics of post-decision confidence processing and to disentangle the processes underlying confidence judgements and decision making, we applied a tactile discrimination task during functional magnetic resonance imaging (fMRI). Our results revealed that reaction times to post-decision confidence depend on the level of confidence, suggesting that post-decision confidence in a perceptual choice is not processed in parallel to perceptual decision making. Moreover, we demonstrated by the parametric analysis of fMRI data that post-decisionally modelled confidence processing can be distinguished from processes related to decision making through anatomical location and through the pattern of neural activity. In contrast to perceptual decision making, post-decision confidence appears to be strictly allocated to a prefrontal network of brain regions, primarily the anterior and dorsolateral prefrontal cortex, areas that have been related to metacognition. Moreover, the processes underlying decision making and post-decision confidence may share recruitment of the dorsolateral prefrontal cortex, although the former probably has distinct functions with regard to processing of perceptual choices and post-decision confidence. Thus, this is the first fMRI study to disentangle the processes underlying post-decision confidence and decision making on behavioural, neuroanatomical, and neurofunctional levels. With regard to the temporal evolution of post-decision confidence, results of the present study provide strong support for the most recent theoretical models of human perceptual decision making, and thus provide implications for investigating confidence in perceptual paradigms.

International audience; Neuroimaging group analyses are used to relate inter-subject signal differences observed in brain imaging with behavioral or genetic variables and to assess risks factors of brain diseases. The lack of stability and of sensitivity of current voxel-based analysis schemes may however lead to non-reproducible results. We introduce a new approach to overcome the limitations of standard methods, in which active voxels are detected according to a consensus on several random parcellations of the brain images, while a permutation test controls the false positive risk. Both on synthetic and real data, this approach shows higher sensitivity, better accuracy and higher reproducibility than state-of-the-art methods. In a neuroimaging-genetic application, we find that it succeeds in detecting a significant association between a genetic variant next to the COMT gene and the BOLD signal in the left thalamus for a functional Magnetic Resonance Imaging contrast associated with incorrect responses of the subjects from a Stop Signal Task protocol.

Experimental MEG source imaging studies have typically been carried out with either a spherically symmetric head model or a single-shell boundary-element (BEM) model that is shaped according to the inner skull surface. The concepts and comparisons behind these simplified models have led to misunderstandings regarding the role of skull and scalp in MEG. In this work, we assess the forward-model errors due to different skull/scalp approximations and due to differences and errors in model geometries. We built five anatomical models of a volunteer using a set of T1-weighted MR scans and three common toolboxes. Three of the models represented typical models in experimental MEG, one was manually constructed, and one contained a major segmentation error at the skull base. For these anatomical models, we built forward models using four simplified approaches and a three-shell BEM approach that has been used as reference in previous studies. Our reference model contained in addition the skull fine-structure (spongy bone). We computed signal topographies for cortically constrained sources in the left hemisphere and compared the topographies using relative error and correlation metrics. The results show that the spongy bone has a minimal effect on MEG topographies, and thus the skull approximation of the three-shell model is justified. The three-shell model performed best, followed by the corrected-sphere and single-shell models, whereas the local-spheres and single-sphere models were clearly worse. The three-shell model was the most robust against the introduced segmentation error. In contrast to earlier claims, there was no noteworthy difference in the computation times between the realistically-shaped and sphere-based models, and the manual effort of building a three-shell model and a simplified model is comparable. We thus recommend the realistically-shaped three-shell model for experimental MEG work. In cases where this is not possible, we recommend a realistically-shaped corrected-sphere or single-shell model.

Intraoperative tissue deformation, known as brain shift, decreases the benefit of using preoperative images to guide neurosurgery. Non-rigid registration of preoperative magnetic resonance (MR) to intraoperative ultrasound (iUS) has been proposed as a means to compensate for brain shift. We focus on the initial registration from MR to predurotomy iUS. We present a method that builds on previous work to address the need for accuracy and generality of MR-iUS registration algorithms in multi-site clinical data. High-dimensional texture attributes were used instead of image intensities for image registration and the standard difference-based attribute matching was replaced with correlation-based attribute matching. A strategy that deals explicitly with the large field-of-view mismatch between MR and iUS images was proposed. Key parameters were optimized across independent MR-iUS brain tumor datasets acquired at 3 institutions, with a total of 43 tumor patients and 758 reference landmarks for evaluating the accuracy of the proposed algorithm. Despite differences in imaging protocols, patient demographics and landmark distributions, the algorithm is able to reduce landmark errors prior to registration in three data sets (5.37±4.27, 4.18±1.97 and 6.18±3.38 mm, respectively) to a consistently low level (2.28±0.71, 2.08±0.37 and 2.24±0.78 mm, respectively). This algorithm was tested against 15 other algorithms and it is competitive with the state-of-the-art on multiple datasets. We show that the algorithm has one of the lowest errors in all datasets (accuracy), and this is achieved while sticking to a fixed set of parameters for multi-site data (generality). In contrast, other algorithms/tools of similar performance need per-dataset parameter tuning (high accuracy but lower generality), and those that stick to fixed parameters have larger errors or inconsistent performance (generality but not the top accuracy). Landmark errors were further characterized according to brain regions and tumor types, a topic so far missing in the literature., National Institute of Health (Grants P41-EB015898, P41-EB015902 and R01-NS049251)

Several independent studies have demonstrated that small amounts of in-scanner motion systematically bias estimates of resting-state functional connectivity. This confound is of particular importance for studies of neurodevelopment in youth because motion is strongly related to subject age during this period. Critically, the effects of motion on connectivity mimic major findings in neurodevelopmental research, specifically an age-related strengthening of distant connections and weakening of short-range connections. Here, in a sample of 780 subjects ages 8-22, we re-evaluate patterns of change in functional connectivity during adolescent development after rigorously controlling for the confounding influences of motion at both the subject and group levels. We find that motion artifact inflates both overall estimates of age-related change as well as specific distance-related changes in connectivity. When motion is more fully accounted for, the prevalence of age-related change as well as the strength of distance-related effects is substantially reduced. However, age-related changes remain highly significant. In contrast, motion artifact tends to obscure age-related changes in connectivity associated with segregation of functional brain modules; improved preprocessing techniques allow greater sensitivity to detect increased within-module connectivity occurring with development. Finally, we show that subject's age can still be accurately estimated from the multivariate pattern of functional connectivity even while controlling for motion. Taken together, these results indicate that while motion artifact has a marked and heterogeneous impact on estimates of connectivity change during adolescence, functional connectivity remains a valuable phenotype for the study of neurodevelopment.

Item does not contain fulltext General linear model (GLM) is a standard and widely used fMRI analysis tool. It enables the detection of hypothesis-driven brain activations. In contrast, Independent Component Analysis (ICA) is a powerful technique, which enables the detection of data-driven spatially independent networks. Hybrid approaches that combine and take advantage of GLM and ICA have been proposed. Yet the choice of the best method is still a challenge, considering that the techniques may yield slightly different results regarding the number of brain regions involved in a task. A poor statistical power or the deviance from the predicted hemodynamic response functions is possible cause for GLM failures in extracting some activations picked by ICA. However, there might be another explanation for different results obtained with GLM and ICA approaches, such as networks cancelation. In this paper, we propose a new supplementary method that can give more insight into the functional data as well as help to clarify inconsistencies between the results of studies using GLM and ICA. We introduce a contributive sources analysis (CSA), which provides a measure of the number and the strength of the neural networks that significantly contribute to brain activation. CSA, applied to fMRI data of anti-saccades, enabled us to verify whether the brain regions involved in the task are dominated by a single network or serve as key nodes for particular networks interaction. Moreover, when applying CSA to the atlas-defined regions-of-interest, results indicated that activity of the parieto-medial temporal network was suppressed by the eye field network and the default mode network. Thus, this effect of networks cancelation explains the absence of parieto-medial temporal activation within the GLM results. Together, those findings indicate that brain activations are a result of complex network interactions. Applying CSA appears to be a useful tool to reveal additional findings outside the scope of the "fixed-model" GLM and data-driven ICA approaches.