Open Access

Regional in vivo transit time measurements of aortic pulse wave velocity in mice with high-field CMR at 17.6 Tesla

  • Marco Parczyk1Email author,
  • Volker Herold1,
  • Gert Klug2,
  • Wolfgang R Bauer2,
  • Eberhard Rommel1 and
  • Peter M Jakob1
Journal of Cardiovascular Magnetic Resonance201012:72

https://doi.org/10.1186/1532-429X-12-72

Received: 12 August 2010

Accepted: 6 December 2010

Published: 6 December 2010

Abstract

Background

Transgenic mouse models are increasingly used to study the pathophysiology of human cardiovascular diseases. The aortic pulse wave velocity (PWV) is an indirect measure for vascular stiffness and a marker for cardiovascular risk.

Results

This study presents a cardiovascular magnetic resonance (CMR) transit time (TT) method that allows the determination of the PWV in the descending murine aorta by analyzing blood flow waveforms. Systolic flow pulses were recorded with a temporal resolution of 1 ms applying phase velocity encoding. In a first step, the CMR method was validated by pressure waveform measurements on a pulsatile elastic vessel phantom. In a second step, the CMR method was applied to measure PWVs in a group of five eight-month-old apolipoprotein E deficient (ApoE(-/-)) mice and an age matched group of four C57Bl/6J mice. The ApoE(-/-) group had a higher mean PWV (PWV = 3.0 ± 0.6 m/s) than the C57Bl/6J group (PWV = 2.4 ± 0.4 m/s). The difference was statistically significant (p = 0.014).

Conclusions

The findings of this study demonstrate that high field CMR is applicable to non-invasively determine and distinguish PWVs in the arterial system of healthy and diseased groups of mice.

Background

Cardiovascular diseases are among the most common causes of death in industrialized countries. Risk factors include increased age, male sex, diabetes, hypertension, and lipoprotein abnormalities. The aorta provides at least 60 to 70% of systemic compliance [1]. Reduced elasticity and compliance of the aorta are etiologic in cardiovascular diseases such as atherosclerosis and therefore serve as early indicators of asymptomatic atherosclerotic lesions [2, 3]. The velocity of pressure and flow pulses travelling down an elastic vessel termed the pulse wave velocity (PWV) increases with arterial stiffness [4]. The PWV is a direct measure of arterial stiffness [3] and serves as an independent predictor for cardiovascular risk and mortality [58] in many cases of CVD, including atherosclerosis [9]. In this study, a high field cardiovascular magnetic resonance (CMR) protocol using the transit time (TT) method [10] was developed and tested for its capability to distinguish groups of healthy and atherosclerotic mice by means of the PWVs in the descending aortas.

Originally, the transit time method has been used to determine the PWV in humans using several invasive methods [11, 12] and non-invasive methods such as ultrasound [13, 14] and CMR [15] before it was applied to smaller mammals such as mice. To this day new CMR methods have been refined by many workgroups [16, 17] and they constitute a comprehensive technique for the characterization of morphologic and functional arterial systemic parameters in humans.

In the past, the need to determine morphologic and functional parameters of the murine arterial system arose. Genetically engineered phenotypes of mice that are deficient in apolipoprotein E (ApoE(-/-) mice) spontaneously develop severe hyperlipidemia and atherosclerotic lesions at the arterial wall [18]. ApoE(-/-) mice develop all stages of lesions observed during atherogenesis that resemble lesions found in humans [19]. Therefore and for its short generation times, the ApoE(-/-) mouse model became an important model for human atherosclerosis [20].

Physical dimensions in mice are approximately 20 times smaller than in humans. With isoflurane anesthetized mice have high heart rates of 8 to 10 beats per second, whereas healthy humans at rest have heart rates of 0.8 to 1.7 beats per second. The necessities of very high spatial and temporal resolutions constitute the challenges for the determination of the PWV. In CMR of the corresponding animal model, the acquirable spatial resolution is limited by the signal-to-noise ratio (SNR). In the past, CMR systems did not provide a sufficient SNR to allow for accurate PWV measurements in mice. New high-field CMR systems facilitate a very high SNR and allow for examinations of morphologic and also functional parameters of the arterial system of mice [21, 22].

Alternatively to CMR, ultrasound methods have been utilized to non-invasively measure PWVs in mice [2325]. However, these ultrasound methods turned out to be highly observer dependent and in general were limited by acoustical windows in the thorax and by angular offsets in the alignment of the ultrasound transducer and the aorta [26]. The limited spatial resolution of ultrasound is a substantial impediment for morphological studies [27].

The primary objective of this study was to test the applicability of CMR and the TT method to assess the regional PWV of the descending murine aorta. The MR measurements were performed on a CMR system with a main magnetic field of 17.6 T. The accuracy and reproducibility of the CMR method were validated on a vessel phantom made of poly(vinyl alcohol) cryogel (PVA-C). In addition, it was tested whether the precision of the CMR method is sufficient to differentiate the mean PWVs of groups of atherosclerotic ApoE(-/-) and healthy wild type mice.

Methods

The Transit Time Method

