Simulation-based optimization and experimental comparison of intracranial T2-weighted DANTE-SPACE vessel wall imaging at 3T and 7T.

PURPOSE: T2-weighted DANTE-SPACE (Delay Alternating with Nutation for Tailored Excitation - Sampling Perfection with Application optimized Contrasts using different flip angle Evolution) sequences facilitate non-invasive intracranial vessel wall imaging at 7T through simultaneous suppression of blood and CSF. However, the achieved vessel wall delineation depends closely on the selected sequence parameters, and little information is available about the performance of the sequence using more widely available 3T MRI. Therefore, in this paper a comprehensive DANTE-SPACE simulation framework is used for the optimization and quantitative comparison of T2-weighted DANTE-SPACE at both 7T and 3T. METHODS: Simulations are used to propose optimized sequence parameters at both 3T and 7T. At 7T, an additional protocol which uses a parallel transmission (pTx) shim during the DANTE preparation for improved suppression of inflowing blood is also proposed. Data at both field strengths using optimized and literature protocols are acquired and quantitatively compared in six healthy volunteers. RESULTS: At 7T, more vessel wall signal can be retained while still achieving sufficient CSF suppression by using fewer DANTE pulses than described in previous implementations. The use of a pTx shim during DANTE at 7T provides a modest further improvement to the inner vessel wall delineation. At 3T, aggressive DANTE preparation is required to achieve CSF suppression, resulting in reduced vessel wall signal. As a result, the achievable vessel wall definition at 3T is around half that of 7T. CONCLUSION: Simulation-based optimization of DANTE parameters facilitates improved T2-weighted DANTE-SPACE contrasts at 7T. The improved vessel definition of T2-weighted DANTE-SPACE at 7T makes DANTE preparation more suitable for T2-weighted VWI at 7T than at 3T.

Accelerated 3D multi-channel B 1 + mapping at 7 T for the brain and heart.

PURPOSE: To acquire accurate volumetric multi-channel B 1 + $$ {\mathrm{B}}_1^{+} $$ maps in under 14 s whole-brain or 23 heartbeats whole-heart for parallel transmit (pTx) applications at 7 T. THEORY AND METHODS: We evaluate the combination of three recently proposed techniques. The acquisition of multi-channel transmit array B 1 + $$ {\mathrm{B}}_1^{+} $$ maps is accelerated using transmit low rank (TxLR) with absolute B 1 + $$ {\mathrm{B}}_1^{+} $$ mapping (Sandwich) acquired in a B 1 + $$ {\mathrm{B}}_1^{+} $$ time-interleaved acquisition of modes (B1TIAMO) fashion. Simulations using synthetic body images derived from Sim4Life were used to test the achievable acceleration for small scan matrices of 24 × 24. Next, we evaluated the method by retrospectively undersampling a fully sampled B 1 + $$ {\mathrm{B}}_1^{+} $$ library of nine subjects in the brain. Finally, Cartesian undersampled phantom and in vivo images were acquired in both the brain of three subjects (8Tx/32 receive [Rx]) and the heart of another three subjects (8Tx/8Rx) at 7 T. RESULTS: Simulation and in vivo results show that volumetric multi-channel B 1 + $$ {\mathrm{B}}_1^{+} $$ maps can be acquired using acceleration factors of 4 in the body, reducing the acquisition time to within 23 heartbeats, which was previously not possible. In silico heart simulations demonstrated a RMS error to the fully sampled native resolution ground truth of 4.2° when combined in first-order circularly polarized mode (mean flip angle 66°) at an acceleration factor of 4. The 14 s 3D B 1 + $$ {\mathrm{B}}_1^{+} $$ maps acquired in the brain have a RMS error of 1.9° to the fully sampled (mean flip angle 86°). CONCLUSION: The proposed method is demonstrated as a fast pTx calibration technique in the brain and a promising method for pTx calibration in the body.

An extended phase graph-based framework for DANTE-SPACE simulations including physiological, temporal, and spatial variations.

