Whole-brain resting-state connectivity underlying impaired inhibitory control during early versus longer-term abstinence in cocaine addiction.

Lapses in inhibitory control have been linked to relapse in human drug addiction. Evidence suggests differences in inhibitory control depending on abstinence duration, but the underlying neural mechanisms remain unknown. We hypothesized that early abstinence (2-5 days) would be characterized by the strongest impairments of inhibitory control and most wide-spread deviations in resting-state functional connectivity of brain networks, while longer-term abstinence (>30 days) would be characterized by weaker impairments as compared to healthy controls. In this laboratory-based cross-sectional study, we compared individuals with Cocaine Use Disorder (iCUD) during early (cocaine urine-positive: N = 19, iCUD+; 32% female; mean age: 46.8 years) and longer-term abstinence (cocaine urine-negative: N = 29, iCUD-; 15% female; mean age: 46.6 years) to healthy controls (N = 33; 24% female; mean age: 40.9 years). We compared the groups on inhibitory control performance (Stop-Signal Task) and, using a whole-brain graph theory analysis (638 region parcellation) of functional magnetic resonance imaging (fMRI) data, we tested for group differences in resting-state brain function (local/global efficiency). We characterized how resting-state brain function was associated with inhibitory control performance within iCUD. Inhibitory control performance was worst in the early abstinence group, and intermediate in the longer-term abstinence group, as compared to the healthy control group (P < 0.01). More recent use of cocaine (CUD+ > CUD- > healthy controls) was characterized by decreased efficiency in fronto-temporal and subcortical networks (primarily in the salience, semantic, and basal ganglia networks) and increased efficiency in visual networks. Importantly, a similar functional connectivity pattern characterized impaired inhibitory control performance within iCUD (all brain analyses P < 0.05, FWE-corrected). Together, we demonstrated that a similar pattern of systematic and widespread deviations in resting-state brain efficiency, extending beyond the networks commonly investigated in human drug addiction, is linked to both abstinence duration and inhibitory control deficits in iCUD.

Age differences in the functional architecture of the human brain.

The intrinsic functional organization of the brain changes into older adulthood. Age differences are observed at multiple spatial scales, from global reductions in modularity and segregation of distributed brain systems, to network-specific patterns of dedifferentiation. Whether dedifferentiation reflects an inevitable, global shift in brain function with age, circumscribed, experience-dependent changes, or both, is uncertain. We employed a multimethod strategy to interrogate dedifferentiation at multiple spatial scales. Multi-echo (ME) resting-state fMRI was collected in younger (n = 181) and older (n = 120) healthy adults. Cortical parcellation sensitive to individual variation was implemented for precision functional mapping of each participant while preserving group-level parcel and network labels. ME-fMRI processing and gradient mapping identified global and macroscale network differences. Multivariate functional connectivity methods tested for microscale, edge-level differences. Older adults had lower BOLD signal dimensionality, consistent with global network dedifferentiation. Gradients were largely age-invariant. Edge-level analyses revealed discrete, network-specific dedifferentiation patterns in older adults. Visual and somatosensory regions were more integrated within the functional connectome; default and frontoparietal control network regions showed greater connectivity; and the dorsal attention network was more integrated with heteromodal regions. These findings highlight the importance of multiscale, multimethod approaches to characterize the architecture of functional brain aging.

Conservative and disruptive modes of adolescent change in human brain functional connectivity.

Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: "conservative" and "disruptive." Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman's correlation between edgewise baseline FC (at 14 y, [Formula: see text]) and adolescent change in FC ([Formula: see text]), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas.

Lateralisation of subcortical functional connectivity during and after general anaesthesia.

