Thermal stimulus task fMRI in the cervical spinal cord at 7 Tesla.

Although functional magnetic resonance imaging (fMRI) is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0 . Increasing field strength enables higher spatial resolution and improved sensitivity to blood oxygenation level-dependent (BOLD) signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D echo-planar imaging (EPI) protocols, which differ in sensitivity to spatial and temporal B0 inhomogeneity. The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. All protocols used 3 mm slice thickness. For each protocol, the BOLD response to 13 10-s noxious thermal stimuli applied to the right thumb was acquired in a 10-min fMRI run. Image quality, temporal signal to noise ratio (SNR), and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols. Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed. Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment. PRACTITIONER POINTS: First stimulus task fMRI results in the spinal cord at 7 T. Single-shot 0.75 mm 2D EPI produced the highest mean z-statistic. Multi-shot 0.60 mm 2D EPI provided the best-localized activation and least distortion.

A Low-Intensity Transcranial Focused Ultrasound Parameter Exploration Study of the Ventral Capsule/Ventral Striatum.

OBJECTIVES: Deep brain stimulation (DBS) of the ventral capsule/ventral striatum (VC/VS) is effective for treatment-resistant obsessive-compulsive disorder (OCD); however, DBS is associated with neurosurgical risks. Transcranial focused ultrasound (tFUS) is a newer form of noninvasive (ie, nonsurgical) stimulation that can modulate deeper regions, such as the VC/VS. tFUS parameters have just begun to be studied and have often not been compared in the same participants. We explored the effects of three VC/VS tFUS protocols and an entorhinal cortex (ErC) tFUS session on the VC/VS and cortico-striato-thalamo-cortical circuit (CSTC) in healthy individuals for later application to patients with OCD. MATERIALS AND METHODS: Twelve individuals participated in a total of 48 sessions of tFUS in this exploratory multisite, within-subject parameter study. We collected resting-state, reward task, and arterial spin-labeled (ASL) magnetic resonance imaging scans before and after ErC tFUS and three VC/VS tFUS sessions with different pulse repetition frequencies (PRFs), pulse widths (PWs), and duty cycles (DCs). RESULTS: VC/VS protocol A (PRF = 10 Hz, PW = 5 ms, 5% DC) was associated with increased putamen activation during a reward task (p = 0.003), and increased VC/VS resting-state functional connectivity (rsFC) with the anterior cingulate cortex (p = 0.022) and orbitofrontal cortex (p = 0.004). VC/VS protocol C (PRF = 125 Hz, PW = 4 ms, 50% DC) was associated with decreased VC/VS rsFC with the putamen (p = 0.017), and increased VC/VS rsFC with the globus pallidus (p = 0.008). VC/VS protocol B (PRF = 125 Hz, PW = 0.4 ms, 5% DC) was not associated with changes in task-related CSTC activation or rsFC. None of the protocols affected CSTC ASL perfusion. CONCLUSIONS: This study began to explore the multidimensional parameter space of an emerging form of noninvasive brain stimulation, tFUS. Our preliminary findings in a small sample suggest that VC/VS tFUS should continue to be investigated for future noninvasive treatment of OCD.

Empirical transmit field bias correction of T1w/T2w myelin maps.

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Impact of autocalibration method on accelerated EPI of the cervical spinal cord at 7 T.

