Associations of quantitative susceptibility mapping with cortical atrophy and brain connectome in Alzheimer's disease: A multi-parametric study

Aberrant susceptibility due to iron level abnormality and brain network disconnections are observed in Alzheimer's disease (AD), with disrupted iron homeostasis hypothesized to be linked to AD pathology and neuronal loss. However, whether associations exist between abnormal quantitative susceptibility mapping (QSM), brain atrophy, and altered brain connectome in AD remains unclear. Based on multi-parametric brain imaging data from 30 AD patients and 26 healthy controls enrolled at the China-Japan Friendship Hospital, we investigated the abnormality of the QSM signal and volumetric measure across 246 brain regions in AD patients. The structural and functional connectomes were constructed based on diffusion MRI tractography and functional connectivity, respectively. The network topology was quantified using graph theory analyses. We identified seven brain regions with both reduced cortical thickness and abnormal QSM (p < 0.05) in AD, including the right superior frontal gyrus, left superior temporal gyrus, right fusiform gyrus, left superior parietal lobule, right superior parietal lobule, left inferior parietal lobule, and left precuneus. Correlations between cortical thickness and network topology computed across patients in the AD group resulted in statistically significant correlations in five of these regions, with higher correlations in functional compared to structural topology. We computed the correlation between network topological metrics, QSM value and cortical thickness across regions at both individual and group-averaged levels, resulting in a measure we call spatial correlations. We found a decrease in the spatial correlation of QSM and the global efficiency of the structural network in AD patients at the individual level. These findings may provide insights into the complex relationships among QSM, brain atrophy, and brain connectome in AD.

Immune-related lncRNAs signature and radiomics signature predict the prognosis and immune microenvironment of glioblastoma multiforme

BACKGROUND

Glioblastoma multiforme (GBM) is the most common primary malignant brain tumor in adults. This study aimed to construct immune-related long non-coding RNAs (lncRNAs) signature and radiomics signature to probe the prognosis and immune infiltration of GBM patients.

METHODS

We downloaded GBM RNA-seq data and clinical information from The Cancer Genome Atlas (TCGA) project database, and MRI data were obtained from The Cancer Imaging Archive (TCIA). Then, we conducted a cox regression analysis to establish the immune-related lncRNAs signature and radiomics signature. Afterward, we employed a gene set enrichment analysis (GSEA) to explore the biological processes and pathways. Besides, we used CIBERSORT to estimate the abundance of tumor-infiltrating immune cells (TIICs). Furthermore, we investigated the relationship between the immune-related lncRNAs signature, radiomics signature and immune checkpoint genes. Finally, we constructed a multifactors prognostic model and compared it with the clinical prognostic model.

RESULTS

We identified four immune-related lncRNAs and two radiomics features, which show the ability to stratify patients into high-risk and low-risk groups with significantly different survival rates. The risk score curves and Kaplan-Meier curves confirmed that the immune-related lncRNAs signature and radiomics signature were a novel independent prognostic factor in GBM patients. The GSEA suggested that the immune-related lncRNAs signature were involved in L1 cell adhesion molecular (L1CAM) interactions and the radiomics signature were involved signaling by Robo receptors. Besides, the two signatures was associated with the infiltration of immune cells. Furthermore, they were linked with the expression of critical immune genes and could predict immunotherapy's clinical response. Finally, the area under the curve (AUC) (0.890,0.887) and C-index (0.737,0.817) of the multifactors prognostic model were greater than those of the clinical prognostic model in both the training and validation sets, indicated significantly improved discrimination.

CONCLUSIONS

We identified the immune-related lncRNAs signature and tradiomics signature that can predict the outcomes, immune cell infiltration, and immunotherapy response in patients with GBM.

Diffusion tensor imaging of optic neuropathies: a narrative review

Background and Objective

Diffusion tensor imaging (DTI) has been implemented in a breadth of scientific investigations of optic neuropathies, though it has yet to be fully adopted for diagnosis or prognosis. This is potentially due to a lack of standardization and weak replication of results. The aim of this investigation was to review DTI results from studies specific to three distinct optic neuropathies in order to probe its current clinical utility.

Methods