During each systole the left ventricle of the heart ejects one stroke volume of blood into the aortic root, hence, generating a pulse beat. The pulse beat that in fact is a local increase of blood pressure and flow velocity propagates down the aorta with a velocity that is named the pulse wave velocity (PWV).

The principle to measure arterial PWV with the TT method is as follows. At least two measurement sites, I and II in Figure 1a, are to be selected on the aorta at a mutual distance, Δz. The transit time, Δt, which the pulse needs to travel from location I to II is to be measured (see Figure 1b). The PWV is calculated as PWV = Δz/Δt. This procedure is called the two-point TT method.
Figure 1

In vivo PWV measurements in the descending murine aorta. Figure 1a shows a reference scan of the murine aorta and the positioning of velocity measurement slices (I: thoracic aorta, II: abdominal aorta). Arrows indicate the regions of velocity measurements (two-point TT method). Figure 1b depicts representative flow velocity waveforms of an ApoE(-/-) mouse and the regression lines that determine the pulse onset time. Only the late diastole and the early systole were recorded. This particular measurement yields a PWV = Δz/Δt = 2.6 m/s.

The generalization of the two-point TT method is the multi-point TT method [2830]. It measures the transit times of the pulse waves in multiple locations along the path of wave propagation. The PWV is the constant of proportionality between the distances and the corresponding transit times.

Figure 1b and 2a reveal that flow and pressure pulse curves show a gradual rather than an instantaneous increase. However, a distinctive feature, the foot of the pulse wave, is necessary to determine the transit time. The foot of the pulse wave is defined as the intersection point of a regression line fitted to the velocity or pressure values before the pulse and another line fitted to the early increasing values.
Figure 2

Pressure wave measurements on the vessel phantom. Figure 2a shows a representative time course of the pressure wave recoded on the vessel phantom (sample rate = 100 ms-1). Figure 2b illustrates the transit of the pressure wave. In order to depict the PWV as the slope of the linear regression line, the independent variable, i.e., the pulse onset time is exceptionally plotted along the abscissa. The diagram shows one of six series of measurements that were performed at different installation heights (here +2 mm) of the vessel phantom relative to the MR iso-center.

Phantom Experiments

The accuracy of the MR method was tested on an elastic vessel phantom that was connected to a custom built flow and pressure pulse generator. The phantom allowed for the determination of PWVs by CMR measurements of flow waveforms and for reference measurements of pressure waveforms utilizing a pressure catheter that was inserted in the outlet end of the vessel (Figure 3). The PWVs obtained by the two different measurement methods were compared to each other in order to give evidence of the accuracy and precision of the MR method.
Figure 3

Simplified scheme and photograph of the pulsatile elastic vessel phantom. The phantom was used to investigate the accuracy and the precision of the developed MR TT method.

The phantom material had to have appropriate physical properties, e.g., mechanical strength, uniformity, long-term stability, and a modulus of elasticity resembling animal tissue. Poly(vinyl alcohol) cryogel (PVA-C), discovered independently by Peppas et al. and Nambu et al. [31, 32], has high breaking strengths and physiologic elastic moduli ranging from 0.1 to 1.0 MPa [33] and was, therefore, used to build the phantom vessel. Aqueous poly(vinyl alcohol) solution (15% w/w) was injected into an acrylic mold and crosslinked to form the gel by repeated freezing and thawing. The elastic modulus of the gel depends amongst others on the number of the freeze-thaw cycles, on the concentration of the aqueous solution, and essentially on the thaw rate [33, 34].

The vessel phantom had an inner diameter of 6 mm and a wall thickness of 0.25 mm. Its length was approximately 8 cm. The vessel was dimensioned according to Eq.[1], the Moens-Korteweg equation [35], to have a PWV in the physiological range of 2.0 to 6.5 m/s.
P W V = E h ρ d
(1)

E is Young's modulus of elasticity, h and d are the vessel wall thickness and the diameter, and ρ is the fluid density. The phantom was stored in a water bath to stabilize for at least 14 days prior to measurements because PVA-C is subject to an aging process that can affect its elasticity [34].

A homemade pressure and flow pulse generator was connected to the inlet end of the vessel phantom. A function generator (F34, Interstate Electronics Corporation, Anaheim, CA, USA) controlled the pulse generator. The function generator was set to generate pressure pulses every 1.8 s. Therewith, residual waves reflected at impedance mismatches inside the vessel phantom setup could decay completely between pulses. The function generator also triggered the data acquisitions of the MR system and of the setup used to record the pressure waves inside the vessel phantom.

The setup to record pressure waves comprised a pressure catheter, a signal amplifier, and a storage oscilloscope. The catheter was compounded of a blunted canula (TSK-Supra, Ebhardt-Söhne GmbH, Geislingen, Germany) and a Statham P23XL pressure transducer (Viggo-Spectramed, Inc., Oxnard, CA, USA). The signal amplifier was an MBS STAT ZAK (ZAK Psychologische und Physiologische Instrumente GmbH, Simbach, Germany). The storage oscilloscope was a TDS 3032 (Tektronix Inc., Beaverton, OR, USA).