PURPOSE: The delay alternating with nutation for tailored excitation (DANTE)-sampling perfection with application-optimized contrasts (SPACE) sequence facilitates 3D intracranial vessel wall imaging with simultaneous suppression of blood and CSF. However, the achieved image contrast depends closely on the selected sequence parameters, and the clinical use of the sequence is limited in vivo by observed signal variations in the vessel wall, CSF, and blood. This paper introduces a comprehensive DANTE-SPACE simulation framework, with the aim of providing a better understanding of the underlying contrast mechanisms and facilitating improved parameter selection and contrast optimization. METHODS: An extended phase graph formalism was developed for efficient spin ensemble simulation of the DANTE-SPACE sequence. Physiological processes such as pulsatile flow velocity variation, varying flow directions, intravoxel velocity variation, diffusion, and B 1 + $$ {\mathrm{B}}_1^{+} $$ effects were included in the framework to represent the mechanisms behind the achieved signal levels accurately. RESULTS: Intravoxel velocity variation improved temporal stability and robustness against small velocity changes. Time-varying pulsatile velocity variation affected CSF simulations, introducing periods of near-zero velocity and partial rephasing. Inclusion of diffusion effects was found to substantially reduce the CSF signal. Blood flow trajectory variations had minor effects, but B 1 + $$ {\mathrm{B}}_1^{+} $$ differences along the trajectory reduced DANTE efficiency in low- B 1 + $$ {\mathrm{B}}_1^{+} $$ areas. Introducing low-velocity pulsatility of both CSF and vessel wall helped explain the in vivo observed signal heterogeneity in both tissue types. CONCLUSION: The presented simulation framework facilitates a more comprehensive optimization of DANTE-SPACE sequence parameters. Furthermore, the simulation framework helps to explain observed contrasts in acquired data.

Head-and-neck multichannel B1+ mapping and RF shimming of the carotid arteries using a 7T parallel-transmit head coil.

PURPOSE: Neurovascular MRI suffers from a rapid drop in B1 + into the neck when using transmit head coils at 7 T. One solution to improving B1 + magnitude in the major feeding arteries in the neck is to use custom RF shims on parallel-transmit head coils. However, calculating such shims requires robust multichannel B1 + maps in both the head and the neck, which is challenging due to low RF penetration into the neck, limited dynamic range of multichannel B1 + mapping techniques, and B0 sensitivity. We therefore sought a robust, large-dynamic-range, parallel-transmit field mapping protocol and tested whether RF shimming can improve carotid artery B1 + magnitude in practice. METHODS: A pipeline is presented that combines B1 + mapping data acquired using circularly polarized (CP) and CP2-mode RF shims at multiple voltages. The pipeline was evaluated by comparing the predicted and measured B1 + for multiple random transmit shims, and by assessing the ability of RF shimming to increase B1 + in the carotid arteries. RESULTS: The proposed method achieved good agreement between predicted and measured B1 + in both the head and the neck. The B1 + magnitude in the carotid arteries can be increased by 43% using tailored RF shims or by 37% using universal RF shims, while also improving the RF homogeneity compared with CP mode. CONCLUSION: B1 + in the neck can be increased using RF shims calculated from multichannel B1 + maps in both the head and the neck. This can be achieved using universal phase-only RF shims, facilitating easy implementation in existing sequences.

Reduced Left Atrial Rotational Flow Is Independently Associated With Embolic Brain Infarcts.