We reviewed the DTI literature specific to primary open-angle glaucoma (POAG), optic neuritis (ON), and traumatic optic neuropathy (TON) by systematically searching the PubMed database on March 1st, 2023. Four distinct DTI metrics are considered: fractional anisotropy (FA), along with mean diffusivity (MD, axial diffusivity (AD), and radial diffusivity (RD). Results from within-group, between-group, and correlational studies were thoroughly assessed.

Key Content and Findings

POAG studies most consistently report a decrease in FA, especially in the optic radiations, followed in prevalence by an increase in RD and then MD, whilst AD yields conflicting results between studies. It is notable that there is not an equal distribution of investigated DTI metrics, with FA utilized the most, followed by MD, RD, and AD. Studies of ON are similar in that the most consistent findings are specific to FA, RD, and MD. These results are specific to the optic nerve and radiation since only one study measured the intermediary regions. More studies are needed to assess the effect that ON has on the tracts of the visual system. Finally, only three studies assessing DTI of TON have been performed to date, displaying low to moderate replicability of results. To improve the level of agreement between studies assessing each optic neuropathy, an increased level of standardization is recommended.

Conclusions

Both POAG and ON studies have yielded some prevalent DTI findings, both for contrast and correlation-based assessments. Although the clinical need is high for TON, considering the limitations of the current diagnostic tools, too few studies exist to make confident conclusions. Future use of standardized and longitudinal DTI, along with the foreseen methodological and technical improvements, is warranted to effectively study optic neuropathies.

Deep learning, data ramping, and uncertainty estimation for detecting artifacts in large, imbalanced databases of MRI images

Magnetic resonance imaging (MRI) is increasingly being used to delineate morphological changes underlying neurological disorders. Successfully detecting these changes depends on the MRI data quality. Unfortunately, image artifacts frequently compromise the MRI utility, making it critical to screen the data. Currently, quality assessment requires visual inspection, a time-consuming process that suffers from inter-rater variability. Automated methods to detect MRI artifacts could improve the efficiency of the process. Such automated methods have achieved high accuracy using small datasets, with balanced proportions of MRI data with and without artifacts. With the current trend towards big data in neuroimaging, there is a need for automated methods that achieve accurate detection in large and imbalanced datasets. Deep learning (DL) is the ideal MRI artifact detection algorithm for large neuroimaging databases. However, the inference generated by DL does not commonly include a measure of uncertainty. Here, we present the first stochastic DL algorithm to generate automated, high-performing MRI artifact detection implemented on a large and imbalanced neuroimaging database. We implemented Monte Carlo dropout in a 3D AlexNet to generate probabilities and epistemic uncertainties. We then developed a method to handle class imbalance, namely data-ramping to transfer the learning by extending the dataset size and the proportion of the artifact-free data instances. We used a 34,800 scans (98% clean) dataset. At baseline, we obtained 89.3% testing accuracy (F1 = 0.230). Following the transfer learning (with data-ramping), we obtained 94.9% testing accuracy (F1 = 0.357) outperforming focal cross-entropy (92.9% testing accuracy, F1 = 0.304) incorporated for comparison at handling class imbalance. By implementing epistemic uncertainties, we improved the testing accuracy to 99.5% (F1 = 0.834), outperforming the results obtained in previous comparable studies. In addition, we estimated aleatoric uncertainties by incorporating random flips to the MRI volumes, and demonstrated that aleatoric uncertainty can be implemented as part of the pipeline. The methods we introduce enhance the efficiency of managing large databases and the exclusion of artifact images from big data analyses.

Comparison of MRI radiomics-based machine learning survival models in predicting prognosis of glioblastoma multiforme

Objective

To compare the performance of radiomics-based machine learning survival models in predicting the prognosis of glioblastoma multiforme (GBM) patients.

Methods

131 GBM patients were included in our study. The traditional Cox proportional-hazards (CoxPH) model and four machine learning models (SurvivalTree, Random survival forest (RSF), DeepSurv, DeepHit) were constructed, and the performance of the five models was evaluated using the C-index.

Results

After the screening, 1792 radiomics features were obtained. Seven radiomics features with the strongest relationship with prognosis were obtained following the application of the least absolute shrinkage and selection operator (LASSO) regression. The CoxPH model demonstrated that age (HR = 1.576, p = 0.037), Karnofsky performance status (KPS) score (HR = 1.890, p = 0.006), radiomics risk score (HR = 3.497, p = 0.001), and radiomics risk level (HR = 1.572, p = 0.043) were associated with poorer prognosis. The DeepSurv model performed the best among the five models, obtaining C-index of 0.882 and 0.732 for the training and test set, respectively. The performances of the other four models were lower: CoxPH (0.663 training set / 0.635 test set), SurvivalTree (0.702/0.655), RSF (0.735/0.667), DeepHit (0.608/0.560).

Conclusion

This study confirmed the superior performance of deep learning algorithms based on radiomics relative to the traditional method in predicting the overall survival of GBM patients; specifically, the DeepSurv model showed the best predictive ability.

Evaluation of whole-brain oxygen metabolism in Alzheimer's disease using QSM and quantitative BOLD

OBJECTIVE

The objective of this study was to evaluate the whole-brain pattern of oxygen extraction fraction (OEF), cerebral blood flow (CBF), and cerebral metabolic rate of oxygen consumption (CMRO2) perturbation in Alzheimer's disease (AD) and investigate the relationship between regional cerebral oxygen metabolism and global cognition.

METHODS

Twenty-six AD patients and 25 age-matched healthy controls (HC) were prospectively recruited in this study. Mini-Mental State Examination (MMSE) was used to evaluate cognitive status. We applied the QQ-CCTV algorithm which combines quantitative susceptibility mapping and quantitative blood oxygen level-dependent models (QQ) for OEF calculation. CBF map was computed from arterial spin labeling and CMRO2 was generated based on Fick's principle. Whole-brain and regional OEF, CBF, and CMRO2 analyses were performed. The associations between these measures in substructures of deep brain gray matter and MMSE scores were assessed.

RESULTS

Whole brain voxel-wise analysis showed that CBF and CMRO2 values significantly decreased in AD predominantly in the bilateral angular gyrus, precuneus gyrus and parieto-temporal regions. Regional analysis showed that CBF value decreased in the bilateral caudal hippocampus and left rostral hippocampus and CMRO2 value decreased in left caudal and rostral hippocampus in AD patients. Considering all subjects in the AD and HC groups combined, the mean CBF and CMRO2 values in the bilateral hippocampus positively correlated with the MMSE score.

CONCLUSION

CMRO2 mapping with the QQ-CCTV method - which is readily available in MR systems for clinical practice - can be a potential biomarker for AD. In addition, CMRO2 in the hippocampus may be a useful tool for monitoring cognitive impairment.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

Effects of repetitive transcranial magnetic stimulation and their underlying neural mechanisms evaluated with magnetic resonance imaging-based brain connectivity network analyses

Repetitive transcranial magnetic stimulation (rTMS) is a noninvasive brain modulation and rehabilitation technique used in patients with neuropsychiatric diseases. rTMS can structurally remodel or functionally induce activities of specific cortical regions and has developed to an important therapeutic method in such patients. Magnetic resonance imaging (MRI) provides brain data that can be used as an explanation tool for the neural mechanisms underlying rTMS effects; brain alterations related to different functions or structures may be reflected in changes in the interaction and influence of brain connections within intrinsic specific networks. In this review, we discuss the technical details of rTMS and the biological interpretation of brain networks identified with MRI analyses, comprehensively summarize the neurobiological effects in rTMS-modulated individuals, and elaborate on changes in the brain network in patients with various neuropsychiatric diseases receiving rehabilitation treatment with rTMS. We conclude that brain connectivity network analysis based on MRI can reflect alterations in functional and structural connectivity networks comprising adjacent and separated brain regions related to stimulation sites, thus reflecting the occurrence of intrinsic functional integration and neuroplasticity. Therefore, MRI is a valuable tool for understanding the neural mechanisms of rTMS and practically tailoring treatment plans for patients with neuropsychiatric diseases.

A consensus protocol for functional connectivity analysis in the rat brain

Task-free functional connectivity in animal models provides an experimental framework to examine connectivity phenomena under controlled conditions and allows for comparisons with data modalities collected under invasive or terminal procedures. Currently, animal acquisitions are performed with varying protocols and analyses that hamper result comparison and integration. Here we introduce StandardRat, a consensus rat functional magnetic resonance imaging acquisition protocol tested across 20 centers. To develop this protocol with optimized acquisition and processing parameters, we initially aggregated 65 functional imaging datasets acquired from rats across 46 centers. We developed a reproducible pipeline for analyzing rat data acquired with diverse protocols and determined experimental and processing parameters associated with the robust detection of functional connectivity across centers. We show that the standardized protocol enhances biologically plausible functional connectivity patterns relative to previous acquisitions. The protocol and processing pipeline described here is openly shared with the neuroimaging community to promote interoperability and cooperation toward tackling the most important challenges in neuroscience.

Altered coupling of resting-state cerebral blood flow and functional connectivity in Meige syndrome

Introduction

Meige syndrome (MS) is an adult-onset segmental dystonia disease, mainly manifested as blepharospasm and involuntary movement caused by dystonic dysfunction of the oromandibular muscles. The changes of brain activity, perfusion and neurovascular coupling in patients with Meige syndrome are hitherto unknown.

Methods

Twenty-five MS patients and thirty age- and sex-matched healthy controls (HC) were prospectively recruited in this study. All the participants underwent resting-state arterial spin labeling and blood oxygen level-dependent examinations on a 3.0 T MR scanner. The measurement of neurovascular coupling was calculated using cerebral blood flow (CBF)-functional connectivity strength (FCS) correlations across the voxels of whole gray matter. Also, voxel-wised analyses of CBF, FCS, and CBF/FCS ratio images between MS and HC were conducted. Additionally, CBF and FCS values were compared between these two groups in selected motion-related brain regions.

Results

MS patients showed increased whole gray matter CBF-FCS coupling relative to HC (t = 2.262, p = 0.028). In addition, MS patients showed significantly increased CBF value in middle frontal gyrus and bilateral precentral gyrus.

Conclusion

The abnormal elevated neurovascular coupling of MS may indicate a compensated blood perfusion in motor-related brain regions and reorganized the balance between neuronal activity and brain blood supply. Our results provide a new insight into the neural mechanism underlying MS from the perspective of neurovascular coupling and cerebral perfusion.

Application of pseudocontinuous arterial spin labeling perfusion imaging in children with autism spectrum disorders

Introduction

Pseudocontinuous Arterial Spin Labeling (pCASL) perfusion imaging allows non-invasive quantification of regional cerebral blood flow (CBF) as part of a multimodal magnetic resonance imaging (MRI) protocol. This study aimed to compare regional CBF in autism spectrum disorders (ASD) individuals with their age-matched typically developing (TD) children using pCASL perfusion imaging.

Materials and methods

This cross-sectional study enrolled 17 individuals with ASD and 13 TD children. All participants underwent pCASL examination on a 3.0 T MRI scanner. Children in two groups were assessed for clinical characteristics and developmental profiles using Autism Behavior Checklist (ABC) and Gesell development diagnosis scale (GDDS), respectively. We compared CBF in different cerebral regions of ASD and TD children. We also assessed the association between CBF and clinical characteristics/developmental profile.

Results

Compared with TD children, individuals with ASD demonstrated a reduction in CBF in the left frontal lobe, the bilateral parietal lobes, and the bilateral temporal lobes. Within the ASD group, CBF was significantly higher in the right parietal lobe than in the left side. Correlation analysis of behavior characteristics and CBF in different regions showed a positive correlation between body and object domain scores on the ABC and CBF of the bilateral occipital lobes, and separately, between language domain scores and CBF of the left frontal lobe. The score of the social and self-help domain was negatively correlated with the CBF of the left frontal lobe, the left parietal lobe, and the left temporal lobe.

Conclusion

Cerebral blood flow was found to be negatively correlated with scores in the social and self-help domain, and positively correlated with those in the body and object domain, indicating that CBF values are a potential MRI-based biomarker of disease severity in ASD patients. The findings may provide novel insight into the pathophysiological mechanisms of ASD.

Resting-state functional MRI of the visual system for characterization of optic neuropathy

Optic neuropathy refers to disease of the optic nerve and can result in loss of visual acuity and/or visual field defects. Combining findings from multiple fMRI modalities can offer valuable information for characterizing and managing optic neuropathies. In this article, we review a subset of resting-state functional magnetic resonance imaging (RS-fMRI) studies of optic neuropathies. We consider glaucoma, acute optic neuritis (ON), discuss traumatic optic neuropathy (TON), and explore consistency between findings from RS and visually driven fMRI studies. Consistent with visually driven studies, glaucoma studies at rest also indicated reduced activation in the visual cortex and dorsal visual stream. RS-fMRI further reported varying levels of functional connectivity in the ventral stream depending on disease severity. ON patients show alterations within the visual cortex in both fMRI techniques. Particularly, higher-than-normal RS activity is observed in the acute phase and decreases as the disease progresses. A similar pattern is observed in the visual cortex of TON-like, open globe injury (OGI), patients. Additionally, visually driven and RS-fMRI studies of ON patients show recovery of brain activity in the visual cortex. RS-fMRI suggests recovery of signals in higher-tier visual areas MT and LOC as well. Finally, RS-fMRI has not yet been applied to TON, although reviewing OGI studies suggests that it is feasible. Future RS-fMRI studies of optic neuropathies could prioritize studying the fine scale RS activity of brain areas that visually driven studies have identified. We suggest that a more systematic longitudinal comparison of optic neuropathies with advanced fMRI would provide improved diagnostic and prognostic information.

Visually driven functional MRI techniques for characterization of optic neuropathy

Optic neuropathies are conditions that cause disease to the optic nerve, and can result in loss of visual acuity and/or visual field defects. An improved understanding of how these conditions affect the entire visual system is warranted, to better predict and/or restore the visual loss. In this article, we review visually-driven functional magnetic resonance imaging (fMRI) studies of optic neuropathies, including glaucoma and optic neuritis (ON); we also discuss traumatic optic neuropathy (TON). Optic neuropathy-related vision loss results in fMRI deficit within the visual cortex, and is often strongly correlated with clinical severity measures. Using predominantly flickering checkerboard stimuli, glaucoma studies indicated retinotopic-specific cortical alteration with more prominent deficits in advanced than in early glaucoma. Some glaucoma studies indicate a reorganized visual cortex. ON studies have indicated that the impacted cortical areas are briefly hyperactive. For ON, brain deficits are greater in the acute stages of the disease, followed by (near) normalization of responses of the LGN, visual cortex, and the dorsal visual stream, but not the ventral extrastriate cortex. Visually-driven fMRI is sensitive, at least in ON, in discriminating patients from controls, as well as the affected eye from the fellow eye within patients. The use of a greater variety of stimuli beyond checkerboards (e.g., visual motion and object recognition) in recent ON studies is encouraging, and needs to continue to disentangle the results in terms of change over time. Finally, visually-driven fMRI has not yet been applied in TON, although preliminary efforts suggest it may be feasible. Future fMRI studies of optic neuropathies should consider using more complex visual stimuli, and inter-regional analysis methods including functional connectivity. We suggest that a more systematic longitudinal comparison of optic neuropathies with advanced fMRI would provide improved diagnostic and prognostic information.

Evidence suggesting common mechanisms underlie contralateral and ipsilateral negative BOLD responses in human visual cortex

The task-evoked positive BOLD response (PBR) to a unilateral visual hemi-field stimulation is often accompanied by robust and sustained contralateral as well as ipsilateral negative BOLD responses (NBRs) in the visual cortex. The signal characteristics and the neural and/or vascular mechanisms that underlie these two types of NBRs are not completely understood. In this paper, we investigated the properties of these two types of NBRs. We first demonstrated the linearity of both NBRs with respect to stimulus duration. Next, we showed that the hemodynamic response functions (HRFs) of the two NBRs were similar to each other, but significantly different from that of the PBR. Moreover, the subject-wise expressions of the two NBRs were tightly coupled to the degree that the correlation between the two NBRs was significantly higher than the correlation between each NBR and the PBR. However, the activation patterns of the two NBRs did not show a high level of interhemispheric spatial similarity, and the functional connectivity between them was not different than the interhemispheric functional connectivity between the NBRs and PBR. Finally, while attention did modulate both NBRs, the attention-related changes in their HRFs were similar. Our findings suggest that the two NBRs might be generated through common neural and/or vascular mechanisms involving distal/deep brain regions that project to the two hemispheres.

Toward next-generation primate neuroscience: A collaboration-based strategic plan for integrative neuroimaging

Open science initiatives are creating opportunities to increase research coordination and impact in nonhuman primate (NHP) imaging. The PRIMatE Data and Resource Exchange community recently developed a collaboration-based strategic plan to advance NHP imaging as an integrative approach for multiscale neuroscience.

Optical Imaging-Based Guidance of Viral Microinjections and Insertion of a Laminar Electrophysiology Probe Into a Predetermined Barrel in Mouse Area S1BF

Wide-field Optical Imaging of Intrinsic Signals (OI-IS; Grinvald et al., 1986) is a method for imaging functional brain hemodynamic responses, mainly used to image activity from the surface of the cerebral cortex. It localizes small functional modules – such as cortical columns – with great spatial resolution and spatial specificity relative to the site of increases in neuronal activity. OI-IS is capable of imaging responses either through an intact or thinned skull or following a craniotomy. Therefore, it is minimally invasive, which makes it ideal for survival experiments. Here we describe OI-IS-based methods for guiding microinjections of optogenetics viral vectors in proximity to small functional modules (S1 barrels) of the cerebral cortex and for guiding the insertion of electrodes for electrophysiological recording into such modules. We validate our proposed methods by tissue processing of the cerebral barrel field area, revealing the track of the electrode in a predetermined barrel. In addition, we demonstrate the use of optical imaging to visualize the spatial extent of the optogenetics photostimulation, making it possible to estimate one of the two variables that conjointly determine which region of the brain is stimulated. Lastly, we demonstrate the use of OI-IS at high-magnification for imaging the upper recording contacts of a laminar probe, making it possible to estimate the insertion depth of all contacts relative to the surface of the cortex. These methods support the precise positioning of microinjections and recording electrodes, thus overcoming the variability in the spatial position of fine-scale functional modules.

Dexmedetomidine - Commonly Used in Functional Imaging Studies - Increases Susceptibility to Seizures in Rats But Not in Wild Type Mice

Functional MRI (fMRI) utilizes changes in metabolic and hemodynamic signals to indirectly infer the underlying local changes in neuronal activity. To investigate the mechanisms of fMRI responses, spontaneous fluctuations, and functional connectivity in the resting-state, it is important to pursue fMRI in animal models. Animal studies commonly use dexmedetomidine sedation. It has been demonstrated that potent sensory stimuli administered under dexmedetomidine are prone to inducing seizures in Sprague-Dawley (SD) rats. Here we combined optical imaging of intrinsic signals and cerebral blood flow with neurophysiological recordings to measure responses in rat area S1FL to electrical forepaw stimulation administered at 8 Hz. We show that the increased susceptibility to seizures starts no later than 1 h and ends no sooner than 3 h after initiating a continuous administration of dexmedetomidine. By administering different combinations of anesthetic and sedative agents, we demonstrate that dexmedetomidine is the sole agent necessary for the increased susceptibility to seizures. The increased susceptibility to seizures prevails under a combination of 0.3–0.5% isoflurane and dexmedetomidine anesthesia. The blood-oxygenation and cerebral blood flow responses to seizures induced by forepaw stimulation have a higher amplitude and a larger spatial extent relative to physiological responses to the same stimuli. The epileptic activity and the associated blood oxygenation and cerebral blood flow responses stretched beyond the stimulation period. We observed seizures in response to forepaw stimulation with 1–2 mA pulses administered at 8 Hz. In contrast, responses to stimuli administered at 4 Hz were seizure-free. We demonstrate that such seizures are generated not only in SD rats but also in Long-Evans rats, but not in C57BL6 mice stimulated with similar potent stimuli under dexmedetomidine sedation. We conclude that high-amplitude hemodynamic functional imaging responses evoked by peripheral stimulation in rats sedated with dexmedetomidine are possibly due to the induction of epileptic activity. Therefore, caution should be practiced in experiments that combine the administration of potent stimuli with dexmedetomidine sedation. We propose stimulation paradigms that elicit seizure-free, well detectable neurophysiological and hemodynamic responses in rats. We further conclude that the increased susceptibility to seizures under dexmedetomidine sedation is species dependent.

Accelerating the Evolution of Nonhuman Primate Neuroimaging

Nonhuman primate neuroimaging is on the cusp of a transformation, much in the same way its human counterpart was in 2010, when the Human Connectome Project was launched to accelerate progress. Inspired by an open data-sharing initiative, the global community recently met and, in this article, breaks through obstacles to define its ambitions.

MRI Based Brain-Specific 3D-Printed Model Aligned to Stereotactic Space for Registering Histology to MRI

Studies that seek to predict the brain microstructure based on MRI require precise alignment of processed brain histology slices to the corresponding 3D MRI data. However, achieving such alignment is a challenging problem, due to tissue distortions and the different contrasts seen in MRI and the processed tissue. Here we present a pipeline for aligning a histology volume to the MRI data of the tissue and to a stereotaxic brain atlas. To this end, we segment the volume of the brain from ex-vivo MRI data, align the MRI data obtained in the native space to an MRI stereotaxic template and create a 3D printed model (a mold or cradle) that precisely fits the brain. The pipeline then makes it possible to create grooves in the 3D model, for guiding blades for cutting slabs of tissue. Placing the brain in the brain-specific 3D printed model aligns the tissue to the MRI data by default. Aligning the MRI data to an MRI stereotaxic template makes it possible to section histology slices parallel to the standard stereotaxic axes of the atlas. This facilitates comparisons to other MRI contrasts and to images of processed tissue aligned to the standard space, while maintaining the high-resolution of the tissue images along the standard stereotaxic plane. Guiding the positioning of the grooves according to species-specific anatomical information from the co-registered atlas facilitates region-specific histology. The pipeline we introduce can be used to create brain-specific sectioning models for a variety of species, including humans, primates, and rodents. To demonstrate the generalizability of the pipeline across species, we show models generated for macaques and rats.

A rat model of somatosensory-evoked reflex seizures induced by peripheral stimulation

OBJECTIVE

We introduce a novel animal model of somatosensory stimulation-induced reflex seizures which generates focal seizures without causing damage to the brain.

METHODS

Specifically, we electrically stimulated digits or forepaws of adult rats sedated with dexmedetomidine while imaging cerebral blood volume and recording neurophysiological activity in cortical area S1FL. For the recordings, we either inserted a linear probe into the D3 digit representation or we performed surface electrocorticography (ECoG) recordings.

RESULTS

Peripheral stimulation of a digit or the forepaw elicited seizures that were followed by a refractory period with decreased neuronal activity, or another seizure or normal response. LFP amplitudes in response to electrical pulses during the seizures (0.28 ± 0.03 mV) were higher than during normal evoked responses (0.25 ± 0.05 mV) and refractory periods (0.2 ± 0.08 mV). Seizures generated during the stimulation period showed prolonged after-discharges that were sustained for 20.9 ± 1.9 s following the cessation of the stimulus. High-frequency oscillations were observed prior to and during the seizures, with amplitudes higher than those associated with normal evoked responses. The seizures were initially focal. Optical imaging of the cerebral blood volume response showed that they propagated from the onset zone to adjacent cortical areas, beyond the S1FL representation of the stimulated digit or forepaw. The spatial extent during seizures was on average 1.74 times larger during the stimulation and 4.1 times following its cessation relative to normal evoked responses. Seizures were recorded not only by probes inserted into cortex but also with ECoG arrays (24.1 ± 5.8 seizures per rat) placed over the dura matter, indicating that the seizures were not induced by damage caused by inserting the probes to the cortex. Stimulation of the forepaw elicited more seizures (18.8 ± 8.5 seizures per rat) than stimulation of a digit (1.7 ± 0.7). Unlike rats sedated with dexmedetomidine, rats anesthetized with urethane showed no seizures, indicating that the seizures may depend on the use of the mild sedative dexmedetomidine.

SIGNIFICANCE

Our proposed animal model generates seizures induced by electrical sensory stimulation free of artifacts and brain damage. It can be used for studying the mechanisms underlying the generation and propagation of reflex seizures and for evaluating antiepileptic drugs.

Using Deep Learning Algorithms to Automatically Identify the Brain MRI Contrast: Implications for Managing Large Databases

Neuroimaging science has seen a recent explosion in dataset size driving the need to develop database management with efficient processing pipelines. Multi-center neuroimaging databases consistently receive magnetic resonance imaging (MRI) data with unlabeled or incorrectly labeled contrast. There is a need to automatically identify the contrast of MRI scans to save database-managing facilities valuable resources spent by trained technicians required for visual inspection. We developed a deep learning (DL) algorithm with convolution neural network architecture to automatically infer the contrast of MRI scans based on the image intensity of multiple slices. For comparison, we developed a random forest (RF) algorithm to automatically infer the contrast of MRI scans based on acquisition parameters. The DL algorithm was able to automatically identify the MRI contrast of an unseen dataset with <0.2% error rate. The RF algorithm was able to identify the MRI contrast of the same dataset with 1.74% error rate. Our analysis showed that reduced dataset sizes caused the DL algorithm to lose generalizability. Finally, we developed a confidence measure, which made it possible to detect, with 100% specificity, all MRI volumes that were misclassified by the DL algorithm. This confidence measure can be used to alert the user on the need to inspect the small fraction of MRI volumes that are prone to misclassification. Our study introduces a practical solution for automatically identifying the MRI contrast. Furthermore, it demonstrates the powerful combination of convolution neural networks and DL for analyzing large MRI datasets.

An Open Resource for Non-human Primate Imaging

Non-human primate neuroimaging is a rapidly growing area of research that promises to transform and scale translational and cross-species comparative neuroscience. Unfortunately, the technological and methodological advances of the past two decades have outpaced the accrual of data, which is particularly challenging given the relatively few centers that have the necessary facilities and capabilities. The PRIMatE Data Exchange (PRIME-DE) addresses this challenge by aggregating independently acquired non-human primate magnetic resonance imaging (MRI) datasets and openly sharing them via the International Neuroimaging Data-sharing Initiative (INDI). Here, we present the rationale, design, and procedures for the PRIME-DE consortium, as well as the initial release, consisting of 25 independent data collections aggregated across 22 sites (total = 217 non-human primates). We also outline the unique pitfalls and challenges that should be considered in the analysis of non-human primate MRI datasets, including providing automated quality assessment of the contributed datasets.

Partial correlation analysis reveals abnormal retinotopically organized functional connectivity of visual areas in amblyopia

Amblyopia is a prevalent developmental visual disorder of childhood that typically persists in adults. Due to altered visual experience during critical periods of youth, the structure and function of adult visual cortex is abnormal. In addition to substantial deficits shown with task-based fMRI, previous studies have used resting state measures to demonstrate altered long-range connectivity in amblyopia. This is the first study in amblyopia to analyze connectivity between regions of interest that are smaller than a single cortical area and to apply partial correlation analysis to reduce network effects. We specifically assess short-range connectivity between retinotopically defined regions of interest within the occipital lobe of 8 subjects with amblyopia and 7 subjects with normal vision (aged 19–45). The representations of visual areas V1, V2, and V3 within each of the four quadrants of visual space were further subdivided into three regions based on maps of visual field eccentricity. Connectivity between pairs of all nine regions of interest in each quadrant was tested via correlation and partial correlation for both groups. Only the tests of partial correlation, i.e., correlation between time courses of two regions following the regression of time courses from all other regions, yielded significant differences between resting state functional connectivity in amblyopic and normal subjects. Subjects with amblyopia showed significantly higher partial correlation between para-foveal and more eccentric representations within V1, and this effect associated with poor acuity of the worse eye. In addition, we observed reduced correlation in amblyopic subjects between isoeccentricity regions in V1 and V2, and separately, between such regions in V2 and V3. We conclude that partial correlation-based connectivity is altered in an eccentricity-dependent pattern in visual field maps of amblyopic patients. Moreover, results are consistent with known clinical and psychophysical vision loss. More broadly, this provides evidence that abnormal cortical adaptations to disease may be better isolated with tests of partial correlation connectivity than with the regular correlation techniques that are currently widely used.

•

Cortical functional connectivity abnormalities exist in amblyopia at a scale finer than previously reported.

•

Connectivity changes within primary visual cortex are consistent with known loss of function.

•

Connectivity changes between visual areas are consistent with concept of deafferentation.

•

Partial correlation differentiates patients from controls, whereas correlation does not.

Performance of machine learning methods in diagnosing Parkinson's disease based on dysphonia measures

Parkinson's disease (PD) is a widespread degenerative syndrome that affects the nervous system. Its early appearing symptoms include tremor, rigidity, and vocal impairment (dysphonia). Consequently, speech indicators are important in the identification of PD based on dysphonic signs. In this regard, computer-aided-diagnosis systems based on machine learning can be useful in assisting clinicians in identifying PD patients. In this work, we evaluate the performance of machine learning based techniques for PD diagnosis based on dysphonia symptoms. Several machine learning techniques were considered and trained with a set of twenty-two voice disorder measurements to classify healthy and PD patients. These machine learning methods included linear discriminant analysis (LDA), k nearest-neighbors (k-NN), naïve Bayes (NB), regression trees (RT), radial basis function neural networks (RBFNN), support vector machine (SVM), and Mahalanobis distance classifier. We evaluated the performance of these methods by means of a tenfold cross validation protocol. Experimental results show that the SVM classifier achieved higher average performance than all other classifiers in terms of overall accuracy, G-mean, and area under the curve of the receiver operating characteristic plot. The SVM classifier achieved higher performance measures than the majority of the other classifiers also in terms of sensitivity, specificity, and F-measure statistics. The LDA, k-NN and RT achieved the highest average precision. The RBFNN method yielded the highest F-measure.; however, it performed poorly in terms of other performance metrics. Finally, t tests were performed to evaluate statistical significance of the results, confirming that the SVM outperformed most of the other classifiers on the majority of performance measures. SVM is a promising method for identifying PD patients based on classification of dysphonia measurements.

Spatial specificity of the functional MRI blood oxygenation response relative to neuronal activity

Previous attempts at characterizing the spatial specificity of the blood oxygenation level dependent functional MRI (BOLD fMRI) response by estimating its point-spread function (PSF) have conventionally relied on retinotopic spatial representations of visual stimuli in area V1. Consequently, their estimates were confounded by the width and scatter of receptive fields of V1 neurons. Here, we circumvent these limits by instead using the inherent cortical spatial organization of ocular dominance columns (ODCs) to determine the PSF for both Gradient Echo (GE) and Spin Echo (SE) BOLD imaging at 7 Tesla. By applying Markov chain Monte Carlo sampling on a probabilistic generative model of imaging ODCs, we quantified the PSFs that best predict the spatial structure and magnitude of differential ODCs' responses. Prior distributions for the ODC model parameters were determined by analyzing published data of cytochrome oxidase patterns from post-mortem histology of human V1 and of neurophysiological ocular dominance indices. The average PSF full-widths at half-maximum obtained from differential ODCs' responses following the removal of voxels influenced by contributions from macroscopic blood vessels were 0.86 mm (SE) and 0.99 mm (GE). Our results provide a quantitative basis for the spatial specificity of BOLD fMRI at ultra-high fields, which can be used for planning and interpretation of high-resolution differential fMRI of fine-scale cortical organizations.

A dataset of multiresolution functional brain parcellations in an elderly population with no or mild cognitive impairment

We present group eight resolutions of brain parcellations for clusters generated from resting-state functional magnetic resonance images for 99 cognitively normal elderly persons and 129 patients with mild cognitive impairment, pooled from four independent datasets. This dataset was generated as part of the following study: Common Effects of Amnestic Mild Cognitive Impairment on Resting-State Connectivity Across Four Independent Studies (Tam et al., 2015) [1]. The brain parcellations have been registered to both symmetric and asymmetric MNI brain templates and generated using a method called bootstrap analysis of stable clusters (BASC) (Bellec et al., 2010) [2]. We present two variants of these parcellations. One variant contains bihemisphereic parcels (4, 6, 12, 22, 33, 65, 111, and 208 total parcels across eight resolutions). The second variant contains spatially connected regions of interest (ROIs) that span only one hemisphere (10, 17, 30, 51, 77, 199, and 322 total ROIs across eight resolutions). We also present maps illustrating functional connectivity differences between patients and controls for four regions of interest (striatum, dorsal prefrontal cortex, middle temporal lobe, and medial frontal cortex). The brain parcels and associated statistical maps have been publicly released as 3D volumes, available in .mnc and .nii file formats on figshare and on Neurovault. Finally, the code used to generate this dataset is available on Github.

Partial Correlation-Based Retinotopically Organized Resting-State Functional Connectivity Within and Between Areas of the Visual Cortex Reflects More Than Cortical Distance

Numerous studies have demonstrated functional magnetic resonance imaging (fMRI)-based resting-state functional connectivity (RSFC) between cortical areas. Recent evidence suggests that synchronous fluctuations in blood oxygenation level-dependent fMRI reflect functional organization at a scale finer than that of visual areas. In this study, we investigated whether RSFCs within and between lower visual areas are retinotopically organized and whether retinotopically organized RSFC merely reflects cortical distance. Subjects underwent retinotopic mapping and separately resting-state fMRI. Visual areas V1, V2, and V3, were subdivided into regions of interest (ROIs) according to quadrants and visual field eccentricity. Functional connectivity (FC) was computed based on Pearson's linear correlation (correlation), and Pearson's linear partial correlation (correlation between two time courses after the time courses from all other regions in the network are regressed out). Within a quadrant, within visual areas, all correlation and nearly all partial correlation FC measures showed statistical significance. Consistently in V1, V2, and to a lesser extent in V3, correlation decreased with increasing eccentricity separation. Consistent with previously reported monkey anatomical connectivity, correlation/partial correlation values between regions from adjacent areas (V1-V2 and V2-V3) were higher than those between nonadjacent areas (V1-V3). Within a quadrant, partial correlation showed consistent significance between regions from two different areas with the same or adjacent eccentricities. Pairs of ROIs with similar eccentricity showed higher correlation/partial correlation than pairs distant in eccentricity. Between dorsal and ventral quadrants, partial correlation between common and adjacent eccentricity regions within a visual area showed statistical significance; this extended to more distant eccentricity regions in V1. Within and between quadrants, correlation decreased approximately linearly with increasing distances separating the tested ROIs. Partial correlation showed a more complex dependence on cortical distance: it decreased exponentially with increasing distance within a quadrant, but was best fit by a quadratic function between quadrants. We conclude that RSFCs within and between lower visual areas are retinotopically organized. Correlation-based FC is nonselectively high across lower visual areas, even between regions that do not share direct anatomical connections. The mechanisms likely involve network effects caused by the dense anatomical connectivity within this network and projections from higher visual areas. FC based on partial correlation, which minimizes network effects, follows expectations based on direct anatomical connections in the monkey visual cortex better than correlation. Last, partial correlation-based retinotopically organized RSFC reflects more than cortical distance effects.

1D Current Source Density (CSD) Estimation in Inverse Theory: A Unified Framework for Higher-Order Spectral Regularization of Quadrature and Expansion-Type CSD Methods

Estimation of current source density (CSD) from the low-frequency part of extracellular electric potential recordings is an unstable linear inverse problem. To make the estimation possible in an experimental setting where recordings are contaminated with noise, it is necessary to stabilize the inversion. Here we present a unified framework for zero- and higher-order singular-value-decomposition (SVD)-based spectral regularization of 1D (linear) CSD estimation from local field potentials. The framework is based on two general approaches commonly employed for solving inverse problems: quadrature and basis function expansion. We first show that both inverse CSD (iCSD) and kernel CSD (kCSD) fall into the category of basis function expansion methods. We then use these general categories to introduce two new estimation methods, quadrature CSD (qCSD), based on discretizing the CSD integral equation with a chosen quadrature rule, and representer CSD (rCSD), an even-determined basis function expansion method that uses the problem's data kernels (representers) as basis functions. To determine the best candidate methods to use in the analysis of experimental data, we compared the different methods on simulations under three regularization schemes (Tikhonov, tSVD, and dSVD), three regularization parameter selection methods (NCP, L-curve, and GCV), and seven different a priori spatial smoothness constraints on the CSD distribution. This resulted in a comparison of 531 estimation schemes. We evaluated the estimation schemes according to their source reconstruction accuracy by testing them using different simulated noise levels, lateral source diameters, and CSD depth profiles. We found that ranking schemes according to the average error over all tested conditions results in a reproducible ranking, where the top schemes are found to perform well in the majority of tested conditions. However, there is no single best estimation scheme that outperforms all others under all tested conditions. The unified framework we propose expands the set of available estimation methods, provides increased flexibility for 1D CSD estimation in noisy experimental conditions, and allows for a meaningful comparison between estimation schemes.

Common Effects of Amnestic Mild Cognitive Impairment on Resting-State Connectivity Across Four Independent Studies

Resting-state functional connectivity is a promising biomarker for Alzheimer's disease. However, previous resting-state functional magnetic resonance imaging studies in Alzheimer's disease and amnestic mild cognitive impairment (aMCI) have shown limited reproducibility as they have had small sample sizes and substantial variation in study protocol. We sought to identify functional brain networks and connections that could consistently discriminate normal aging from aMCI despite variations in scanner manufacturer, imaging protocol, and diagnostic procedure. We therefore combined four datasets collected independently, including 112 healthy controls and 143 patients with aMCI. We systematically tested multiple brain connections for associations with aMCI using a weighted average routinely used in meta-analyses. The largest effects involved the superior medial frontal cortex (including the anterior cingulate), dorsomedial prefrontal cortex, striatum, and middle temporal lobe. Compared with controls, patients with aMCI exhibited significantly decreased connectivity between default mode network nodes and between regions of the cortico-striatal-thalamic loop. Despite the heterogeneity of methods among the four datasets, we identified common aMCI-related connectivity changes with small to medium effect sizes and sample size estimates recommending a minimum of 140 to upwards of 600 total subjects to achieve adequate statistical power in the context of a multisite study with 5-10 scanning sites and about 10 subjects per group and per site. If our findings can be replicated and associated with other established biomarkers of Alzheimer's disease (e.g., amyloid and tau quantification), then these functional connections may be promising candidate biomarkers for Alzheimer's disease.

Laminar Distribution of Phase-Amplitude Coupling of Spontaneous Current Sources and Sinks

Although resting-state functional connectivity is a commonly used neuroimaging paradigm, the underlying mechanisms remain unknown. Thalamo-cortical and cortico-cortical circuits generate oscillations at different frequencies during spontaneous activity. However, it remains unclear how the various rhythms interact and whether their interactions are lamina-specific. Here we investigated intra- and inter-laminar spontaneous phase-amplitude coupling (PAC). We recorded local-field potentials using laminar probes inserted in the forelimb representation of rat area S1. We then computed time-series of frequency-band- and lamina-specific current source density (CSD), and PACs of CSD for all possible pairs of the classical frequency bands in the range of 1–150 Hz. We observed both intra- and inter-laminar spontaneous PAC. Of 18 possible combinations, 12 showed PAC, with the highest measures of interaction obtained for the pairs of the theta/gamma and delta/gamma bands. Intra- and inter-laminar PACs involving layers 2/3–5a were higher than those involving layer 6. Current sinks (sources) in the delta band were associated with increased (decreased) amplitudes of high-frequency signals in the beta to fast gamma bands throughout layers 2/3–6. Spontaneous sinks (sources) of the theta and alpha bands in layers 2/3–4 were on average linked to dipoles completed by sources (sinks) in layer 6, associated with high (low) amplitudes of the beta to fast-gamma bands in the entire cortical column. Our findings show that during spontaneous activity, delta, theta, and alpha oscillations are associated with periodic excitability, which for the theta and alpha bands is lamina-dependent. They further emphasize the differences between the function of layer 6 and that of the superficial layers, and the role of layer 6 in controlling activity in those layers. Our study links theories on the involvement of PAC in resting-state functional connectivity with previous work that revealed lamina-specific anatomical thalamo-cortico-cortical connections.

BRAIN initiative: fast and parallel solver for real-time monitoring of the eddy current in the brain for TMS applications

This paper introduce a fast and efficient solver for simulating the induced (eddy) current distribution in the brain during transcranial magnetic stimulation procedure. This solver has been integrated with MRI and neuronavigation software to accurately model the electromagnetic field and show eddy current in the head almost in real-time. To examine the performance of the proposed technique, we used a 3D anatomically accurate MRI model of the 25 year old female subject.

Variability of the coupling of blood flow and oxygen metabolism responses in the brain: a problem for interpreting BOLD studies but potentially a new window on the underlying neural activity

Recent studies from our group and others using quantitative fMRI methods have found that variations of the coupling ratio of blood flow (CBF) and oxygen metabolism (CMRO2) responses to a stimulus have a strong effect on the BOLD response. Across a number of studies an empirical pattern is emerging in the way CBF and CMRO2 changes are coupled to neural activation: if the stimulus is modulated to create a stronger response (e.g., increasing stimulus contrast), CBF is modulated more than CMRO2; on the other hand, if the brain state is altered such that the response to the same stimulus is increased (e.g., modulating attention, adaptation, or excitability), CMRO2 is modulated more than CBF. Because CBF and CMRO2 changes conflict in producing BOLD signal changes, this finding has an important implication for conventional BOLD-fMRI studies: the BOLD response exaggerates the effects of stimulus variation but is only weakly sensitive to modulations of the brain state that alter the response to a standard stimulus. A speculative hypothesis is that variability of the coupling ratio of the CBF and CMRO2 responses reflects different proportions of inhibitory and excitatory evoked activity, potentially providing a new window on neural activity in the human brain.

Quantification of the impact of a confounding variable on functional connectivity confirms anti-correlated networks in the resting-state

The effect of regressing out the global average signal (GAS) in resting state fMRI data has become a concern for interpreting functional connectivity analyses. It is not clear whether the reported anti-correlations between the Default Mode and the Dorsal Attention Networks are intrinsic to the brain, or are artificially created by regressing out the GAS. Here we introduce a concept, Impact of the Global Average on Functional Connectivity (IGAFC), for quantifying the sensitivity of seed-based correlation analyses to the regression of the GAS. This voxel-wise IGAFC index is defined as the product of two correlation coefficients: the correlation between the GAS and the fMRI time course of a voxel, times the correlation between the GAS and the seed time course. This definition enables the calculation of a threshold at which the impact of regressing-out the GAS would be large enough to introduce spurious negative correlations. It also yields a post-hoc impact correction procedure via thresholding, which eliminates spurious correlations introduced by regressing out the GAS. In addition, we introduce an Artificial Negative Correlation Index (ANCI), defined as the absolute difference between the IGAFC index and the impact threshold. The ANCI allows a graded confidence scale for ranking voxels according to their likelihood of showing artificial correlations. By applying this method, we observed regions in the Default Mode and Dorsal Attention Networks that were anti-correlated. These findings confirm that the previously reported negative correlations between the Dorsal Attention and Default Mode Networks are intrinsic to the brain and not the result of statistical manipulations. Our proposed quantification of the impact that a confound may have on functional connectivity can be generalized to global effect estimators other than the GAS. It can be readily applied to other confounds, such as systemic physiological or head movement interferences, in order to quantify their impact on functional connectivity in the resting state.

White matter abnormalities and structural hippocampal disconnections in amnestic mild cognitive impairment and Alzheimer's disease

The purpose of this project was to evaluate white matter degeneration and its impact on hippocampal structural connectivity in patients with amnestic mild cognitive impairment, non-amnestic mild cognitive impairment and Alzheimer’s disease. We estimated white matter fractional anisotropy, mean diffusivity and hippocampal structural connectivity in two independent cohorts. The ADNI cohort included 108 subjects [25 cognitively normal, 21 amnestic mild cognitive impairment, 47 non-amnestic mild cognitive impairment and 15 Alzheimer’s disease]. A second cohort included 34 subjects [15 cognitively normal and 19 amnestic mild cognitive impairment] recruited in Montreal. All subjects underwent clinical and neuropsychological assessment in addition to diffusion and T1 MRI. Individual fractional anisotropy and mean diffusivity maps were generated using FSL-DTIfit. In addition, hippocampal structural connectivity maps expressing the probability of connectivity between the hippocampus and cortex were generated using a pipeline based on FSL-probtrackX. Voxel-based group comparison statistics of fractional anisotropy, mean diffusivity and hippocampal structural connectivity were estimated using Tract-Based Spatial Statistics. The proportion of abnormal to total white matter volume was estimated using the total volume of the white matter skeleton. We found that in both cohorts, amnestic mild cognitive impairment patients had 27-29% white matter volume showing higher mean diffusivity but no significant fractional anisotropy abnormalities. No fractional anisotropy or mean diffusivity differences were observed between non-amnestic mild cognitive impairment patients and cognitively normal subjects. Alzheimer’s disease patients had 66.3% of normalized white matter volume with increased mean diffusivity and 54.3% of the white matter had reduced fractional anisotropy. Reduced structural connectivity was found in the hippocampal connections to temporal, inferior parietal, posterior cingulate and frontal regions only in the Alzheimer’s group. The severity of white matter degeneration appears to be higher in advanced clinical stages, supporting the construct that these abnormalities are part of the pathophysiological processes of Alzheimer’s disease.

Coupling of cerebral blood flow and oxygen metabolism is conserved for chromatic and luminance stimuli in human visual cortex

The ratio of the changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) during brain activation is a critical determinant of the magnitude of the blood oxygenation level dependent (BOLD) response measured with functional magnetic resonance imaging (fMRI). Cytochrome oxidase (CO), a key component of oxidative metabolism in the mitochondria, is non-uniformly distributed in visual area V1 in distinct blob and interblob regions, suggesting significant spatial variation in the capacity for oxygen metabolism. The goal of this study was to test whether CBF/CMRO(2) coupling differed when these subpopulations of neurons were preferentially stimulated, using chromatic and luminance stimuli to preferentially stimulate either the blob or interblob regions. A dual-echo spiral arterial spin labeling (ASL) technique was used to measure CBF and BOLD responses simultaneously in 7 healthy human subjects. When the stimulus contrast levels were adjusted to evoke similar CBF responses (mean 65.4% ± 19.0% and 64.6% ± 19.9%, respectively for chromatic and luminance contrast), the BOLD responses were remarkably similar (1.57% ± 0.39% and 1.59% ± 0.35%) for both types of stimuli. We conclude that CBF-CMRO(2) coupling is conserved for the chromatic and luminance stimuli used, suggesting a consistent coupling for blob and inter-blob neuronal populations despite the difference in CO concentration.