BACKGROUND: Arousal and awareness are two important components of consciousness states. Functional neuroimaging has furthered our understanding of cortical and thalamocortical mechanisms of awareness. Investigating the relationship between subcortical functional connectivity and arousal has been challenging owing to the relatively small size of brainstem structures and thalamic nuclei, and their depth in the brain. METHODS: Resting state functional MRI scans of 72 healthy volunteers were acquired before, during, 1 h after, and 1 day after sevoflurane general anaesthesia. Functional connectivity of subcortical regions of interest vs whole brain and homotopic functional connectivity for assessment of left-right symmetry analyses of both cortical and subcortical regions of interest were performed. Both analyses used high resolution atlases generated from deep brain stimulation applications. RESULTS: Functional connectivity in subcortical loci within the thalamus and of the ascending reticular activating system was sharply restricted under anaesthesia, featuring a general lateralisation of connectivity. Similarly, left-right homology was sharply reduced under anaesthesia. Subcortical bilateral functional connectivity was not fully restored after emergence from anaesthesia, although greater restoration was seen between ascending reticular activating system loci and specific thalamic nuclei thought to be involved in promoting and maintaining arousal. Functional connectivity was fully restored to baseline by the following day. CONCLUSIONS: Functional connectivity in the subcortex is sharply restricted and lateralised under general anaesthesia. This restriction may play a part in loss and return of consciousness. CLINICAL TRIAL REGISTRATION: NCT02275026.

Advances in resting state fMRI acquisitions for functional connectomics.

Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.

Sex-specific impact of prenatal androgens on social brain default mode subsystems.

Early-onset neurodevelopmental conditions (e.g., autism) affect males more frequently than females. Androgens may play a role in this male-bias by sex-differentially impacting early prenatal brain development, particularly neural circuits that later develop specialized roles in social cognition. Here, we find that increasing prenatal testosterone in humans is associated with later reduction of functional connectivity between social brain default mode (DMN) subsystems in adolescent males, but has no effect in females. Since testosterone can work directly via the androgen receptor (AR) or indirectly via the estrogen receptor through aromatase conversion to estradiol, we further examined how a potent non-aromatizable androgen, dihydrotestosterone (DHT), acts via the AR to influence gene expression in human neural stem cells (hNSC)-particularly for genes of high-relevance for DMN circuitry. DHT dysregulates a number of genes enriched for syndromic causes of autism and intellectual disability and for genes that in later development are expressed in anatomical patterns that highly correspond to the cortical midline DMN subsystem. DMN-related and DHT-affected genes (e.g., MEF2C) are involved in a number of synaptic processes, many of which impact excitation-inhibition balance. Androgens have male-specific prenatal influence over social brain circuitry in humans and may be relevant towards explaining some component of male-bias in early-onset neurodevelopmental conditions.

Ridding fMRI data of motion-related influences: Removal of signals with distinct spatial and physical bases in multiecho data.

"Functional connectivity" techniques are commonplace tools for studying brain organization. A critical element of these analyses is to distinguish variance due to neurobiological signals from variance due to nonneurobiological signals. Multiecho fMRI techniques are a promising means for making such distinctions based on signal decay properties. Here, we report that multiecho fMRI techniques enable excellent removal of certain kinds of artifactual variance, namely, spatially focal artifacts due to motion. By removing these artifacts, multiecho techniques reveal frequent, large-amplitude blood oxygen level-dependent (BOLD) signal changes present across all gray matter that are also linked to motion. These whole-brain BOLD signals could reflect widespread neural processes or other processes, such as alterations in blood partial pressure of carbon dioxide (pCO2) due to ventilation changes. By acquiring multiecho data while monitoring breathing, we demonstrate that whole-brain BOLD signals in the resting state are often caused by changes in breathing that co-occur with head motion. These widespread respiratory fMRI signals cannot be isolated from neurobiological signals by multiecho techniques because they occur via the same BOLD mechanism. Respiratory signals must therefore be removed by some other technique to isolate neurobiological covariance in fMRI time series. Several methods for removing global artifacts are demonstrated and compared, and were found to yield fMRI time series essentially free of motion-related influences. These results identify two kinds of motion-associated fMRI variance, with different physical mechanisms and spatial profiles, each of which strongly and differentially influences functional connectivity patterns. Distance-dependent patterns in covariance are nearly entirely attributable to non-BOLD artifacts.

