Donate Help Contact The AHA Sign In Home
American Heart Association
Circulation Research
Search: search_blue_button Advanced Search
Circulation Research. 2001;88:340-346

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bauer, W. R.
Right arrow Articles by Ertl, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bauer, W. R.
Right arrow Articles by Ertl, G.
Right arrowPubmed/NCBI databases
Medline Plus Health Information
*MRI Scans
Related Collections
Right arrow CT and MRI
Right arrow Coronary circulation
(Circulation Research. 2001;88:340.)
© 2001 American Heart Association, Inc.


Integrative Physiology

Fast High-Resolution Magnetic Resonance Imaging Demonstrates Fractality of Myocardial Perfusion in Microscopic Dimensions

Wolfgang R. Bauer, Karl-Heinz Hiller, Paolo Galuppo, Stefan Neubauer, Julian Köpke, Axel Haase, Christiane Waller, Georg Ertl

From the Medizinische Universitätsklinik, Universität Würzburg (W.R.B., K.-H.H., P.G., C.W., G.E.), Germany; Department of Cardiovascular Medicine (S.B.), John Redcliffe Hospital, Oxford University, United Kingdom; Institut für Röntgendiagnostik (J.K.), Klinikum Mannheim, Universität Heidelberg, Germany; and Physikalisches Institut (EPV), Universität Würzburg (A.H.), Germany.

Correspondence to Wolfgang R. Bauer, Medizinische Universitätsklinik Würzburg, Josef Schneider Strasse 2, 97080 Würzburg. E-mail w.bauer{at}medizin.uni-wuerzburg.de


*    Abstract
up arrowTop
*Abstract
down arrowIntroduction
down arrowMaterials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Abstract—The fractal nature of heterogeneity of myocardial blood flow and its implications for the healthy and diseased heart is not yet understood. The main hindrance for investigation of blood flow heterogeneity and its role in physiology and pathophysiology is that conventional methods for determination of myocardial perfusion have severe limitations concerning temporal and spatial resolution and invasiveness. In isolated rat hearts, we developed a nuclear magnetic resonance technique that does not depend on contrast agents and in which the apparent longitudinal relaxation time is made perfusion sensitive by selective preparation of the imaging slice. This perfusion-sensitive relaxation time is determined within 40 seconds as a map with a high spatial in-plane resolution of 140x140 µm2 and a thickness of 1.5 mm. Perfusion imaging was validated with the established microsphere technique. Additionally, the congruence between perfusion-sensitive T1 maps and first-pass perfusion imaging was demonstrated. As an application of high-resolution perfusion imaging, fractal analysis of the spatial distribution of perfusion was performed. We were able to demonstrate that the fractality of this distribution exists even in microscopic dimensions. Vasodilation by nitroglycerin modulated the fractal pattern of perfusion, and the decrease of the fractal dimension indicated a shift toward homogeneity. This implies that parameters of the fractal distribution depend on the microvascular tone rather than on anatomic preformations; ie, fractality is a functional characteristic of perfusion.


Key Words: magnetic resonance imaging • perfusion • microcirculation • fractality


*    Introduction
up arrowTop
up arrowAbstract
*Introduction
down arrowMaterials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
The understanding of the local regulation and coordination of myocardial blood flow is an important goal in cardiac and integrative physiology. It is well established that myocardial blood flow shows considerable spatial heterogeneity.1 2 Quantitative analysis demonstrated3 that after division of the heart into N sections, the relative dispersion RD=SD/mean value of flow of these sections follows a power law,

(1)
where D is denoted as the fractal dimension, which takes values of {approx}1.1 to 1.3. A gaussian distribution of regional flow would reveal a value of 1.5 and a homogeneous flow of 1. This implies that heterogeneity of blood flow is neither organized in a purely random way, nor is it uniform, and that the spatial correlation of regional flow is long ranging. The origin of this pattern of flow heterogeneity is not understood. Some theories relate this to the branching of the vessel tree, ie, anatomic preformations; others postulate that this pattern of heterogeneity is a result of adaptive coordination of perfusion between myocardial segments; ie, fractality is a functional characteristic of perfusion.

Although the understanding of myocardial flow heterogeneity is of paramount interest in cardiac and integrative physiology, progress in this field of research is hampered by the limitations of available methods. The "gold standard" method is the application of microspheres,3 the deposition of which in tissue is proportional to perfusion. However, this is an invasive procedure, the microspheres itself affect perfusion, the spatial resolution is limited to pieces of {approx}0.2 g, the preparation of tissue implies systematic errors, and there are only a few states of perfusion assessable because of the restricted number of labels of microspheres. Quantitative magnetic resonance techniques that depend on contrast agent determine perfusion from signal-time curves4 during the first pass of the contrast agent; however, repetitive measurements that are necessary to increase precision or to investigate different perfusion states are hampered by residual contrast agent.

All methods mentioned have a rather low spatial resolution; ie, at the present it is not clear whether the fractal principle of flow heterogeneity also holds in microscopic dimensions. Therefore, our aim was to develop and validate a fast high-resolution nuclear magnetic resonance (NMR) perfusion imaging technique that does not require contrast agents. With this technique, spatial heterogeneity of perfusion should become visible on high-resolution images, and our next aim was to perform quantitative analysis of this heterogeneity in pilot studies.


