High-sensitivity TMS/fMRI of the Human Motor Cortex Using a Dedicated Multichannel MR Coil

PURPOSE

To validate a novel setup for concurrent TMS/fMRI in the human motor cortex based on a dedicated, ultra-thin, multichannel receive MR coil positioned between scalp and TMS system providing greatly enhanced sensitivity compared to the standard birdcage coil setting.

METHODS

A combined TMS/fMRI design was applied over the primary motor cortex based on 1Hz stimulation with stimulation levels of 80%, 90%, 100%, and 110% of the individual active motor threshold, respectively. Due to the use of a multichannel receive coil we were able to use multiband-accelerated (MB=2) EPI sequences for the acquisition of functional images. Data were analysed with SPM12 and BOLD-weighted signal intensity time courses were extracted in each subject from two local maxima (individual functional finger tapping localiser, fixed MNI coordinate of the hand knob) next to the hand area of the primary motor cortex (M1) and from the global maximum.

RESULTS

We report excellent image quality without noticeable signal dropouts or image distortions. Parameter estimates in the three peak voxels showed monotonically ascending activation levels over increasing stimulation intensities. Across all subjects, mean BOLD signal changes for 80%, 90%, 100%, 110% of the individual active motor threshold were 0.43%, 0.63%, 1.01%, 2.01% next to the individual functional finger tapping maximum, 0.73%, 0.91%, 1.34%, 2.21% next to the MNI-defined hand knob and 0.88%, 1.09%, 1.65%, 2.77% for the global maximum, respectively.

CONCLUSION

Our results show that the new setup for concurrent TMS/fMRI experiments using a dedicated MR coil array allows for high-sensitivity fMRI particularly at the site of stimulation. Contrary to the standard birdcage approach, the results also demonstrate that the new coil can be successfully used for multiband-accelerated EPI acquisition. The gain in flexibility due to the new coil can be easily combined with neuronavigation within the MR scanner to allow for accurate targeting in TMS/fMRI experiments.

Normalised time-to-peak-distribution curves correlate with cerebral white matter hyperintensities - Could this improve early diagnosis?

Parameter-free assessment of the time-to-peak (TTP) histogram, termed 'TTP-distribution curve' (TDC), of dynamic susceptibility contrast-enhanced magnetic resonance imaging (DSC-MRI) was introduced as a robust method to evaluate cerebral perfusion. TDC-assessment works fully automatically without the need of an arterial input function, thereby providing full comparability between different measurements. In the investigated sample of 106 patients, a strong dependency of TDC on the hemodynamic state of cerebral microvessels and the arterio-venous bolus-transit time [Formula: see text] was demonstrated. Accordingly, TDC-derived [Formula: see text] was 3.3-3.7 s for control patients and 4.4 s for cerebral small vessel disease patients. Measurements of associated bolus spread velocities ν and accelerations [Formula: see text] additionally revealed a direct effect from spin-spin relaxation time T₂-weighted white matter hyperintensity volume, considered to indicate microangiopathy in cerebral small vessel disease, on the TDC-measurements. This strongly supports the prevailing hypothesis that cerebral small vessel disease directly influences DSC-measurements, where the degree could be estimated from an analysis of TDC. While this may be used to correct DSC-parameters for undesirable effects from cerebral small vessel disease, it could also serve to potentially identify patients at risk for cerebral small vessel disease at an early stage, since a subset of patients without yet significant WHM-volume, but clearly altered hemodynamics in TDC-measurements, was identified in this study.

Multi-turn multi-gap transmission line resonators - Concept, design and first implementation at 4.7T and 7T