Resting-state functional connectivity in early postanaesthesia recovery is characterised by globally reduced anticorrelations.

BACKGROUND: A growing body of literature addresses the possible long-term cognitive effects of anaesthetics, but no study has delineated the normal trajectory of neural recovery attributable to anaesthesia alone in adults. We obtained resting-state functional MRI scans on 72 healthy human volunteers between ages 40 and 80 (median: 59) yr before, during, and after general anaesthesia with sevoflurane, in the absence of surgery, as part of a larger study on cognitive function postanaesthesia. METHODS: Region-of-interest analysis, independent component analysis, and seed-to-voxel analysis were used to characterise resting-state functional connectivity and to differentiate between correlated and anticorrelated connectivity before, during, and after general anaesthesia. RESULTS: Whilst positively correlated functional connectivity remained essentially unchanged across these perianaesthetic states, anticorrelated functional connectivity decreased globally by 35% 1 h after emergence from general anaesthesia compared with baseline, as seen by the region-of-interest analysis. This decrease corresponded to a consistent reduction in expression of canonical resting-state networks, as seen by independent component analysis. All measures returned to baseline 1 day later. CONCLUSIONS: The normal perianaesthesia trajectory of resting-state connectivity in healthy adults is characterised by a transient global reduction in anticorrelated activity shortly after emergence from anaesthesia that returns to baseline by the following day. CLINICAL TRIAL REGISTRATION: NCT02275026.

The Integration of Functional Brain Activity from Adolescence to Adulthood.

Age-related changes in human functional neuroanatomy are poorly understood. This is partly due to the limits of interpretation of standard fMRI. These limits relate to age-related variation in noise levels in data from different subjects, and the common use of standard adult brain parcellations for developmental studies. Here we used an emerging MRI approach called multiecho (ME)-fMRI to characterize functional brain changes with age. ME-fMRI acquires blood oxygenation level-dependent (BOLD) signals while also quantifying susceptibility-weighted transverse relaxation time (T2*) signal decay. This approach newly enables reliable detection of BOLD signal components at the subject level as opposed to solely at the group-average level. In turn, it supports more robust characterization of the variability in functional brain organization across individuals. We hypothesized that BOLD components in the resting state are not stable with age, and would decrease in number from adolescence to adulthood. This runs counter to the current assumptions in neurodevelopmental analyses of brain connectivity that the number of BOLD signal components is a random effect. From resting-state ME-fMRI of 51 healthy subjects of both sexes, between 8.3 and 46.2 years of age, we found a highly significant (r = -0.55, p ≪ 0.001) exponential decrease in the number of BOLD components with age. The number of BOLD components were halved from adolescence to the fifth decade of life, stabilizing in middle adulthood. The regions driving this change were dorsolateral prefrontal cortices, parietal cortex, and cerebellum. The functional network of these regions centered on the cerebellum. We conclude that an age-related decrease in BOLD component number concurs with the hypothesis of neurodevelopmental integration of functional brain activity. We show evidence that the cerebellum may play a key role in this process.SIGNIFICANCE STATEMENT Human brain development is ongoing from childhood to at least 30 years of age. Functional MRI (fMRI) is key for characterizing changes in brain function that accompany development. However, developmental fMRI studies have relied on reference maps of adult brain organization in the analysis of data from younger subjects. This approach may limit the characterization of functional activity patterns that are particular to children and adolescents. Here we used an emerging fMRI approach called multi-echo fMRI that is not susceptible to such biases when analyzing the variation in functional brain organization over development. We hypothesized an integration of the components of brain activity over development, and found that the number of components decreases exponentially, halving from 8 to 35 years of age. The brain regions most affected underlie executive function and coordination. In summary, we show major changes in the organization and integration of functional networks over development into adulthood, with both methodological and neurobiological implications for future lifespan and disease studies on brain connectivity.