PURPOSE: The spinal cord contains sensorimotor neural circuits of scientific and clinical interest. However, spinal cord functional MRI (fMRI) is significantly more technically demanding than brain fMRI, due primarily to its proximity to the lungs. Accelerated echo-planar imaging (EPI) at 7 T is particularly vulnerable to k-space phase inconsistencies induced by motion or B0 fluctuation, during either autocalibration signal (ACS) or time-series acquisition. For 7 T brain fMRI, sensitivity to motion and B0 fluctuation can be reduced using a re-ordered segmented EPI ACS based on the fast low-angle excitation echo-planar technique (FLEET). However, respiration-induced B0 fluctuations (exceeding 100 Hz at C7) are greater, and fewer k-space lines per slice are required for cervical spinal cord fMRI at 7 T, necessitating a separate evaluation of ACS methods. METHODS: We compared 24-line single-shot EPI with 48-line two-shot segmented EPI, two-shot FLEET, and gradient echo (GRE)-based ACS acquisition methods, performed under various physiological conditions, in terms of temporal signal-to-noise ratio and prevalence of artifacts in generalized autocalibrating partially parallel acquisition (GRAPPA)-accelerated EPI of the cervical spinal cord at 7 T. RESULTS: Segmented EPI and FLEET ACS produce images with nearly identical patterns of severe image artifacts. GRE and single-shot EPI ACS consistently produce images free from significant artifacts, and temporal signal-to-noise ratio is significantly greater for GRE ACS, particularly in lower slices where through-slice dephasing is most severe. CONCLUSIONS: GRE and single-shot EPI-ACS acquisition methods, which are robust to respiration-induced phase errors between k-space segments, produce images with fewer and less severe artifacts than either FLEET or conventionally segmented EPI for accelerated EPI of the cervical spinal cord at 7 T.

Comparison of multicenter MRI protocols for visualizing the spinal cord gray matter.

PURPOSE: Spinal cord gray-matter imaging is valuable for a number of applications, but remains challenging. The purpose of this work was to compare various MRI protocols at 1.5 T, 3 T, and 7 T for visualizing the gray matter. METHODS: In vivo data of the cervical spinal cord were collected from nine different imaging centers. Data processing consisted of automatically segmenting the spinal cord and its gray matter and co-registering back-to-back scans. We computed the SNR using two methods (SNR_single using a single scan and SNR_diff using the difference between back-to-back scans) and the white/gray matter contrast-to-noise ratio per unit time. Synthetic phantom data were generated to evaluate the metrics performance. Experienced radiologists qualitatively scored the images. We ran the same processing on an open-access multicenter data set of the spinal cord MRI (N = 267 participants). RESULTS: Qualitative assessments indicated comparable image quality for 3T and 7T scans. Spatial resolution was higher at higher field strength, and image quality at 1.5 T was found to be moderate to low. The proposed quantitative metrics were found to be robust to underlying changes to the SNR and contrast; however, the SNR_single method lacked accuracy when there were excessive partial-volume effects. CONCLUSION: We propose quality assessment criteria and metrics for gray-matter visualization and apply them to different protocols. The proposed criteria and metrics, the analyzed protocols, and our open-source code can serve as a benchmark for future optimization of spinal cord gray-matter imaging protocols.

Angiogenic gene networks are dysregulated in opioid use disorder: evidence from multi-omics and imaging of postmortem human brain.

Opioid use disorder (OUD) is a public health crisis in the U.S. that causes over 50 thousand deaths annually due to overdose. Using next-generation RNA sequencing and proteomics techniques, we identified 394 differentially expressed (DE) coding and long noncoding (lnc) RNAs as well as 213 DE proteins in Brodmann Area 9 of OUD subjects. The RNA and protein changes converged on pro-angiogenic gene networks and cytokine signaling pathways. Four genes (LGALS3, SLC2A1, PCLD1, and VAMP1) were dysregulated in both RNA and protein. Dissecting these DE genes and networks, we found cell type-specific effects with enrichment in astrocyte, endothelial, and microglia correlated genes. Weighted-genome correlation network analysis (WGCNA) revealed cell-type correlated networks including an astrocytic/endothelial/microglia network involved in angiogenic cytokine signaling as well as a neuronal network involved in synaptic vesicle formation. In addition, using ex vivo magnetic resonance imaging, we identified increased vascularization in postmortem brains from a subset of subjects with OUD. This is the first study integrating dysregulation of angiogenic gene networks in OUD with qualitative imaging evidence of hypervascularization in postmortem brain. Understanding the neurovascular effects of OUD is critical in this time of widespread opioid use.

Detailed mapping of human habenula resting-state functional connectivity.