*    Materials and Methods
up arrowTop
up arrowAbstract
up arrowIntroduction
*Materials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Experimental Preparations
All experiments with animals were conducted in accordance with the European regulation on care and use of laboratory animals.

Isolated Heart Model
Isolated hearts of male Wistar rats (body weight 300 to 350 g) were perfused (Langendorff mode) with oxygenated Krebs-Henseleit buffer (37°C) in a perfusion system that was adapted to a NMR-microscopy system.5 Coronary flow was measured by an ultrasonic flowmeter (T106, Transonic Systems Inc), and left ventricular pressure was measured by a balloon that was inserted in the left ventricle and connected to a Statham P23XL pressure transducer (Gould Instruments). Left ventricular end-diastolic pressure (LVEDP) was adjusted to 5 mm Hg by the balloon volume at the beginning of the experiment.

Induction of Coronary Artery Stenosis
To obtain hearts with a stenotic coronary artery, rats were intubated under ether anesthesia, and after thoracotomy the heart was exposed. A suture including a probe of 300 µm was positioned 2 to 3 mm from the origin of the left coronary artery. After ligation the probe was removed, the heart replaced, and the thorax closed. Hearts were excised after 2 weeks and treated as described above.

NMR Imaging
Pulse Triggering
Because the timing and the mode of NMR pulse triggering by the left ventricular pressure were described in detail elsewhere,5 only essentials are repeated. NMR imaging was triggered in early diastole. Cine NMR imaging had revealed a maximum transversal displacement during systole of 1.5 to 2 mm in the short-axis view. Motion was absent in diastole; ie, NMR imaging triggered at the beginning of diastole was not affected by artifacts.

Noncontrast Perfusion Imaging
Non–contrast-agent techniques make the apparent T1 of tissue perfusion sensitive by different labeling of inflowing arterial spins and stationary spins in the imaging slice. To achieve this, we performed, in a vertical high-field system (Bruker AMX 500, 11.75 T), a slice-selective spin inversion6 (thickness 3 mm, rectangular shaped) in the short-axis view 3 mm below the valvular plane of the heart. This was followed by a series of 16 2-dimensional FLASH images7 8 (TE=1 ms, TR=2.5 ms, flip angle=3°, pixel size in plane=140x140 µm2) within this slice. The imaging slice had a thickness of 1.5 mm (gaussian shaped) and was located symmetrically within the inversion slice; ie, >95% of the gaussian-shaped imaging profile was covered symmetrically by the inversion slice. It was verified in static fluid phantoms that the thickness of the inversion slice was the minimum value for which T1 without slice-selective preparation and T1 with slice-selective preparation were identical. This procedure provided high perfusion sensitivity and excluded systematic errors due to insufficient inversion of spins in the imaging slice. Imaging was performed in the short-axis view, because epicardial flow is directed from base to apex; ie, myocardium in the selected slice is mainly perfused by equilibrium spins.

A T1 map was obtained by a single exponential fit for each pixel from the time course of the 16 signal intensities.7 Four acquisitions were performed to gain the final T1 map. The T1 value of a pixel was rejected when the average relative deviation of the fitting curve and the time course of the signal intensities within this pixel were >5%. Those pixels appeared black on the T1 map. The high-resolution T1 map was obtained in <40 seconds. In our setup, the magnetization of inflowing spins was that of thermomagnetic equilibrium; ie, the apparent T1 relaxation in the imaging slice was accelerated. This acceleration depended on perfusion and mixing of intra- and extravascular spins, which were predominantly located in the capillary region. Estimation of the intracapillary-extracapillary water exchange rate6 and taking into account the volume distribution of coronary vessel segments and the geometry of the experimental setup implied a linear dependence between perfusion and relaxation rate,6

(2)
where {lambda} is the tissue/perfusate partition coefficient of water.

Error Analysis and Accuracy of T1 Measurement
Spatial deviations of T1 in myocardium arise as a superposition of spatial heterogeneity of (1) morphological and physiological parameters, (2) the random error, and (3) systematic spatial deviations due to the NMR system. The latter 2 were estimated by measuring T1 and its SD for stationary buffer, which was added contrast agent (Gd–diethylenetriamine penta-acetic acid [DTPA]), to cover a range of T1=1 to T1=3.5 seconds. There was no systematic spatial deviation of T1. The relative spatial SDs of the T1 maps were <2%.

It was demonstrated previously in detail5 that almost no motion of the heart was present in diastole in which 2-dimensional FLASH imaging was performed. This implies that T1 measurements were not affected by cardiac motion.

Validation of Perfusion Imaging
Experimental Protocols
The partition coefficient {lambda} of Equation 2Up was determined from a group of 8 isolated, buffer-perfused hearts. After beating 40 minutes at a perfusion pressure of 100 mm Hg, hearts were removed and the wet weight of left and right ventricular wall was measured. Hearts were placed in an oven (55°C) for 48 hours, and the partition coefficient was determined according to the following:

In 2 groups, the NMR perfusion imaging technique was validated with the established microspheres technique. Maps of the perfusion-sensitive relaxation time T1 were obtained in the short-axis view at different perfusion levels, and microspheres were injected simultaneously. In one group (n=5), perfusion was varied by reduction of perfusion pressure (to 100, 80, and 70 mm Hg in 4 hearts and to 100, 80, and 50 mm Hg in 1 heart). To induce more heterogeneous variations of perfusion, we elevated in another group (n=12) left LVEDP from 5 to 60 and thereafter to 80 mm Hg at a constant coronary perfusion pressure of 100 mm Hg.