Evaluation and calibration of functional network modeling methods based on known anatomical connections

Recent studies have identified large scale brain networks based on the spatio-temporal structure of spontaneous fluctuations in resting-state fMRI data. It is expected that functional connectivity based on resting-state data is reflective of - but not identical to - the underlying anatomical connectivity. However, which functional connectivity analysis methods reliably predict the network structure remains unclear. Here we tested and compared network connectivity analysis methods by applying them to fMRI resting-state time-series obtained from the human visual cortex. The methods evaluated here are those previously tested against simulated data in Smith et al. (Neuroimage, 2011). To this end, we defined regions within retinotopic visual areas V1, V2, and V3 according to their eccentricity in the visual field, delineating central, intermediate, and peripheral eccentricity regions of interest (ROIs). These ROIs served as nodes in the models we study. We based our evaluation on the "ground-truth", thoroughly studied retinotopically-organized anatomical connectivity in the monkey visual cortex. For each evaluated method, we computed the fractional rate of detecting connections known to exist ("c-sensitivity"), while using a threshold of the 95th percentile of the distribution of interaction magnitudes of those connections not expected to exist. Under optimal conditions - including session duration of 68min, a relatively small network consisting of 9 nodes and artifact-free regression of the global effect - each of the top methods predicted the expected connections with 67-85% c-sensitivity. Correlation methods, including Correlation (Corr; 85%), Regularized Inverse Covariance (ICOV; 84%) and Partial Correlation (PCorr; 81%) performed best, followed by Patel's Kappa (80%), Bayesian Network method PC (BayesNet; 77%), General Synchronization measures (67-77%), and Coherence (CohB; 74%). With decreased session duration, these top methods saw decreases in c-sensitivities, achieving 59-76% for 17min sessions. With a short resting-state fMRI scan of 8.5min, none of the methods predicted the real network well, with Corr (65%) performing best. With increased complexity of the network from 9 to 36 nodes, multivariate methods including PCorr and BayesNet saw a decrease in performance. Artifact-free regression of the global effect increased the c-sensitivity of the top-performing methods. In an overall evaluation across all tests we performed, correlation methods (Corr, ICOV, and PCorr), Patel's Kappa, and BayesNet method PC set themselves somewhat above all other methods. We propose that data-based calibration based on known anatomical connections be integrated into future network studies, in order to maximize sensitivity and reduce false positives.