The habenula (Hb) inhibits dopaminergic reward signaling in response to negative outcomes and has been linked to numerous functional domains relevant to mental health, including reward prediction, motivation, and aversion processing. Despite its important neuroscientific and clinical implications, however, the human Hb remains poorly understood due to its small size and the associated technical hurdles to in vivo functional magnetic resonance imaging (fMRI) investigation. Using high-resolution 3 T fMRI data from 68 healthy young adults acquired through the Human Connectome Project, we developed a rigorous approach for mapping the whole-brain resting-state functional connectivity of the human Hb. Our study combined an optimized strategy for defining subject-level connectivity seeds to maximize Hb blood-oxygen-level-dependent (BOLD) signal sensitivity with high-quality surface-based alignment for robust functional localization and cortical sensitivity. We identified significant positive Hb connectivity with: (i) conserved brainstem targets, including the dopaminergic ventral tegmental area, serotonergic raphe nuclei, and periaqueductal gray; (ii) subcortical structures related to reward and motor function, including the nucleus accumbens, dorsal striatum, pallidum, thalamus, and cerebellum; and (iii) cortical areas associated with the Salience Network and early sensory processing, including the dorsal anterior cingulate, anterior insula, and primary visual and auditory cortices. Hb connectivity was strongly biased towards task-positive brain regions, with weak or negative connectivity observed throughout the task-negative Default Mode Network. Our study provides a detailed characterization of Hb resting-state functional connectivity in healthy young adults, demonstrating both the feasibility and clinical potential of studying the human Hb using high-resolution 3 T fMRI.

Incorporating non-linear alignment and multi-compartmental modeling for improved human optic nerve diffusion imaging.

In vivo human optic nerve diffusion magnetic resonance imaging (dMRI) is technically challenging with two outstanding issues not yet well addressed: (i) non-linear optic nerve movement, independent of head motion, and (ii) effect from partial-volumed cerebrospinal fluid or interstitial fluid such as in edema. In this work, we developed a non-linear optic nerve registration algorithm for improved volume alignment in axial high resolution optic nerve dMRI. During eyes-closed dMRI data acquisition, optic nerve dMRI measurements by diffusion tensor imaging (DTI) with and without free water elimination (FWE), and by diffusion basis spectrum imaging (DBSI), as well as optic nerve motion, were characterized in healthy adults at various locations along the posterior-to-anterior dimension. Optic nerve DTI results showed consistent trends in microstructural parametric measurements along the posterior-to-anterior direction of the entire intraorbital optic nerve, while the anterior portion of the intraorbital optic nerve exhibited the largest spatial displacement. Multi-compartmental dMRI modeling, such as DTI with FWE or DBSI, was less subject to spatially dependent biases in diffusivity and anisotropy measurements in the optic nerve which corresponded to similar spatial distributions of the estimated fraction of isotropic diffusion components. DBSI results derived from our clinically feasible (∼10 min) optic nerve dMRI protocol in this study are consistent with those from small animal studies, which provides the basis for evaluating the utility of multi-compartmental dMRI modeling in characterizing coexisting pathophysiology in human optic neuropathies.

Formalin tissue fixation biases myelin-sensitive MRI.

PURPOSE: Chemical fixatives such as formalin form cross-links between proteins and affect the relaxation times and diffusion properties of tissue. These fixation-induced changes likely also affect myelin density measurements produced by quantitative magnetization transfer and myelin water imaging. In this work, we evaluate these myelin-sensitive MRI methods for fixation-induced biases. METHODS: We perform quantitative magnetization transfer, myelin water imaging, and deuterium oxide-exchanged zero TE imaging on unfixed human spinal cord tissue at 9.4 Tesla and repeat these measurements after 1 day and 31 days of formalin fixation. RESULTS: The quantitative magnetization-transfer bound pool fraction increased by 30.7% ± 21.1% after 1 day of fixation and by 42.6% ± 33.9% after 31 days of fixation. Myelin water fraction increased by 39.7% ± 15.5% and 37.0% ± 15.9% at these same time points, and mean T2 of the myelin water pool nearly doubled. Reference-normalized deuterium oxide-exchanged zero TE signal intensity increased by 8.17% ± 6.03% after 31 days of fixation but did not change significantly after 1 day of fixation. After fixation, specimen cross-sectional area decreased by approximately 5%; after correction for shrinkage, changes in deuterium oxide-exchanged zero TE intensity were nearly eliminated. CONCLUSION: Bound pool fraction and myelin water fraction are significantly increased by formalin fixation, whereas deuterium oxide-exchanged zero TE intensity is minimally affected. Changes in quantitative magnetization transfer and myelin water imaging may be due in part to delamination and formation of vacuoles in the myelin sheath. Deuterium oxide-exchanged signal intensity may be altered by fixation-induced changes in myelin lipid solid-state 1 H T1 . We urge caution in the comparison of these measurements across subjects or specimens in different states, especially unfixed versus fixed tissue.