BACKGROUND: Up to 25% of embolic strokes occur in individuals without atrial fibrillation (AF) or other identifiable mechanisms. OBJECTIVES: This study aims to assess whether left atrial (LA) blood flow characteristics are associated with embolic brain infarcts, independently of AF. METHODS: The authors recruited 134 patients: 44 with a history of ischemic stroke and 90 with no history of stroke but CHA2DS2VASc score ≥1. Cardiac magnetic resonance (CMR) evaluated cardiac function and LA 4-dimensional flow parameters, including velocity and vorticity (a measure of rotational flow), and brain magnetic resonance imaging (MRI) was performed to detect large noncortical or cortical infarcts (LNCCIs) (likely embolic), or nonembolic lacunar infarcts. RESULTS: Patients (41% female; age 70 ± 9 years) had moderate stroke risk (median CHA2DS2VASc = 3, Q1-Q3: 2-4). Sixty-eight (51%) had diagnosed AF, of whom 58 (43%) were in AF during CMR. Thirty-nine (29%) had ≥1 LNCCI, 20 (15%) had ≥1 lacunar infarct without LNCCI, and 75 (56%) had no infarct. Lower LA vorticity was significantly associated with prevalent LNCCIs after adjustment for AF during CMR, history of AF, CHA2DS2VASc score, LA emptying fraction, LA indexed maximum volume, left ventricular ejection fraction, and indexed left ventricular mass (OR: 2.06 [95% CI: 1.08-3.92 per SD]; P = 0.027). By contrast, LA flow peak velocity was not significantly associated with LNCCIs (P = 0.21). No LA parameter was associated with lacunar infarcts (all P > 0.05). CONCLUSIONS: Reduced LA flow vorticity is significantly and independently associated with embolic brain infarcts. Imaging LA flow characteristics may aid identification of individuals who would benefit from anticoagulation for embolic stroke prevention, regardless of heart rhythm.

Rapid 3D absolute B1 + mapping using a sandwiched train presaturated TurboFLASH sequence at 7 T for the brain and heart.

PURPOSE: To shorten the acquisition time of magnetization-prepared absolute transmit field (B1 + ) mapping known as presaturation TurboFLASH, or satTFL, to enable single breath-hold whole-heart 3D B1 + mapping. METHODS: SatTFL is modified to remove the delay between the reference and prepared images (typically 5 T1 ), with matching transmit configurations for excitation and preparation RF pulses. The new method, called Sandwich, is evaluated as a 3D sequence, measuring whole-brain and gated whole-heart B1 + maps in a single breath-hold. We evaluate the sensitivity to B1 + and T1 using numerical Bloch, extended phase graph, and Monte Carlo simulations. Phantom and in vivo images were acquired in both the brain and heart using an 8-channel transmit 7 Tesla MRI system to support the simulations. A segmented satTFL with a short readout train was used as a reference. RESULTS: The method significantly reduces acquisition times of 3D measurements from 360 s to 20 s, in the brain, while simultaneously reducing bias in the measured B1 + due to T1 and magnetization history. The mean coefficient of variation was reduced by 81% for T1 s of 0.5-3 s compared to conventional satTFL. In vivo, the reproducibility coefficient for flip angles in the range 0-130° was 4.5° for satTFL and 4.7° for our scheme, significantly smaller than for a short TR satTFL sequence, which was 12°. The 3D sequence measured B1 + maps of the whole thorax in 26 heartbeats. CONCLUSION: Our adaptations enable faster B1 + mapping, with minimal T1 sensitivity and lower sensitivity to magnetization history, enabling single breath-hold whole-heart absolute B1 + mapping.

External Hardware and Sensors, for Improved MRI.

Complex engineered systems are often equipped with suites of sensors and ancillary devices that monitor their performance and maintenance needs. MRI scanners are no different in this regard. Some of the ancillary devices available to support MRI equipment, the ones of particular interest here, have the distinction of actually participating in the image acquisition process itself. Most commonly, such devices are used to monitor physiological motion or variations in the scanner's imaging fields, allowing the imaging and/or reconstruction process to adapt as imaging conditions change. "Classic" examples include electrocardiography (ECG) leads and respiratory bellows to monitor cardiac and respiratory motion, which have been standard equipment in scan rooms since the early days of MRI. Since then, many additional sensors and devices have been proposed to support MRI acquisitions. The main physical properties that they measure may be primarily "mechanical" (eg acceleration, speed, and torque), "acoustic" (sound and ultrasound), "optical" (light and infrared), or "electromagnetic" in nature. A review of these ancillary devices, as currently available in clinical and research settings, is presented here. In our opinion, these devices are not in competition with each other: as long as they provide useful and unique information, do not interfere with each other and are not prohibitively cumbersome to use, they might find their proper place in future suites of sensors. In time, MRI acquisitions will likely include a plurality of complementary signals. A little like the microbiome that provides genetic diversity to organisms, these devices can provide signal diversity to MRI acquisitions and enrich measurements. Machine-learning (ML) algorithms are well suited at combining diverse input signals toward coherent outputs, and they could make use of all such information toward improved MRI capabilities. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Optimization of undersampling parameters for 3D intracranial compressed sensing MR angiography at 7 T.