A novel design scheme for monolithic transmission line resonators (TLRs) is presented - the multi-turn multi-gap TLR (MTMG-TLR) design. The MTMG-TLR design enables the construction of TLRs with multiple turns and multiple gaps. This presents an additional degree of freedom in tuning self-resonant TLRs, as their resonance frequency is fully determined by the coil geometry (e.g. diameter, number of turns, conductor width, etc.). The novel design is evaluated at 4.7T and 7T by simulations and experiments, where it is demonstrated that MTMG-TLRs can be used for MRI, and that the B1 distribution of MTMG-TLRs strongly depends on the number and distribution of turns. A comparison to conventional loop coils revealed that the B1 performance of MTMG-TLRs is comparable to a loop coil with the same mean diameter; however, lower 10g SAR values were found for MTMG-TLRs. The MTMG-TLR design is expected to bring most benefits at high static field, where it allows for independent size and frequency selection, which cannot be achieved with standard TLR design. However, it also enables more accurate geometric optimization at low static field. Thereby, the MTMG-TLR design preserves the intrinsic advantages of TLRs, i.e. mechanical flexibility, high SAR efficiency, mass production, and coil miniaturization.

Simultaneous and interleaved acquisition of NMR signals from different nuclei with a clinical MRI scanner

PURPOSE

Modification of a clinical MRI scanner to enable simultaneous or rapid interleaved acquisition of signals from two different nuclei.

METHODS

A device was developed to modify the local oscillator signal fed to the receive channel(s) of an MRI console. This enables external modification of the frequency at which the receiver is sensitive and rapid switching between different frequencies. Use of the device was demonstrated with interleaved and simultaneous ³¹ P and ¹ H spectroscopic acquisitions, and with interleaved ³¹ P and ¹ H imaging.

RESULTS

Signal amplitudes and signal-to-noise ratios were found to be unchanged for the modified system, compared with data acquired with the MRI system in the standard configuration.

CONCLUSION

Interleaved and simultaneous ¹ H and ³¹ P signal acquisition was successfully demonstrated with a clinical MRI scanner, with only minor modification of the RF architecture. While demonstrated with ³¹ P, the modification is applicable to any detectable nucleus without further modification, enabling a wide range of simultaneous and interleaved experiments to be performed within a clinical setting. Magn Reson Med 76:1636-1641, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Sex differences in the functional connectivity of the amygdalae in association with cortisol

Human amygdalae are involved in various behavioral functions such as affective and stress processing. For these behavioral functions, as well as for psychophysiological arousal including cortisol release, sex differences are reported. Here, we assessed cortisol levels and resting-state functional connectivity (rsFC) of left and right amygdalae in 81 healthy participants (42 women) to investigate potential modulation of amygdala rsFC by sex and cortisol concentration. Our analyses revealed that rsFC of the left amygdala significantly differed between women and men: Women showed stronger rsFC than men between the left amygdala and left middle temporal gyrus, inferior frontal gyrus, postcentral gyrus and hippocampus, regions involved in face processing, inner-speech, fear and pain processing. No stronger connections were detected for men and no sex difference emerged for right amygdala rsFC. Also, an interaction of sex and cortisol appeared: In women, cortisol was negatively associated with rsFC of the amygdalae with striatal regions, mid-orbital frontal gyrus, anterior cingulate gyrus, middle and superior frontal gyri, supplementary motor area and the parietal-occipital sulcus. Contrarily in men, positive associations of cortisol with rsFC of the left amygdala and these structures were observed. Functional decoding analyses revealed an association of the amygdalae and these regions with emotion, reward and memory processing, as well as action execution. Our results suggest that functional connectivity of the amygdalae as well as the regulatory effect of cortisol on brain networks differs between women and men. These sex-differences and the mediating and sex-dependent effect of cortisol on brain communication systems should be taken into account in affective and stress-related neuroimaging research. Thus, more studies including both sexes are required.

Big Data Approaches for the Analysis of Large-Scale fMRI Data Using Apache Spark and GPU Processing: A Demonstration on Resting-State fMRI Data from the Human Connectome Project

Technologies for scalable analysis of very large datasets have emerged in the domain of internet computing, but are still rarely used in neuroimaging despite the existence of data and research questions in need of efficient computation tools especially in fMRI. In this work, we present software tools for the application of Apache Spark and Graphics Processing Units (GPUs) to neuroimaging datasets, in particular providing distributed file input for 4D NIfTI fMRI datasets in Scala for use in an Apache Spark environment. Examples for using this Big Data platform in graph analysis of fMRI datasets are shown to illustrate how processing pipelines employing it can be developed. With more tools for the convenient integration of neuroimaging file formats and typical processing steps, big data technologies could find wider endorsement in the community, leading to a range of potentially useful applications especially in view of the current collaborative creation of a wealth of large data repositories including thousands of individual fMRI datasets.