Reproducibility of myelin content-based human habenula segmentation at 3 Tesla.

In vivo morphological study of the human habenula, a pair of small epithalamic nuclei adjacent to the dorsomedial thalamus, has recently gained significant interest for its role in reward and aversion processing. However, segmenting the habenula from in vivo magnetic resonance imaging (MRI) is challenging due to the habenula's small size and low anatomical contrast. Although manual and semi-automated habenula segmentation methods have been reported, the test-retest reproducibility of the segmented habenula volume and the consistency of the boundaries of habenula segmentation have not been investigated. In this study, we evaluated the intra- and inter-site reproducibility of in vivo human habenula segmentation from 3T MRI (0.7-0.8 mm isotropic resolution) using our previously proposed semi-automated myelin contrast-based method and its fully-automated version, as well as a previously published manual geometry-based method. The habenula segmentation using our semi-automated method showed consistent boundary definition (high Dice coefficient, low mean distance, and moderate Hausdorff distance) and reproducible volume measurement (low coefficient of variation). Furthermore, the habenula boundary in our semi-automated segmentation from 3T MRI agreed well with that in the manual segmentation from 7T MRI (0.5 mm isotropic resolution) of the same subjects. Overall, our proposed semi-automated habenula segmentation showed reliable and reproducible habenula localization, while its fully-automated version offers an efficient way for large sample analysis.

7 Tesla 22-channel wrap-around coil array for cervical spinal cord and brainstem imaging.

PURPOSE: Increased signal-to-noise ratio and blood oxygenation level-dependent sensitivity at 7 Tesla (T) have the potential to enable high-resolution imaging of the human cervical spinal cord and brainstem. We propose a new two-panel radiofrequency coil design for these regions to fully exploit the advantages of ultra-high field. METHODS: A two-panel array, containing four transmit/receive and 18 receive-only elements fully encircling the head and neck, was constructed following simulations demonstrating the B1+ and specific absorption rate (SAR) benefits of two-panel over one-panel arrays. This array was compared with a previously reported posterior-only array and tested for safety using a phantom. Its anatomical, functional, and diffusion MRI performance was demonstrated in vivo. RESULTS: The two-panel array produced more uniform B1+ across the brainstem and cervical spinal cord without compromising SAR, and achieved 70% greater receive sensitivity than the posterior-only array. The two-panel design enabled acceleration of R = 2 × 2 in two dimensions or R = 3 in a single dimension. High quality in vivo anatomical, functional, and diffusion images of the human cervical spinal cord and brainstem were acquired. CONCLUSION: We have designed and constructed a wrap-around coil array with excellent performance for cervical spinal cord and brainstem MRI at 7T, which enables simultaneous human cervical spinal cord and brainstem functional MRI. Magn Reson Med 78:1623-1634, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Human habenula segmentation using myelin content.

The habenula consists of a pair of small epithalamic nuclei located adjacent to the dorsomedial thalamus. Despite increasing interest in imaging the habenula due to its critical role in mediating subcortical reward circuitry, in vivo neuroimaging research targeting the human habenula has been limited by its small size and low anatomical contrast. In this work, we have developed an objective semi-automated habenula segmentation scheme consisting of histogram-based thresholding, region growing, geometric constraints, and partial volume estimation steps. This segmentation scheme was designed around in vivo 3 T myelin-sensitive images, generated by taking the ratio of high-resolution T1w over T2w images. Due to the high myelin content of the habenula, the contrast-to-noise ratio with the thalamus in the in vivo 3T myelin-sensitive images was significantly higher than the T1w or T2w images alone. In addition, in vivo 7 T myelin-sensitive images (T1w over T2*w ratio images) and ex vivo proton density-weighted images, along with histological evidence from the literature, strongly corroborated the in vivo 3 T habenula myelin contrast used in the proposed segmentation scheme. The proposed segmentation scheme represents a step toward a scalable approach for objective segmentation of the habenula suitable for both morphological evaluation and habenula seed region selection in functional and diffusion MRI applications.