Global and system-specific resting-state fMRI fluctuations are uncorrelated: principal component analysis reveals anti-correlated networks

The influence of the global average signal (GAS) on functional-magnetic resonance imaging (fMRI)-based resting-state functional connectivity is a matter of ongoing debate. The global average fluctuations increase the correlation between functional systems beyond the correlation that reflects their specific functional connectivity. Hence, removal of the GAS is a common practice for facilitating the observation of network-specific functional connectivity. This strategy relies on the implicit assumption of a linear-additive model according to which global fluctuations, irrespective of their origin, and network-specific fluctuations are super-positioned. However, removal of the GAS introduces spurious negative correlations between functional systems, bringing into question the validity of previous findings of negative correlations between fluctuations in the default-mode and the task-positive networks. Here we present an alternative method for estimating global fluctuations, immune to the complications associated with the GAS. Principal components analysis was applied to resting-state fMRI time-series. A global-signal effect estimator was defined as the principal component (PC) that correlated best with the GAS. The mean correlation coefficient between our proposed PC-based global effect estimator and the GAS was 0.97±0.05, demonstrating that our estimator successfully approximated the GAS. In 66 out of 68 runs, the PC that showed the highest correlation with the GAS was the first PC. Since PCs are orthogonal, our method provides an estimator of the global fluctuations, which is uncorrelated to the remaining, network-specific fluctuations. Moreover, unlike the regression of the GAS, the regression of the PC-based global effect estimator does not introduce spurious anti-correlations beyond the decrease in seed-based correlation values allowed by the assumed additive model. After regressing this PC-based estimator out of the original time-series, we observed robust anti-correlations between resting-state fluctuations in the default-mode and the task-positive networks. We conclude that resting-state global fluctuations and network-specific fluctuations are uncorrelated, supporting a Resting-State Linear-Additive Model. In addition, we conclude that the network-specific resting-state fluctuations of the default-mode and task-positive networks show artifact-free anti-correlations.