Local hypoperfusion was investigated in hearts (n=5) with a stenotic coronary artery. Because the spatial resolution of the microsphere technique is lower than the hypoperfused area in the imaging slice, we compared perfusion-sensitive T1 mapping with first-pass imaging after application of contrast agent. First, 3-dimensional NMR angiography5 documented that the affected coronary artery was not occluded. After T1 mapping without contrast agent, a bolus of 0.2 mL Gd-DTPA (1.5 mmol/l) was injected into the Windkessel of the perfusion system, which was located just above the aorta of the heart. A series of T1-weighted diastolically triggered FLASH images (TE=1 ms, TR=2.5 ms, spatial resolution=140x140 µm2) followed.

Perfusion Measurement by Colored Microspheres
The microsphere technique and its adaptation to the NMR microscope were recently described by us.9 10 The number of microspheres (Dye-Trak; Triton Technology Inc; mean diameter 15 µm) of 1 color was {approx}30 000. Microspheres were flushed into the aorta with Krebs-Henseleit buffer containing 0.002% Tween 80 at 0.5 mL per minute over 10 minutes. After the experiments, hearts were frozen in liquid nitrogen. Transverse slices 4 to 6 mm in thickness were cut from ventricular base to apex. Within the first slice, which definitely included the imaging slice, the right and left ventricular wall were prepared. In 6 hearts of the group in which LVEDP was elevated, the left ventricular wall was further subdivided into endocardial and epicardial layers (ratio of wall thickness, {approx}1:2). After digestion of tissue and solvating the dye from the microspheres, mathematical decomposition of the composite absorption spectrum (200 to 820 nm) revealed the absorbance of a single color. Perfusion in the left and right ventricular wall of the slice was determined from absorbance of the sample/absorbance of all heart samplesxcoronary flow (flow from flowmeter+0.5 mL/min) per sample weight.

Data Analysis
In those hearts in which non–contrast-agent perfusion imaging was correlated with the microsphere technique, the right and left ventricular walls were delineated manually in T1 maps of each perfusion level. In the 6 hearts in which endo- and epicardial perfusion was analyzed, the endo- and epicardial layers were defined in the same way as described in the microsphere section. The spatial mean T1 of these structures was obtained automatically from averaging the T1 values of each pixel within the demarcated area. Variations of perfusion were determined from the mean T1 values according to Equation 2Up, which relates perfusion to tissue volume, whereas it is normally related to the mass. Because specific weight of myocardium is {approx}1 g/cm3, NMR perfusion data could be directly converted. The variations of perfusion obtained from T1 measurements were compared with the corresponding data from microsphere measurements by linear regression analysis and Bland-Altman plots (Figure 2Down).



View larger version (30K):
[in this window]
[in a new window]
 
Figure 2. Regression analysis (top) and Bland-Altman plots (bottom) of variations of perfusion obtained from microspheres {Delta}Pms and T1 data {Delta}PT1. Data are presented for the group in which perfusion pressure varied and for the group in which LVEDP was elevated (m indicates slope; b, offset of the regression line).

In those hearts in which perfusion-sensitive T1 mapping and first-pass imaging were compared, a stripe of 2x20 pixels ranging from a normal to a hypoperfused region was analyzed. Signal-to-noise ratio was too low to analyze signal intensity–time curves in each single pixel. Therefore, the 40-pixel stripe was divided into 10 (2x2-pixel) squares in which signal intensity–time curves were fitted to a {Gamma}-variate function. From this function the perfusion index of each square was determined as the slope, as follows: peak intensity/time of peak intensity (arbitrary units/number of heartbeats). This perfusion index was correlated with the mean 1/T1 of each square.

Fractal Analysis of Perfusion Heterogeneity
Equation 2Up shows that the analysis of heterogeneity of perfusion according to relation 1 is equivalent with that of the relaxation rate 1/T1. For fractal analysis, we therefore considered 4 nonoverlapping squares containing 16x16 pixels that were inserted into T1 maps of left ventricular myocardium. In each square, the relative dispersion RD of 1/T1 was obtained for N=16x16=256 pixels and thereafter for N=128, 64, 32, 16, 8, and 4 blocks of neighboring pixels (aggregates). Linear regression analysis of logRD(N) versus logN followed (Figure 4Down). This procedure was performed for each square under various flow conditions. We investigated 5 hearts, in each of which coronary flow was passively altered by adjustment of perfusion pressure at 100, 80, and 50 mm Hg. To investigate the effect of active vasodynamics on the heterogeneity of perfusion, we investigated another group of 5 hearts before, during, and after continuous infusion of nitroglycerin (0.5 mg/min) at constant perfusion pressure (100 mm Hg).



View larger version (52K):
[in this window]
[in a new window]
 
Figure 4. Fractal analysis; gray-scale T1 map of representative heart. A square of 16x16 pixels in the left ventricular wall is delineated and magnified (top right panel). The linear relation between relaxation rate 1/T1 and perfusion (Equation 2Up) implies a linear relation between the relative dispersion (RD) of both parameters. Hence, the fractal dimension D (Equation 1Up) can be determined from linear regression analysis of logRD1/T1 (N) vs logN, where N is the number of blocks of neighboring pixels (aggregates). Note that the power law of relation 1 holds up to microscopic dimensions (lower panel; r indicates correlation coefficient).