Resting-state functional connectivity of the human habenula in healthy individuals: Associations with subclinical depression.

INTRODUCTION: The habenula (Hb) is postulated to play a critical role in reward and aversion processing across species, including humans, and has been increasingly implicated in depression. However, technical constraints have limited in vivo investigation of the human Hb, and its function remains poorly characterized. We sought to overcome these challenges by examining the whole-brain resting-state functional connectivity of the Hb and its possible relationship to depressive symptomatology using the high-resolution WU-Minn Human Connectome Project (HCP) dataset. METHODS: Anatomical and resting-state functional MRI data from 50 healthy subjects with low or high subclinical depression scores (n = 25 each) were analyzed. Using novel semi-automated segmentation and optimization techniques, we generated individual-specific Hb seeds and calculated whole-brain functional connectivity for the entire cohort and the contrast of high vs. low depression groups. RESULTS: In the entire cohort, the Hb exhibited significant connectivity with key brainstem structures (i.e., ventral tegmental area, substantia nigra, pons) as well as the anterior and posterior cingulate cortices, precuneus, thalamus, and sensorimotor cortex. Multiple regions showed differential Hb connectivity based on subclinical depression scores, including the amygdala, insula, and prefrontal, mid-cingulate, and entorhinal cortices. CONCLUSIONS: Hb connectivity findings converged on areas associated with salience processing, sensorimotor systems, and the default mode network. We also detected substantial Hb-brainstem connectivity, consistent with prior histological and animal research. High and low subclinical depression groups exhibited differences in Hb connectivity with multiple regions previously linked to depression, suggesting the relationship between these structures as a potential target for future research and treatment. Hum Brain Mapp 37:2369-2384, 2016. © 2016 Wiley Periodicals, Inc.

A SEmi-Adiabatic matched-phase spin echo (SEAMS) PINS pulse-pair for B1 -insensitive simultaneous multislice imaging.

PURPOSE: Simultaneous multislice (SMS) imaging is a powerful technique that can reduce image acquisition time for anatomical, functional, and diffusion weighted magnetic resonance imaging. At higher magnetic fields, such as 7 Tesla, increased radiofrequency (RF) field inhomogeneity, power deposition, and changes in relaxation parameters make SMS spin echo imaging challenging. We designed an adiabatic 180° Power Independent of Number of Slices (PINS) pulse and a matched-phase 90° PINS pulse to generate a SEmi-Adiabatic Matched-phase Spin echo (SEAMS) PINS sequence to address these issues. METHODS: We used the adiabatic Shinnar Le-Roux (SLR) algorithm to generate a 180° pulse. The SLR polynomials for the 180° pulse were then used to create a matched-phase 90° pulse. The pulses were sub-sampled to produce a SEAMS PINS pulse-pair and the performance of this pulse-pair was validated in phantoms and in vivo. RESULTS: Simulations as well as phantom and in vivo results, demonstrate multislice capability and improved B1 -insensitivity of the SEAMS PINS pulse-pair when operating at RF amplitudes of up to 40% above adiabatic threshold. CONCLUSION: The SEAMS PINS approach presented here achieves multislice spin echo profiles with improved B1 -insensitivity when compared with a conventional spin echo.

High-Resolution Mapping of Myeloarchitecture In Vivo: Localization of Auditory Areas in the Human Brain.

The precise delineation of auditory areas in vivo remains problematic. Histological analysis of postmortem tissue indicates that the relation of areal borders to macroanatomical landmarks is variable across subjects. Furthermore, functional parcellation schemes based on measures of, for example, frequency preference (tonotopy) remain controversial. Here, we propose a 7 Tesla magnetic resonance imaging method that enables the anatomical delineation of auditory cortical areas in vivo and in individual brains, through the high-resolution visualization (0.6 × 0.6 × 0.6 mm(3)) of intracortical anatomical contrast related to myelin. The approach combines the acquisition and analysis of images with multiple MR contrasts (T1, T2*, and proton density). Compared with previous methods, the proposed solution is feasible at high fields and time efficient, which allows collecting myelin-related and functional images within the same measurement session. Our results show that a data-driven analysis of cortical depth-dependent profiles of anatomical contrast allows identifying a most densely myelinated cortical region on the medial Heschl's gyrus. Analyses of functional responses show that this region includes neuronal populations with typical primary functional properties (single tonotopic gradient and narrow frequency tuning), thus indicating that it may correspond to the human homolog of monkey A1.