PURPOSE: 3D time-of-flight MRA can accurately visualize the intracranial vasculature but is limited by long acquisition times. Compressed sensing reconstruction can be used to substantially accelerate acquisitions. The quality of those reconstructions depends on the undersampling patterns used. In this work, we optimize sets of undersampling parameters for various acceleration factors of Cartesian 3D time-of-flight MRA. METHODS: Fully sampled datasets, acquired at 7 Tesla, were retrospectively undersampled using variable-density Poisson disk sampling with various autocalibration region sizes, polynomial orders, and acceleration factors. The accuracy of reconstructions from the different undersampled datasets was assessed using the vessel-masked structural similarity index. Identified optimal undersampling parameters were then evaluated in additional prospectively undersampled datasets. Compressed sensing reconstruction parameters were chosen based on a preliminary reconstruction parameter optimization. RESULTS: For all acceleration factors, using a fully sampled calibration area of 12 × 12 k-space lines and a polynomial order of 2 resulted in the highest image quality. The importance of parameter optimization of the sampling was found to increase for higher acceleration factors. The results were consistent across resolutions and regions of interest with vessels of varying sizes and tortuosity. The number of visible small vessels increased by 7.0% and 14.2% when compared to standard parameters for acceleration factors of 7.2 and 15, respectively. CONCLUSION: The image quality of compressed sensing time-of-flight MRA can be improved by appropriate choice of undersampling parameters. The optimized sets of parameters are independent of the acceleration factor and enable a larger number of vessels to be visualized.

The impact of atrial fibrillation and stroke risk factors on left atrial blood flow characteristics.

AIMS: Altered left atrial (LA) blood flow characteristics account for an increase in cardioembolic stroke risk in atrial fibrillation (AF). Here, we aimed to assess whether exposure to stroke risk factors is sufficient to alter LA blood flow even in the presence of sinus rhythm (SR). METHODS AND RESULTS: We investigated 95 individuals: 37 patients with persistent AF, who were studied before and after cardioversion [Group 1; median CHA2DS2-VASc = 2.0 (1.5-3.5)]; 35 individuals with no history of AF but similar stroke risk to Group 1 [Group 2; median CHA2DS2-VASc = 3.0 (2.0-4.0)]; and 23 low-risk individuals in SR [Group 3; median CHA2DS2-VASc = 0.0 (0.0-0.0)]. Cardiac function and LA flow characteristics were evaluated using cardiac magnetic resonance. Before cardioversion, Group 1 displayed impaired left ventricular (LV) and LA function, reduced LA flow velocities and vorticity, and a higher normalized vortex volume (all P < 0.001 vs. Groups 2 and 3). After restoration of SR at ≥4-week post-cardioversion, LV systolic function and LA flow parameters improved significantly (all P < 0.001 vs. pre-cardioversion) and were no longer different from those in Group 2. However, in the presence of SR, LA flow peak and mean velocity, and vorticity were lower in Groups 1 and 2 vs. Group 3 (all P < 0.01), and were associated with impaired LA emptying fraction (LAEF) and LV diastolic dysfunction. CONCLUSION: Patients at moderate-to-high stroke risk display altered LA flow characteristics in SR in association with an LA myopathic phenotype and LV diastolic dysfunction, regardless of a history of AF.

Accelerated calibrationless parallel transmit mapping using joint transmit and receive low-rank tensor completion.

PURPOSE: To evaluate an algorithm for calibrationless parallel imaging to reconstruct undersampled parallel transmit field maps for the body and brain. METHODS: Using a combination of synthetic data and in vivo measurements from brain and body, 3 different approaches to a joint transmit and receive low-rank tensor completion algorithm are evaluated. These methods included: 1) virtual coils using the product of receive and transmit sensitivities, 2) joint-receiver coils that enforces a low rank structure across receive coils of all transmit modes, and 3) transmit low rank that uses a low rank structure for both receive and transmit modes simultaneously. The performance of each is investigated for different noise levels and different acceleration rates on an 8-channel parallel transmit 7 Tesla system. RESULTS: The virtual coils method broke down with increasing noise levels or acceleration rates greater than 2, producing normalized RMS error greater than 0.1. The joint receiver coils method worked well up to acceleration factors of 4, beyond which the normalized RMS error exceeded 0.1. Transmit low rank enabled an eightfold acceleration, with most normalized RMS errors remaining below 0.1. CONCLUSION: This work demonstrates that undersampling factors of up to eightfold are feasible for transmit array mapping and can be reconstructed using calibrationless parallel imaging methods.