The outlet end of the vessel was connected to a reservoir filled with an aqueous solution of 1.75 mM CuSO4. The height of the reservoir was adjusted in order that the vessel did not bloat or collapse. Two plastic rods supported the vessel on one side to prevent it from swinging in the transversal direction but allowed it to expand during the transit of the pulse.

The multi-point TT method was used to determine the PWV of the vessel phantom by pressure measurements. The measurement positions of the pressure catheter were determined at its protruding end with a calliper gauge. The pressure data were analyzed in MATLAB (The Mathworks, Inc., Natick, MA, USA) to determine the onset times of the pressure waves for each measurement position of the pressure catheter. The measurement positions were plotted versus the corresponding pressure pulse onset times. The slope of the linear regression yielded the PWV (see Figure 2b). The multi-point TT measurements were performed at six different installation heights of the vessel phantom that ranged from -13 mm to 12 mm from the MR iso-center. The PWV at the MR iso-center was determined by linear regression through the PWV and installation height values.

The MR system that was used for the phantom and in vivo measurements was a Bruker AVANCE 750 spectrometer (Bruker Biospin, Rheinstetten, Germany) with a vertical main magnetic field of 17.6 T and a bore size of 89 mm. The self-shielded gradient insert, a Bruker Micro 2.5, had a maximum gradient strength of 993 mT/m and an inner diameter of 40 mm. The homebuilt transversal electromagnetic mode radio frequency (RF) resonator had an accessible inner diameter of 25 mm.

The MR sequence applied phase velocity encoding to record the time courses of the flow velocities. The sequence was based on a two-dimensional FLASH sequence with incorporated velocity compensating gradients in all three gradient directions (Figure 4). To encode through-slice flow velocities, bipolar velocity encoding gradients were superposed onto the velocity compensating gradients in the slice direction. Three flow encoding steps were applied to each scan. The velocity-encoding window of the MR sequence was set to ± 0.20 m/s to accommodate the maximal flow velocities. The flow encoding gradients required the maximal available gradient power.
Figure 4

Scheme of the flow encoding MR method used to measure the PWV in the descending murine aorta. A fivefold repetition of the frames loop with interjacent incrementation of the variable delay produced an effective temporal resolution of 1000 time frames per second.

Transit times of the pulse waves are in the range of a few milliseconds, therefore, a temporal resolution of 1 ms was necessary to be able to obtain the PWV from the time courses of arterial blood flow velocities. The intrinsic temporal resolution of the MR sequence, equivalent to its repetition time, was 5 ms. Hence the time courses of the flow velocities had to be sampled in an interleaved fashion. The sequence was initiated five times, every time with an additional delay of 1 ms between the trigger signal and the initiation of the sequence. The resulting five data streams were interleaved in the post processing to generate a recording of flow velocities with a temporal resolution of 1 ms. The detailed timing and loop structure of the MR sequence are visualized in Figure 4.

The spatial in-plane resolution was 147 × 147 μm2 and the slice thickness was 1 mm. The field of view (FOV) was 22 × 22 mm2. The echo time was 1.6 ms. The Gauss RF excitation pulse had a length of 200 μs and an excitation angle of 20°. Signal averaging was not applied.

The two-point and multi-point TT methods were used to determine the PWV of the vessel phantom by MR measurements. A set of two-dimensional Fast Low Angle Shot (FLASH) experiments was leading the measurement protocol to localize the measurement positions on the vessel phantom. The MR imaging slices were positioned perpendicularly to the vessel at positions ±2.5 mm, ± 5.0 mm, ± 7.5 mm, and ± 10.0 mm from the magnets iso-center. A perpendicular slice orientation was crucial to avoid displacement artifacts caused by fluid flowing out of the imaging slice.

MR data were processed with a custom written routine using MATLAB. Velocity information was computed for pixels inside the vessel lumen by fitting a line to the phase values as a function of the first moments of the velocity encoding gradients (phase difference method). For every time frame the velocities were averaged over the luminal area.

64 pressure measurements were performed before 80 MR measurements. Differences in the mean PWV values obtained on the vessel phantom by the two measurement methods were tested by a two-sided t-test for unpaired data points for the hypothesis of equality against the alternative of differing values. The data points were not paired because many parameters, such as static pressure, temperature, and acoustic noise could not be held constant during the transition from pressure to MR measurements due to limited access and the strong magnetic field inside the magnet. The hypothesis was accepted for p-values ≥ 0.05. PWV values are given as mean ± standard error (SE).

In Vivo Experiments

The proposed MR method was designed to utilize the two-point TT method in vivo in the descending murine aorta. Nine MR measurements were performed on a group of five female eight-month-old ApoE(-/-) mice and eight measurements on a group of four age- and sex-matched C57Bl/6J mice. The ApoE(-/-) mice were fed a western type diet (TD 88137, Harlan Laboratories, Inc., Indianapolis, IN, USA) 10 weeks prior to MR measurements. During the MR examinations, the mice were anesthetized with an isoflurane inhalation (1.5 - 2.0 Vol.%) in O2 (2 L/min) applied by means of a nose cone. Mice were placed vertically (head up) in the RF resonator. Due to the small diameters of the gradient insert and the RF resonator, their body temperature could be kept constant at 37°C by adjusting the temperature of the gradient insert temperature control unit.