Temporally-independent functional modes of spontaneous brain activity.

Resting-state functional magnetic resonance imaging has become a powerful tool for the study of functional networks in the brain. Even "at rest," the brain's different functional networks spontaneously fluctuate in their activity level; each network's spatial extent can therefore be mapped by finding temporal correlations between its different subregions. Current correlation-based approaches measure the average functional connectivity between regions, but this average is less meaningful for regions that are part of multiple networks; one ideally wants a network model that explicitly allows overlap, for example, allowing a region's activity pattern to reflect one network's activity some of the time, and another network's activity at other times. However, even those approaches that do allow overlap have often maximized mutual spatial independence, which may be suboptimal if distinct networks have significant overlap. In this work, we identify functionally distinct networks by virtue of their temporal independence, taking advantage of the additional temporal richness available via improvements in functional magnetic resonance imaging sampling rate. We identify multiple "temporal functional modes," including several that subdivide the default-mode network (and the regions anticorrelated with it) into several functionally distinct, spatially overlapping, networks, each with its own pattern of correlations and anticorrelations. These functionally distinct modes of spontaneous brain activity are, in general, quite different from resting-state networks previously reported, and may have greater biological interpretability.

ICA-based artefact removal and accelerated fMRI acquisition for improved resting state network imaging.

The identification of resting state networks (RSNs) and the quantification of their functional connectivity in resting-state fMRI (rfMRI) are seriously hindered by the presence of artefacts, many of which overlap spatially or spectrally with RSNs. Moreover, recent developments in fMRI acquisition yield data with higher spatial and temporal resolutions, but may increase artefacts both spatially and/or temporally. Hence the correct identification and removal of non-neural fluctuations is crucial, especially in accelerated acquisitions. In this paper we investigate the effectiveness of three data-driven cleaning procedures, compare standard against higher (spatial and temporal) resolution accelerated fMRI acquisitions, and investigate the combined effect of different acquisitions and different cleanup approaches. We applied single-subject independent component analysis (ICA), followed by automatic component classification with FMRIB's ICA-based X-noiseifier (FIX) to identify artefactual components. We then compared two first-level (within-subject) cleaning approaches for removing those artefacts and motion-related fluctuations from the data. The effectiveness of the cleaning procedures was assessed using time series (amplitude and spectra), network matrix and spatial map analyses. For time series and network analyses we also tested the effect of a second-level cleaning (informed by group-level analysis). Comparing these approaches, the preferable balance between noise removal and signal loss was achieved by regressing out of the data the full space of motion-related fluctuations and only the unique variance of the artefactual ICA components. Using similar analyses, we also investigated the effects of different cleaning approaches on data from different acquisition sequences. With the optimal cleaning procedures, functional connectivity results from accelerated data were statistically comparable or significantly better than the standard (unaccelerated) acquisition, and, crucially, with higher spatial and temporal resolution. Moreover, we were able to perform higher dimensionality ICA decompositions with the accelerated data, which is very valuable for detailed network analyses.

Estimation of the CSA-ODF using Bayesian compressed sensing of multi-shell HARDI.