Statistics
Group data are expressed as mean±SD and were analyzed according to the Wilcoxon-Mann-Whitney U test or Scheffe F test when >2 groups were compared. Differences were considered as significant when P<0.05.


*    Results
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
*Results
down arrowDiscussion
down arrowReferences
 
Validation of Perfusion Imaging
The partition coefficient was determined as {lambda}=0.83±0.003, and there was no difference between left and right ventricular wall.

Variations of myocardial perfusion were directly reflected by myocardial T1 (Figure 1Down). When coronary flow was varied by reduction of perfusion pressure, microspheres measured initial perfusion of the imaged myocardium as 13±1.3 mL/g per minute, and almost equal decrements were obtained in the right and left ventricular wall. When variations of perfusion of the imaged region were obtained from T1 mapping according to Equation 2Up, a very close correlation was found with microsphere data (Figure 2Up), which was also present when left and right ventricular wall were analyzed separately as linear regression analysis demonstrated, as follows: r=0.96, slope=1.1, and offset=0.15 mL/g per minute for the left, and r=0.91, slope=0.92, and offset=-0.2 mL/g per minute for the right ventricular wall.



View larger version (80K):
[in this window]
[in a new window]
 
Figure 1. Color-coded T1 maps of isolated rat hearts after slice-selective spin inversion in short-axis view. Indicated are the left ventricle (LV) with the balloon inside; the right ventricle (RV); and cross sections of some coronary vessels (CV), which could be identified by their low T1 due to intraluminal perfusion. In the top image series, constant perfusion pressure (100 mm Hg) was maintained, and LVEDP was elevated by inflating the left ventricular balloon. Left ventricular developed pressure was 75, 50, and 45 mm Hg, and coronary flow was 15, 12, and 10 mL/min at the respective LVEDP levels. The increasing size of the left ventricular balloon compresses the right ventricle. The elevating LVEDP prolongs T1 predominantly in the left ventricular subendocardial layer, whereas the epicardial layer and right ventricular wall are less affected. This reflects the respective shifts of perfusion ({Delta}P), which are shown in the diagram at left. There is a close correlation between microsphere ({Delta}Pms) and NMR ({Delta}PT1) data. In the bottom image series, the heart was perfused at 3 different levels of coronary perfusion pressure (100, 80, and 50 mm Hg). Heart rate was 310 to 320 bpm, and left ventricular developed pressure was 75, 60, and 46 mm Hg at the respective pressure levels. T1 increases with decreasing coronary flow (CF). Note in both series the heterogeneous distribution of T1, which reflects that of perfusion.

In hearts in which LVEDP was elevated, microspheres demonstrated significant different decrements of perfusion in the right and left ventricular wall at each pressure step, 2.2±1.8 versus 3.3±1.3 mL/g per minute and 1.0±0.36 versus 1.6±0.7 mL/g per minute. Again, a close correlation was found between perfusion obtained from T1 mapping according to Equation 2Up and microsphere data (Figure 2Up). At LVEDP=5 mm Hg, perfusion was higher in the endocardial layer than in the epicardial one (endocardial-epicardial difference=2.3±0.3 mL/g per minute [NMR] and 2.5±0.2 mL/g per minute [microspheres]). The increments of LVEDP induced significant different decrements (Figure 1Up) of endo- and epicardial perfusion (endocardium/epicardium), as follows: 5->60 mm Hg, 3.8±0.8/2.5±0.4 mL/g per minute (NMR) and 4.0±0.9/2.7±0.5 mL/g per minute (microspheres), and 60->80 mm Hg, 3.5±0.6/1.8±0.3 mL/g per minute (NMR) and 3.4±0.5/1.9±0.6 mL/g per minute (microspheres). A close correlation was found between variations of perfusion determined by T1 mapping and microspheres, as follows: r=0.90, slope=0.94, and offset=0.2 mL/g per minute for the epicardium, and r=0.88, slope=0.90, and offset=0.13 mL/g per minute for the endocardium.

Evident congruence was found between the hypoperfused regions obtained from T1 mapping and first-pass imaging in the hearts with coronary artery stenosis (Figure 3Down). This was also confirmed by the similarity of the spatial dependence of 1/T1 and the perfusion index in the stripes (Figure 3Down), which is reflected in the close correlation of both parameters (mean correlation coefficient 0.92, range 0.89 to 0.95).



View larger version (38K):
[in this window]
[in a new window]
 
Figure 3. Perfusion-sensitive T1 map and first-pass perfusion maps (horizontal black lines in top right panels are orientation markers) of a heart with a stenotic left coronary artery. Three corresponding hypoperfused regions are indicated (top panels, 1 through 3). Even details of the T1 map in hypoperfused regions such as the arc-like structure in region 3 become visible by first-pass imaging. The stripe of 10 2x2-pixel squares was positioned along the x direction from normally to hypoperfused myocardium in region 3. The parallel courses of 1/T1, which is proportional to perfusion (Equation 2Up), and the perfusion index obtained from first-pass imaging are evident. Both parameters demonstrate a close correlation (lower right panel: m indicates slope; b, offset of the regression line; and r, correlation coefficient). Regarding the correlation plot, the data points (70, 0.805 s-1) (fourth square) and (70, 0.807 s-1) (sixth square) appear as 1 point.

