The current study presents a new protocol for local pulse wave velocity (PWV) measurement using dynamic MR sequences, which have a high temporal resolution (TR < 6 ms). MR images were obtained at two positions along the common carotid artery, separated by a distance of 5 cm. In each phase of a MR series, carotid region was automatically extracted and then its area distension waveform could be obtained. Sixteen volunteers with no symptoms of cardiovascular diseases were studied. For local PWV estimation, three delay estimation principles were tested and produced the following values: intersecting tangents method (M1): 4.72 ± 1.40 m/s, second derivative method (M2): 4.94 ± 1.68 m/s and cross-correlation method (M3): 5.03 ± 1.17 m/s. The cross-correlation method showed a relative high reliability as its least standard deviation.
Heartbeat generates pressure waves traveling along the aorta and main pulmonary artery to the peripheral branches. The velocity of a pressure wave, which is called the pulse wave velocity (PWV), is considered to be an important surrogate marker for biomechanical properties of vessels [
In recent years, different techniques have been developed to assess the PWV noninvasively. It is commonly estimated by dividing a known distance between two measurement sites by the propagation time of the pressure or flow wave. Usually, the PWV is measured over a long pathway, often the carotid-femoral arterial trajectory.
Commercially available PWV measuring devices can directly record pulse waves by means of transducers, such as mechanotransducers or applanation tonometers, calculate the PWV, which is a global PWV value, and provide a general arterial stiffness evaluation [
The global PWV that we have described is an average PWV over a long segment composed of arteries with different mechanical characteristics [
On the other side, magnetic resonance imaging (MRI) has also been proposed to determine PWV. It has the advantage of determining the central aortic PWV globally and eliminating errors of vessel segment length measurements [
To our knowledge, few literatures describe a “gold standard” method for local PWV assessment.
In the present study, we propose a new method to derive the local carotid PWV using gated dynamic MRI sequences with a high temporal resolution (<6 ms) for normal and healthy persons. In contrast to other existent methods, cross-sectional area distension waveforms were analyzed. These waveforms are composed by about one hundred sample points, so that pulse propagation time can be calculated more accurately.
Sixteen volunteers with no symptoms of cardiovascular disease were studied (11 male and 5 female, mean age 26.1 ± 10.2 years; systolic blood pressure: 116.8 ± 11.1 mmHg; diastolic blood pressure: 81.6 ± 9.7 mmHg). This population with a low mean age allows us to study and approach the standard value of local PWV in carotid.
All acquisitions were performed on two 1.5T MRI units (Signa® HDxt 1.5T imaging system, GE Healthcare, Milwaukee, WI, USA), one in CIMA-UTC (Compiègne, France) and the other in Necker Hospital (Paris, France). Meanwhile, we used 16-channel head-neck-spine volumetric coils (GE Signa® HD Head-Neck-Spine array coils, GE Healthcare) as reception coils.
Images were acquired through a Fast Cardiac Gating Gradient Recalled Echo Sequence (FastCard GRE, GE Healthcare) which is an ECG-gated, 2D gradient echo sequence that acquires multiple phases of a cardiac cycle at a single slice (repetition time TR = 6 ms, echo time TE = 1.7 ms, matrix size = 256 × 256, field of view (FOV) = 16 cm, flip angle = 50˚, slice thickness (Tck) = 5 mm, views per segment (VPS) = 1, Arrhythmia Rejection Window (ARW) = 10%, trigger delay minimum = 10 ms). As a result of time of flight (TOF) effect, this sequence ensures an ideal contrast between flowing blood and static tissues. TR was so short that fluid T2 decay could be ignored [
During each phase of a cardiac cycle, a row of data space would be filled successively. The number of phases was calculated based on the heart rate, minimum TR, trigger window and views per segment. Data acquisition is launched after a trigger delay of 10 ms. This delay is applied similarly on both slices, so its effect on delay estimation was neglected. Irregular R-R intervals are rejected by ARW. An entire image (256 × 256) filled by 256 rows requires 256 heartbeats.
PWV could be determined by the following equation:
where L is the distance between two slices and Δt is the transit time.
In order to measure local PWV, images were obtained at two perpendicular slice positions along the common carotid artery, separated by a distance (L) of 4 to 5 cm. The superior slice locates near the bifurcation (
In each image of MR series, carotid region was extracted using a previously developed algorithm [
To estimate pulse wave velocity, different characteristics of pulse wave have been analyzed in other studies, such as arterial pressure [
As no standard value of local PWV exists, different transit time estimation methods have been used and compared.
A crucial point called “foot of waveform” is considered as the reference part of pulse wave. The transit time is determined by the delay between the two feet of the pulse waves of superior and inferior slices. But early wave reflection hinders the identification of the foot of waveform [
Two principal algorithms have been used to identify the foot of pulse waves: intersecting tangents method and second derivative maximum method [
1) Intersecting tangents method (denoted as M1)
In this method, the foot of the waveform is identified by the intersection of an upstroke line with a horizontal line through the minimal value of the pre-systolic portion. The upstroke line is defined as the tangent to the maximum slope point in systolic upstroke [
2) Second derivative maximum method (denoted as M2)
An assumption is made that the foot of the waveform corresponds the maximum value in artery cross-sectional area acceleration, which corresponds to the second derivative maximum of the distension waveform.
When the second derivative curve was obtained, we notice that it was very sensitive to interferences by the original one. Hence, a triangular moving average filter has been applied to smooth the first and the second derivative waveforms [
This method is also called cross-correlation method (denoted as M3). The two distension waveforms were normalized using the maximal value, then the inferior waveform has been moved horizontally, and an enclosed area between the two curves was determined. When the area is equal to 0, the horizontal displacement is considered as the transit time.
In a first step, the amplitude of systolic parts has been normalized and scaled up into an interval between 0
and 1. Besides, phase index of waveforms were kept and not modified during the normalization. For the determination of correlation area, we only considered the normalized values in the interval between 0.2 and 0.4 in order to avoid the confounding influence of wave reflection (
To approach a standard value of local PWV in carotid artery for healthy persons, all mentioned algorithms were
applied to all subjects. Carotid arteries in both two sides were considered, and no sign shows that there exists a correlation between the PWV values acquired from different sides. As a result, they may be considered as the values measured from independent subjects.
The estimated PWV obtained from those methods were listed in
With the results in
Some errors may originate from algorithms. In M2, the shape of second derivative is too sensitive and has small variations. This prevents the determination of the characteristic point even with the application of a triangular moving average filter.
With the cross-correlation method, only the portions between 0.2 and 0.4 of normalized waveform have been considered. That permitted to eliminate irregular shapes at the beginning part of the waveforms, and to minimize reduce the influence of reflected waves [
Subject | M1 | M2 | M3 | Distance | TR | |||
---|---|---|---|---|---|---|---|---|
Right | Left | Right | Left | Right | Left | (mm) | (ms) | |
1 | 6.14 | 3.58 | 4.74 | 5.18 | 4.27 | 3.86 | 64.00 | 5.66 |
2 | 3.88 | - | 6.41 | - | 5.88 | - | 50.49 | 5.66 |
3 | 6.52 | 3.14 | 5.00 | 2.73 | 4.79 | 2.58 | 50.49 | 5.66 |
4 | 4.16 | 7.62 | 2.28 | 3.84 | 5.69 | 5.75 | 51.22 | 5.64 |
5 | 3.71 | 3.39 | 4.62 | 9.24 | 5.92 | 3.46 | 51.95 | 5.66 |
6 | 5.90 | 4.49 | 5.62 | 8.99 | 5.82 | 4.56 | 50.61 | 5.66 |
7 | 3.99 | - | 4.14 | - | 4.60 | - | 51.22 | 5.66 |
8 | - | 5.90 | - | 4.56 | - | 6.60 | 51.22 | 5.66 |
9 | 3.91 | 3.84 | 3.50 | 3.25 | 4.73 | 5.90 | 51.22 | 5.66 |
10 | 4.63 | 3.68 | 4.64 | 5.97 | 4.76 | 3.19 | 50.72 | 6.12 |
11 | 6.32 | 7.84 | 7.22 | 6.19 | 4.33 | 6.98 | 48.59 | 5.65 |
12 | 5.29 | 3.55 | 4.13 | 3.75 | 6.32 | 5.62 | 50.16 | 6.12 |
13 | 5.50 | 5.03 | 4.03 | 4.03 | 6.63 | 5.80 | 49.76 | 5.65 |
14 | 6.58 | 4.21 | 5.80 | 5.80 | 5.36 | 4.24 | 49.39 | 6.12 |
15 | 2.78 | - | 5.90 | - | 5.40 | - | 50.25 | 6.12 |
16 | 3.65 | 3.04 | 3.06 | 3.75 | 2.58 | 5.20 | 52.10 | 6.12 |
Mean PWV (m/s) | 4.72 | 4.94 | 5.03 | |||||
Coefficient of Variance | 30% | 34% | 23% |
In this part, the reproducibility of segmentation has been studied. For each subject, we have applied five times of segmentation process to extract the right carotid region (15 subjects in all) with different sizes of determination windows (used to frame the carotid). Then we observed their influences to the PWV. The results (mean value and standard deviation) are shown in
In this study, we present a new method to assess the local PWV in normal carotid artery using dynamic MRI. This method is totally noninvasive and is suited to studies of large population. Compared with other MRI based methods, distension waveforms obtained by the proposed technique has a much higher temporal resolution (6 ms). This allows us to calculate characteristic time-points (“foot”) more accurately. In addition, the variation of carotid during a complete cardiac cycle could be observed, with the FastCard gating method, although a tiny delay (10 ms) exists before each data acquisition period.
In this context, few attempts of local PWV measurement in carotid exist at present, and the present study provides seem to introduce a new alternative of “gold standard” in this field.
Hermeling et al. and Meinders et al. used multiple M-line ultrasound system to obtain artery diameter distension waveform simultaneously at multiple positions along the common carotid artery. The former got 14 waveforms spaced over 17 mm while the latter applied 16 echo lines spanning a total width of 15.86 mm. These two studies have found estimated PWV of 3.5 ± 0.7 m/s with 4 subjects [
We do not have a method to determine the true normal local carotid PWV, or a “gold standard”, therefore the absolute accuracy could not be evaluated. However, we have executed three different transit time calculating algorithms in order to test the feasibility of the local estimation and to have an average values range.
We presented a new protocol for local PWV estimation in normal common carotid artery, using dynamic MR sequence having high temporal resolution. Dynamic images were used to obtain two distension waveforms which illustrate the variation of cross-sectional area of the carotid artery.
With the obtained waveforms, we effectuated the delay estimation using three different algorithms.
Cross-correlation method (M3) showed an estimated PWV result of 5.03 ± 1.40 m/s (CV = 23%) with a relative high reliability, while the results of the other two methods, intersecting tangents (M1) and second derivative method (M2), are 4.72 ± 1.68 and 4.94 ± 1.17 m/s (CV = 30% and 34%) respectively.
It would be highly recommended to apply this new protocol in pathologic cases. This would determine the influence of different pathologies on vessel properties. In addition, detections and diagnostics of lesions or pathologic changes could be performed locally and more precisely.