PURPOSE: Diffusion MRI provides important information about the brain white matter structures and has opened new avenues for neuroscience and translational research. However, acquisition time needed for advanced applications can still be a challenge in clinical settings. There is consequently a need to accelerate diffusion MRI acquisitions. METHODS: A multi-task Bayesian compressive sensing (MT-BCS) framework is proposed to directly estimate the constant solid angle orientation distribution function (CSA-ODF) from under-sampled (i.e., accelerated image acquisition) multi-shell high angular resolution diffusion imaging (HARDI) datasets, and accurately recover HARDI data at higher resolution in q-space. The proposed MT-BCS approach exploits the spatial redundancy of the data by modeling the statistical relationships within groups (clusters) of diffusion signal. This framework also provides uncertainty estimates of the computed CSA-ODF and diffusion signal, directly computed from the compressive measurements. Experiments validating the proposed framework are performed using realistic multi-shell synthetic images and in vivo multi-shell high angular resolution HARDI datasets. RESULTS: Results indicate a practical reduction in the number of required diffusion volumes (q-space samples) by at least a factor of four to estimate the CSA-ODF from multi-shell data. CONCLUSION: This work presents, for the first time, a multi-task Bayesian compressive sensing approach to simultaneously estimate the full posterior of the CSA-ODF and diffusion-weighted volumes from multi-shell HARDI acquisitions. It demonstrates improvement of the quality of acquired datasets by means of CS de-noising, and accurate estimation of the CSA-ODF, as well as enables a reduction in the acquisition time by a factor of two to four, especially when "staggered" q-space sampling schemes are used. The proposed MT-BCS framework can naturally be combined with parallel MR imaging to further accelerate HARDI acquisitions.

Study protocol: The Whitehall II imaging sub-study.

BACKGROUND: The Whitehall II (WHII) study of British civil servants provides a unique source of longitudinal data to investigate key factors hypothesized to affect brain health and cognitive ageing. This paper introduces the multi-modal magnetic resonance imaging (MRI) protocol and cognitive assessment designed to investigate brain health in a random sample of 800 members of the WHII study. METHODS/DESIGN: A total of 6035 civil servants participated in the WHII Phase 11 clinical examination in 2012-2013. A random sample of these participants was included in a sub-study comprising an MRI brain scan, a detailed clinical and cognitive assessment, and collection of blood and buccal mucosal samples for the characterisation of immune function and associated measures. Data collection for this sub-study started in 2012 and will be completed by 2016. The participants, for whom social and health records have been collected since 1985, were between 60-85 years of age at the time the MRI study started. Here, we describe the pre-specified clinical and cognitive assessment protocols, the state-of-the-art MRI sequences and latest pipelines for analyses of this sub-study. DISCUSSION: The integration of cutting-edge MRI techniques, clinical and cognitive tests in combination with retrospective data on social, behavioural and biological variables during the preceding 25 years from a well-established longitudinal epidemiological study (WHII cohort) will provide a unique opportunity to examine brain structure and function in relation to age-related diseases and the modifiable and non-modifiable factors affecting resilience against and vulnerability to adverse brain changes.

Apparent diffusion coefficient values predict response to brachytherapy in bulky cervical cancer.

BACKGROUND: Diffusion-weighted magnetic resonance imaging (DWI) provides a measurement of tumor cellularity. We evaluated the potential of apparent diffusion coefficient (ADC) values obtained from post-external beam radiation therapy (EBRT) DWI and prior to brachytherapy (BT) to predict for complete metabolic response (CMR) in bulky cervical cancer. METHODS: Clinical and DWI (b value = 500 s/mm2) data were obtained from patients undergoing interstitial BT with high-risk clinical target volumes (HR-CTVs) > 30 cc. Volumes were contoured on co-registered T2 weighted images and 90th percentile ADC values were calculated. Patients were stratified by CMR (defined by PET-CT at three months post-BT). Relation of CMR with 90th percentile ADC values and other clinical factors (International Federation of Gynecology and Obstetrics (FIGO) stage, histology, tumor and HR-CTV size, pre-treatment hemoglobin, and age) was assessed both in univariate and multivariate logistic regression analyses. Youden's J statistic was used to identify a threshold value. RESULTS: Among 45 patients, twenty-eight (62%) achieved a CMR. On univariate analysis for CMR, only 90th percentile ADC value was significant (p = 0.029) while other imaging and clinical factors were not. Borderline significant factors were HR-CTV size (p = 0.054) and number of chemotherapy cycles (p = 0.078). On multivariate analysis 90th percentile ADC (p < 0.0001) and HR-CTV size (p < 0.003) were highly significant. Patients with 90th percentile ADC values above 2.10 × 10- 3 mm2/s were 5.33 (95% CI, 1.35-24.4) times more likely to achieve CMR. CONCLUSIONS: Clinical DWI may serve to risk-stratify patients undergoing interstitial BT for bulky cervical cancer.