Mapping the organization of axis of motion selective features in human area MT using high-field fMRI

Functional magnetic resonance imaging (fMRI) at high magnetic fields has made it possible to investigate the columnar organization of the human brain in vivo with high degrees of accuracy and sensitivity. Until now, these results have been limited to the organization principles of early visual cortex (V1). While the middle temporal area (MT) has been the first identified extra-striate visual area shown to exhibit a columnar organization in monkeys, evidence of MT's columnar response properties and topographic layout in humans has remained elusive. Research using various approaches suggests similar response properties as in monkeys but failed to provide direct evidence for direction or axis of motion selectivity in human area MT. By combining state of the art pulse sequence design, high spatial resolution in all three dimensions (0.8 mm isotropic), optimized coil design, ultrahigh field magnets (7 Tesla) and novel high resolution cortical grid sampling analysis tools, we provide the first direct evidence for large-scale axis of motion selective feature organization in human area MT closely matching predictions from topographic columnar-level simulations.

Modeling and analysis of mechanisms underlying fMRI-based decoding of information conveyed in cortical columns

Multivariate machine learning algorithms applied to human functional MRI (fMRI) data can decode information conveyed by cortical columns, despite the voxel-size being large relative to the width of columns. Several mechanisms have been proposed to underlie decoding of stimulus orientation or the stimulated eye. These include: (I) aliasing of high spatial-frequency components, including the main frequency component of the columnar organization, (II) contributions from local irregularities in the columnar organization, (III) contributions from large-scale non-columnar organizations, (IV) functionally selective veins with biased draining regions, and (V) complex spatio-temporal filtering of neuronal activity by fMRI voxels. Here we sought to assess the plausibility of two of the suggested mechanisms: (I) aliasing and (II) local irregularities, using a naive model of BOLD as blurring and MRI voxel sampling. To this end, we formulated a mathematical model that encompasses both the processes of imaging ocular dominance (OD) columns and the subsequent linear classification analysis. Through numerical simulations of the model, we evaluated the distribution of functional differential contrasts that can be expected when considering the pattern of cortical columns, the hemodynamic point spread function, the voxel size, and the noise. We found that with data acquisition parameters used at 3 Tesla, sub-voxel supra-Nyquist frequencies, including frequencies near the main frequency of the OD organization (0.5 cycles per mm), cannot contribute to the differential contrast. The differential functional contrast of local origin is dominated by low-amplitude contributions from low frequencies, associated with irregularities of the cortical pattern. Realizations of the model with parameters that reflected best-case scenario and the reported BOLD point-spread at 3 Tesla (3.5mm) predicted decoding performances lower than those that have been previously obtained at this magnetic field strength. We conclude that low frequency components that underlie local irregularities in the columnar organization are likely to play a role in decoding. We further expect that fMRI-based decoding relies, in part, on signal contributions from large-scale, non-columnar functional organizations, and from complex spatio-temporal filtering of neuronal activity by fMRI voxels, involving biased venous responses. Our model can potentially be used for evaluating and optimizing data-acquisition parameters for decoding information conveyed by cortical columns.