A pressure sensitive pneumatic balloon (Graseby Medical Limited, Watford, United Kingdom) was placed between the inner RF resonator wall and thoraces of mice to detect cardiac trigger and respiratory gating signals. Outside of the gradient insert, a pressure transducer (24PCEFA6 D, Honeywell S&C, Golden Valley, MN, USA) transformed the pressure signal from the balloon into an electrical signal that was amplified and processed in real-time by a homebuilt unit. Thus, electrical interference between the trigger signal and gradient fields oscillating in the same frequency domain was avoided.

All experimental procedures were in accordance with institutional and internationally recognized guidelines and were approved by the Regierung von Unterfranken (Government of Lower Franconia, Germany). The reference number of the permit of the animal experiments is 55.2-2531.01-19/07.

Two imaging slices were positioned perpendicularly to the thoracic and the abdominal aorta (as illustrated in Figure 1a). Reasonable in vivo experiment times only allowed for the two-point TT method to be applied. Slice separations were maximized in order to minimize the errors in the measured PWVs. Again, a perpendicular slice orientation was of crucial importance. To localize the descending aorta, a set of two-dimensional FLASH experiments was leading the measurement protocol.

The velocity-encoding window was set to ± 1.66 m/s to accommodate the maximal blood flow velocities. The RF excitation pulse had an excitation angle of 40°. The mice had a heart period of approximately 115 ms. A time window of 40 ms was sufficient to sample the late diastole and early systole. All other parameters of the measurement protocol were set as in the measurements on the vessel phantom. The total image acquisition times ranged from 15 to 20 min.

MR data were processed using the MATLAB routine that was also used for the analysis of the data acquired on the vessel phantom. The curve progression in the velocity-time diagrams showed sharp bends between the sections before and during the pulse. This allowed for a semi automatic selection of the fit ranges that defined the onset times of the flow pulse. The fit range before the pulse was selected manually. It included data points on an approximate horizontal line. The beginning of the pulse fit region was set automatically as a series of three consecutive data points with velocity values at least two standard deviations above the extrapolated fit line of the section before the pulse. The end point was selected manually as the last data point on the first straight portion of the flow pulse section. The program AMIRA (Visage Imaging, Inc., San Diego, CA, USA) was applied to gauge the distance between the two measurement locations. Therefore, straight-line segments were drawn along the luminal midline in a longitudinal reference scan of the aorta.

Isoflurane causes respiratory depression and, therefore, anesthetized mice develop a gasping breathing pattern with a period of approximately 1 s to 1.5 s. The duration of each respiratory movement is approximately 0.4 s. The MR data were acquired in the intermediate movement-free time intervals. From the periodicity of previous respirations, the homebuilt heart-triggering/breath-gating unit anticipates the next phase of respiratory motion and stops the MR data acquisition before the respiratory motion sets in. Data acquisition is continued approximately 0.4 s after the detected onset of breathing motion. Occasionally, the breathing pattern was not precisely periodic. Then data was acquired during the respiratory movement and motional artifacts in the MR images were the result. Sporadic non-periodic respiration was observed during all measurements, but showed no adverse motional artifacts in most of the measurements. In some of the motional artifacts, the signal amplitude and phase information of extra-luminal pixels were shifted into the vessel lumen. Those pixels were excluded from the calculation of the flow velocities when the anatomical structures that those pixels belonged to could be identified visually.

Differences in the mean PWV values of the two animal groups were tested by a two sided t-test for the hypothesis of equality against the alternative of a higher value for the ApoE(-/-) group. The hypothesis was rejected in favor of the alternative for p-values < 0.05. In this work all in vivo PWV values are given as mean ± standard deviation (SD).

Results

Phantom Experiments

On the vessel phantom, a reference PWV = 3.31 ± 0.18 m/s was obtained by the pressure catheter measurements. The multi-point MR measurements determined a PWV = 3.32 ± 0.18 m/s. The PWV values of the pressure and MR measurements are not different with statistical significance (p = 0.999). The parameters and results of the validation measurements are summarized in Table 1. The PWV values of the two-point MR measurements agree within their standard errors. A representative pressure wave in the vessel phantom and the multi-point TT method (explained in the method section) are shown in Figure 2.
Table 1

Parameters and results of the validation PWV measurements.

pressure TT measurements

number of data points

measurement positions

PWV

64

multi-point

3.31 ± 0.18 m/s*

flow velocity TT method

number of data points

measurement positions

PWV

80

multi-point

3.32 ± 0.18 m/s*

20

two-point at ± 10.0 mm

3.30 ± 0.20 m/s

20

two-point at ± 7.5 mm

3.39 ± 0.49 m/s

20

two-point at ± 5.0 mm

3.17 ± 0.29 m/s

20

two-point at ± 2.5 mm

3.44 ± 0.94 m/s

PWV values are given as mean ± standard error. *The mean values of the pressure and MR measurements do not show a statistically significant difference (p = 0.99).

In Vivo Experiments