Analysis of Heterogeneity of Perfusion
The perfusion-sensitive T1 maps appeared very heterogeneous (Figure 1Up), and the relative dispersion of the maps was in the range of 10% to 34%, which was far above the relative error of T1 determination of 2%. Hence, the heterogeneity of T1 reflected that of perfusion.

The correlation analysis of relative dispersion of perfusion versus sample size (logRD versus logN) demonstrated a linear logarithmic dependence of both parameters (Figures 4Up and 5Down) with correlation coefficients ranging between 0.90 and 0.98 (mean 0.96). In the hearts subjected to varying perfusion pressure, the mean of the fractal dimension D was 1.118±0.007 (100 mm Hg), 1.119±0.017 (80 mm Hg), and 1.120±0.016 (50 mm Hg).



View larger version (40K):
[in this window]
[in a new window]
 
Figure 5. Influence of nitroglycerin on heterogeneity of perfusion at constant perfusion pressure of 100 mm Hg. A 16x16-pixel 1/T1 map (top panel) of left ventricular myocardium is shown before, during, and after nitroglycerin infusion. Bottom, Graphs of corresponding fractal analysis within the map (see also Figure 4Up and text). 1/T1 maps demonstrate also that perfusion increases during nitroglycerin infusion, and it declines below baseline level after cessation of the drug (rebound effect). The fractal dimension reveals the opposite behavior; ie, heterogeneity of perfusion is modulated by nitroglycerin.

Nitroglycerin infusion increased coronary flow from 16.2±4.3 mL/min to19.0±4.1 mL/min, and after cessation of the drug, coronary flow declined to 15.3±3.9 mL/min. The corresponding time course of the mean fractal dimension revealed a decrease from 1.124±0.02 to 1.117±0.013 and thereafter a significant increase to 1.154±0.01 (Figure 5Up). The 1/T1 maps in Figure 5Up reveal that onset and cessation of nitroglycerin do not simply shift local perfusion to higher or lower values, but rearrange spatial heterogeneity of perfusion.


*    Discussion
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
up arrowResults
*Discussion
down arrowReferences
 
Development of Perfusion Imaging
There exist qualitative11 12 and quantitative13 14 15 16 17 non–contrast-agent perfusion imaging techniques. Most of the latter techniques work with active labeling of the inflowing spins. This requires spatial separation of tissue and its associated supplying artery; ie, this approach is suitable for the brain, the kidney, or the isolated perfused heart. However, in the intact heart the complex topography of coronary arteries, left ventricle, and myocardium hampers the active spin labeling. Another obstacle is that many sensitive parameters of spin labeling have to be considered.

In contrast to active spin labeling, we6 and others18 exploited the inflow of equilibrium spins to accelerate longitudinal relaxation by perfusion after selective inversion. This acceleration is due to inflow and subsequent transcapillary mixing of equilibrium spins with stationary spins. The dynamics of the latter are much faster than T1,6 which implies the linear dependence of relaxation rate on perfusion (Equation 2Up). In contrast to us, Reeders et al18 performed T1-weighted spin-echo imaging after selective inversion. A repetition time of {approx}12 seconds and 64 phase-encoding steps imply an acquisition time of 13 minutes, whereas we performed direct T1 mapping from 16 complete FLASH images during one relaxation course (snapshot FLASH); ie, 4 acquisitions provided a T1 map in <40 seconds. This technique combines high accuracy concerning determination of T1 with a high speed of the imaging procedure.

Potential delay of the inflow effect and, hence, underestimation of perfusion due to arterial spins in the inversion slice was minimized by short-axis imaging, because epicardial flow is directed from base to apex. The contribution of intramyocardial arteries in the inversion slice such as vessels directing blood from epi- to endocardium is negligible as well, given that their volume fraction is <1% of tissue volume. This implies that perfusate inside these arteries is exchanged in <200 ms for perfusion values >3 mL/g per minute. Because we determine relaxation time from an observation period of (16x1 image per heart cycle) >3 seconds, the initial delay is negligible.

The validation of perfusion imaging with the microsphere technique is hampered in part by the fact that the imaging slice is not quite identical with the slice prepared for microsphere measurements. This is due to obvious technical restrictions of the preparation and to the fact that the microsphere technique demands more myocardial mass9 than that of the imaging slice to provide reliable data. On the other hand, the close correlation of microsphere and NMR perfusion data suggests that the difference of the slices is of minor importance. Furthermore, the elevated endocardial perfusion obtained from NMR experiments under baseline conditions (LVEDP=5 mm Hg, perfusion pressure=100 mm Hg) is in agreement with our microsphere data and data from literature.19

It is obvious that the maximum spatial resolution on which the perfusion imaging technique could be validated was that of the microsphere technique, which we already reached in our preparations. The fact that our mathematical model of determining perfusion from T1 is independent from the sample size,6 and the accuracy of the T1 measurements on the scale of 140x140x1500 µm3 makes it evident that perfusion imaging is also valid in these dimensions. The congruence between the perfusion-sensitive T1 map and signal enhancement in the first-pass experiments in hearts with coronary artery stenosis confirm the validity of the non–contrast-agent technique.