Brain-behavior patterns define a dimensional biotype in medication-naïve adults with attention-deficit hyperactivity disorder.

BACKGROUND: Childhood-onset attention-deficit hyperactivity disorder (ADHD) in adults is clinically heterogeneous and commonly presents with different patterns of cognitive deficits. It is unclear if this clinical heterogeneity expresses a dimensional or categorical difference in ADHD. METHODS: We first studied differences in functional connectivity in multi-echo resting-state functional magnetic resonance imaging (rs-fMRI) acquired from 80 medication-naïve adults with ADHD and 123 matched healthy controls. We then used canonical correlation analysis (CCA) to identify latent relationships between symptoms and patterns of altered functional connectivity (dimensional biotype) in patients. Clustering methods were implemented to test if the individual associations between resting-state brain connectivity and symptoms reflected a non-overlapping categorical biotype. RESULTS: Adults with ADHD showed stronger functional connectivity compared to healthy controls, predominantly between the default-mode, cingulo-opercular and subcortical networks. CCA identified a single mode of brain-symptom co-variation, corresponding to an ADHD dimensional biotype. This dimensional biotype is characterized by a unique combination of altered connectivity correlating with symptoms of hyperactivity-impulsivity, inattention, and intelligence. Clustering analyses did not support the existence of distinct categorical biotypes of adult ADHD. CONCLUSIONS: Overall, our data advance a novel finding that the reduced functional segregation between default-mode and cognitive control networks supports a clinically important dimensional biotype of childhood-onset adult ADHD. Despite the heterogeneity of its presentation, our work suggests that childhood-onset adult ADHD is a single disorder characterized by dimensional brain-symptom mediators.

Delineating the Trajectory of Cognitive Recovery From General Anesthesia in Older Adults: Design and Rationale of the TORIE (Trajectory of Recovery in the Elderly) Project.

BACKGROUND: Mechanistic aspects of cognitive recovery after anesthesia and surgery are not yet well characterized, but may be vital to distinguishing the contributions of anesthesia and surgery in cognitive complications common in the elderly such as delirium and postoperative cognitive dysfunction. This article describes the aims and methodological approach to the ongoing study, Trajectory of Recovery in the Elderly (TORIE), which focuses on the trajectory of cognitive recovery from general anesthesia. METHODS: The study design employs cognitive testing coupled with neuroimaging techniques such as functional magnetic resonance imaging, diffusion tensor imaging, and arterial spin labeling to characterize cognitive recovery from anesthesia and its biological correlates. Applying these techniques to a cohort of age-specified healthy volunteers 40-80 years of age, who are exposed to general anesthesia alone, in the absence of surgery, will assess cognitive and functional neural network recovery after anesthesia. Imaging data are acquired before, during, and immediately after anesthesia, as well as 1 and 7 days after. Detailed cognitive data are captured at the same time points as well as 30 days after anesthesia, and brief cognitive assessments are repeated at 6 and 12 months after anesthesia. RESULTS: The study is underway. Our primary hypothesis is that older adults may require significantly longer to achieve cognitive recovery, measured by Postoperative Quality of Recovery Scale cognitive domain, than younger adults in the immediate postanesthesia period, but all will fully recover to baseline levels within 30 days of anesthesia exposure. Imaging data will address systems neuroscience correlates of cognitive recovery from general anesthesia. CONCLUSIONS: The data acquired in this project will have both clinical and theoretical relevance regardless of the outcome by delineating the mechanism behind short-term recovery across the adult age lifespan, which will have major implications for our understanding of the effects of anesthetic drugs.

Multi-echo fMRI: A review of applications in fMRI denoising and analysis of BOLD signals.