Identification of Voxels Confounded by Venous Signals Using Resting-State fMRI Functional Connectivity Graph Community Identification

Identifying venous voxels in fMRI datasets is important to increase the specificity of fMRI analyses to microvasculature in the vicinity of the neural processes triggering the BOLD response. This is, however, difficult to achieve in particular in typical studies where magnitude images of BOLD EPI are the only data available. In this study, voxelwise functional connectivity graphs were computed on minimally preprocessed low TR (333 ms) multiband resting-state fMRI data, using both high positive and negative correlations to define edges between nodes (voxels). A high correlation threshold for binarization ensures that most edges in the resulting sparse graph reflect the high coherence of signals in medium to large veins. Graph clustering based on the optimization of modularity was then employed to identify clusters of coherent voxels in this graph, and all clusters of 50 or more voxels were then interpreted as corresponding to medium to large veins. Indeed, a comparison with SWI reveals that 75.6±5.9% of voxels within these large clusters overlap with veins visible in the SWI image or lie outside the brain parenchyma. Some of the remaining differences between the two modalities can be explained by imperfect alignment or geometric distortions between the two images. Overall, the graph clustering based method for identifying venous voxels has a high specificity as well as the additional advantages of being computed in the same voxel grid as the fMRI dataset itself and not needing any additional data beyond what is usually acquired (and exported) in standard fMRI experiments.

Simultaneous and interleaved acquisition of NMR signals from different nuclei with a clinical MRI scanner

Purpose

Modification of a clinical MRI scanner to enable simultaneous or rapid interleaved acquisition of signals from two different nuclei.

Methods

A device was developed to modify the local oscillator signal fed to the receive channel(s) of an MRI console. This enables external modification of the frequency at which the receiver is sensitive and rapid switching between different frequencies. Use of the device was demonstrated with interleaved and simultaneous ³¹P and ¹H spectroscopic acquisitions, and with interleaved ³¹P and ¹H imaging.

Results

Signal amplitudes and signal‐to‐noise ratios were found to be unchanged for the modified system, compared with data acquired with the MRI system in the standard configuration.

Conclusion

Interleaved and simultaneous ¹H and ³¹P signal acquisition was successfully demonstrated with a clinical MRI scanner, with only minor modification of the RF architecture. While demonstrated with ³¹P, the modification is applicable to any detectable nucleus without further modification, enabling a wide range of simultaneous and interleaved experiments to be performed within a clinical setting. Magn Reson Med 76:1636–1641, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

A functional polymorphism in the prodynorphin gene affects cognitive flexibility and brain activation during reversal learning

Whether the opioid system plays a role in the ability to flexibly adapt behavior is still unclear. We used fMRI to investigate the effect of a nucleotide tandem repeat (68-bp VNTR) functional polymorphism of the prodynorphin (PDYN) gene on cerebral activation during a reversal learning task in which participants had to flexibly adapt stimulus-response associations. Past studies suggested that alleles with 3 or 4 repeats (HH genotype) of this polymorphism are associated with higher levels of dynorphin peptides than alleles with 1 or 2 repeats (LL genotype). On the behavioral level, the HH group made more perseverative errors than the LL group. On the neural level, the HH group demonstrated less engagement of left orbitofrontal cortex (lOFC) and cortico-striatal circuitry, and lower effective connectivity of lOFC with anterior midcingulate cortex and anterior insula/ventrolateral prefrontal cortex during reversal learning and processing negative feedback. This points to a lower ability of the HH genotype to monitor or adapt to changes in reward contingencies. These findings provide first evidence that dynorphins may contribute to individual differences in reversal learning, and that considering the opioid system may shed new light on the neurochemical correlates of decision-making and behavioral regulation.

Dynamic PCr and pH imaging of human calf muscles during exercise and recovery using (31) P gradient-Echo MRI at 7 Tesla