In contrast to techniques with active spin labeling, the method of perfusion mapping by selective preparation of the imaging slice is not restricted to the isolated rat heart but can be adopted to the intact rat20 21 and with some modifications also to humans, as pilot studies demonstrate.22 Another advantage compared with active spin labeling is that our technique is independent from the T1 of perfusate; ie, the strong field dependence of this parameter is excluded. The T1 of myocardium in low field systems is in the range of 1 second, ie, only a factor of {approx}2 to 2.5 smaller than at 11.7 T. These facts imply that the field dependence of T1 should not hamper the transfer to low field systems. A real challenge for the transfer to the intact animal is that baseline perfusion (perfusion pressure {approx}100 mm Hg) is much lower than in the isolated saline perfused heart (in rats, 3.5 to 4 versus 12 to 15 mL/g per minute), which implies that variations of regional blood flow are much smaller as well. However, preliminary data in intact rats at 7 T with coronary artery stenosis demonstrate that under resting conditions hypoperfused regions can be detected with a spatial in-plane resolution of {approx}0.5x0.5 mm2.23 In humans, perfusion is even lower, and an increase of perfusion from baseline to maximum (0.75 versus 3.5 to 4 mL/g per minute) results in a 5% decrease of T1. This implies that the accuracy of T1 determination should be 2% in the segments under consideration to detect perfusion abnormalities.

Analysis of Heterogeneity of Perfusion
Our correlation analysis of relative dispersion versus sample size demonstrates that heterogeneity of perfusion can be described by the power law of Equation 1Up (Figures 4Up and 5Up). This implies that heterogeneity of perfusion reveals a fractal behavior in microscopic dimensions that has not been demonstrated yet. The origin of this fractality still remains a matter of discussion. The branching of the vessel tree, which itself follows a self-similar pattern,24 may be responsible as was mentioned by Van Beek et al.25 Our experiments demonstrate that this fractality is not a static characteristic of perfusion. Instead, the data obtained from hearts subjected to nitroglycerin infusion (Figure 5Up) demonstrate that vasodilation modulates the pattern of flow heterogeneity. After cessation of infusion, this pattern differs from the initial state. The fractal dimensions in these experiments demonstrate a shift toward a more uniform distribution of flow during nitroglycerin infusion. In contrast, when perfusion was varied by adjustment of perfusion pressure, almost no variations of the fractal dimension were observed. These observations speak in favor of the dependence of flow heterogeneity on the microvascular tone rather than on anatomic preformations.

The fractal analysis as indicated above and the physiological experiments we performed have to be considered as the first steps for our future goals. The way data are obtained from our method allows the investigation of more sophisticated problems and the application of complex mathematical methods, eg, multifractal analysis. Because the perfusion-sensitive T1 map is obtained in <40 seconds, we achieve a high temporal resolution that allows us to study the dynamics of heterogeneity in future projects. The perfusion imaging technique can be combined with metabolic26 and oxygenation22 sensitive imaging, which will provide further insight into the origin of heterogeneity of myocardial blood flow.

Conclusion
We have demonstrated that fast high-resolution NMR imaging of cardiac perfusion is feasible without contrast agent. The basis for this technique is the selective labeling of spins in the imaging slice, which made T1 perfusion sensitive. We have demonstrated that the fractal behavior of perfusion is maintained in microscopic dimensions and that this fractality is most likely a result of the microvascular tone.


*    Acknowledgments
 
This work was supported by Grant SFB 355 and Graduiertenkolleg "NMR" Grant HA 1232/8-1, Deutsche Forschungsgesellschaft, Forschungsfonds des Klinikum Mannheim.


*    Footnotes
 
Original received July 6, 2000; revision received December 14, 2000; accepted December 15, 2000.


*    References
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
up arrowResults
up arrowDiscussion
*References
 
1. Yipintsoi T, Dobbs WA, Scanlon PD, Knopp TJ, Bassingthwaighte JB. Regional distribution of diffusible tracers and carbonized microspheres in the left ventricle of isolated dog hearts. Circ Res. 1973;33:573–587.[Abstract/Free Full Text]

2. King RB, Bassingthwaighte JB, Hales JRS, Rowell LB. Stability of heterogeneity of myocardial flow in normal awake baboons. Circ Res. 1985;57:285–295.[Abstract/Free Full Text]

3. Bassingthwaighte JB, King RB, Roger SA. Fractal nature of myocardial blood flow heterogeneity. Circ Res. 1989;65:578–590.[Abstract/Free Full Text]

4. Wilke N, Jerosch-Herold M, Stillman AE, Kroll K, Tsekos N, Merkle H, Parrish T, Hu X, Wang Y, Bassinghtwaighte J, Bache RJ, Ugurbil K. Concepts of myocardial perfusion imaging in magnetic resonance imaging. Magn Reson Q. 1994;10:249–286.[Medline] [Order article via Infotrieve]

5. Bauer WR, Hiller K-H, Roder F, Neubauer S, Fuchs A, Grosse Boes C, Lutz R, Gaudron P, Hu K, Haase A, Ertl G. Investigation of coronary vessels in microscopic dimensions by two- and three-dimensional NMR-microscopic imaging in the isolated rat heart. Circulation. 1995;92:968–977.[Abstract/Free Full Text]

6. Bauer WR, Roder F, Hiller K-H, Han H, Fröhlich S, Rommel E, Haase A, Ertl G. The effect of perfusion on T1 after slice selective spin inversion in the isolated cardioplegic rat heart: measurement of a lower bound of intracapillary-extravascular water proton exchange rate. Magn Res Med. 1997;38:917–923.[Medline] [Order article via Infotrieve]

7. Deichmann R, Haase A. Quantification of T1 values by snapshot FLASH NMR imaging. J Magn Reson. 1992;96:608–612.

