Abstract
Cardiovascular magnetic resonance imaging (CMR) has become an indispensable clinical tool for the assessment of morphology, function and structure of the heart muscle. By exploiting quantification of the effective transverse relaxation time () CMR also affords myocardial tissue characterization and probing of cardiac physiology, both being in the focus of ongoing research. These developments are fueled by the move to ultrahigh magnetic field strengths, which permits enhanced sensitivity and spatial resolution that help to overcome limitations of current clinical MR systems with the goal to contribute to a better understanding of myocardial (patho)physiology in vivo. In this context, the aim of this report is to introduce myocardial mapping at ultrahigh magnetic fields as a promising technique to non-invasively assess myocardial (patho)physiology. For this purpose the basic principles of assessment, the biophysical mechanisms determining and (pre)clinical applications of myocardial mapping are presented. Technological challenges and solutions for sensitized CMR at ultrahigh magnetic field strengths are discussed followed by a review of acquisition techniques and post-processing approaches. Preliminary results derived from myocardial mapping in healthy subjects and cardiac patients at 7.0 T are presented. A concluding section discusses remaining questions and challenges and provides an outlook on future developments and potential clinical applications.
Introduction
sensitized cardiovascular magnetic resonance
Myocardial tissue characterization plays an important role in the diagnosis and treatment of cardiac diseases. Thanks to its soft tissue contrast and versatility, cardiovascular magnetic resonance imaging (CMR) has become a vital clinical tool for diagnosis and for guiding therapy of cardiac diseases [1, 2]. CMR can provide morphologic and functional information as well as insights into microstructural changes of the heart muscle [2]. Quantitative mapping of MR relaxation times which govern the MR signal evolution offers the potential of non-invasive myocardial tissue characterization without the need of exogenous contrast agents. Mapping of the effective transverse relaxation time is the subject of intense clinical interest in CMR. By exploiting the blood oxygenation level-dependent (BOLD) effect [3], sensitized CMR has been proposed as a means of assessing myocardial tissue oxygenation and perfusion. mapping has been shown to be capable of detecting myocardial ischemia caused by coronary artery stenosis [4], to reveal myocardial perfusion deficits under pharmacological stress [5–10], to study endothelial function [11] or to assess breathing maneuver-dependent oxygenation changes in the myocardium [12–15]. Preclinical studies have also demonstrated the potential of mapping to detect structural changes in the infarcted heart muscle and even to distinguish between focal and diffuse fibrosis [16–18]. In clinical application mapping is the method of choice for quantification of myocardial iron content, an essential parameter for guiding therapy in patients with myocardial iron overload [19–23].
The linear increase of susceptibility effects with magnetic field strength together with the availability of ultrahigh field (B0 ≥ 7.0 T) whole body human MR systems has fueled explorations into myocardial mapping at 7.0 T. In this context, the aim of this report is to introduce the biophysical background of as a promising MR biomarker, present challenges and technical solutions for myocardial assessment at ultrahigh magnetic field strengths, discuss its merits and current limitations, as well as to show early applications in healthy volunteers and cardiac patients along with providing a look beyond the horizon.
Biophysics of the effective transverse relaxation time
The fundamental principle behind relaxation is the loss of phase coherence of an ensemble of spins contained within a volume of interest or voxel after a radio frequency (RF) excitation. Unlike T1 relaxation which is based on spin-lattice interactions or T2 relaxation which is caused by spin-spin interactions both being inherent properties of tissues in a magnetic field, relaxation includes a tissue inherent part as well as contributions from external magnetic field perturbations [24]. These magnetic field inhomogeneities influence the effective transversal MR relaxation time [25, 26]. is defined as:
with T2 being the transverse relaxation time and representing magnetic susceptibility related contributions [27].
The most common way of -weighted imaging is gradient recalled echo (GRE) imaging. The MR signal magnitude Sm(θ) created by a spoiled GRE pulse sequence is:
with S0 representing the spin density, TR the repetition time, TE the echo time defined by the time between MR signal excitation and MR signal readout [28], T1 and are tissue specific longitudinal and effective transversal relaxation time constants and θ is the flip angle about which the magnetization is deflected by the excitation RF pulse. If TR and T1 are being kept constant Equation (2) can be simplified to:
Exploiting this relationship can be estimated by acquiring a series of images at different echo times TE followed by an exponential fit of the measured signal intensity vs. the echo time TE. This is commonly realized by using multi echo gradient echo (MEGRE) pulse sequences, which employ a series of dephasing and refocusing gradients to quickly acquire a series of sensitized images at several echo times as illustrated in Figure 1. weighted MRI is most sensitive to field perturbations when TE is equal to [29]. Exponential fitting of the signal decay can be done either for each voxel individually or for the mean signal within a region of interest. Single voxel fitting is more prone to noise but provides spatially resolved information in the form of relaxation maps (Figure 1). Besides mono-exponential fitting also multi-exponential fitting can be applied, if multiple signal compartments with different relaxation times are expected within an imaging voxel.
Here S1, S2, Sn represent the relative volume fractions of the different compartments with their corresponding effective transverse relaxation times 1, 2, and n.
Figure 1
relaxation is blood oxygenation level dependent and provides a functional MR contrast which serves as the basis of functional brain mapping [3, 26]. The effect results from a change of the magnetic susceptibility of hemoglobin (Hb) depending on its oxygenation state. Oxygenated hemoglobin is diamagnetic and has minor effect on magnetic field homogeneity. Deoxygenated hemoglobin in contrast is paramagnetic and causes magnetic field perturbations on a microscopic level resulting in spin dephasing and signal loss. -weighted MRI is sensitive to changes in the amount of deoxygenated Hb (deoxy Hb) per tissue volume element (voxel). changes and corresponding signal attenuation in -weighted MR images can hence result from a change in hemoglobin oxygenation or a change of the tissue blood volume fraction. The discovery of the BOLD phenomenon led to the development of functional MRI for mapping of human brain function, but also inspired research into BOLD imaging and mapping of the heart [9, 30].
sensitized imaging and mapping are widely assumed to provide a surrogate of oxygenation. Yet the factors impacting the transverse relaxation rate other than oxygenation are numerous including macroscopic magnetic field inhomogeneities, blood volume fraction and hematocrit [31]. Considering a biologic tissue with a specific blood volume fraction BVf, a hematocrit Hct, and a local blood oxygen saturation So2, can be modeled as:
with γ|ΔBother| representing additional magnetic field inhomogeneities such as macroscopic field changes [32, 33] and Δχ0 = 3.318 ppm being the magnetic susceptibility difference of fully oxygenated and fully deoxygenated hemoglobin (in SI units) [34]. When tissue blood volume fraction, hematocrit level and macroscopic B0 contributions are known and echo times are greater than a characteristic time Equation (5) can be employed to non-invasively estimate tissue oxygenation using MRI [33]. It should be noted that a reduction in the tissue blood volume fraction can result in a increase which could be misinterpreted as an oxygenation increase and hence result in premature conclusions if the effect of blood volume fraction is not taken into account [35]. If all the parameters are considered correctly, can serve as a non-invasive means to probe physiology in vivo. It should be noted that T2 changes, e.g., caused by alterations in tissue water content or distribution are also reflected in and hence should be considered as potential confounders.
Benefits of myocardial mapping at higher magnetic field strengths
The magnetization M of a material in response to an applied magnetic field is given by its magnetic susceptibility χ and the strength of the applied magnetic field H:
This relationship results in a linear increase of magnetic field perturbations induced by microscopic susceptibility changes—the main driving force behind decay—when moving to higher magnetic fields. The effect has been confirmed for myocardial ( = 1/) in vivo rendering mapping at ultrahigh magnetic fields (B0 ≥ 7.0 T) particularly appealing [36] (Figure 2). The enhanced susceptibility effects at 7.0 T may be useful to extend the dynamic range of the sensitivity for monitoring changes and to lower their detection level. Another advantage of performing weighted imaging and mapping at ultrahigh magnetic field strengths (UHF) is that the signal-to-noise ratio (SNR) gain achieved at higher fields can be used to improve the spatial resolution [37, 38]. This reduction in voxel sizes lowers the impact of macroscopic magnetic field gradients on intra-voxel dephasing and hence which otherwise can be a concern especially in the vicinity of strong susceptibility transitions. Transitioning to higher magnetic field strengths runs the additional benefit that the in-phase inter-echo time governed by the fat-water phase shift between the water and main fat peak of about 3.5 ppm is reduced from approximately 4.5 ms (223 Hz) at 1.5 T to 0.96 ms (1,043 Hz) at 7.0 T. This enables rapid acquisition of multiple echoes with different sensitization and facilitates high spatio-temporally resolved myocardial CINE mapping of the human heart [39]. Taking advantage of this technique, mapping at ultrahigh magnetic fields has been suggested as a means to probe myocardial physiology and to advance myocardial tissue characterization.
Figure 2
Challenges and technical solutions for cardiac MRI at ultrahigh magnetic fields
Enabling radio frequency antenna technology
Imaging the heart—a deep-lying target region surrounded by the lung within the large volume of the thorax—at ultrahigh magnetic field strengths poses a severe challenge due to the short wavelength of the proton resonance frequency in tissue (λmyocardium ≈ 12 cm at 7.0 T). As a result, non-uniformities in the transmission field () can cause shading, massive signal drop-off or signal void in the images up to non-diagnostic image quality. These constraints were reported being a concern in CMR at 3.0 T [40] and were to be expected to pose a major obstacle for CMR at UHF.
A plethora of reports have evolved during the past years introducing technical innovations in RF antenna design to overcome non-uniform transmission fields. Local transceiver (TX/RX) and multi-channel transmission arrays in conjunction with multi-channel local receive arrays have been suggested as possible solutions. Eminent developments put building blocks to use consisting of stripline elements [41–45], electrical dipoles [45–51], dielectric resonant antennas [52], slot antennas [53], and loop elements [54–59]. Rigid, flexible and modular configurations have been exploited. Irrespective of the building block technology, a trend toward higher numbers of transmit and receive elements can be observed with the purpose to advance anatomic coverage [47, 54–59] and to add degrees of freedom for transmission field shaping [60].
Figure 3 compiles developments of loop element based transceiver configurations optimized for CMR at 7.0 T. A 4-channel TX/RX [55] (Figure 3A) and an 8-channel TX/RX [58] (Figure 3B) one-dimensional array were reported and extended to a 16-channel two-dimensional design [56] (Figure 3C). A modular 32-channel TX/RX [59] (Figure 3D) array further exploited the two-dimensional building block layout.
Figure 3
Electric dipoles hold the benefit of a linearly polarized current pattern with the RF energy being directed perpendicular to the dipole along the Poynting vector to the subject. As a consequence, the excitation field is symmetrical and uniform and comes with ample depth penetration [48] which renders electric dipoles particularly promising for MR of the upper torso and the heart. This property formed the starting point for explorations into electric dipole configurations [46–48, 50, 51]. Due to their length straight dipole elements are unfavorable if not unfeasible for high density multi-dimensional transceiver coil arrays [48]. To address this issue, the fractionated dipole concept splits the dipole's legs into segments interconnected by capacitors or inductors to achieve dipole shortening. Reduced SAR levels, moderate coupling and homogeneous have been reported for prostate imaging using an eight-element array consisting of fractionated dipoles [51]. A combined 16-channel loop-dipole transceiver array exploiting the fractionated dipole antenna design provided cardiac images acquired at 7.0 T exhibiting high SNR and transmit efficiency [50]. As an alternative, shortening of the effective antenna length can be achieved for a bow tie shaped λ/2-dipole antenna by immersing it in D2O. Following this achievement electric dipole configurations optimized for UHF-CMR have been reported using 8 or 16 bow tie antenna building blocks [47] (Figures 3E,F). As a result of these research efforts, dedicated RF antenna arrays are now available which facilitate cardiac MRI at 7.0 T with rather uniform signal intensities across the heart.
In summary, explorations into enabling RF antenna technology underlined the benefits of many-channel high-density arrays for UHF-CMR.
Ancillary devices for cardiac synchronization
Imaging the heart requires synchronization of the data acquisition with the cardiac cycle. Magneto hydrodynamic (MHD) effects severely disturb the electrocardiogram (ECG) [61–63] commonly applied for cardiac triggering and gating at clinical field strengths [64–66] (Figure 4). Distortions of the ECG's S-T segment are caused by the increased MHD impact during systolic aortic flow [67]. The S-T elevation might be mis-interpreted as an R-wave. Consequently, image quality is impaired due to the mis-detected onset of a cardiac cycle. The propensity to MHD effects is pronounced at ultrahigh magnetic field strengths [42, 68, 69]. An MR-stethoscope has been proposed as an alternative to ECG gating and triggering putting acoustic signals to use which have been reported to be immune to interferences with electromagnetic fields. With this practical solution the first heart tone of the phonocardiogram is detected, which marks the onset of the acoustic cardiac cycle. The minor latency between the onset of the electrophysiological cardiac activity and the onset of the acoustic cardiac activity allows prospectively triggered and retrospectively gated acquisitions. Acoustic triggering can hence be used with all pulse sequences that support ECG triggering without the need of sequence adjustments. Reliable trigger information has been demonstrated when using acoustic cardiac triggering and gating for UHF-CMR (Figure 4) [61, 68, 70]. Further alternatives for cardiac synchronization include post-processing of the ECG signal to reduce MHD induced distortions of the ECG trace [71, 72]. Wideband radar, magnetic field probes or optical systems have been proposed for physiological monitoring [73–75].
Figure 4
imaging techniques
sensitized imaging and mapping is commonly performed employing gradient echo imaging independent of magnetic field strength. To decrease acquisition times multi echo gradient echo techniques (MEGRE) acquiring multiple echoes after each RF excitation instead of only one echo per repetition time TR are recommended for fast mapping (Figure 5A, top). For myocardial mapping cardiac triggered segmented acquisitions are commonly performed in end-expiratory breath hold conditions to avoid respiratory and cardiac motion and to reduce related macroscopic B0 field fluctuations.
Figure 5
The used echo times should be adapted to sufficiently cover the decay. As the contributing fat and water signal are oscillating at different frequencies mapping algorithms must either account for or compensate the varying signal intensity from fat and water. Acquiring -weighted images at times when fat and water are equally contributing to the MR signal (in-phase) is the simplest approach to achieve this goal. At 7.0 T this is the case for echo times being a multiple of about 0.96 ms due to the chemical shift between the main fat and water peaks of approximately 3.5 ppm (1,043 Hz). Acquisition of echoes with an inter-echo spacing of ~1 ms constitutes a challenge due to gradient amplitude and rise time limitations, especially when large acquisition matrix sizes are used. Alternatively, interleaved acquisitions can be performed by distributing the acquisition of neighboring echoes across multiple excitations (Figure 5A, middle). This approach permits low inter-echo spacing even for large matrix sizes but results in longer scan times since more than one TR is required to acquire a full decay series. While mapping at clinical field strengths is limited to single cardiac phase acquisitions, CINE mapping covering the entire cardiac cycle is feasible at UHF [39]. This advanced capability is facilitated by two benefits of UHF-MR. First, due to transversal relaxation time shortening at ultrahigh magnetic fields, TE can be limited to a range of approximately TE = 0 ms to TE = 15 ms to properly sample the myocardial decay. This approach is beneficial for reducing the duration of the gradient echo trains compared to lower magnetic field strengths enabling breath held multi echo CINE acquisitions. Second, the reduced in-phase echo spacing permits acquisition of a sufficient number of echoes needed to cover the signal decay and to provide an appropriate number of data points for signal fitting. Interleaving of echo times can be combined with distributing the acquisition across multiple breath-holds to ease gradient constraints and limit breath-hold durations for each acquisition (Figure 5A, bottom). All described acquisition strategies are capable of producing maps of similar fidelity as illustrated in Figure 5B for a homogenous MR phantom resembling the relaxation properties of human myocardium. To reduce the effect of macroscopic magnetic field gradients on spin dephasing and , small voxel sizes are preferable. Of course this preference has to be balanced with SNR requirements for accurate mapping which can be challenging particularly at lower magnetic field strengths. Figure 5B compares the effect of slice thickness on . While maps acquired with slice thicknesses of 6 mm or above show intra-voxel dephasing and decrease pronounced at the phantom interfaces, this effect is mitigated for slice thicknesses of 4 mm or less resulting in a more uniform map. Employing the described multi-breath hold CINE technique at 7.0 T, CINE mapping with more than 20 cardiac phases is feasible which allows monitoring of myocardial across the cardiac cycle (Figure 6).
Figure 6
In contrast to gradient echo imaging Rapid Acquisition with Relaxation Enhancement (RARE) imaging is rather insensitive to B0 inhomogeneities, provides images free of distortion due to the use of RF refocused echoes and inherently suppresses blood signal. Cardiac RARE imaging at 7.0 T has been shown to be feasible [76]. These results—in conjunction with the challenges and opportunities of myocardial mapping at UHF—build the starting point for explorations into RARE based mapping of the heart at 7.0 T. An evolution time τ inserted after the excitation RF pulse accrues an additional phase to the magnetization that reflects the effect [77] (Figure 7A). As a consequence, the Carr-Purcell-Meiboom-Gill condition cannot be met and destructive interferences between odd and even echo groups configuring the signal in RARE imaging impair the image [78, 79]. Measures to account for this effect include displaced RARE [77, 80], avoiding interferences by discarding one of both echo groups, resulting in an SNR loss of factor two. Split-echo variants hold the benefit that the full available signal is maintained [81] (Figure 7A). A series of weighted images derived from RARE imaging using evolution times τ ranging from 2 to 14 ms is displayed in Figure 7B and demonstrates that the geometric integrity of the RARE images is maintained over the range of applied weighting. The corresponding map is shown in Figure 7C. MEGRE imaging results are shown for comparison and exhibited less myocardium to blood contrast due to the bright-blood characteristic of the technique. Consequently, the delineation of the myocardium in the corresponding map (Figure 7C) is more challenging compared to the RARE based map. The average effective transversal relaxation time derived from RARE imaging compares well to values previously reported for MEGRE approaches [39]. The concerns of RF power deposition and RF non-uniformity of RARE imaging were offset in this preliminary study. Split-echo RARE hence holds the potential to provide an alternative for mapping free of geometric distortion and with high blood myocardium contrast.
Figure 7
Assessment and control of main magnetic field homogeneity
The complex MR signal S (r, t) obtained by a gradient echo technique at a location r and time t after signal excitation is given by:
Ŝ(t) is the magnitude signal, i is the imaginary unit and ϕ is the signal phase. The signal phase ϕ can be written as a function of spatial location and time [82]:
with γ being the gyromagnetic constant of the nucleus (for 1H γ = 2.675 × 108 rad/s/T) and ΔB(r) representing local magnetic field deviations with regard to the main magnetic field strength. ϕ0(r) represents a constant receiver phase offset, while the time-dependent phase component γ·ΔB(r)·t is dominated by the deviation from the static magnetic field and evolves linearly over time [83]. Assuming there are no other external sources of dephasing such as, e.g., motion or flow, the signal phase ϕ serves as a direct measure of deviations from the main magnetic field B0.
Equation (1) can be approximated as:
if a linear B0 gradient within a voxel is assumed [32]. This assumption is justified for typical voxel sizes used in cardiac mapping at 7.0 T with an in-plane spatial resolution of about 1 mm and a slice thickness of 2–4 mm [39, 84]. |ΔB| in Equation (9) includes microscopic magnetic field perturbations resulting from microstructural changes, blood oxygenation changes, iron accumulation etc. which are of diagnostic interest as well as macroscopic field changes, e.g., due to magnet imperfections or strong susceptibility transitions at air tissue interfaces. This dependency highlights the need to monitor and possibly compensate B0 inhomogeneities when performing assessment to make sure that decay reflects microscopic susceptibility changes rather than macroscopic field perturbations.
By acquiring images at two echo times (TE), the local magnetic field variations can be calculated at each voxel by making use of Equation (8):
where ΔB(r) is given in Tesla. This procedure is referred to as B0 mapping. It is also very common to represent local magnetic field variations in Hz by means of off-center frequency maps:
B0 shimming describes the process of adjusting the main magnetic field B0 to improve macroscopic field homogeneity. Active shimming refers to the adjustment of the main magnetic field by making use of dedicated shim coils thereby creating compensatory magnetic fields up to the 5th order of spherical harmonics [85]. Despite the existence of 5th order shim systems, most scanners provide only shim coils up to 2nd order. Active shimming options include fixed shim current settings or shimming modes based on an individual B0 map acquired for a specific subject. For cardiac B0 shimming a cardiac triggered field map acquisition is recommended to avoid motion artifacts. Usually a shim volume of interest is defined covering the target anatomy. Ideally the target volume should cover a small region of interest like the heart, to achieve satisfying field homogeneity even when a limited order of shim coils is available. Figure 8 compares B0 homogeneity in the heart of a healthy subject at 7.0 T after applying a volume shim (Figure 8A, top), which is focused only on the heart and after applying a global shim (Figure 8A, bottom) which includes the entire field of view. Volume selective shimming was found to lead to a significant improvement in macroscopic B0 homogeneity vs. global B0 shimming (Figure 8B) [39].
Figure 8
Due to increased susceptibility effects, magnetic field inhomogeneities are pronounced at higher magnetic field strengths [36] (compare Equation 6). This is often a concern for weighted imaging or mapping at high and ultrahigh magnetic fields. However, by employing dedicated shimming approaches, a mean peak-to-peak off-resonance frequency across the human heart of 80 Hz [39] was reported at 7.0 T. For the left ventricle a B0 peak-to-peak difference of approximately 65 Hz was observed at 7.0 T after volume selective shimming [39]. These results were obtained with a second order shim system so that the same level of B0 uniformity reported here should be achievable with the current and next generation of 7.0 T scanners. The frequency shift across the heart obtained at 7.0 T compares well with previous 3.0 T studies which reported a peak-to-peak off-resonance variation of (176 ± 30) Hz over the left ventricle and (121 ± 31) Hz over the right ventricle with the use of localized linear and second-order shimming [86]. The use of an enhanced locally optimized shim algorithm, which is tailored to the geometry of the heart, afforded a reduction of the peak-to-peak frequency variation over the heart from 235 to 86 Hz at 3.0 T [87]. Another study showed a peak-to-peak off-resonance of (71 ± 14) Hz for short axis views acquired at 1.5 T [88] using global shimming. While dedicated shim routines revealed competitive results at 7.0 T vs. lower magnetic fields strength, the feasibility of using these approaches in a clinical setting remains to be investigated.
The achievement of macroscopic B0 homogeneity across the heart at 7.0 T which is competitive with that obtained at lower magnetic field strengths provides encouragement to pursue susceptibility-based myocardial mapping at ultrahigh fields. Yet, with the arrival of CINE mapping techniques enabled by 7.0 T [39], also temporal B0 fluctuations across the cardiac cycle and their implications on need to be considered for a meaningful interpretation of dynamic -weighted acquisitions. Temporal B0 variation across the cardiac cycle was reported to be negligible at 1.5 T [89], but due to the increase of susceptibility effects at ultrahigh fields further investigations of this potential confounder were required at 7.0 T. It should be noted that -weighted contrast is determined by magnetic field gradients rather than absolute magnetic field strength, thus it is essential to investigate the change of these gradients over the cardiac cycle. This was done at 7.0 T by assessing macroscopic B0 gradients across the cardiac cycle in the heart of healthy volunteers together with high temporal and spatial resolution maps [84, 90]. -weighted series of short-axis views were acquired using a MEGRE CINE approach (Figure 9A, top) for cardiac phase resolved B0 and mapping. Temporally-resolved B0 maps of the heart were calculated (Figure 9A, middle). Macroscopic intra-voxel field gradients were determined for each cardiac phase and their fluctuations were analyzed across the cardiac cycle. The septal in-plane gradients were found to be significantly larger compared to through-plane gradients within a voxel [with a mean in-plane field dispersion of (2.5 ± 0.2) Hz/mm against the a mean through-plane field dispersion of (0.4 ± 0.1) Hz/mm] [91].
Figure 9
In order to evaluate how these B0 gradients affect measurements, the B0 gradient-induced change of represented as Δ was estimated [84, 91]. Figure 9B shows the plot of mean septal , intra-voxel B0 gradients and estimated gradient-induced Δ over the cardiac cycle averaged over a group of healthy subjects. The mean septal per cardiac phase was found to vary over the cardiac cycle in a range of approximately 23% of the total mean for all phases. Yet, the temporal range of mean Δ induced by the calculated intra-voxel macroscopic B0 gradients represented only a 5% change of total mean . The remaining 18% were suggested to reflect microscopic B0 gradient changes (potentially caused by physiological events) rather than macroscopic field inhomogeneities [84, 91].
In summary it can be concluded that if careful shimming is applied, macroscopic magnetic field inhomogeneities are not of concern for myocardial mapping even at a magnetic field strength as high as 7.0 T. Also dynamic B0 fluctuations across the cardiac cycle can be considered negligible in the ventricular septum. These findings represent an essential prerequisite for meaningful interpretation of myocardial and its dynamics across the cardiac cycle.
Data post-processing
The effective transverse relaxation time can be estimated by fitting an exponential function to a series of gradient echo images with different weighting, i.e., different echo times (compare Equations 3, 4; Figure 1). A common way of calculating such a fit by avoiding non-linear fitting procedures, is to calculate the natural logarithm of the acquired signal intensities and apply a least squares linear fit to the resulting data. This procedure is fast and produces the best solution to represent the linearized data in a least squares sense. Further to this, non-linear fitting approaches which can be applied directly to the measured data are available. The algorithms employ non-linear models in conjunction with optimization algorithms like Levenberg-Marquardt or Simplex [92]. Iterative non-linear fitting algorithms can provide improved fitting accuracy, but can be prone to careful initialization. Fast linear fitting can be used to initialize non-linear fitting algorithms leading to increased robustness and faster convergence. Irrespective of the kind of employed fitting procedure, care should be taken for voxels with intensities being at the noise level and hence potentially deteriorating the fit quality. Low SNR is a common problem particularly for myocardial mapping where acquisition times are often limited by the tolerable breath hold time. Such voxels should either be excluded from fitting—a procedure referred to as truncation—or included in the fit model as a constant noise term.
Most commercial MR systems support exponential fitting algorithms as part of the systems' software, but using customized fitting routines is often beneficial for research. First, it is often unclear what model or fitting approach is used by commercial software and how good the fit quality was, i.e., how well the fit describes the measured data. Measures like the coefficient of determination R2 or the standard deviation of the fit [93] should be used to evaluate the reliability of the results. Taking fitting results for granted without considering fit quality may lead to wrong results and eventually wrong conclusions. Second, tailored fitting procedures offer the freedom to select the most appropriate fit model, optimization approach, truncation threshold, etc. for the particular research question. Depending on the kind of application, it may for example make sense to use a bi-exponential or multi-exponential model instead of a mono-exponential approach.
In contrast to qualitative signal intensity images, relaxation maps like maps offer the advantage of providing quantitative, comparable results. Yet, effects like B0 inhomogeneities or signal noise can impair the assessment of and lead to wrong results. Dedicated B0 shimming approaches help to mitigate the impact of macroscopic magnetic field inhomogeneities. A reduction in voxel size can further reduce the influence of B0 gradients on (Figure 5B). Yet, reducing voxel size results in an SNR loss, which can induce poor fit quality offsetting the benefit of the smaller voxel size. While a reduction in SNR might be counteracted by signal averaging and increasing acquisition times in static acquisition situations (e.g., MRI of the brain), it is often not feasible in cardiac applications, where acquisitions need to be synchronized with the cardiac cycle and where it is common to utilize breath held conditions constraining the viable window of data acquisition to few seconds. This issue is even further pronounced in patients suffering from cardiac diseases and for acquisitions at high spatial or temporal resolution such as CINE mapping.
Image de-noising presents a viable solution to address this constraint and allows the use of small voxel sizes while still achieving acceptable fit quality. Powerful de-noising approaches like non-local means filtering [94] are readily available and can greatly improve SNR with minimal loss of information. De-noising of the maps is not recommended, because low fit quality of fitting noisy data can result in large errors that might even be enlarged or spread out by filtering. Also algorithms that estimate the noise level from the provided data will fail if presented with maps. Filtering of the magnitude images prior to fitting instead represents a robust way of improving fit quality and has been shown to increase fitting accuracy and precision [84, 95, 96] without the risk of introducing large errors. Figure 10 illustrates an example of how mapping can benefit from noise filtering. Here a 30% reduction of fit standard deviation was achieved by noise filtering. For applications where SNR is limited such as myocardial mapping, image de-noising approaches provide a good solution to improve mapping results.
Figure 10
Insights from in vivo human myocardial mapping at ultrahigh fields
The technological and methodological developments in ultrahigh field CMR outlined above, permit for the first time the in vivo assessment of temporal myocardial changes across the cardiac cycle. Initial studies have applied these advances to gain first insights from using this technique in healthy volunteers and patients suffering from cardiovascular diseases to investigate their feasibility and potential [39, 84].
CINE myocardial mapping in healthy subjects
A first study systematically investigating the temporal changes of myocardial across the cardiac cycle in healthy subjects at 7.0 T was published in 2016 [84]. The authors analyzed the time course of myocardial throughout the cardiac cycle along with basic myocardial morphology, i.e., ventricular septal wall thickness and inner left ventricular radius as potential confounders of . The results demonstrated that myocardial obtained correlates linearly with the myocardial wall thickness [84], (Figure 11). The same study also showed that in the ventricular septum changes periodically across the cardiac cycle [84]. It increases in systole—the part of the cardiac cycle when the ventricles contract—and decreases in diastole—the part of the cardiac cycle when the heart relaxes and refills with blood (Figure 12A). The mean systole to diastole ratio was found to be approximately 1.1. Despite the numerous factors affecting such as blood volume fraction, hematocrit etc., myocardial is still often regarded as surrogate for tissue oxygenation. Interpreting to reflect tissue oxygenation, the observed systolic increase would imply an increase in left myocardial oxygenation during systole. This is contrary to physiological knowledge. Instead, changes in myocardial blood volume fraction induced by variations in blood pressure and resulting myocardial wall stress are believed to be responsible for the observed cyclic changes [84]. The contraction of the heart muscle compresses the intramyocardial vasculature such that inflow of arterial blood is interrupted while deoxygenated blood is squeezed out of the myocardium toward the venous coronary sinus [97–100] (Figure 12B). The major systolic decrease in blood volume fraction in the myocardium reduces the amount of deoxygenated hemoglobin per tissue volume, thereby increasing—instead of lowering— during systole. Previous studies of skeletal muscle have also linked changes to alterations in tissue pH and resulting changes of the tissue water content and distribution after exercise [101, 102]. These studies have examined baseline and post-exercise conditions, which are difficult to compare with the heart which is constantly exercising. Still, T2 changes driven by tissue water content and distribution changes should be considered as a potential source of changes also in the heart. The hypothesis, that the observed periodic changes could be induced by macroscopic B0 field variations induced by changes in bulk morphology between systole and diastole, was carefully investigated but not confirmed. Both, in silico magneto static simulations and in vivo temporally resolved B0 mapping, showed negligible impact of cardiac morphology on the macroscopic B0 field in the ventricular septum and hence [84].
Figure 11
Figure 12
Myocardial mapping in patients with cardiovascular diseases
Besides the application of myocardial mapping at ultrahigh magnetic fields in healthy volunteers, first investigations were carried out to explore the potential of the technique to distinguish between healthy and pathologic myocardium. These early UHF-CMR studies focused on mapping in patients with hypertrophic cardiomyopathy (HCM). HCM is the most common inherited cardiac disease affecting about 0.2–0.5% of the general population [103, 104]. The disease is characterized by an increase in myocardial wall thickness related to myocyte hypertrophy, microstructural changes like myocardial disarray, fibrosis and microvascular dysfunction. HCM patients often remain asymptomatic, but the disease can have a severe outcome in a subgroup of patients where it may cause heart failure and sudden unexpected cardiac death (SCD) in any age group. SCD and has been reported to affect about 6% of HCM patients within a mean follow up time of (8 ± 7) years [105]. This renders risk stratification vital for HCM patients. CMR plays an important role in the diagnosis and prognosis of HCM [106]. While a number of SCD risk factors in HCM have been identified such as degree of hypertrophy or presence of fibrosis the task remains challenging [107]. Consequently, basic research efforts and clinical science activities are required to better characterize HCM patient populations and to direct appropriate therapies to those at risk.
Based on the structural and physiologic changes, differences in myocardial were hypothesized in HCM patients compared to healthy controls. This hypothesis was investigated using high spatiotemporal resolution mapping at 7.0 T (Figure 13). It was found that septal is significantly increased in HCM with mean septal being (17.5 ± 1.4) ms in a cohort of HCM patients compared to (13.7 ± 1.1) ms in a group of gender, age and body mass index matched healthy controls. While variations of myocardial across the cardiac cycle have been attributed to changes in myocardial blood volume fraction rather than changes in tissue oxygenation [84], two main factors are assumed to cause the observed overall increase in HCM. Improved tissue oxygenation in the diseased myocardium in the case of HCM is unlikely. Instead, T2 has been reported to be elevated in HCM [108] related to presence of inflammation and edema. A T2 increase would also result in increased as seen from Equation (1). Further to this, reduced myocardial perfusion and ischemia are common in HCM [109], effectively reducing the tissue blood volume fraction resulting in a increase as suggested by Equation (5). These conditions are also associated with a higher risk for a poor outcome in HCM patients [110]. With this in mind it is fair to conclude that myocardial mapping could be beneficial for a better understanding of cardiac (patho)physiology in vivo with the ultimate goal to support risk stratification in HCM.
Figure 13
Conclusion and future directions
The progress in myocardial mapping at ultrahigh magnetic fields is promising [111–113]. Yet, there are still a number of questions to be answered and the clinical benefit remains to be carefully investigated. This requires further efforts to tackle unsolved problems and unmet needs standing in the way en route to broader clinical studies. For example, the relatively long breath hold times required for the acquisition of high spatiotemporal resolution maps constitute a challenge particularly in cardiac patients. Free breathing acquisition techniques could offset this constraint permitting broader application and full 3D heart coverage. This would also help to further investigate the effect of through-plane motion. Acquisition approaches like simultaneous multi-slice excitation can be used to reduce scan times while multi-channel transmit systems can be employed to balance excitation field homogeneity and RF power deposition constraints [114, 115].
Based on the multifaceted contributions of physiological parameters on research will not stop at just mapping myocardial . Tailored acquisition schemes, data post-processing, analysis and interpretation will allow exploiting the wealth of information encoded into . For example high spatial resolution mapping facilitated by ultrahigh magnetic field strengths might be beneficial to gain a better insight into the myocardial microstructure in vivo with the ultimate goal to visualize myocardial fibers or to examine their helical angulation, since the susceptibility effects depend on the orientation of blood filled capillaries with regard to the external magnetic field [116]. Myocardial fiber tracking using mapping holds the promise to be less sensitive to bulk motion than diffusion-weighted MR of the myocardium [117, 118]. The increased susceptibility contrast available at 7.0 T could be exploited to quantitatively study iron accumulation in the heart with high sensitivity and spatial resolution superior to what can be achieved at 1.5 and 3.0 T. This requires the determination of norm values for healthy myocardial at 7.0 T as a mandatory precursor to broader clinical studies.
At the same time small animal studies employing cardiac disease models can provide a valuable contribution to understanding the underlying biophysical principles and (patho)physiological contrast mechanisms governing . Unlike human studies they offer the unique possibility to directly compare in vivo findings by MRI with ex vivo histology, the gold standard for tissue characterization. Recent such studies indicate that might provide not only an alternative for detection of both replacement and diffuse fibrosis without the need for exogenous contrast agents, but also has potential to distinguish the two by means of relaxation time changes induced by the presence of collagen and other fibrotic elements in the extracellular matrix [17, 18]. This could provide new diagnostic means to a large group of patients excluded from contrast agent injections due to renal insufficiencies. A study employing a mouse myocardial ischemia/reperfusion model has provided first insights into in vivo quantification of changes in the mouse myocardium in relation to tissue damage [16]. Local decrease of was found in the infarct zone and associated with deposition of collagen. The authors describe that varies dynamically during infarct development suggesting that it may be used to discriminate between acute and chronic infarctions. Taken together, by concordance the findings between human studies and cardiac disease models of small rodents will provide stronger evidence for fundamental understandings of myocyte biology, and cardiac performance with the goal to provide a more accurate diagnosis and risk stratification. Thanks to the sensitivity gain at 7.0 T the spatial fidelity feasible for T2* mapping in humans approaches the relative anatomical spatial resolution—in terms of number of voxels with respect to anatomy—demonstrated for cardiac imaging in animal models [119, 120]. This achievement is translatable into opportunities for discovery and translational research.
The ability to probe for changes in myocardial tissue oxygenation using sensitized imaging/mapping offers the potential to address some of the spatial and temporal resolution constraints of conventional first pass perfusion imaging and holds the promise to obviate the need for exogenous contrast agents. Since microscopic susceptibility increases with field strength, thus making oxygenation sensitivity due to ischemic (patho) physiology more pronounced, mapping at 7.0 T might be beneficial to address some of the BOLD sensitivity constraints reported for the assessment of regional myocardial oxygenation changes in the presence of coronary artery stenosis [121] or for the characterization of vasodilator-induced changes of myocardial oxygenation at 1.5 T and at 3.0 T [10].
The pace of discovery is exciting and a powerful motivator to transfer the lessons learned from mapping research at 7.0 T into the clinical scenario. These efforts are fueled by the quest for advancing the capabilities of quantitative MRI and the wish to overcome the need of exogenous contrast agent injections. The requirements of mapping at 7.0 T are likely to pave the way for further advances in MR technology and MR systems design. With appropriate multi transmit systems offering more than 16 transmit channels each providing at least 4 kW peak power, an optimistically-inclined scientist might envision the implementation of high density transceiver arrays with 64 and more elements with the ultimate goal to break ground for a many element upper torso or even a body coil array. This vision continues to motivate new research on integrated multi-channel transmission systems [122], on novel RF pulse design, on RF coil design together with explorations into ideal current patterns yielding optimal signal-to-noise-ratio for UHF-CMR [123]. Perhaps another development is the move toward myocardial mapping using reduced field of views zoomed into the target anatomy enabled by spatially selective excitation techniques which put the capabilities of parallel transmission technology to good use. With more than 45,000 examinations already performed at 7.0 T, the reasons for employing UHF-MR in translational research and for moving UHF-MR into clinical applications are more compelling than ever. This provides strong motivation to put further weight behind pushing the solution of the many remaining problems. As an important step toward this goal a system manufacturer has recently filed for FDA clearing for the clinical use of a 7.0 T system. With this development we can expect more pioneering research institutions, university hospitals and large clinics to become early adopters of CMR at 7.0 T and start harvesting knowledge and know-how that will benefit clinical applications.
Statements
Ethics statement
In vivo studies of which data is presented in this work were carried out in accordance with the guidelines of the local ethical committee (registration number DE/CA73/5550/09, Landesamt für Arbeitsschutz, Gesundheitsschutz und technische Sicherheit, Berlin, Germany) with written informed consent from all subjects in compliance with the local institutional review board guidelines. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocols were approved by the local ethical committee.
Author contributions
TH and TN wrote the manuscript with help from KP, MK, and TS.
Funding
This work was supported (in part, TH and TN) by the DZHK (German Centre for Cardiovascular Research, partner site Berlin, BER 601) and by the BMBF (Federal Ministry of Education and Research, Berlin, Germany, FKZ 81Z6100161). TN received support by the BMBF (Federal Ministry of Education and Research, Berlin, Germany, FKZ 01QE1501B) and the EUROSTARS program (E! 9340 hearRT-4-EU).
Acknowledgments
The authors wish to acknowledge the team at the Berlin Ultrahigh Field Facility (B.U.F.F.) at the Max-Delbrueck Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany; Jeanette Schulz-Menger from the working group for Cardiac Magnetic Resonance, Charite', Berlin, Germany; Peter Kellman (National Institutes of Health, NHLBI, Laboratory of Cardiac Energetics, Bethesda, USA); who kindly contributed examples of their pioneering work or other valuable assistance.
Conflict of interest
TN is founder and CEO of MRI.TOOLS GmbH, Berlin, Germany. The reviewer AIS and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review. The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
von Knobelsdorff-BrenkenhoffFSchulz-MengerJ. Role of cardiovascular magnetic resonance in the guidelines of the European Society of Cardiology. J Cardiovasc Magn Reson. (2016) 18:6. 10.1186/s12968-016-0225-6
2.
CapturGManistyCMoonJC. Cardiac MRI evaluation of myocardial disease. Heart (2016) 102:1429–35. 10.1136/heartjnl-2015-309077
3.
OgawaSMenonRSTankDWKimSGMerkleHEllermannJMet al. Functional brain mapping by blood oxygenation level-dependent contrast magnetic resonance imaging. A comparison of signal characteristics with a biophysical model. Biophys J. (1993) 64:803–12. 10.1016/S0006-3495(93)81441-3
4.
WackerCMHartlepAWPflegerSSchadLRErtlGBauerWR. Susceptibility-sensitive magnetic resonance imaging detects human myocardium supplied by a stenotic coronary artery without a contrast agent. J Am Coll Cardiol. (2003) 41:834–40. 10.1016/S0735-1097(02)02931-5
5.
FriedrichMNiendorfTSchulz-MengerJ. Blood oxygen level-dependent magnetic resonance imaging in patients with stress-induced angina. ACC Curr J Rev. (2004) 13:30. 10.1016/j.accreview.2004.03.013
6.
JahnkeCGebkerRMankaRSchnackenburgBFleckEPaetschI. Navigator-gated 3D blood oxygen level–dependent CMR at 3.0-T for detection of stress-induced myocardial ischemic reactions. JACC Cardiovasc Imaging (2010) 3:375–84. 10.1016/j.jcmg.2009.12.008
7.
KaramitsosTDLeccisottiLArnoldJRRecio-MayoralABhamra-ArizaPHowellsRKet al. Relationship between regional myocardial oxygenation and perfusion in patients with coronary artery disease insights from cardiovascular magnetic resonance and positron emission tomography. Circ Cardiovasc Imaging (2010) 3:32–40. 10.1161/CIRCIMAGING.109.860148
8.
TsaftarisSATangRZhouXLiDDharmakumarR. Ischemic extent as a biomarker for characterizing severity of coronary artery stenosis with blood oxygen-sensitive MRI. J Magn Reson Imaging (2012) 35:1338–48. 10.1002/jmri.23577
9.
WackerCMBockMHartlepAWBeckGvan KaickGErtlGet al. Changes in myocardial oxygenation and perfusion under pharmacological stress with dipyridamole: assessment using T*2 and T1 measurements. Magn Reson Med. (1999) 41:686–95.
10.
VohringerMFlewittJAGreenJDDharmakumarRWangJJr.TybergJVet al. Oxygenation-sensitive CMR for assessing vasodilator-induced changes of myocardial oxygenation. J Cardiovasc Magn Reson. (2010) 12:20. 10.1186/1532-429X-12-20
11.
UtzWJordanJNiendorfTStoffelsMLuftFCDietzRet al. Blood oxygen level-dependent MRI of tissue oxygenation: relation to endothelium-dependent and endothelium-independent blood flow changes. Arterioscler Thromb Vasc Biol. (2005) 25:1408–13. 10.1161/01.ATV.0000170131.13683.d7
12.
GuenschDPFischerKFlewittJAFriedrichMG. Impact of intermittent apnea on myocardial tissue oxygenation–a study using oxygenation-sensitive cardiovascular magnetic resonance. PLoS ONE (2013) 8:e53282. 10.1371/journal.pone.0053282
13.
GuenschDPFischerKFlewittJAFriedrichMG. Myocardial oxygenation is maintained during hypoxia when combined with apnea–a cardiovascular MR study. Physiol. Rep. (2013) 1:e00098. 10.1002/phy2.98
14.
GuenschDPFischerKFlewittJAYuJLukicRFriedrichJAet al. Breathing manoeuvre-dependent changes in myocardial oxygenation in healthy humans. Eur Heart J Cardiovasc Imaging (2014) 15:409–14. 10.1093/ehjci/jet171
15.
FischerKGuenschDPShieNLebelJFriedrichMG. Breathing maneuvers as a vasoactive stimulus for detecting inducible myocardial ischemia - an experimental cardiovascular magnetic resonance study. PLoS ONE (2016) 11:e0164524. 10.1371/journal.pone.0164524
16.
AguorENArslanFvan de KolkCWNederhoffMGDoevendansPAvan EchteldCJet al. Quantitative T2* assessment of acute and chronic myocardial ischemia/reperfusion injury in mice. Magn Reson Mater Phys Biol Med. (2012) 25:369–79. 10.1007/s10334-012-0304-0
17.
de JongSZwanenburgJJVisserFvan der NagelRvan RijenHVVosMAet al. Direct detection of myocardial fibrosis by MRI. J Mol Cell Cardiol. (2011) 51:974–9. 10.1016/j.yjmcc.2011.08.024
18.
van NieropBJBaxNANelissenJLArslanFMotaalAGde GraafLet al. Assessment of myocardial fibrosis in mice using a T2*-weighted 3D radial magnetic resonance imaging sequence. PLoS ONE (2015) 10:e0129899. 10.1371/journal.pone.0129899
19.
HeT. Cardiovascular magnetic resonance T2* for tissue iron assessment in the heart. Quant Imaging Med Surg. (2014) 4:407–12. 10.3978/j.issn.2223-4292.2014.10.05
20.
AndersonLHoldenSDavisBPrescottECharrierCBunceNet al. Cardiovascular T 2-star(T 2*) magnetic resonance for the early diagnosis of myocardial iron overload. Eur Heart J. (2001) 22:2171–9. 10.1053/euhj.2001.2822
21.
MavrogeniS. Evaluation of myocardial iron overload using magnetic resonance imaging. Blood Transfusion (2009) 7:183–7. 10.2450/2008.0063-08
22.
CarpenterJ-PHeTKirkPRoughtonMAndersonLJde NoronhaSVet al. On T2* magnetic resonance and cardiac iron. Circulation (2011) 123:1519–28. 10.1161/CIRCULATIONAHA.110.007641
23.
FriedrichMGKaramitsosTD. Oxygenation-sensitive cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2013) 15:43. 10.1186/1532-429X-15-43
24.
HaackeEMBrownRWThompsonMRVenkatesanR editors. Magnetization, relaxation and the bloch equation. In: Magnetic Resonance Imaging - Physical Principles and Sequence Design. 1st ed. New York, NY: John Wiley & Sons (1999). p. 51–64.
25.
ZhaoJMClingmanCSNärväinenMJKauppinenRAvan ZijlP. Oxygenation and hematocrit dependence of transverse relaxation rates of blood at 3T. Magn. Reson. Med. (2007) 58:592–7. 10.1002/mrm.21342
26.
OgawaSLeeTMKayARTankDW. Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc Natl Acad Sci USA. (1990) 87:9868–72. 10.1073/pnas.87.24.9868
27.
BrownRWChengY-CNHaackeEMThompsonMRVenkatesanR. Magnetic Resonance Imaging: Physical Principles and Sequence Design. New York, NY: John Wiley & Sons (2014).
28.
BernsteinMAKingKFZhouXJ. Handbook of MRI Pulse Sequences. New York, NY: Elsevier (2004).
29.
OgawaSTankDWMenonREllermannJMKimSGMerkleHet al. Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging. Proc Natl Acad Sci USA. (1992) 89:5951–5. 10.1073/pnas.89.13.5951
30.
BauerWRNadlerWBockMSchadLRWackerCHartlepAet al. Theory of the BOLD effect in the capillary region: an analytical approach for the determination of T2 in the capillary network of myocardium. Magn Reson Med. (1999) 41:51–62.
31.
GuenschDPNadeshalingamGFischerKStalderAFFriedrichMG. The impact of hematocrit on oxygenation-sensitive cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2016) 18:42. 10.1186/s12968-016-0262-1
32.
YablonskiyDAHaackeEM. Theory of NMR signal behavior in magnetically inhomogeneous tissues: the static dephasing regime. Magn Reson Med. (1994) 32:749–63. 10.1002/mrm.1910320610
33.
ChristenTLemassonBPannetierNFarionRSegebarthCRémyCet al. Evaluation of a quantitative blood oxygenation level-dependent (qBOLD) approach to map local blood oxygen saturation. NMR Biomed. (2011) 24:393–403. 10.1002/nbm.1603
34.
SpeesWMYablonskiyDAOswoodMCAckermanJJ. Water proton MR properties of human blood at 1.5 tesla: magnetic susceptibility, T1, T2, T2*, and non-lorentzian signal behavior. Magn Reson Med. (2001) 45:533–42. 10.1002/mrm.1072
35.
NiendorfTPohlmannAArakelyanKFlemmingBCantowKHentschelJet al. How bold is blood oxygenation level-dependent (BOLD) magnetic resonance imaging of the kidney? Opportunities, challenges and future directions. Acta Physiol. (2015) 213:19–38. 10.1111/apha.12393
36.
MeloniAHezelFPositanoVKeilbergPPepeALombardiMet al. Detailing magnetic field strength dependence and segmental artifact distribution of myocardial effective transverse relaxation rate at 1.5, 3.0, and 7.0 T. Magn Reson Med. (2014) 71:2224–30. 10.1002/mrm.24856
37.
OcaliOAtalarE. Ultimate intrinsic signal-to-noise ratio in MRI. Magn Reson Med. (1998) 39:462–73. 10.1002/mrm.1910390317
38.
UgurbilKAdrianyGAndersenPChenWGarwoodMGruetterRet al. Ultrahigh field magnetic resonance imaging and spectroscopy. Magn Reson Imaging (2003) 21:1263–81. 10.1016/j.mri.2003.08.027
39.
HezelFThalhammerCWaicziesSSchulz-MengerJNiendorfT. High spatial resolution and temporally resolved mapping of normal human myocardium at 7.0 tesla: an ultrahigh field magnetic resonance feasibility study. PLoS ONE (2012) 7:e52324. 10.1371/journal.pone.0052324
40.
MuellerAKouwenhovenMNaehleCPGiesekeJStrachKWillinekWAet al. Dual-source radiofrequency transmission with patient-adaptive local radiofrequency shimming for 3.0-T cardiac mr imaging: initial experience. Radiology (2012) 263:77–85. 10.1148/radiol.11110347
41.
VaughanJTSnyderCJDelaBarreLJBolanPJTianJBolingerLet al. Whole-body imaging at 7T: preliminary results. Magn Reson Med. (2009) 61:244–8. 10.1002/mrm.21751
42.
SnyderCJDelaBarreLMetzgerGJvan de MoortelePFAkgunCUgurbilKet al. Initial results of cardiac imaging at 7 Tesla. Magn Reson Med. (2009) 61:517–24. 10.1002/mrm.21895
43.
MaderwaldSOrzadaSSchaeferLCBitzAKBroteIKraffOet al editors. 7T Human in vivo cardiac imaging with an 8-channel transmit/receive array. In: Proceedings of the International Society for Magnetic Resonance Medicine.Honolulu, HI (2009).
44.
van den BergenBKlompDWRaaijmakersAJde CastroCABoerVOKroezeHet al. Uniform prostate imaging and spectroscopy at 7 T: comparison between a microstrip array and an endorectal coil. NMR Biomed. (2011) 24:358–65. 10.1002/nbm.1599
45.
IpekORaaijmakersAJKlompDWLagendijkJJLuijtenPRvan den BergCA. Characterization of transceive surface element designs for 7 tesla magnetic resonance imaging of the prostate: radiative antenna and microstrip. Phys Med Biol. (2012) 57:343–55. 10.1088/0031-9155/57/2/343
46.
WinterLÖzerdemCHoffmannWSantoroDMüllerAWaicziesHet al. Design and evaluation of a hybrid radiofrequency applicator for magnetic resonance imaging and RF induced hyperthermia: electromagnetic field simulations up to 14.0 tesla and proof-of-concept at 7.0 tesla. PLoS ONE (2013) 8:e61661. 10.1371/journal.pone.0061661
47.
OezerdemCWinterLGraesslAPaulKElsAWeinbergerOet al. 16-channel bow tie antenna transceiver array for cardiac MR at 7.0 tesla. Magn Reson Med. (2016) 75:2553–65. 10.1002/mrm.25840
48.
RaaijmakersAJIpekOKlompDWPossanziniCHarveyPRLagendijkJJet al. Design of a radiative surface coil array element at 7 T: the single-side adapted dipole antenna. Magn Reson Med. (2011) 66:1488–97. 10.1002/mrm.22886
49.
IpekORaaijmakersAJLagendijkJJLuijtenPRvan den BergCA. Intersubject local SAR variation for 7T prostate MR imaging with an eight-channel single-side adapted dipole antenna array. Magn Reson Med. (2014) 71:1559–67. 10.1002/mrm.24794
50.
ErturkMARaaijmakersAJAdrianyGUgurbilKMetzgerGJ. A 16-channel combined loop-dipole transceiver array for 7 Tesla body MRI. Magn Reson Med. (2017) 77:884–94. 10.1002/mrm.26153
51.
RaaijmakersAJEItaliaanderMVoogtIJLuijtenPRHoogduinJMKlompDWJet al. The fractionated dipole antenna: a new antenna for body imaging at 7 Tesla. Magn Reson Med. (2016) 75:1366–74. 10.1002/mrm.25596
52.
AussenhoferSAWebbAG. An eight-channel transmit/receive array of TE01 mode high permittivity ceramic resonators for human imaging at 7T. J Magn Reson. (2014) 243:122–9. 10.1016/j.jmr.2014.04.001
53.
AlonLDenizCMBrownRSodicksonDKCollinsCM editors. A slot antenna concept for high fidelity body imaging at ultra high field. In: Proceedings of the International Society for Magnetic Resonance Medicine.Singapore (2016).
54.
VersluisMJTsekosNSmithNBWebbAG. Simple RF design for human functional and morphological cardiac imaging at 7tesla. J Magn Reson. (2009) 200:161–6. 10.1016/j.jmr.2009.06.014
55.
DieringerMARenzWLindelTSeifertFFrauenrathTvon Knobelsdorff-BrenkenhoffFet al. Design and application of a four-channel transmit/receive surface coil for functional cardiac imaging at 7T. J Magn Reson Imaging (2011) 33:736–41. 10.1002/jmri.22451
56.
ThalhammerCRenzWWinterLHezelFRiegerJPfeifferHet al. Two-Dimensional sixteen channel transmit/receive coil array for cardiac MRI at 7.0 T: design, evaluation, and application. J Magn Reson Imaging (2012) 36:847–57. 10.1002/jmri.23724
57.
WinterLKellmanPRenzWGräßlAHezelFThalhammerCet al. Comparison of three multichannel transmit/receive radiofrequency coil configurations for anatomic and functional cardiac MRI at 7.0T: implications for clinical imaging. Eur Radiol. (2012) 22:2211–20. 10.1007/s00330-012-2487-1
58.
GraesslAWinterLThalhammerCRenzWKellmanPMartinCet al. Design, evaluation and application of an eight channel transmit/receive coil array for cardiac MRI at 7.0 T. Eur J. Radiol. (2013) 82:752–9. 10.1016/j.ejrad.2011.08.002
59.
GraesslARenzWHezelFDieringerMAWinterLOezerdemCet al. Modular 32-channel transceiver coil array for cardiac MRI at 7.0T. Magn Reson Med. (2014) 72:276–90. 10.1002/mrm.24903
60.
AdrianyGRitterJVaughanJTUgurbilKMoortelePF editors. Experimental verification of enhanced B1 shim performance with a Z-encoding RF coil array at 7 Tesla. In: Proceedings of the International Society for Magnetic Resonance Medicine.Stockholm (2010).
61.
FrauenrathTHezelFRenzWde GeyerTdODieringerMKnobelsdorff-BrenkenhoffFVet al. Acoustic cardiac triggering: a practical solution for synchronization and gating of cardiovascular magnetic resonance at 7 Tesla. J Cardiovasc Magn Reson. (2010) 12:67. 10.1186/1532-429X-12-67
62.
TogawaTOkaiOOshimaM. Observation of blood flow E.M.F. in externally applied strong magnetic field by surface electrodes. Med Biol Eng. (1967) 5:169–70. 10.1007/BF02474505
63.
StuberMBotnarRMFischerSELamerichsRSminkJHarveyPet al. Preliminary report on in vivo coronary MRA at 3 Tesla in humans. Magn Reson Med. (2002) 48:425–9. 10.1002/mrm.10240
64.
LanzerPBartaCBotvinickEHWiesendangerHUModinGHigginsCB. ECG-synchronized cardiac MR imaging: method and evaluation. Radiology (1985) 155:681–6. 10.1148/radiology.155.3.4001369
65.
FischerSEWicklineSALorenzCH. Novel real-time R-wave detection algorithm based on the vectorcardiogram for accurate gated magnetic resonance acquisitions. Magn Reson Med. (1999) 42:361–70.
66.
ChiaJMFischerSEWicklineSALorenzCH. Performance of QRS detection for cardiac magnetic resonance imaging with a novel vectorcardiographic triggering method. J Magn Reson Imaging (2000) 12:678–88. 10.1002/1522-2586(200011)12:5<678::aid-jmri4>3.0.co;2-5
67.
BeckerMFrauenrathTHezelFKrombachGAKremerUKoppersBet al. Comparison of left ventricular function assessment using phonocardiogram- and electrocardiogram-triggered 2D SSFP CINE MR imaging at 1.5 T and 3.0 T. Eur Radiol. (2010) 20:1344–55. 10.1007/s00330-009-1676-z
68.
FrauenrathTHezelFHeinrichsUKozerkeSUttingJFKobMet al. Feasibility of cardiac gating free of interference with electro-magnetic fields at 1.5 Tesla, 3.0 Tesla and 7.0 Tesla using an MR-stethoscope. Invest Radiol. (2009) 44:539–47. 10.1097/RLI.0b013e3181b4c15e
69.
BrandtsAWestenbergJJVersluisMJKroftLJSmithNBWebbAGet al. Quantitative assessment of left ventricular function in humans at 7 T. Magn Reson Med. (2010) 64:1471–7. 10.1002/mrm.22529
70.
FrauenrathTNiendorfTMalteK. Acoustic method for synchronization of magnetic resonance imaging (MRI). Acta Acust. United Acust. (2008) 94:148–55. 10.3813/AAA.918017
71.
KrugJWRoseGCliffordGDOsterJ. ECG-based gating in ultra high field cardiovascular magnetic resonance using an independent component analysis approach. J Cardiovasc Magn Reson. (2013) 15:104. 10.1186/1532-429X-15-104
72.
GregoryTSSchmidtEJZhangSHHoTse ZT. 3DQRS: a method to obtain reliable QRS complex detection within high field MRI using 12-lead electrocardiogram traces. Magn Reson Med. (2014) 71:1374–80. 10.1002/mrm.25078
73.
KoschOThielFDi ClementeFSHeinMSeifertF editors. Monitoring of human cardio-pulmonary activity by multi-channel UWB-radar. Antennas and Propagation in Wireless Communications (APWC), 2011 IEEE-APS Topical Conference on.Torino: IEEE (2011).
74.
GrossSVionnetLKasperLDietrichBEPruessmannKP. Physiology recording with magnetic field probes for fMRI denoising. Neuroimage (in press). [Epub ahead of print]. 10.1016/j.neuroimage.2017.01.022
75.
MaclarenJAksoyMBammerR. Contact-free physiological monitoring using a markerless optical system. Magn Reson Med. (2015) 74:571–7. 10.1002/mrm.25781
76.
FuchsKHezelFWinterLOezerdemCGraesslADieringerMet al. editors. Feasibility of cardiac fast spin echo imaging at 7.0 T using a two-dimensional 16 channel array of bowtie transceivers. In: Proceedings of the International Society for Magnetic Resonance Medicine.Salt Lake City (2013).
77.
NorrisDGBoernertPReeseTLeibfritzD. On the application of ultra-fast rare experiments. Magn Reson Med. (1992) 27:142–64. 10.1002/mrm.1910270114
78.
CarrHYPurcellEM. Effects of diffusion on free precession in nuclear magnetic resonance experiments. Phys Rev. (1954) 94:630–8. 10.1103/PhysRev.94.630
79.
MeiboomSGillD. Modified spin-echo method for measuring nuclear relaxation times. Rev Sci Instrum. (1958) 29:688–91. 10.1063/1.1716296
80.
HeinrichsUUttingJFFrauenrathTHezelFKrombachGAHodeniusMAet al. Myocardial T2* mapping free of distortion using susceptibility-weighted fast spin-echo imaging: a feasibility study at 1.5 T and 3.0 T. Magn Reson Med. (2009) 62:822–8. 10.1002/mrm.22054
81.
SchickF. SPLICE: sub-second diffusion-sensitive MR imaging using a modified fast spin-echo acquisition mode. Magn Reson Med. (1997) 38:638–44. 10.1002/mrm.1910380418
82.
HaackeEBrownRThompsonMVenkatesanRBernsteinMKingKet al. Magnetic Resonance Imaging - Physical Principles and Sequence Design. New York, NY: John Wiley and Sons (1999).
83.
RobinsonSSchodlHTrattnigS. A method for unwrapping highly wrapped multi-echo phase images at very high field: UMPIRE. Magn Reson Med. (2013) 72:80–92. 10.1002/mrm.24897
84.
HuelnhagenTHezelFSerradas DuarteTPohlmannAOezerdemCFlemmingBet al. Myocardial effective transverse relaxation time T2* correlates with left ventricular wall thickness: a 7.0 T MRI study. Magn Reson Med. (2016) 77:2381–9. 10.1002/mrm.26312
85.
JafferFAWenHBalabanRSWolffSD. A method to improve the BO homogeneity of the heart in vivo. Magn Reson Med. (1996) 36:375–83. 10.1002/mrm.1910360308
86.
ScharMKozerkeSFischerSEBoesigerP. Cardiac SSFP imaging at 3 Tesla. Magn Reson Med. (2004) 51:799–806. 10.1002/mrm.20024
87.
ScharMVonkenEJStuberM. Simultaneous B(0)- and B(1)+-map acquisition for fast localized shim, frequency, and RF power determination in the heart at 3 T. Magn Reson Med. (2010) 63:419–26. 10.1002/mrm.22234
88.
ReederSBFaraneshAZBoxermannJLMcVeighER. In vivo measurement of and field inhomogeneity maps in the human heart at 1.5 T. Magn Reson Med. (1998) 39:988–98. 10.1002/mrm.1910390617
89.
ShahSKellmanPGreiserAWealePJZuehlsdorffSJerecicR editors. Rapid fieldmap estimation for cardiac shimming. In: Proceedings of the International Society for Magnetic Resonance Medicine 17. Honolulu, HI (2009).
90.
Serradas DuarteTHuelnhagenTNiendorfT editors. Assessment of myocardial B0 over the cardiac cycle at 7.0T: implications for susceptibility-based cardiac MR techniques. In: Proceedings of the International Society for Magnetic Resonance Medicine 24. Singapore (2016).
91.
Serradas DuarteT. Detailing Myocardial B0 across the Cardiac Cycle at UHF: B0 Assessment and Implications for Susceptibility-based CMR Techniques. Lisbon: Universidade Nova de Lisboa (2016).
92.
KelleyCT. Iterative Methods for Optimization. Philadelphia: Siam (1999).
93.
SandinoCMKellmanPAraiAEHansenMSXueH. Myocardial T2* mapping: influence of noise on accuracy and precision. J Cardiovasc Magn Reson. (2015) 17:7. 10.1186/s12968-015-0115-3
94.
ManjonJVCoupePMarti-BonmatiLCollinsDLRoblesM. Adaptive non-local means denoising of MR images with spatially varying noise levels. J Magn Reson Imaging (2010) 31:192–203. 10.1002/jmri.22003
95.
FengYHeTFengMCarpenterJPGreiserAXinXet al. Improved pixel-by-pixel MRI R2* relaxometry by nonlocal means. Magn Reson Med. (2014) 72:260–8. 10.1002/mrm.24914
96.
HuelnhagenTPohlmannANiendorfT editors. Improving T2* mapping accuracy by spatially adaptive non-local means noise filtering. In: Proceedings of the International Society for Magnetic Resonance Medicine 23. Toronto, ON (2015).
97.
GuytonACHallJE editors. Muscle blood flow and cardiac output during exercise; the coronary circulation and ischemic heart disease. In: Guyton and Hall Textbook of Medical Physiology. 10th ed. Philadelphia: Saunders (2000). p. 223–34.
98.
SchmidtRFLangFHeckmannM. Herzmechanik. Physiologie des Menschen: Mit Pathophysiologie. 31st ed. Heidelberg: Springer (2010). p. 539–64.
99.
DeussenA. Herzstoffwechsel und Kroronardurchblutung. In: Schmidt RF, Lang F, Heckmann M, editors. Physiologie des Menschen: Mit Pathophysiologie. 31st ed. Heidelberg: Springer (2010). p. 565–71.
100.
JuddRMLevyBI. Effects of barium-induced cardiac contraction on large-and small-vessel intramyocardial blood volume. Circ Res. (1991) 68:217–25. 10.1161/01.RES.68.1.217
101.
FleckensteinJLCanbyRCParkeyRWPeshockRM. Acute effects of exercise on MR imaging of skeletal muscle in normal volunteers. AJR Am J Roentgenol. (1988) 151:231–7. 10.2214/ajr.151.2.231
102.
SchmidAISchewzowKFiedlerGBGoluchSLaistlerEWolztMet al. Exercising calf muscle T2* changes correlate with pH, PCr recovery and maximum oxidative phosphorylation. NMR Biomed. (2014) 27:553–60. 10.1002/nbm.3092
103.
MozaffarianDBenjaminEJGoASArnettDKBlahaMJCushmanMet al. Heart disease and stroke statistics—2016 update a report from the american heart association. Circulation (2015) 133:e38–360. 10.1161/CIR.0000000000000350
104.
KramerCMAppelbaumEDesaiMYDesvigne-NickensPDiMarcoJPFriedrichMGet al. Hypertrophic Cardiomyopathy Registry: the rationale and design of an international, observational study of hypertrophic cardiomyopathy. Am Heart J. (2015) 170:223–30. 10.1016/j.ahj.2015.05.013
105.
MaronBJOlivottoISpiritoPCaseySABellonePGohmanTEet al. Epidemiology of hypertrophic cardiomyopathy–related death revisited in a large non-referral-based patient population. Circulation (2000) 102:858–64. 10.1161/01.CIR.102.8.858
106.
NoureldinRALiuSNacifMSJudgeDPHalushkaMKAbrahamTPet al. The diagnosis of hypertrophic cardiomyopathy by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2012) 14:17. 10.1186/1532-429X-14-17
107.
SteriotisAKSharmaS. Risk stratification in hypertrophic cardiomyopathy. Eur Cardiol Rev. (2015) 10:31–6. 10.15420/ecr.2015.10.01.31
108.
Abdel-AtyHCockerMStrohmOFilipchukNFriedrichMG. Abnormalities in T2-weighted cardiovascular magnetic resonance images of hypertrophic cardiomyopathy: regional distribution and relation to late gadolinium enhancement and severity of hypertrophy. J Magn Reson Imaging (2008) 28:242–5. 10.1002/jmri.21381
109.
JohanssonBMörnerSWaldenströmAStålP. Myocardial capillary supply is limited in hypertrophic cardiomyopathy: a morphological analysis. Int J Cardiol. (2008) 126:252–7. 10.1016/j.ijcard.2007.04.003
110.
CecchiFOlivottoIGistriRLorenzoniRChiriattiGCamiciPG. Coronary microvascular dysfunction and prognosis in hypertrophic cardiomyopathy. New Engl J. Med. (2003) 349:1027–35. 10.1056/NEJMoa025050
111.
NiendorfTGraesslAThalhammerCDieringerMAKrausOSantoroDet al. Progress and promises of human cardiac magnetic resonance at ultrahigh fields: a physics perspective. J Magn Reson. (2013) 229:208–22. 10.1016/j.jmr.2012.11.015
112.
NiendorfTSodicksonDKKrombachGASchulz-MengerJ. Toward cardiovascular MRI at 7 T: clinical needs, technical solutions and research promises. Eur Radiol. (2010) 20:2806–16. 10.1007/s00330-010-1902-8
113.
NiendorfTPaulKOezerdemCGraesslAKlixSHuelnhagenTet al. W(h)ither human cardiac and body magnetic resonance at ultrahigh fields? Technical advances, practical considerations, applications, and clinical opportunities. NMR Biomed. (2016) 29:1173–97. 10.1002/nbm.3268
114.
SchmitterSWuXUğurbilKde MoorteleV. Design of parallel transmission radiofrequency pulses robust against respiration in cardiac MRI at 7 Tesla. Magn Reson Med. (2015) 74:1291–305. 10.1002/mrm.25512
115.
SchmitterSMoellerSWuXAuerbachEJMetzgerGJVan de MoortelePFet al. Simultaneous multislice imaging in dynamic cardiac MRI at 7T using parallel transmission. Magn Reson Med. (2016) 77:1010–20. 10.1002/mrm.26180.
116.
ReichenbachJRHaackeEM. High-resolution BOLD venographic imaging: a window into brain function. NMR Biomed. (2001) 14:453–67. 10.1002/nbm.722
117.
ReeseTGWeisskoffRMSmithRNRosenBRDinsmoreREWedeenVJ. Imaging myocardial fiber architecture in vivo with magnetic resonance. Magn Reson Med. (1995) 34:786–91. 10.1002/mrm.1910340603
118.
WuMTTsengWYSuMYLiuCPChiouKRWedeenVJet al. Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infarction: correlation with viability and wall motion. Circulation (2006) 114:1036–45. 10.1161/CIRCULATIONAHA.105.545863
119.
WagenhausBPohlmannADieringerMAElsAWaicziesHWaicziesSet al. Functional and morphological cardiac magnetic resonance imaging of mice using a cryogenic quadrature radiofrequency coil. PLoS ONE (2012) 7:e42383. 10.1371/journal.pone.0042383
120.
NiendorfTPohlmannAReimannHMWaicziesHPeperEHuelnhagenTet al. Advancing cardiovascular, neurovascular, and renal magnetic resonance imaging in small rodents using cryogenic radiofrequency coil technology. Front Pharmacol. (2015) 6:255. 10.3389/fphar.2015.00255
121.
DharmakumarRArumanaJMTangRHarrisKZhangZLiD. Assessment of regional myocardial oxygenation changes in the presence of coronary artery stenosis with balanced SSFP imaging at 3.0 T: theory and experimental evaluation in canines. J Magn Reson Imaging (2008) 27:1037–45. 10.1002/jmri.21345
122.
PouloLAlonLDenizCHaefnerRSodicksonDStoeckelBet al. A 32-Channel parallel exciter/amplifier transmit system for 7T imaging. In: ISMRM, editor. Proceedings of the International Society for Magnetic Resonance Medicine 19. Montréal (2011). p. 1867.
123.
LattanziRSodicksonDK. Ideal current patterns yielding optimal signal-to-noise ratio and specific absorption rate in magnetic resonance imaging: computational methods and physical insights. Magn Reson Med. (2012) 68:286–304. 10.1002/mrm.23198
Summary
Keywords
magnetic resonance, MRI, ultrahigh field, magnetic susceptibility, MR technology, cardiac physiology, cardiovascular imaging, myocardial tissue characterization
Citation
Huelnhagen T, Paul K, Ku M-C, Serradas Duarte T and Niendorf T (2017) Myocardial T2* Mapping with Ultrahigh Field Magnetic Resonance: Physics and Frontier Applications. Front. Phys. 5:22. doi: 10.3389/fphy.2017.00022
Received
24 March 2017
Accepted
30 May 2017
Published
14 June 2017
Volume
5 - 2017
Edited by
Ewald Moser, Medical University of Vienna, Austria
Reviewed by
Bernhard Gruber, University Medical Center Utrecht, Netherlands; Albrecht Ingo Schmid, Medical University of Vienna, Austria
Updates
Copyright
© 2017 Huelnhagen, Paul, Ku, Serradas Duarte and Niendorf.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thoralf Niendorf thoralf.niendorf@mdc-berlin.de
This article was submitted to Biomedical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.