PURPOSE

Simultaneous acquisition of spatially resolved (31) P-MRI data for evaluation of muscle specific energy metabolism, i.e., PCr and pH kinetics.

METHODS

A three-dimensional (3D) gradient-echo sequence for multiple frequency-selective excitations of the PCr and Pi signals in an interleaved sampling scheme was developed and tested at 7 Tesla (T). The pH values were derived from the chemical shift-induced phase difference between the resonances. The achieved spatial resolution was ∼2 mL with image acquisition time below 6 s. Ten healthy volunteers were studied performing plantar flexions during the delay between (31) P-MRI acquisitions, yielding a temporal resolution of 9-10 s.

RESULTS

Signal from anatomically matched regions of interest had sufficient signal-to-noise ratio to allow single-acquisition PCr and pH quantification. The Pi signal was clearly detected in voxels of actively exercising muscles. The PCr depletions were in gastrocnemius 42 ± 14% (medialis), 48 ± 17% (lateralis) and in soleus 20 ± 11%. The end exercise pH values were 6.74 ± 0.18 and 6.65 ± 0.27 for gastrocnemius medialis and lateralis, respectively, and 6.96 ± 0.12 for soleus muscle.

CONCLUSION

Simultaneous acquisition of PCr and Pi images with high temporal resolution, suitable for measuring PCr and pH kinetics in exercise-recovery experiments, was demonstrated at 7T. This study presents a fast alternative to MRS for quantifying energy metabolism of posterior muscle groups of the lower leg. Magn Reson Med 75:2324-2331, 2016. © 2015 Wiley Periodicals, Inc.

fMRI measurements of amygdala activation are confounded by stimulus correlated signal fluctuation in nearby veins draining distant brain regions

Imaging the amygdala with functional MRI is confounded by multiple averse factors, notably signal dropouts due to magnetic inhomogeneity and low signal-to-noise ratio, making it difficult to obtain consistent activation patterns in this region. However, even when consistent signal changes are identified, they are likely to be due to nearby vessels, most notably the basal vein of rosenthal (BVR). Using an accelerated fMRI sequence with a high temporal resolution (TR = 333 ms) combined with susceptibility-weighted imaging, we show how signal changes in the amygdala region can be related to a venous origin. This finding is confirmed here in both a conventional fMRI dataset (TR = 2000 ms) as well as in information of meta-analyses, implying that “amygdala activations” reported in typical fMRI studies are likely confounded by signals originating in the BVR rather than in the amygdala itself, thus raising concerns about many conclusions on the functioning of the amygdala that rely on fMRI evidence alone.

Reduced default mode network suppression during a working memory task in remitted major depression

Insufficient default mode network (DMN) suppression was linked to increased rumination in symptomatic Major Depressive Disorder (MDD). Since rumination is known to predict relapse and a more severe course of MDD, we hypothesized that similar DMN alterations might also exist during full remission of MDD (rMDD), a condition known to be associated with increased relapse rates specifically in patients with adolescent onset. Within a cross-sectional functional magnetic resonance imaging study activation and functional connectivity (FC) were investigated in 120 adults comprising 78 drug-free rMDD patients with adolescent- (n = 42) and adult-onset (n = 36) as well as 42 healthy controls (HC), while performing the n-back task. Compared to HC, rMDD patients showed diminished DMN deactivation with strongest differences in the anterior-medial prefrontal cortex (amPFC), which was further linked to increased rumination response style. On a brain systems level, rMDD patients showed an increased FC between the amPFC and the dorsolateral prefrontal cortex, which constitutes a key region of the antagonistic working-memory network. Both whole-brain analyses revealed significant differences between adolescent-onset rMDD patients and HC, while adult-onset rMDD patients showed no significant effects. Results of this study demonstrate that reduced DMN suppression exists even after full recovery of depressive symptoms, which appears to be specifically pronounced in adolescent-onset MDD patients. Our results encourage the investigation of DMN suppression as a putative predictor of relapse in clinical trials, which might eventually lead to important implications for antidepressant maintenance treatment.