Left atrial 4D flow cardiovascular magnetic resonance: a reproducibility study in sinus rhythm and atrial fibrillation.

BACKGROUND: Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) allows sophisticated quantification of left atrial (LA) blood flow, and could yield novel biomarkers of propensity for intra-cardiac thrombus formation and embolic stroke. As reproducibility is critically important to diagnostic performance, we systematically investigated technical and temporal variation of LA 4D flow in atrial fibrillation (AF) and sinus rhythm (SR). METHODS: Eighty-six subjects (SR, n = 64; AF, n = 22) with wide-ranging stroke risk (CHA2DS2VASc 0-6) underwent LA 4D flow assessment of peak and mean velocity, vorticity, vortex volume, and stasis. Eighty-five (99%) underwent a second acquisition within the same session, and 74 (86%) also returned at 30 (27-35) days for an interval scan. We assessed variability attributable to manual contouring (intra- and inter-observer), and subject repositioning and reacquisition of data, both within the same session (same-day scan-rescan), and over time (interval scan). Within-subject coefficients of variation (CV) and bootstrapped 95% CIs were calculated and compared. RESULTS: Same-day scan-rescan CVs were 6% for peak velocity, 5% for mean velocity, 7% for vorticity, 9% for vortex volume, and 10% for stasis, and were similar between SR and AF subjects (all p > 0.05). Interval-scan variability was similar to same-day scan-rescan variability for peak velocity, vorticity, and vortex volume (all p > 0.05), and higher for stasis and mean velocity (interval scan CVs of 14% and 8%, respectively, both p < 0.05). Longitudinal changes in heart rate and blood pressure at the interval scan in the same subjects were associated with significantly higher variability for LA stasis (p = 0.024), but not for the remaining flow parameters (all p > 0.05). SR subjects showed significantly greater interval-scan variability than AF patients for mean velocity, vortex volume, and stasis (all p < 0.05), but not peak velocity or vorticity (both p > 0.05). CONCLUSIONS: LA peak velocity and vorticity are the most reproducible and temporally stable novel LA 4D flow biomarkers, and are robust to changes in heart rate, blood pressure, and differences in heart rhythm.

Assessment of radio-frequency heating of a parallel transmit coil in a phantom using multi-echo proton resonance frequency shift thermometry.

We propose a workflow for validating parallel transmission (pTx) radio-frequency (RF) magnetic field heating patterns using Proton-Resonance Frequency shift (PRF)-based MR thermometry. Electromagnetic (EM) and thermal simulations of a 7 T 8-channel dipole coil were done using commercially available software (Sim4Life) to assess RF heating. The fabrication method for a phantom with electrical properties matched to human tissue is also described, along with methods for its electrical and thermal characterisation. Energy was deposited to specific transmit channels, whilst acquiring 3D PRF data using a pair of interleaved RF shim transmit modes. A multi-echo readout and pre-scan stabilisation protocol were used for increased sensitivity and to correct for measurement-to-measurement instabilities. The electrical properties of the phantom were found to be within 10% of the intended values. Adoption of a 14-min stabilisation scan gave sufficient suppression of any evolving background spatial variation in the B0 field to achieve <0.001 °C/mm thermometry drift over 10 min of subsequent scanning. Using two RF shim transmit modes enabled full phantom coverage and combining multiple echo times enabled a 13-54% improvement in the RMSE sensitivity to temperature changes. Combining multiple echoes reduced the peak RMSE by 45% and visually reduced measurement-to-measurement instabilities. A reference fibre optic probe showed temperature deviations from the PRF-estimated temperature to be smaller than 0.5 °C. Given the importance of RF safety in pTx applications, this workflow enables accurate validation of RF heating simulations with minimal additional hardware requirements.