Figure 1 shows representative time courses of flow velocities, which were recorded on an ApoE(-/-) mouse. The measured PWV values (Figure 5) ranged from 1.9 m/s in a C57Bl/6J mouse to 3.8 m/s in an ApoE(-/-) mouse. The average values of each animal group are PWV = 3.0 ± 0.6 m/s for ApoE(-/-) mice and PWV = 2.4 ± 0.4 m/s for C57Bl/6J mice. The mean value of the ApoE(-/-) group was higher than that of the C57Bl/6J group with statistical significance (p = 0.014). The change of the heart periods of examined mice was smaller than 5 ms during experiments. The experiment parameters and the results for both animal groups are summarized in Table 2.
Figure 5

PWV values in descending murine aortas. Figure 5a shows the individually measured PWVs of ApoE(-/-) and C57Bl/6J mice plotted versus the body masses. Vertically aligned datapoints belong to the same animal. Figure 5b shows a box-and-whiskers plot of the PWV values of the groups of eight-month-old female ApoE(-/-) and C57Bl/6J mice. The whiskers indicate the minimum and maximum values in each animal group. The boxes indicate the lower and upper quartiles; the band near the middle of each box represents the median. The ApoE(-/-) mice received a western type diet 10 weeks prior to the measurements. The PWV values are stated as mean ± standard deviation. The number of measurements on each group is indicated by n.

Table 2

Parameters and results of the in vivo MR PWV measurements.

parameter (mean ± standard deviation)

ApoE (-/-)group

C57Bl/6J group

age of animals (mo.)

8

8

number of animals

5

4

number of measurements

9

8

body mass (g)

28.4 ± 4.5

30.4 ±0.9

anesthetic

  

   isofluorane (%)

1.8 ± 0.2

1.7 ± 0.2

   in O2 (l/min)

1.5 ± 0.0

1.8 ± 0.4

avg. heart period (ms)

109 ± 9

112 ± 4

change in heart period of single experiments (ms)

2.0 ± 1.6

2.2 ± 1.5

distance between measurement locations (mm)

12.2 ± 1.1

9.5 ± 1.4

transit time between measurement locations (mm)

4.1 ± 0.9

4.2 ± 1.1

PWV* (m/s)

3.0 ± 0.6

2.4 ± 0.4

*The mean PWV of the ApoE(-/-) group is higher than the mean PWV of the C57Bl/6J group with statistical significance (p = 0.014).

Discussion

A pulsatile elastic vessel phantom with a physiological PWV was developed. PWV values of the vessel phantom were investigated by MR and pressure measurements. The mean PWV determined by the MR multi-point TT method showed no statistically significant difference to the reference value. The developed MR TT method is accurate. The agreement of the MR two-point TT measurements at the different measurement positions on the vessel phantom indicates that the MR method delivers reproducible PWV values. These measurements also indicate that the separation of the imaging slices shall be maximized in order to reduce the uncertainty in the measured PWV. A statement about the in vivo precision of the MR method in a best-case scenario cannot be made from the results of the phantom measurements, because the standard deviation of the phantom measurements is higher than that of the in vivo measurements. Different error mechanisms, such as residual transversal vessel wall oscillations [36], water droplets forming on the outside of the phantom vessel wall, and a jitter in the pressure pulses reduced the precision of the phantom measurements.

The results of this study show that the PWV is measurable in vivo in the descending murine aorta using a two-point TT CMR method and an CMR system with a main magnetic field strength of 17.6 T. The measurements are made possible by the high SNR intrinsic to high-field CMR and the large filling factor of the used RF resonator. The interleaved acquisition scheme and the pneumatic triggering technique (explained in the methods section) provide the necessary temporal resolution. Heart periods are sufficiently stable throughout the experiments, due to reliable anaesthesia and body temperature control. Flow compensation and the short echo time of the MR sequence greatly reduce motional artifacts. Additionally, the short echo time alleviates susceptibility artifacts, which are usually caused by boundaries of tissues with differing magnetic susceptibilities in high-field CMR.

Wave reflections, mainly generated by the micro-vascular bed, i.e. the arterioles, will change the shape of the flow wave and will affect the estimation of PWVs [37, 38]. It is critical for the determination of the PWV of the forward travelling pulse wave that only the reflection-free part of the wave is analyzed. Assuming a hypothetical PWV of 5 m/s and a distance of 20 mm between the abdominal measurement position and the microvascular bed, the reflected pulse wave needs 8 ms to return to the abdominal aorta. In this study, at most the early 7 ms of the flow velocity upstroke were used for the calculation of the PWV, therewith, wave reflections did not impose an impediment.

The PWV measured on the group of eight-month-old ApoE(-/-) mice was 25% higher than the PWV in the age-matched wild-type group. The in vivo experiments showed, that the precision of the measurement method suffices to distinguish between the PWVs of different animal groups. The difference is statistically significant (p = 0.014). The PWV values measured in mice in this study are lower than values previously stated in literature but agree in general with most publications [25, 39, 40]. The deviations between the findings of this study and literature might be caused by systematic deviations inherent to the different methods of measurement, as the different results of ultrasound and pressure measurements on alike animal groups by Hartley and Wang suggest. Deviations to the findings can also be caused by the differences in age of the animals, because different animal strains were used, or because the animals were positioned vertically in the MR system. The PWVs measured in vivo in this study highly deviate from values found in an MR study that analyzed the part of the pulse wave that contains wave reflections [38].