Localized semi-LASER dynamic (31)P magnetic resonance spectroscopy of the soleus during and following exercise at 7 T

OBJECTIVES

This study demonstrates the applicability of semi-LASER localized dynamic (31)P MRS to deeper lying areas of the exercising human soleus muscle (SOL). The effect of accurate localization and high temporal resolution on data specificity is investigated.

MATERIALS AND METHODS

To achieve high signal-to-noise ratio (SNR) at a temporal resolution of 6 s, a custom-built human calf coil array was used at 7T. The kinetics of phosphocreatine (PCr) and intracellular pH were quantified separately in SOL and gastrocnemius medialis (GM) muscle of nine volunteers, during rest, plantar flexion exercise, and recovery.

RESULTS

The average SNR of PCr at rest was [Formula: see text] in SOL ([Formula: see text] in GM). End exercise PCr depletion in SOL ([Formula: see text] %) was far lower than in GM ([Formula: see text] %). The pH in SOL increased rapidly and, in contrast to GM, remained elevated until the end of exercise.

CONCLUSION

(31)P MRS in single-shots every 6 s localized in the deeper-lying SOL enabled quantification of PCr recovery times at low depletions and of fast pH changes, like the initial rise. Both high temporal resolution and accurate spatial localization improve specificity of Pi and, thus, pH quantification by avoiding multiple, and potentially indistinguishable sources for changing the Pi peak shape.

Unexplained high BMD in DXA-scanned patients is generalized throughout the skeleton and characterized by thicker cortical and trabecular bone

Unexplained high bone mineral density (BMD) is a rare condition and the mechanisms responsible are yet to be described in detail. The aim of the study was to identify patients with unexplained high BMD from a local DXA database and compare their radiological phenotype with an age- and a gender-matched group of population-based controls. We defined high BMD as a DXA Z-score ≥ + 2.5 at the total hip and lumbar spine. We characterized the findings as "unexplained" if no osteodegenerative changes, bone metabolic disease, or arthritis at the hip or lumbar spine was observed. All participants were investigated with high-resolution peripheral quantitative computed tomography (HR-pQCT), QCT, DXA, fasting blood samples, a 24-h urine sample, and questionnaires. The DXA database contained data on 25,118 patients. Initially, 138 (0.55%) potential participants with high BMD were identified, and during the study ten additional cases were identified from new DXA scans. Sixty-seven patients accepted to participate in the study, and among these we identified 15 women and one man with unexplained high BMD. These 15 women had higher BMD throughout the skeleton relative to controls, similar area/volume at the hip and the distal extremities, a higher number of trabeculae, which was thicker than in the controls, and a higher finite element estimated bone strength. The 15 women were heavier and had a higher fat mass then controls. We conclude that patients with unexplained high BMD have a generalized high BMD phenotype throughout their skeleton, which is characterized with a denser microarchitecture.

Disrupted effective connectivity between the amygdala and orbitofrontal cortex in social anxiety disorder during emotion discrimination revealed by dynamic causal modeling for FMRI

Social anxiety disorder (SAD) is characterized by over-reactivity of fear-related circuits in social or performance situations and associated with marked social impairment. We used dynamic causal modeling (DCM), a method to evaluate effective connectivity, to test our hypothesis that SAD patients would exhibit dysfunctions in the amygdala-prefrontal emotion regulation network. Thirteen unmedicated SAD patients and 13 matched healthy controls performed a series of facial emotion and object discrimination tasks while undergoing fMRI. The emotion-processing network was identified by a task-related contrast and motivated the selection of the right amygdala, OFC, and DLPFC for DCM analysis. Bayesian model averaging for DCM revealed abnormal connectivity between the OFC and the amygdala in SAD patients. In healthy controls, this network represents a negative feedback loop. In patients, however, positive connectivity from OFC to amygdala was observed, indicating an excitatory connection. As we did not observe a group difference of the modulatory influence of the FACE condition on the OFC to amygdala connection, we assume a context-independent reduction of prefrontal control over amygdalar activation in SAD patients. Using DCM, it was possible to highlight not only the neuronal dysfunction of isolated brain regions, but also the dysbalance of a distributed functional network.