In recent years the field of fMRI research has enjoyed expanded technical abilities related to resolution, as well as use across many fields of brain research. At the same time, the field has also dealt with uncertainty related to many known and unknown effects of artifact in fMRI data. In this review we discuss an emerging fMRI technology, called multi-echo (ME)-fMRI, which focuses on improving the fidelity and interpretability of fMRI. Where the essential problem of standard single-echo fMRI is the indeterminacy of sources of signals, whether BOLD or artifact, this is not the case for ME-fMRI. By acquiring multiple echo images per slice, the ME approach allows T2* decay to be modeled at every voxel at every time point. Since BOLD signals arise by changes in T2* over time, an fMRI experiment sampling the T2* signal decay can be analyzed to distinguish BOLD from artifact signal constituents. While the ME approach has a long history of use in theoretical and validation studies, modern MRI systems enable whole-brain multi-echo fMRI at high resolution. This review covers recent multi-echo fMRI acquisition methods, and the analysis steps for this data to make fMRI at once more principled, straightforward, and powerful. After a brief overview of history and theory, T2* modeling and applications will be discussed. These applications include T2* mapping and combining echoes from ME data to increase BOLD contrast and mitigate dropout artifacts. Next, the modeling of fMRI signal changes to detect signal origins in BOLD-related T2* versus artifact-related S0 changes will be reviewed. A focus is on the use of ME-fMRI data to extract and classify components from spatial ICA, called multi-echo ICA (ME-ICA). After describing how ME-fMRI and ME-ICA lead to a general model for analysis of fMRI signals, applications in animal and human imaging will be discussed. Applications include removing motion artifacts in resting state data at subject and group level. New imaging methods such as multi-band multi-echo fMRI and imaging at 7T are demonstrated throughout the review, and a practical analysis pipeline is described. The review culminates with evidence from recent studies of major boosts in statistical power from using multi-echo fMRI for detecting activation and connectivity in healthy individuals and patients with neuropsychiatric disease. In conclusion, the review shows evidence that the multi-echo approach expands the range of experiments that is practicable using fMRI. These findings suggest a compelling future role of the multi-echo approach in subject-level and clinical fMRI.

Improving effect size estimation and statistical power with multi-echo fMRI and its impact on understanding the neural systems supporting mentalizing.

Functional magnetic resonance imaging (fMRI) research is routinely criticized for being statistically underpowered due to characteristically small sample sizes and much larger sample sizes are being increasingly recommended. Additionally, various sources of artifact inherent in fMRI data can have detrimental impact on effect size estimates and statistical power. Here we show how specific removal of non-BOLD artifacts can improve effect size estimation and statistical power in task-fMRI contexts, with particular application to the social-cognitive domain of mentalizing/theory of mind. Non-BOLD variability identification and removal is achieved in a biophysical and statistically principled manner by combining multi-echo fMRI acquisition and independent components analysis (ME-ICA). Without smoothing, group-level effect size estimates on two different mentalizing tasks were enhanced by ME-ICA at a median rate of 24% in regions canonically associated with mentalizing, while much more substantial boosts (40-149%) were observed in non-canonical cerebellar areas. Effect size boosting occurs via reduction of non-BOLD noise at the subject-level and consequent reductions in between-subject variance at the group-level. Smoothing can attenuate ME-ICA-related effect size improvements in certain circumstances. Power simulations demonstrate that ME-ICA-related effect size enhancements enable much higher-powered studies at traditional sample sizes. Cerebellar effects observed after applying ME-ICA may be unobservable with conventional imaging at traditional sample sizes. Thus, ME-ICA allows for principled design-agnostic non-BOLD artifact removal that can substantially improve effect size estimates and statistical power in task-fMRI contexts. ME-ICA could mitigate some issues regarding statistical power in fMRI studies and enable novel discovery of aspects of brain organization that are currently under-appreciated and not well understood.

Integrated strategy for improving functional connectivity mapping using multiecho fMRI.