An investigation into the minimum number of tissue groups required for 7T in-silico parallel transmit electromagnetic safety simulations in the human head.

PURPOSE: Safety limits for the permitted specific absorption rate (SAR) place restrictions on pulse sequence design, especially at ultrahigh fields (≥ 7 tesla). Due to intersubject variability, the SAR is usually conservatively estimated based on standard human models that include an applied safety margin to ensure safe operation. One approach to reducing the restrictions is to create more accurate subject-specific models from their segmented MR images. This study uses electromagnetic simulations to investigate the minimum number of tissue groups required to accurately determine SAR in the human head. METHODS: Tissue types from a fully characterized electromagnetic human model with 47 tissue types in the head and neck region were grouped into different tissue clusters based on the conductivities, permittivities, and mass densities of the tissues. Electromagnetic simulations of the head model inside a parallel transmit head coil at 7 tesla were used to determine the minimum number of required tissue clusters to accurately determine the subject-specific SAR. The identified tissue clusters were then evaluated using 2 additional well-characterized electromagnetic human models. RESULTS: A minimum of 4-clusters-plus-air was found to be required for accurate SAR estimation. These tissue clusters are centered around gray matter, fat, cortical bone, and cerebrospinal fluid. For all 3 simulated models, the parallel transmit maximum 10g SAR was consistently determined to within an error of <12% relative to the full 47-tissue model. CONCLUSION: A minimum of 4-clusters-plus-air are required to produce accurate personalized SAR simulations of the human head when using parallel transmit at 7 tesla.

Motion correction methods for MRS: experts' consensus recommendations.

Long acquisition times due to intrinsically low signal-to-noise ratio and the need for highly homogeneous B0 field make MRS particularly susceptible to motion or scanner instability compared with MRI. Motion-induced changes in both localization and shimming (ie B0 homogeneity) degrade MRS data quality. To mitigate the effects of motion three approaches can be employed: (1) subject immobilization, (2) retrospective correction, and (3) prospective real-time correction using internal and/or external tracking methods. Prospective real-time correction methods can simultaneously update localization and the B0 field to improve MRS data quality. While localization errors can be corrected with both internal (navigators) and external (optical camera, NMR probes) tracking methods, the B0 field correction requires internal navigator methods to measure the B0 field inside the imaged volume and the possibility to update the scanner shim hardware in real time. Internal and external tracking can rapidly update the MRS localization with submillimeter and subdegree precision, while scanner frequency and first-order shims of scanner hardware can be updated by internal methods every sequence repetition. These approaches are most well developed for neuroimaging, for which rigid transformation is primarily applicable. Real-time correction greatly improves the stability of MRS acquisition and quantification, as shown in clinical studies on subjects prone to motion, including children and patients with movement disorders, enabling robust measurement of metabolite signals including those with low concentrations, such as gamma-aminobutyric acid and glutathione. Thus, motion correction is recommended for MRS users and calls for tighter integration and wider availability of such methods by MR scanner manufacturers.

Navigator-based reacquisition and estimation of motion-corrupted data: Application to multi-echo spin echo for carotid wall MRI.

PURPOSE: To assess whether artifacts in multi-slice multi-echo spin echo neck imaging, thought to be caused by brief motion events such as swallowing, can be corrected by reacquiring corrupted central k-space data and estimating the remainder with parallel imaging. METHODS: A single phase-encode line (ky = 0, phase-encode direction anteroposterior) navigator echo was used to identify motion-corrupted data and guide the online reacquisition. If motion corruption was detected in the 7 central k-space lines, they were replaced with reacquired data. Subsequently, GRAPPA reconstruction was trained on the updated central portion of k-space and then used to estimate the remaining motion-corrupted k-space data from surrounding uncorrupted data. Similar compressed sensing-based approaches have been used previously to compensate for respiration in cardiac imaging. The g-factor noise amplification was calculated for the parallel imaging reconstruction of data acquired with a 10-channel neck coil. The method was assessed in scans with 9 volunteers and 12 patients. RESULTS: The g-factor analysis showed that GRAPPA reconstruction of 2 adjacent motion-corrupted lines causes high noise amplification; therefore, the number of 2-line estimations should be limited. In volunteer scans, median ghosting reduction of 24% was achieved with 2 adjacent motion-corrupted lines correction, and image quality was improved in 2 patient scans that had motion corruption close to the center of k-space. CONCLUSION: Motion-corrupted echo-trains can be identified with a navigator echo. Combined reacquisition and parallel imaging estimation reduced motion artifacts in multi-slice MESE when there were brief motion events, especially when motion corruption was close to the center of k-space.