Neuronal correlates of spontaneous fluctuations in fMRI signals in monkey visual cortex: Implications for functional connectivity at rest

Recent studies have demonstrated large amplitude spontaneous fluctuations in functional-MRI (fMRI) signals in humans in the resting state. Importantly, these spontaneous fluctuations in blood-oxygenation-level-dependent (BOLD) signal are often synchronized over distant parts of the brain, a phenomenon termed functional-connectivity. Functional-connectivity is widely assumed to reflect interregional coherence of fluctuations in activity of the underlying neuronal networks. Despite the large body of human imaging literature on spontaneous activity and functional-connectivity in the resting state, the link to underlying neural activity remains tenuous. Through simultaneous fMRI and intracortical neurophysiological recording, we demonstrate correlation between slow fluctuations in BOLD signals and concurrent fluctuations in the underlying locally measured neuronal activity. This correlation varied with time-lag of BOLD relative to neuronal activity, resembling a traditional hemodynamic response function with peaks at approximately 6 s lag of BOLD signal. The correlations were reliably detected when the neuronal signal consisted of either the spiking rate of a small group of neurons, or relative power changes in the multi-unit activity band, and particularly in the local field potential gamma band. Analysis of correlation between the voxel-by-voxel fMRI time-series and the neuronal activity measured within one cortical site showed patterns of correlation that slowly traversed cortex. BOLD fluctuations in widespread areas in visual cortex of both hemispheres were significantly correlated with neuronal activity from a single recording site in V1. To the extent that our V1 findings can be generalized to other cortical areas, fMRI-based functional-connectivity between remote regions in the resting state can be linked to synchronization of slow fluctuations in the underlying neuronal signals.