Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome.

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Two-dimensional XD-GRASP provides better image quality than conventional 2D cardiac cine MRI for patients who cannot suspend respiration.

OBJECTIVES: Residual respiratory motion degrades image quality in conventional cardiac cine MRI (CCMRI). We evaluated whether a free-breathing (FB) radial imaging CCMRI sequence with compressed sensing reconstruction [extradimensional (e.g. cardiac and respiratory phases) golden-angle radial sparse parallel, or XD-GRASP] could provide better image quality than a conventional Cartesian breath-held (BH) sequence in an unselected population of patients undergoing clinical CCMRI. MATERIALS AND METHODS: One hundred one patients who underwent BH and FB imaging in a midventricular short-axis plane at a matching location were included. Visual and quantitative image analysis was performed by two blinded experienced readers, using a five-point qualitative scale to score overall image quality and visual signal-to-noise ratio (SNR) grade, with measures of noise and sharpness. End-diastolic and end-systolic left ventricular areas were also measured and compared for both BH and FB images. RESULTS: Image quality was generally better with the BH cines (overall quality grade for BH vs FB images 4 vs 2.9, p < 0.001; noise 0.06 vs 0.08 p < 0.001; SNR grade 4.1 vs 3, p < 0.001), except for sharpness (p = 0.48). There were no significant differences between BH and FB images regarding end-diastolic or end-systolic areas (p = 0.35 and p = 0.12). Eighteen of the 101 patients had poor BH image quality (grade 1 or 2). In this subgroup, the quality of the FB images was better (p = 0.0032), as was the SNR grade (p = 0.003), but there were no significant differences regarding noise and sharpness (p = 0.45 and p = 0.47). CONCLUSION: Although FB XD-GRASP CCMRI was visually inferior to conventional BH CCMRI in general, it provided improved image quality in the subgroup of patients with respiratory-motion-induced artifacts on BH images.

High resolution CBV assessment with PEAK-EPI: k-t-undersampling and reconstruction in echo planar imaging.

PURPOSE: Achieving higher spatial resolution and improved brain coverage while mitigating in-plane susceptibility artifacts in the assessment of perfusion parameters, such as cerebral blood volume, in echo planar imaging (EPI)-based dynamic susceptibility contrast weighted cerebral perfusion measurements. METHODS: PEAK-EPI, an EPI sequence with interleaved readout trajectories and three different strategies for autocalibration-signal acquisition (inplace, dynamic extra and extra) is presented. Performance of each approach is analyzed in vivo based on flip angle variation induced dynamics, assessing temporal fidelity, temporal SNR and g-factors. All approaches are compared with conventional GRAPPA reconstructions. PEAK-EPI with inplace autocalibration-signal at R = 5 is then compared with the standard clinical EPI protocol in six patients, using two half-dose dynamic susceptibility contrast weighted cerebral perfusion measurements per subject. RESULTS: PEAK-EPI acquisition facilitates a substantial increase of spatial resolution at a higher number of slices per TR and provides improved SNR compared to conventional GRAPPA. High dependency of the resulting reconstruction quality on the type of autocalibration-signal acquisition is observed. PEAK-EPI with inplace autocalibration-signal achieves high temporal fidelity and initial feasibility is shown. CONCLUSION: The obtained high resolution cerebral blood volume maps reveal more detailed information than in corresponding standard EPI measurements and facilitate detailed delineation of tumorous tissue. Magn Reson Med 77:2153-2166, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Multi-cycle Reconstruction of Cardiac MRI for the Analysis of Inter-ventricular Septum Motion During Free Breathing.

Small variations in left-ventricular preload due to respiration produce measurable changes in cardiac function in normal subjects. We show that this mechanism is altered in patients with reduced ejection fraction (EF), hypertrophy, or volume-loaded right ventricle (RV). We propose a multi-dimensional retrospective image reconstruction, based on an adaptive, soft classification of data into respiratory and cardiac phases, to study these effects.

Use of self-gated radial cardiovascular magnetic resonance to detect and classify arrhythmias (atrial fibrillation and premature ventricular contraction).