Scattering matrix imaging pulse design for real-time respiration and cardiac motion monitoring.

PURPOSE: The scattering matrix (S-matrix) of a parallel transmit (pTx) coil is sensitive to physiological motion but requires additional monitoring RF pulses to be measured. In this work, we present and evaluate pTx RF pulse designs that simultaneously excite for imaging and measure the S-matrix to generate real-time motion signals without prolonging the image sequence. THEORY AND METHODS: Three pTx waveforms for measuring the S-matrix were identified and superimposed onto the imaging excitation RF pulses: (1) time division multiplexing, (2) frequency division multiplexing, and (3) code division multiplexing. These 3 methods were evaluated in healthy volunteers for scattering sensitivity and image artefacts. The S-matrix and real-time motion signals were calculated on the image calculation environment of the MR scanner. Prospective cardiac triggers were identified in early systole as a high rate of change of the cardiac motion signal. Monitoring accuracy was compared against electrocardiogram or the imaged diaphragm position. RESULTS: All 3 monitoring approaches measure the S-matrix during image excitation with quality correlated to input power. No image artefacts were observed for frequency multiplexing, and low energy artefacts were observed in the other methods. The accuracy of the achieved prospective cardiac gating was 15 ± 16 ms for breath hold and 24 ± 17 ms during free breathing. The diaphragm position prediction accuracy was 1.3 ± 0.9 mm. In all volunteers, good quality cine images were acquired for breath hold scans and dual gated CINEs were demonstrated. CONCLUSION: The S-matrix can be measured during image excitation to generate real-time cardiac and respiratory motion signals for prospective gating. No artefacts are introduced when frequency division multiplexing is used.

Left Ventricular Flow Analysis.

BACKGROUND: Cardiac remodeling, after a myocardial insult, often causes progression to heart failure. The relationship between alterations in left ventricular blood flow, including kinetic energy (KE), and remodeling is uncertain. We hypothesized that increasing derangements in left ventricular blood flow would relate to (1) conventional cardiac remodeling markers, (2) increased levels of biochemical remodeling markers, (3) altered cardiac energetics, and (4) worsening patient symptoms and functional capacity. Methods Thirty-four dilated cardiomyopathy patients, 30 ischemic cardiomyopathy patients, and 36 controls underwent magnetic resonance including 4-dimensional flow, BNP (brain-type natriuretic peptide) measurement, functional capacity assessment (6-minute walk test), and symptom quantification. A subgroup of dilated cardiomyopathy and control subjects underwent cardiac energetic assessment. Left ventricular flow was separated into 4 components: direct flow, retained inflow, delayed ejection flow, and residual volume. Average KE throughout the cardiac cycle was calculated. RESULTS: Patients had reduced direct flow proportion and direct-flow average KE compared with controls ( P<0.0001). The residual volume proportion and residual volume average KE were increased in patients ( P<0.0001). Importantly, in a multiple linear regression model to predict the patient's 6-minute walk test, the independent predictors were age (ß=-0.3015; P=0.019) and direct-flow average KE (ß=0.280, P=0.035; R2 model, 0.466, P=0.002). In contrast, neither ejection fraction nor left ventricular volumes were independently predictive. CONCLUSIONS: This study demonstrates an independent predictive relationship between the direct-flow average KE and a prognostic measure of functional capacity. Intracardiac 4-dimensional flow parameters are novel biomarkers in heart failure and may provide additive value in monitoring new therapies and predicting prognosis.

Test-retest variability of left ventricular 4D flow cardiovascular magnetic resonance measurements in healthy subjects.