Functional connectivity analysis of resting state blood oxygen level-dependent (BOLD) functional MRI is widely used for noninvasively studying brain functional networks. Recent findings have indicated, however, that even small (≤1 mm) amounts of head movement during scanning can disproportionately bias connectivity estimates, despite various preprocessing efforts. Further complications for interregional connectivity estimation from time domain signals include the unaccounted reduction in BOLD degrees of freedom related to sensitivity losses from high subject motion. To address these issues, we describe an integrated strategy for data acquisition, denoising, and connectivity estimation. This strategy builds on our previously published technique combining data acquisition with multiecho (ME) echo planar imaging and analysis with spatial independent component analysis (ICA), called ME-ICA, which distinguishes BOLD (neuronal) and non-BOLD (artifactual) components based on linear echo-time dependence of signals-a characteristic property of BOLD T*2 signal changes. Here we show for 32 control subjects that this method provides a physically principled and nearly operator-independent way of removing complex artifacts such as motion from resting state data. We then describe a robust estimator of functional connectivity based on interregional correlation of BOLD-independent component coefficients. This estimator, called independent components regression, considerably simplifies statistical inference for functional connectivity because degrees of freedom equals the number of independent coefficients. Compared with traditional connectivity estimation methods, the proposed strategy results in fourfold improvements in signal-to-noise ratio, functional connectivity analysis with improved specificity, and valid statistical inference with nominal control of type 1 error in contrasts of connectivity between groups with different levels of subject motion.

Separating slow BOLD from non-BOLD baseline drifts using multi-echo fMRI.

The functional magnetic resonance (fMRI) baseline is known to drift over the course of an experiment and is often attributed to hardware instability. These ultraslow fMRI fluctuations are inseparable from blood oxygenation level dependent (BOLD) changes in standard single echo fMRI and they are therefore typically removed before further analysis in both resting-state and task paradigms. However, some part of these fluctuations may be of neuronal origin, as neural activity can indeed fluctuate at the scale of several minutes or even longer, such as after the administration of drugs or during the ultradian rhythms. Here, we show that it is possible to separate the slow BOLD and non-BOLD drifts automatically using multi-echo fMRI and multi-echo independent components analysis (ME-ICA) denoising by demonstrating the detection of a visual signal evoked from a flickering checkerboard with slowly changing contrast.

Enhanced identification of BOLD-like components with multi-echo simultaneous multi-slice (MESMS) fMRI and multi-echo ICA.

The recent introduction of simultaneous multi-slice (SMS) acquisitions has enabled the acquisition of blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) data with significantly higher temporal sampling rates. In a parallel development, the use of multi-echo fMRI acquisitions in conjunction with a multi-echo independent component analysis (ME-ICA) approach has been introduced as a means to automatically distinguish functionally-related BOLD signal components from signal artifacts, with significant gains in sensitivity, statistical power, and specificity. In this work, we examine the gains that can be achieved with a combined approach in which data obtained with a multi-echo simultaneous multi-slice (MESMS) acquisition are analyzed with ME-ICA. We find that ME-ICA identifies significantly more BOLD-like components in the MESMS data as compared to data acquired with a conventional multi-echo single-slice acquisition. We demonstrate that the improved performance of MESMS derives from both an increase in the number of temporal samples and the enhanced ability to filter out high-frequency artifacts.

Differentiating BOLD and non-BOLD signals in fMRI time series from anesthetized rats using multi-echo EPI at 11.7 T.

The study of spontaneous brain activity using fMRI is central to mapping brain networks. However, current fMRI methodology has limitations in the study of small animal brain organization using ultra-high field fMRI experiments, as imaging artifacts are difficult to control and the relationship between classical neuroanatomy and spontaneous functional BOLD activity is not fully established. Challenges are especially prevalent during the fMRI study of individual rodent brains, which could be instrumental to studies of disease progression and pharmacology. A recent advance in fMRI methodology enables unbiased, accurate, and comprehensive identification of functional BOLD signals by interfacing multi-echo (ME) fMRI acquisition, NMR signal decay analysis, and independent components analysis (ICA), in a procedure called ME-ICA. Here we present a pilot study on the suitability of ME-ICA for ultra high field animal fMRI studies of spontaneous brain activity under anesthesia. ME-ICA applied to 11.7 T fMRI data of rats first showed robust performance in automatic high dimensionality estimation and ICA decomposition, similar to that previously reported for 3.0 T human data. ME sequence optimization for 11.7 T indicated that 3 echoes, 0.5mm isotropic voxel size and TR=3s was adequate for sensitive and specific BOLD signal acquisition. Next, in seeking optimal inhaled isoflurane anesthesia dosage, we report that progressive increase in anesthesia goes with concomitant decrease in statistical complexity of "global" functional activity, as measured by the number of BOLD components, or degrees of freedom (DOF). Finally, BOLD functional connectivity maps for individual rodents at the component level show that spontaneous BOLD activity follows classical neuroanatomy, and seed-based analysis shows plausible cortical-cortical and cortical-subcortical functional interactions.

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.