Wang et al. found no significant differences in PWV values between ApoE(-/-) and C57Bl/6J mice that were younger than 13 months [24]. In this study, a significant difference was observed at eight months of age already. In the current work, ApoE(-/-) mice received a different, a western type diet, which accelerates aortic wall remodelling and increases elastic destruction and thus the PWV [41].

To make a general remark about the field of application of the TT method, it has to be considered that the TT method delivers a regional PWV value, which is averaged over the propagation pathway in between the measurement locations. In early stages of atherosclerosis diminutive vascular lesions are scattered along the aorta, which affect the PWV locally. Therefore, interest of research is also on the local elastic properties of the murine aorta. Our workgroup successfully applied the QA method [17] to mice in another study [22] to determine local PWVs. The QA method uses the volume flow (Q) to cross-sectional area (A) relation during the early systole to calculate the PWV. The QA method is unrivalled in local estimation of the PWV. But since it employs cross-sectional vessel areas, it comprises an additional source of error. The TT method is more robust and requires less user interpretation because it operates without the determination of cross-sectional areas [42].

Conclusions

This study demonstrates the feasibility of the non-invasive determination of the pulse wave velocity by the transit time method and high-field CMR. The measured PWVs are accurate and reproducible in comparison to reference pressure measurements on a vessel phantom and sufficiently precise to distinguish between animal groups. Therefore, the measured PWV can be used as a marker to classify the state of arterial dysfunction in living mice. Due to the short in vivo acquisition times of 15 to 20 minutes, the presented MR method is applicable to enhance studies of morphologic parameters characterizing the murine arterial system [21, 43] by adding the PWV as a functional parameter to form a comprehensive set of parameters.

Declarations

Acknowledgements

This work was funded by the Deutsche Forschungsgemeinschaft in the scope of Sonderforschungsbereich 688 'Mechanisms and imaging of cardiovascular cell-cell interactions' (grant number: SFB 688).

Authors’ Affiliations

(1)
Julius-Maximilians-Universität Würzburg
(2)
Julius-Maximilians-Universität Würzburg, Medizinische Klinik und Poliklinik I