BACKGROUND: Arrhythmia can significantly alter the image quality of cardiovascular magnetic resonance (CMR); automatic detection and sorting of the most frequent types of arrhythmias during the CMR acquisition could potentially improve image quality. New CMR techniques, such as non-Cartesian CMR, can allow self-gating: from cardiac motion-related signal changes, we can detect cardiac cycles without an electrocardiogram. We can further use this data to obtain a surrogate for RR intervals (valley intervals: VV). Our purpose was to evaluate the feasibility of an automated method for classification of non-arrhythmic (NA) (regular cycles) and arrhythmic patients (A) (irregular cycles), and for sorting of common arrhythmia patterns between atrial fibrillation (AF) and premature ventricular contraction (PVC), using the cardiac motion-related signal obtained during self-gated free-breathing radial cardiac cine CMR with compressed sensing reconstruction (XD-GRASP). METHODS: One hundred eleven patients underwent cardiac XD-GRASP CMR between October 2015 and February 2016; 33 were included for retrospective analysis with the proposed method (6 AF, 8 PVC, 19 NA; by recent ECG). We analyzed the VV, using pooled statistics (histograms) and sequential analysis (Poincaré plots), including the median (medVV), the weighted mean (meanVV), the total number of VV values (VVval), and the total range (VVTR) and half range (VVHR) of the cumulative frequency distribution of VV, including the median to half range (medVV/VVHR) and the half range to total range (VVHR/VVTR) ratios. We designed a simple algorithm for using the VV results to differentiate A from NA, and AF from PVC. RESULTS: Between NA and A, meanVV, VVval, VVTR, VVHR, medVV/VVHR and VVHR/VVTR ratios were significantly different (p values = 0.00014, 0.0027, 0.000028, 5×10-9, 0.002, respectively). Between AF and PVC, meanVV, VVval and medVV/VVHR ratio were significantly different (p values = 0.018, 0.007, 0.044, respectively). Using our algorithm, sensitivity, specificity, and accuracy were 93 %, 95 % and 94 % to discriminate between NA and A, and 83 %, 71 %, and 77 % to discriminate between AF and PVC, respectively; areas under the ROC curve were 0.93 and 0.89. CONCLUSIONS: Our study shows we can reliably detect arrhythmias and differentiate AF from PVC, using self-gated cardiac cine XD-GRASP CMR.

A g-factor metric for k-t SENSE and k-t PCA based parallel imaging.

PURPOSE: To propose and validate a g-factor formalism for k-t SENSE, k-t PCA and related k-t methods for assessing SNR and temporal fidelity. METHODS: An analytical gxf -factor formulation in the spatiotemporal frequency domain is derived, enabling assessment of noise and depiction fidelity in both the spatial and frequency domain. Using pseudoreplica analysis of cardiac cine data the gxf -factor description is validated and example data are used to analyze the performance of k-t methods for various parameter settings. RESULTS: Analytical gxf -factor maps were found to agree well with pseudoreplica analysis for 3x, 5x, and 7x k-t SENSE and k-t PCA. While k-t SENSE resulted in lower average gxf values (gx (avg) ) in static regions when compared with k-t PCA, k-t PCA yielded lower gx (avg) values in dynamic regions. Temporal transfer was better preserved with k-t PCA for increasing undersampling factors. CONCLUSION: The proposed gxf -factor and temporal transfer formalism allows assessing noise performance and temporal depiction fidelity of k-t methods including k-t SENSE and k-t PCA. The framework enables quantitative comparison of different k-t methods relative to frame-by-frame parallel imaging reconstruction.

A g-factor metric for k-t-GRAPPA- and PEAK-GRAPPA-based parallel imaging.

PURPOSE: The aim of this work is to derive a theoretical framework for quantitative noise and temporal fidelity analysis of time-resolved k-space-based parallel imaging methods. THEORY: An analytical formalism of noise distribution is derived extending the existing g-factor formulation for nontime-resolved generalized autocalibrating partially parallel acquisition (GRAPPA) to time-resolved k-space-based methods. The noise analysis considers temporal noise correlations and is further accompanied by a temporal filtering analysis. METHODS: All methods are derived and presented for k-t-GRAPPA and PEAK-GRAPPA. A sliding window reconstruction and nontime-resolved GRAPPA are taken as a reference. Statistical validation is based on series of pseudoreplica images. The analysis is demonstrated on a short-axis cardiac CINE dataset. RESULTS: The superior signal-to-noise performance of time-resolved over nontime-resolved parallel imaging methods at the expense of temporal frequency filtering is analytically confirmed. Further, different temporal frequency filter characteristics of k-t-GRAPPA, PEAK-GRAPPA, and sliding window are revealed. CONCLUSION: The proposed analysis of noise behavior and temporal fidelity establishes a theoretical basis for a quantitative evaluation of time-resolved reconstruction methods. Therefore, the presented theory allows for comparison between time-resolved parallel imaging methods and also nontime-resolved methods. Magn Reson Med 74:125-135, 2015. © 2014 Wiley Periodicals, Inc.

The impact of latent confounders in directed network analysis in neuroscience.

In the analysis of neuroscience data, the identification of task-related causal relationships between various areas of the brain gives insights about the network of physiological pathways that are active during the task. One increasingly used approach to identify causal connectivity uses the concept of Granger causality that exploits predictability of activity in one region by past activity in other regions of the brain. Owing to the complexity of the data, selecting components for the analysis of causality as a preprocessing step has to be performed. This includes predetermined-and often arbitrary-exclusion of information. Therefore, the system is confounded by latent sources. In this paper, the effect of latent confounders is demonstrated, and paths of influence among three components are studied. While methods for analysing Granger causality are commonly based on linear vector autoregressive models, the effects of latent confounders are expected to be present also in nonlinear systems. Therefore, all analyses are also performed for a simulated nonlinear system and discussed with regard to applications in neuroscience.