RAD51 deficiency disrupts the corticospinal lateralization of motor control.

Mirror movements are involuntary symmetrical movements of one side of the body that mirror voluntary movements of the other side. Congenital mirror movement disorder is a rare condition characterized by mirror movements that persist throughout adulthood in subjects with no other clinical abnormalities. The affected individuals have mirror movements predominating in the muscles that control the fingers and are unable to perform purely unimanual movements. Congenital mirror movement disorder thus provides a unique paradigm for studying the lateralization of motor control. We conducted a multimodal, controlled study of patients with congenital mirror movements associated with RAD51 haploinsufficiency (n = 7, mean age 33.3 ± 16.8 years) by comparison with age- and gender-matched healthy volunteers (n = 14, mean age 33.9 ± 16.1 years). We showed that patients with congenital mirror movements induced by RAD51 deficiency had: (i) an abnormal decussation of the corticospinal tract; (ii) abnormal interhemispheric inhibition and bilateral cortical activation of primary motor areas during intended unimanual movements; and (iii) an abnormal involvement of the supplementary motor area during both unimanual and bimanual movements. The lateralization of motor control thus requires a fine interplay between interhemispheric communication and corticospinal wiring. This fine interplay determines: (i) the delivery of appropriate motor plans from the supplementary motor area to the primary motor cortex; (ii) the lateralized activation of the primary motor cortex; and (iii) the unilateral transmission of the motor command to the limb involved in the intended movement. Our results also unveil an unexpected function of RAD51 in corticospinal development of the motor system.

Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI.

A central challenge in the fMRI based study of functional connectivity is distinguishing neuronally related signal fluctuations from the effects of motion, physiology, and other nuisance sources. Conventional techniques for removing nuisance effects include modeling of noise time courses based on external measurements followed by temporal filtering. These techniques have limited effectiveness. Previous studies have shown using multi-echo fMRI that neuronally related fluctuations are Blood Oxygen Level Dependent (BOLD) signals that can be characterized in terms of changes in R(2)* and initial signal intensity (S(0)) based on the analysis of echo-time (TE) dependence. We hypothesized that if TE-dependence could be used to differentiate BOLD and non-BOLD signals, non-BOLD signal could be removed to denoise data without conventional noise modeling. To test this hypothesis, whole brain multi-echo data were acquired at 3 TEs and decomposed with Independent Components Analysis (ICA) after spatially concatenating data across space and TE. Components were analyzed for the degree to which their signal changes fit models for R(2)* and S(0) change, and summary scores were developed to characterize each component as BOLD-like or not BOLD-like. These scores clearly differentiated BOLD-like "functional network" components from non BOLD-like components related to motion, pulsatility, and other nuisance effects. Using non BOLD-like component time courses as noise regressors dramatically improved seed-based correlation mapping by reducing the effects of high and low frequency non-BOLD fluctuations. A comparison with seed-based correlation mapping using conventional noise regressors demonstrated the superiority of the proposed technique for both individual and group level seed-based connectivity analysis, especially in mapping subcortical-cortical connectivity. The differentiation of BOLD and non-BOLD components based on TE-dependence was highly robust, which allowed for the identification of BOLD-like components and the removal of non BOLD-like components to be implemented as a fully automated procedure.