References

  1. Stergiopulos N, Segers P, Westerhof N: Use of pulse pressure method for estimating total arterial compliance in vivo. Am J Physiol. 1999, 276: H424-428.PubMedGoogle Scholar
  2. Arnett DK, Evans GW, Riley WA: Arterial stiffness: a new cardiovascular risk factor?. Am J Epidemiol. 1994, 140: 669-682.PubMedGoogle Scholar
  3. Laurent S, Cockcroft J, Van Bortel L, Boutouyrie P, Giannattasio C, Hayoz D, Pannier B, Vlachopoulos C, Wilkinson I, Struijker-Boudier H: Expert consensus document on arterial stiffness: methodological issues and clinical applications. Eur Heart J. 2006, 27: 2588-2605. 10.1093/eurheartj/ehl254.View ArticlePubMedGoogle Scholar
  4. Callaghan FJ, Geddes LA, Babbs CF, Bourland JD: Relationship between pulse-wave velocity and arterial elasticity. Med Biol Eng Comput. 1986, 24: 248-254. 10.1007/BF02441620.View ArticlePubMedGoogle Scholar
  5. Laurent S, Katsahian S, Fassot C, Tropeano AI, Gautier I, Laloux B, Boutouyrie P: Aortic stiffness is an independent predictor of fatal stroke in essential hypertension. Stroke. 2003, 34: 1203-1206. 10.1161/01.STR.0000065428.03209.64.View ArticlePubMedGoogle Scholar
  6. Blacher J, Asmar R, Djane S, London GM, Safar ME: Aortic pulse wave velocity as a marker of cardiovascular risk in hypertensive patients. Hypertension. 1999, 33: 1111-1117.View ArticlePubMedGoogle Scholar
  7. Blacher J, Guerin AP, Pannier B, Marchais SJ, Safar ME, London GM: Impact of aortic stiffness on survival in end-stage renal disease. Circulation. 1999, 99: 2434-2439.View ArticlePubMedGoogle Scholar
  8. Blacher J, Pannier B, Guerin AP, Marchais SJ, Safar ME, London GM: Carotid arterial stiffness as a predictor of cardiovascular and all-cause mortality in end-stage renal disease. Hypertension. 1998, 32: 570-574.View ArticlePubMedGoogle Scholar
  9. Sato M, Katsuki Y, Kanehiro H, Iimura M, Akada Y, Mizota M, Kunihiro Y: Effects of ethyl all-cis-5,8,11,14,17-icosapentaenoate on the physical properties of arterial walls in high cholesterol diet-fed rabbits. j Cardiovasc Pharmacol. 1993, 22: 1-9. 10.1097/00005344-199307000-00001.View ArticlePubMedGoogle Scholar
  10. Bramwell JC, Hill AV: Velocity of transmission of the pulse-wave and elasticity of arteries. The Lancet. 1922, 199: 891-892. 10.1016/S0140-6736(00)95580-6.View ArticleGoogle Scholar
  11. Latham RD, Westerhof N, Sipkema P, Rubal BJ, Reuderink P, Murgo JP: Regional wave travel and reflections along the human aorta: a study with six simultaneous micromanometric pressures. Circulation. 1985, 72: 1257-1269.View ArticlePubMedGoogle Scholar
  12. Kelly R, Hayward C, Avolio A, O'Rourke M: Noninvasive determination of age-related changes in the human arterial pulse. Circulation. 1989, 80: 1652-1659.View ArticlePubMedGoogle Scholar
  13. Benthin M, Dahl P, Ruzicka R, Lindstrom K: Calculation of pulse-wave velocity using cross correlation--effects of reflexes in the arterial tree. Ultrasound Med Biol. 1991, 17: 461-469. 10.1016/0301-5629(91)90182-V.View ArticlePubMedGoogle Scholar
  14. Lehmann ED, Hopkins KD, Gosling RG: Aortic compliance measurements using Doppler ultrasound: in vivo biochemical correlates. Ultrasound Med Biol. 1993, 19: 683-710. 10.1016/0301-5629(93)90087-5.View ArticlePubMedGoogle Scholar
  15. Mohiaddin RH, Firmin DN, Longmore DB: Age-related changes of human aortic flow wave velocity measured noninvasively by magnetic resonance imaging. J Appl Physiol. 1993, 74: 492-497.PubMedGoogle Scholar
  16. Kraft KA, Itskovich VV, Fei DY: Rapid measurement of aortic wave velocity: in vivo evaluation. Magn Reson Med. 2001, 46: 95-102. 10.1002/mrm.1164.View ArticlePubMedGoogle Scholar
  17. Vulliemoz S, Stergiopulos N, Meuli R: Estimation of local aortic elastic properties with MRI. Magn Reson Med. 2002, 47: 649-654. 10.1002/mrm.10100.View ArticlePubMedGoogle Scholar
  18. Reddick RL, Zhang SH, Maeda N: Atherosclerosis in mice lacking apo E. Evaluation of lesional development and progression. Arterioscler Thromb. 1994, 14: 141-147.View ArticlePubMedGoogle Scholar
  19. Nakashima Y, Plump AS, Raines EW, Breslow JL, Ross R: ApoE-deficient mice develop lesions of all phases of atherosclerosis throughout the arterial tree. Arterioscler Thromb. 1994, 14: 133-140.View ArticlePubMedGoogle Scholar
  20. Osada J, Joven J, Maeda N: The value of apolipoprotein E knockout mice for studying the effects of dietary fat and cholesterol on atherogenesis. Curr Opin Lipidol. 2000, 11: 25-29. 10.1097/00041433-200002000-00004.View ArticlePubMedGoogle Scholar
  21. Herold V, Wellen J, Ziener CH, Weber T, Hiller KH, Nordbeck P, Rommel E, Haase A, Bauer WR, Jakob PM, Sarkar SK: In vivo comparison of atherosclerotic plaque progression with vessel wall strain and blood flow velocity in apoE(-/-) mice with MR microscopy at 17.6 T. MAGMA. 2009, 22: 159-166. 10.1007/s10334-008-0160-0.View ArticlePubMedGoogle Scholar
  22. Herold V, Parczyk M, Morchel P, Ziener CH, Klug G, Bauer WR, Rommel E, Jakob PM: In vivo measurement of local aortic pulse-wave velocity in mice with MR microscopy at 17.6 Tesla. Magn Reson Med. 2009, 61: 1293-1299. 10.1002/mrm.21957.View ArticlePubMedGoogle Scholar
  23. Hartley CJ, Taffet GE, Michael LH, Pham TT, Entman ML: Noninvasive determination of pulse-wave velocity in mice. Am J Physiol. 1997, 273: H494-500.PubMedGoogle Scholar
  24. Wang YX: Cardiovascular functional phenotypes and pharmacological responses in apolipoprotein E deficient mice. Neurobiol Aging. 2005, 26: 309-316. 10.1016/j.neurobiolaging.2004.06.018.View ArticlePubMedGoogle Scholar
  25. Reddy AK, Li YH, Pham TT, Ochoa LN, Trevino MT, Hartley CJ, Michael LH, Entman ML, Taffet GE: Measurement of aortic input impedance in mice: effects of age on aortic stiffness. Am J Physiol Heart Circ Physiol. 2003, 285: H1464-1470.View ArticlePubMedGoogle Scholar
  26. Stadler RW, Taylor JA, Lees RS: Comparison of B-mode, M-mode and echo-tracking methods for measurement of the arterial distension waveform. Ultrasound Med Biol. 1997, 23: 879-887. 10.1016/S0301-5629(97)00074-4.View ArticlePubMedGoogle Scholar
  27. Williams R, Needles A, Cherin E, Zhou YQ, Henkelman RM, Adamson SL, Foster FS: Noninvasive ultrasonic measurement of regional and local pulse-wave velocity in mice. Ultrasound Med Biol. 2007, 33: 1368-1375. 10.1016/j.ultrasmedbio.2007.03.012.View ArticlePubMedGoogle Scholar
  28. Dumoulin CL, Doorly DJ, Caro CG: Quantitative measurement of velocity at multiple positions using comb excitation and Fourier velocity encoding. Magn Reson Med. 1993, 29: 44-52. 10.1002/mrm.1910290110.View ArticlePubMedGoogle Scholar
  29. Hardy C, Bolster B, McVeigh M, Adams W, Zerhouni E: One-Dimensional Velocity Technique for NMR Measurement of Aortic Distensibility. Magnetic Resonance in Medicine. 1994, 31: 513-520. 10.1002/mrm.1910310507.PubMed CentralView ArticlePubMedGoogle Scholar
  30. Fielden S, Fornwalt B, Jerosch-Herold M, Eisner R, Stillman A, Ohinski J: A New Method for the Determination of Aortic Pulse Wave Velocity Using Cross-Correlation on 2 D PCMR Velocity Data. Journal of Magnetic Resonance Imaging. 2008, 27: 1382-1397. 10.1002/jmri.21387.View ArticlePubMedGoogle Scholar
  31. Peppas N: Turbidimetric studies of aqueous poly(vinyl alcohol) solutions. Die Makromolekulare Chemie. 1975, 176: 3433-3440. 10.1002/macp.1975.021761125.View ArticleGoogle Scholar
  32. Nambu M: A hydrogel. Book A hydrogel (Editor ed.^eds.). 1982, City, 57/130543:Google Scholar
  33. Chu KC, Rutt BK: Polyvinyl alcohol cryogel: an ideal phantom material for MR studies of arterial flow and elasticity. Magn Reson Med. 1997, 37: 314-319. 10.1002/mrm.1910370230.View ArticlePubMedGoogle Scholar
  34. Surry KJ, Austin HJ, Fenster A, Peters TM: Poly(vinyl alcohol) cryogel phantoms for use in ultrasound and MR imaging. Phys Med Biol. 2004, 49: 5529-5546. 10.1088/0031-9155/49/24/009.View ArticlePubMedGoogle Scholar
  35. Zamir M: The Physics of Pulsatile Flow. 2000, New York, NY, USA: Springer-Verlag New York, IncView ArticleGoogle Scholar
  36. Dineley J, Meagher S, Poepping TL, McDicken WN, Hoskins PR: Design and characterisation of a wall motion phantom. Ultrasound Med Biol. 2006, 32: 1349-1357. 10.1016/j.ultrasmedbio.2006.05.033.View ArticlePubMedGoogle Scholar
  37. Vlachopoulos C, O'Rourke M: Genesis of the normal and abnormal arterial pulse. Curr Probl Cardiol. 2000, 25: 303-367. 10.1067/mcd.2000.104057.View ArticlePubMedGoogle Scholar
  38. Zhao X, Pratt R, J W: Quantification of Aortic Compliance in Mice Using Radial Phase Contrast MRI. Journal of Magnetic Resonance Imaging. 2009, 30: 286-291. 10.1002/jmri.21846.View ArticlePubMedGoogle Scholar
  39. Wang YX, Halks-Miller M, Vergona R, Sullivan ME, Fitch R, Mallari C, Martin-McNulty B, da Cunha V, Freay A, Rubanyi GM, Kauser K: Increased aortic stiffness assessed by pulse wave velocity in apolipoprotein E-deficient mice. Am J Physiol Heart Circ Physiol. 2000, 278: H428-434.PubMedGoogle Scholar
  40. Hartley CJ, Reddy AK, Madala S, Martin-McNulty B, Vergona R, Sullivan ME, Halks-Miller M, Taffet GE, Michael LH, Entman ML, Wang YX: Hemodynamic changes in apolipoprotein E-knockout mice. Am J Physiol Heart Circ Physiol. 2000, 279: H2326-2334.PubMedGoogle Scholar
  41. Zhang SH, Reddick RL, Burkey B, Maeda N: Diet-induced atherosclerosis in mice heterozygous and homozygous for apolipoprotein E gene disruption. J Clin Invest. 1994, 94: 937-945. 10.1172/JCI117460.PubMed CentralView ArticlePubMedGoogle Scholar
  42. Ibrahim EH, Johnson KR, Miller AB, Schaffer JM, White RD: Measuring aortic pulse wave velocity using high-field cardiovascular magnetic resonance: comparison of techniques. J Cardiovasc Magn Reson. 2010, 12: 26-38. 10.1186/1532-429X-12-26.PubMed CentralView ArticleGoogle Scholar
  43. Klug G, Kampf T, Ziener CH, Parczyk M, Bauer E, Herold V, Rommel E, Jakob PM, Bauer WR: Murine atherosclerotic plaque imaging with the USPIO Ferumoxtran-10. Front Biosci. 2009, 14: 2546-2552. 10.2741/3396.View ArticleGoogle Scholar

Copyright

© Parczyk et al; licensee BioMed Central Ltd. 2010

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.