Improved quantification of cerebral hemodynamics using individualized time thresholds for assessment of peak enhancement parameters derived from dynamic susceptibility contrast enhanced magnetic resonance imaging

Purpose

Assessment of cerebral ischemia often employs dynamic susceptibility contrast enhanced magnetic resonance imaging (DSC-MRI) with evaluation of various peak enhancement time parameters. All of these parameters use a single time threshold to judge the maximum tolerable peak enhancement delay that is supposed to reliably differentiate sufficient from critical perfusion. As the validity of this single threshold approach still remains unclear, in this study, (1) the definition of a threshold on an individual patient-basis, nevertheless (2) preserving the comparability of the data, was investigated.

Methods

The histogram of time-to-peak (TTP) values derived from DSC-MRI, the so-called TTP-distribution curve (TDC), was modeled using a double-Gaussian model in 61 patients without severe cerebrovascular disease. Particular model-based z_f-scores were used to describe the arterial, parenchymal and venous bolus-transit phase as time intervals I_a,_p,v. Their durations (delta I_a,p,v), were then considered as maximum TTP-delays of each phase.

Results

Mean-R² for the model-fit was 0.967. Based on the generic z_f-scores the proposed bolus transit phases could be differentiated. The I_p-interval reliably depicted the parenchymal bolus-transit phase with durations of 3.4 s–10.1 s (median = 4.3s), where an increase with age was noted (∼30 ms/year).

Conclusion

Individual threshold-adjustment seems rational since regular bolus-transit durations in brain parenchyma obtained from the TDC overlap considerably with recommended critical TTP-thresholds of 4 s–8 s. The parenchymal transit time derived from the proposed model may be utilized to individually correct TTP-thresholds, thereby potentially improving the detection of critical perfusion.

A novel coil array for combined TMS/fMRI experiments at 3 T

Purpose

To overcome current limitations in combined transcranial magnetic stimulation (TMS) and functional magnetic resonance imaging (fMRI) studies by employing a dedicated coil array design for 3 Tesla.

Methods

The state‐of‐the‐art setup for concurrent TMS/fMRI is to use a large birdcage head coil, with the TMS between the subject's head and the MR coil. This setup has drawbacks in sensitivity, positioning, and available imaging techniques. In this study, an ultraslim 7‐channel receive‐only coil array for 3 T, which can be placed between the subject's head and the TMS, is presented. Interactions between the devices are investigated and the performance of the new setup is evaluated in comparison to the state‐of‐the‐art setup.

Results

MR sensitivity obtained at the depth of the TMS stimulation is increased by a factor of five. Parallel imaging with an acceleration factor of two is feasible with low g‐factors. Possible interactions between TMS and the novel hardware were investigated and were found negligible.

Conclusion

The novel coil array is safe, strongly improves signal‐to‐noise ratio in concurrent TMS/fMRI experiments, enables parallel imaging, and allows for flexible positioning of the TMS on the head while ensuring efficient TMS stimulation due to its ultraslim design. Magn Reson Med 74:1492–1501, 2015. © 2014 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

On the generalizability of resting-state fMRI machine learning classifiers

Machine learning classifiers have become increasingly popular tools to generate single-subject inferences from fMRI data. With this transition from the traditional group level difference investigations to single-subject inference, the application of machine learning methods can be seen as a considerable step forward. Existing studies, however, have given scarce or no information on the generalizability to other subject samples, limiting the use of such published classifiers in other research projects. We conducted a simulation study using publicly available resting-state fMRI data from the 1000 Functional Connectomes and COBRE projects to examine the generalizability of classifiers based on regional homogeneity of resting-state time series. While classification accuracies of up to 0.8 (using sex as the target variable) could be achieved on test datasets drawn from the same study as the training dataset, the generalizability of classifiers to different study samples proved to be limited albeit above chance. This shows that on the one hand a certain amount of generalizability can robustly be expected, but on the other hand this generalizability should not be overestimated. Indeed, this study substantiates the need to include data from several sites in a study investigating machine learning classifiers with the aim of generalizability.