BACKGROUND: Quantification and visualisation of left ventricular (LV) blood flow is afforded by three-dimensional, time resolved phase contrast cardiovascular magnetic resonance (CMR 4D flow). However, few data exist upon the repeatability and variability of these parameters in a healthy population. We aimed to assess the repeatability and variability over time of LV 4D CMR flow measurements. METHODS: Forty five controls underwent CMR 4D flow data acquisition. Of these, 10 underwent a second scan within the same visit (scan-rescan), 25 returned for a second visit (interval scan; median interval 52 days, IQR 28-57 days). The LV-end diastolic volume (EDV) was divided into four flow components: 1) Direct flow: inflow that passes directly to ejection; 2) Retained inflow: inflow that enters and resides within the LV; 3) Delayed ejection flow: starts within the LV and is ejected and 4) Residual volume: blood that resides within the LV for > 2 cardiac cycles. Each flow components' volume was related to the EDV (volume-ratio). The kinetic energy at end-diastole (ED) was measured and divided by the components' volume. RESULTS: The dominant flow component in all 45 controls was the direct flow (volume ratio 38 ± 4%) followed by the residual volume (30 ± 4%), then delayed ejection flow (16 ± 3%) and retained inflow (16 ± 4%). The kinetic energy at ED for each component was direct flow (7.8 ± 3.0 microJ/ml), retained inflow (4.1 ± 2.0 microJ/ml), delayed ejection flow (6.3 ± 2.3 microJ/ml) and the residual volume (1.2 ± 0.5 microJ/ml). The coefficients of variation for the scan-rescan ranged from 2.5%-9.2% for the flow components' volume ratio and between 13.5%-17.7% for the kinetic energy. The interval scan results showed higher coefficients of variation with values from 6.2-16.1% for the flow components' volume ratio and 16.9-29.0% for the kinetic energy of the flow components. CONCLUSION: LV flow components' volume and their associated kinetic energy values are repeatable and stable within a population over time. However, the variability of these measurements in individuals over time is greater than can be attributed to sources of error in the data acquisition and analysis, suggesting that additional physiological factors may influence LV flow measurements.

Differential flow improvements after valve replacements in bicuspid aortic valve disease: a cardiovascular magnetic resonance assessment.

BACKGROUND: Abnormal aortic flow patterns in bicuspid aortic valve disease (BAV) may be partly responsible for the associated aortic dilation. Aortic valve replacement (AVR) may normalize flow patterns and potentially slow the concomitant aortic dilation. We therefore sought to examine differences in flow patterns post AVR. METHODS: Ninety participants underwent 4D flow cardiovascular magnetic resonance: 30 BAV patients with prior AVR (11 mechanical, 10 bioprosthetic, 9 Ross procedure), 30 BAV patients with a native aortic valve and 30 healthy subjects. RESULTS: The majority of subjects with mechanical AVR or Ross showed normal flow pattern (73% and 67% respectively) with near normal rotational flow values (7.2 ± 3.9 and 10.6 ± 10.5 mm2/ms respectively vs 3.8 ± 3.1 mm2/s for healthy subjects; both p > 0.05); and reduced in-plane wall shear stress (0.19 ± 0.13 N/m2 for mechanical AVR vs. 0.40 ± 0.28 N/m2 for native BAV, p < 0.05). In contrast, all subjects with a bioprosthetic AVR had abnormal flow patterns (mainly marked right-handed helical flow), with comparable rotational flow values to native BAV (20.7 ± 8.8 mm2/ms and 26.6 ± 16.6 mm2/ms respectively, p > 0.05), and a similar pattern for wall shear stress. Data before and after AVR (n = 16) supported these findings: mechanical AVR showed a significant reduction in rotational flow (30.4 ± 16.3 → 7.3 ± 4.1 mm2/ms; p < 0.05) and in-plane wall shear stress (0.47 ± 0.20 → 0.20 ± 0.13 N/m2; p < 0.05), whereas these parameters remained similar in the bioprosthetic AVR group. CONCLUSIONS: Abnormal flow patterns in BAV disease tend to normalize after mechanical AVR or Ross procedure, in contrast to the remnant abnormal flow pattern after bioprosthetic AVR. This may in part explain different aortic growth rates post AVR in BAV observed in the literature, but requires confirmation in a prospective study.