8. Nekolla S, Gneiting T, Syha J, Deichmann R, Haase A. T1 maps by K-space reduced snapshot-FLASH MRI. J Comput Assist Tomogr. 1992;16:327–332.[Medline] [Order article via Infotrieve]

9. Hiller K-H, Adami P, Voll S, Roder F, Kowallik P, Bauer WR, Haase A, Ertl G. In vivo colored microspheres in the isolated perfused rat heart for use in NMR. J Mol Cell Cardiol. 1996;28:571–577.[Medline] [Order article via Infotrieve]

10. Hiller K-H, Roder F, Adami P, Voll S, Kowallik P, Haase A, Ertl G, Bauer WR. Study of microcirculation by coloured microspheres and NMR-microscopy in isolated rat heart: effect of ischaemia, endothelin-1, and endothelin-1 antagonist BQ 610. J Mol Cell Cardiol. 1997;29:3115–3122.[Medline] [Order article via Infotrieve]

11. Prasad PV, Burstein D, Edelman RR. MRI evaluation of myocardial perfusion without a contrast agent using magnetization transfer. Magn Res Med. 1993;30:267–270.[Medline] [Order article via Infotrieve]

12. Edelman RR, Siewert B, Darby DG, Thangaraj V, Nobre AC, Mesulam MM, Warach S. Qualitative mapping of cerebral blood flow and functional localization with echo-planar MR imaging and signal targeting with alternating radio frequency. Radiology. 1994;192:513–520.[Abstract/Free Full Text]

13. Williams DS, Grandis DJ, Zhang W, Koretzky AP. Magnetic resonance imaging of perfusion in the isolated rat heart using spin inversion of arterial water. Magn Res Med. 1993;30:361–365.[Medline] [Order article via Infotrieve]

14. Detre JA, Leigh JS, Williams DS, Koretzky AP. Perfusion imaging. Magn Res Med. 1992;23:37–45.[Medline] [Order article via Infotrieve]

15. Zhang W, Williams DS, Detre JA, Koretzky AP. Measurement of brain perfusion by volume-localized NMR spectroscopy using inversion of arterial water spins: accounting for transit time and cross relaxation. Magn Res Med. 1992;25:362–371.[Medline] [Order article via Infotrieve]

16. Williams DS, Detre JA, Leigh JS, Koretzky AP. Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc Natl Acad Sci U S A. 1992;89:213–216.

17. Zhang W, Williams DS, Koretzky AP. Measurement of rat brain perfusion by NMR using spin labeling of arterial water: in vivo determination of the degree of spin labeling. Magn Res Med. 1993;29:416–421.[Medline] [Order article via Infotrieve]

18. Reeders SB, Atalay MK, McVeigh ER, Zerhouni EA, Forder JR. Quantitative cardiac perfusion: a noninvasive spin-labeling method that exploits coronary vessel geometry. Radiology. 1996;200:177–184.[Abstract/Free Full Text]

19. Anderson PG, Bishop SP, Digerness SB. Transmural progression of morphological changes during ischemic contracture and reperfusion in the normal and hypertrophied rat heart. Am J Pathol. 1987;129:152–167.[Abstract]

20. Belle V, Kahler E, Waller C, Rommel E, Voll S, Hiller K-H, Bauer WR, Haase A. In vivo quantitative mapping of cardiac perfusion in rats using a non invasive MR spin-labeling method. J Magn Res Imaging. 1998;6:1240–1245.

21. Waller C, Belle V, Rommel E, Hiller K-H, Voll S, Hu K, Schnackertz KD, Haase A, Ertl G, Bauer WR. Assessment of myocardial perfusion and intracapillary blood volume in intact rats at rest and coronary dilation: combination of a non contrast agent and contrast agent dependent NMR-imaging technique. Radiology. 2000;215:189–197.[Abstract/Free Full Text]

22. Wacker CM, Bock M, Hartlep AW, Beck G, van Kaick G, Ertl G, Bauer WR, Schad LR. Changes in myocardial oxygenation and perfusion under pharmacological stress with dipyridamol: assessment using T*2 and T1 measurements. Magn Reson Med. 1999;41:686–695.[Medline] [Order article via Infotrieve]

23. Waller C, Hiller K-H, Hu K, Nahrendorf M, Voll S, Ertl G, Bauer WR. Chronische Koronarstenosierung führt zur regionalen Mikrozirkulationsstörung: eine MR-Bildgebungsstudie am Rattenherzen. Z Kardiol. 2000;89(suppl 5):678. Abstract.

24. Tanaka A, Mori H, Tanaka A, Mohammed MU, Tanaka Y, Sekka T, Ito K, Shinozaki Y, Hyodo K, Ando M, Umetani K, Tanioka K, Kubota M, Abe S, Handa S, Nakazawa H. Branching patterns of intramural coronary vessels determined by microangiography using synchrotron radiation. Am J Physiol. 1999;276:H2262–H2267.[Abstract/Free Full Text]

25. Van Beek JHGM, Roger SA, Bassingthwaighte JB. Regional myocardial flow heterogeneity explained with fractal networks. Am J Physiol. 1989;257:H1670–H1680.[Abstract/Free Full Text]

26. von Kienlin M, Rosch C, Le Fur Y, Behr W, Roder F, Haase A, Horn M, Illing B, Hu K, Ertl G, Neubauer S. Three-dimensional 31P magnetic resonance spectroscopic imaging of regional high-energy phosphate metabolism in injured rat heart. Magn Reson Med. 1998;39:731–741.[Medline] [Order article via Infotrieve]




This article has been cited by other articles:


Home page
Cardiovasc ResHome page
A. R. Pries and T. W. Secomb
Origins of heterogeneity in tissue perfusion and metabolism
Cardiovasc Res, February 1, 2009; 81(2): 328 - 335.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. Marxen, J. G. Sled, L. X. Yu, C. Paget, and R. M. Henkelman
Comparing microsphere deposition and flow modeling in 3D vascular trees
Am J Physiol Heart Circ Physiol, November 1, 2006; 291(5): H2136 - H2141.
[Abstract] [Full Text] [PDF]


Home page
Physiol. Rev.Home page
N. Westerhof, C. Boer, R. R. Lamberts, and P. Sipkema
Cross-talk between cardiac muscle and coronary vasculature.
Physiol Rev, October 1, 2006; 86(4): 1263 - 1308.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. Nahrendorf, J. U. Streif, K.-H. Hiller, K. Hu, P. Nordbeck, O. Ritter, D. Sosnovik, L. Bauer, S. Neubauer, P. M. Jakob, et al.
Multimodal functional cardiac MRI in creatine kinase-deficient mice reveals subtle abnormalities in myocardial perfusion and mechanics
Am J Physiol Heart Circ Physiol, June 1, 2006; 290(6): H2516 - H2521.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
C. Waller, T. Engelhorn, K.-H. Hiller, G. Heusch, G. Ertl, W. R. Bauer, and R. Schulz
Impaired resting perfusion in viable myocardium distal to chronic coronary stenosis in rats
Am J Physiol Heart Circ Physiol, June 1, 2005; 288(6): H2588 - H2593.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
P. G. Carlier and D. Bertoldi
In vivo functional NMR imaging of resistance artery control
Am J Physiol Heart Circ Physiol, March 1, 2005; 288(3): H1028 - H1036.
[Abstract] [Full Text] [PDF]


Home page
Exp PhysiolHome page
R. S. Kirton, A. J. Taberner, P. M. F. Nielsen, A. A. Young, and D. S. Loiselle
Strain softening behaviour in nonviable rat right-ventricular trabeculae, in the presence and the absence of butanedione monoxime
Exp Physiol, September 1, 2004; 89(5): 593 - 604.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
T. Matsumoto, H. Tachibana, T. Asano, M. Takemoto, Y. Ogasawara, K. Umetani, and F. Kajiya
Pattern differences between distributions of microregional myocardial flows in crystalloid- and blood-perfused rat hearts
Am J Physiol Heart Circ Physiol, April 1, 2004; 286(4): H1331 - H1338.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. Nahrendorf, K.-H. Hiller, A. Greiser, S. Kohler, T. Neuberger, K. Hu, C. Waller, M. Albrecht, S. Neubauer, A. Haase, et al.
Chronic coronary artery stenosis induces impaired function of remote myocardium: MRI and spectroscopy study in rat
Am J Physiol Heart Circ Physiol, December 1, 2003; 285(6): H2712 - H2721.
[Abstract] [Full Text] [PDF]


Home page
JGPHome page
R. Karch, F. Neumann, B. K. Podesser, M. Neumann, P. Szawlowski, and W. Schreiner
Fractal Properties of Perfusion Heterogeneity in Optimized Arterial Trees: A Model Study
J. Gen. Physiol., August 25, 2003; 122(3): 307 - 322.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. Marxen and R. M. Henkelman
Branching tree model with fractal vascular resistance explains fractal perfusion heterogeneity
Am J Physiol Heart Circ Physiol, May 1, 2003; 284(5): H1848 - H1857.
[Abstract] [Full Text] [PDF]


Home page
PhysiologyHome page
U. K. M. Decking
Spatial Heterogeneity in the Heart: Recent Insights and Open Questions
Physiology, December 1, 2002; 17(6): 246 - 250.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
M. Nahrendorf, K.-H. Hiller, D. Theisen, K. Hu, C. Waller, R. Kaiser, A. Haase, G. Ertl, R. Brinkmann, and W. R. Bauer
Effect of Transmyocardial Laser Revascularization on Myocardial Perfusion and Left Ventricular Remodeling after Myocardial Infarction in Rats
Radiology, November 1, 2002; 225(2): 487 - 493.
[Abstract] [Full Text] [PDF]


Home page
JNMHome page
K. P. Schafers, T. J. Spinks, P. G. Camici, P. M. Bloomfield, C. G. Rhodes, M. P. Law, C. S.R. Baker, and O. Rimoldi
Absolute Quantification of Myocardial Blood Flow with H215O and 3-Dimensional PET: An Experimental Validation
J. Nucl. Med., August 1, 2002; 43(8): 1031 - 1040.
[Abstract] [Full Text] [PDF]


Home page
Circ. Res.Home page
R. S. Balaban and A. Arai
Function, Metabolic, and Flow Heterogeneity of the Heart : The View Is Getting Better
Circ. Res., February 16, 2001; 88(3): 265 - 267.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bauer, W. R.
Right arrow Articles by Ertl, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bauer, W. R.
Right arrow Articles by Ertl, G.
Right arrowPubmed/NCBI databases
Medline Plus Health Information
*MRI Scans
Related Collections
Right arrow CT and MRI
Right arrow Coronary circulation