- Open Access
Sequentially based analysis versus image based analysis of Intima Media Thickness in common carotid arteries studies - Do major IMT studies underestimate the true relations for cardio- and cerebrovascular risk?
Cardiovascular Ultrasound volume 6, Article number: 32 (2008)
Image-based B-mode ultrasound has gained popularity in major studies as a non-invasive method of measuring cardio- and cerebrovascular risk factors. However, none of the major studies appears to have paid sufficient attention to the variation in end diastolic wall process. By using sequentially based analyses (SBA) of Intima-Media Thickness (IMT), the general purpose of this study was to show that the current image based (ECG tracked) analysis (IBA) has some major variations and might underestimate the true relations for cardiovascular events and stroke for IMT measurement.
The study group consisted of 2500 healthy male subjects aged between 35 to 55 years. 4 sequences (300 images) were analyzed per subject. 750,000 images were analysed throughout the course of this study.
IBA showed significantly lower mean, maximal, and minimal values for IMT in CCA than for SBA. The correlation analysis between IBA and SBA with the cardio- and cerebrovascular risk factors showed a higher correlation of SBA for all risk factors. The Pearson coefficient was 0.81, p < 0.01, for SBA versus Framingham CHD risk level (FCRL) and 0.49, p = 0.01, for IBA versus FCRL.
IBA did not measure the true maximal values of the IMT in this study. Together with the correlation analysis, this indicates that IBA might underestimate the true relations for IMT and risk factors.
Over the last few years, numerous studies have linked carotid Intima Media Thickness (IMT) and IMT progression with prevalent symptomatic coronary artery disease (CAD) [1–5] and cerebrovascular disease [1, 4–7], as well as with cerebral white matter lesions[8, 9]. Increased IMT as identified by B-mode ultrasound of the extra cranial carotid arteries is a predictor of incident coronary events [1, 10, 11], stroke  and all-cause mortality .
O'Leary and colleagues  reported that the risk of incident stroke or myocardial infarction adjusted for age and sex increased among Cardiovascular Health Study participants. In the Atherosclerosis Risk in Communities study, the addition of IMT to a basic model including traditional risk factors increased the area under the receiver-operator curve for stroke, but this was only statistically significant when combined with markers of peripheral vascular disease . In a head-to-head comparison, CT-identified coronary artery calcium (CAC) discriminated CAD patients from controls better than IMT .
It seems logical to develop additional ultrasound methods to measure the risk of cardio- and cerebrovascular events and stroke more precisely.
The carotid artery seems to be a useful window for cardio- and cerebrovascular risk, with ultrasound able to provide unique information on morphological, hemodynamic and elastic properties. The majority of clinical wall motion studies have derived elasticity parameters from measurements of the wall motion at a single location in vessels. Vessels have included the common carotid artery [16, 17], femoral artery , brachial artery [17, 18], cerebral arteries  and aorta [20, 21].
Various ultrasound techniques have been used to detect and track vessel wall motion – manual as well as computer based measurements. Computational techniques have mainly been based on analysis of the B-mode greyscale images , M-mode  and analyses of the raw RF ultrasound data.
All methods – manual and automatic – have one main part in common; all methods use an end diastolic ECG tracked image . Is this truly the maximal value during the systolic and diastolic time period?
In heart rate regulation, baroreflex plays a very important role, and is commonly related to carotid baroreceptor function. Systemic changes of the blood vessel wall and heart contraction may affect carotid artery pulse wave phase characteristics: carotid artery dilates more or less, decreases or increases the impulsation from baroreceptors. Because of that, heart rate variability changes in the baroreflex band of the frequency spectrum. Murgo et al [22, 23] divide arterial systolic pressure pulsation into three classes, based on the timing and prominence of the secondary wave caused by reflections. Each class has a different restoring force for the IMT. Due to these variations in systolic movement, the point in time of the end diastolic maximal expanse of the IM layer differs from cycle to cycle. Based on this concept and knowledge about the variety of HR variation, a measurement of end diastolic IMT images seems to be very complicated due to the class of arterial pulse pressure. We think that the computer can solve this problem by calculating each image during the end diastolic period. The maximal value of this sequence should be independent from heart cycle regulations.
Most of the major studies about IMT in major journals[10, 24] have not dealt adequately with this fundamental problem of IMT measurement. Due to this, they have more measurement error (primary error) and thus the reported relations are likely an underestimation of the true underlying relation of IMT compared to cardio- and cerebrovascular risk factors.
The aim of this study was to evaluate the correlation of image based IMT measurement to sequentially based analysis and cardio- and cerebrovascular risk factors.
Materials and methods
The study group consists of 2500 healthy male subjects (HS) aged between 35 to 55 years (age: 46.1 ± 5.6 years). All subjects underwent a detailed interview, completed a questionnaire (26) to evaluate their physical activity and food habits, a clinical examination, an ultrasound examination of the carotid tree, and a blood test for the lipid profile.
All subjects were free of cardiovascular disease as determined by medical history. The characteristics of the study population are given in Table 1.
Subjects were advised to avoid caffeine intake for 8 h before the study and to maintain their usual lifestyle (diet, alcohol intake and exercise level) between studies. None of the subjects took any medications.
The research was carried out in accordance with the Declaration of Helsinki (1989) of the World Medical Association, and has been approved by the Freiburg University Research and Ethics Committee.
After at least 15 minutes rest in supine position, each patient was scanned by six experienced ultrasonographic observers.
We used a Toshiba SSA-380A, Japan, ultrasound scanner with high resolution 10 MHz transducer and an aperture of 32 mm. The carotid tree was scanned the transducer positioned longitudinal to the CCA. The lumen was maximised and the position changed until the typical double lines of the intima and adventitia boundaries of the arterial wall could be seen. Both, near and far, boundaries could be seen in the scan.
We stored one sequence for each subject for the lateral and for antero-lateral view for left and right CCA (4 sequences per subject) with the simultaneously recorded ECG at the top of the R-wave on digital disk (Sony DKR-700, Japan) (Figure 1).
Each sequence contains a time period of 3 seconds including at least 2 heart cycles. Each sequence was stored with 25 Hertz, which included 75 images per sequence. Additionally, all examinations were stored on videotape in order assist in the interpretation, if necessary.
After all we used 300 images per subject for the analysis (left and right artery, for every artery a lateral and for antero-lateral view, every sequence contains 75 images, 2 × 2 × 75 = 300). In the course of this study, 750,000 (2 × 2 × 75 × 2500 subjects) images for the sequence and 10,000 (2 × 2 × 2500 subjects) single images were analysed.
Analysis of the images
All measurements and analyses were done by one experienced reader blinded to all clinical information. The reader visually controlled stored sequences of each subject.
Image based analysis
For image based analysis (IBA), the reader chose the image from the digital videotape of the patient during the ultrasound examination with the R-wave in the ECG signal. This is the standardized procedure actually recommended by the recent consensus papers for IMT measurement.
This image was analysed for Image based Analyses of IMT measurement (Figure 1, dotted line). To understand this study, it is very important to emphasise that the image used for Image Based Analysis was part of the sequence later used for sequence analysis.
Sequentially based analysis
The first sequentially based analysis was presented by Selzer et al., in 2001. Based on this idea we used the same digital videotape sequence as described above (Figure 2) for the sequential based analysis (SBA). The same image chosen for the image based analysis was used as the starting point for the sequence analysis. The sequence array contains 75 images (SeqImage0–74). Because of this procedure, we reduced the variability in the measurement procedure due to rescanning, reselecting images and remeasuring processes to nearly zero.
For analysis of the IM layer we used a length of 1 cm from the beginning of the bulb backwards into the CCA like recommended by the recent consensus papers for IMT measurement .
We measured IMT at the far wall of the distal CCA and not at the near wall, because the detecting and trailing edge of the IM layers in the near wall lead to problems in the interpretation of the interfaces .
Results from left and right carotid artery were averaged.
Image based analysis
Mean Intima-Media Thickness (mIMT)
mIMT is defined as the mean value of the total IM layer in a chosen segment.
Maximal Intima-Media Thickness (maxIMT)
maxIMT is defined as the maximal value of the total IM layer in a chosen segment.
Minimal Intima-Media Thickness (minIMT)
minIMT is defined as the maximal value of the minimal total IM layer in a chosen segment.
Sequence based analysis
Maximal Mean Intima-Media Thickness (MmIMT)
MmIMT is defined as the maximal mean value of the total IM layer in a chosen sequence and segment.
Maximal Maximal Intima-Media Thickness (MmaxIMT)
maxIMT is defined as the maximal maximal value of the total IM layer in a chosen sequence and segment.
Minimal Intima-Media Thickness (MminIMT)
MminIMT is defined as the maximal minimal value of the total IM layer in a chosen sequence and segment.
Assessment of Individual Coronary Heart Disease (CHD) Risk Level
For all subjects, the individual CHD risk level was estimated based on the Framingham CHD score sheets presented by Wilson et al. . This scoring system was applied to subjects with an age ranging from 30 to 79 years. The details of the algorithm for determining the risk score are described in the original report. The higher the value of the score, the higher is the risk for a CHD.
All data were calculated using the software program SPSS (Software package for the social sciences) for Windows, Version 12. The arithmetic mean and the standard deviations (SD) were used for descriptive statistics.
Pearson's correlation coefficient was used to measure the strength of the linear relationship between IMT and WI (0.80–1.00: strong association, 0.60–0.79: strong-moderate association, 0.40–0.59: weak-moderate association)
A Bland & Altman analysis was performed to describe reproducibility of the method. The agreement for the reproducibility of WI was accepted if the upper and lower boundaries of the 95% confidence interval for WI covered zero.
Characteristics of subjects
A total of 2500 subjects were examined in the study, including a total number of 750,000 images. Anthropometric data, resting heart rate, and blood pressure for both groups are shown in Table 1.
Characteristics of IMT measurement
In the ultrasound analysis, IMT image based analysis showed significantly lower mean (IBA mean: 0.63 ± 0.18, SBA mean: 0.67 ± 0.19, p < 0.01), maximal (IBA max: 0.77 ± 0.19, SBA max: 0.83 ± 0.18, p < 0.01) and minimal (IBA min: 0.48 ± 0.15, SBA min: 0.51 ± 0.15, p < 0.01) values for the left, right, and combined CCA than for sequential based analysis (Table 2).
The Bland-Altman analysis revealed a difference of 0.04 mm for maximal IMT common carotid artery and a correlation coefficient of 0.86 with systematically lower values for imaged based analysis than for sequential tracing in mean and minimum like the maximum IMT in left, right, and common carotid artery (Table 3).
The correlation analysis between IBA and SBA with the cardio- and cerebrovascular risk factors showed a higher correlation of SBA for all risk factors. The Pearson coefficient was 0.81, p < 0.01 for SBA versus Framingham, and 0.49, p = 0.01 for IBA versus Framingham. The same results were found for the other risk factors as well (Figure 3, Table 4).
The main findings of this study are:
1) Sequential analysis could detect higher values in vivo situation of the maximal IM layers than image based IMT measurement (see Characteristics and Bland-Altman analysis).
2) Sequential analysis was significantly different from image based analysis
3) The results of correlation analyses may lead to the assumption that current IMT studies on image based analysis might underestimate the true relations of individual risk for atherosclerosis (see Correlation analysis).
First it is important to emphasize that this study used the same computer based detecting algorithm for the image based analysis as for the sequential based analysis. Using the same algorithm (25) reduces the errors of detecting the correct layers to a minimum.
The image used for the image analysis (IBA) was part of the sequence used for the sequential analysis. This is important to emphasize. Using this procedure – we used the same detection algorithm, the same reader, and the same video sequence for the analysis of the image based and sequence based analysis – we reduced the variation and reading errors for the IBA and the SBA analysis to a possible minimum. One variation left is the possibility that the differences seen in IMT measures between one base image and a sequence of images might be due to motion during the acquisition of the carotid artery images. The ultrasound observer tried very hard to get a stable scan sequence during the 3 seconds of data recording. Due to our ultrasound procedure we think that this possibility is very small.
The remaining difference can be assumed as true difference of both methods.
In this study we found a significant difference between IBA and SBA. The results of Bland-Altman analysis showed that both methods did not detect significantly the same values of IM layers. The image based analysis (IBA) measures a significantly lower IMT value than the sequence based analysis. What is the reason for the difference?
A detection error can be excluded, since we used an automatic system. The obvious and logical answer is that image based analysis did not measure the maximal value of an Intima Media layer during a heart cycle (Finding 1, 2).
However, in epidemiology, it is not the absolute value that is important, but how it ranks subjects as having high or low risk. For example, the CLAS study had an IMT value of 0.60 mm whereas this was a very high risk group. The ARIC study had mean levels of 0.70 mm in middle aged subjects with lower cardiovascular risk. However, the ability to rank subjects using the method is also very relevant. The differences between the results of both studies might be explained by the different detection methods and by the end diastolic time variations.
The CLAS and ARIC studies showed the possibility of IMT measurement as a parameter for to evaluate cardio- and cerebrovascular risk, but they also showed the lack of the actual IMT measurement protocols.
Is a person with a mean IMT of 0.6 mm a person with a higher risk for cardio- or cerebrovascular disease or not? To our knowledge there is no data available which can answer this question. All major studies differ significantly with regard to the baseline values of their subjects. The main reason for this difference might be the measurement method. The answer cannot be solved by the current protocols for image based analysis due to the large number of variations and end diastolic variation errors.
Using a large number of subjects, as was recently done in major IMT studies, can decrease the variation and end diastolic variation errors of other influences on the measurement. Due to the large number of subjects, the statements of the major studies are correct that IMT is a proven parameter of atherosclerotic processes in arteries.
We added a correlation analysis to our study to evaluate the correlation of IBA and SBA to the Framingham CHD risk level and other proven risk parameters. We found a weak to moderate correlation for IBA and a strong association for the SBA for the risk factors. Figure 3 shows that an increase in IMT is better correlated with a higher Framingham risk level for SBA than for IBA. When looking at all data together, there is a high indication that IBA analysis might underestimate the true relations for cardio- and cerebrovascular risk.
However, in this study, we hoped to show that IBA might underestimate the true relations. This statement has also been shown in histopathological studies [29–31], but the actual situation of image based analysis – manual or automatic-might need to be reconsidered (Finding 3).
In conclusion, this study shows that sequential analysis has a higher correlation to risk factors than image based analysis. We showed that the image based analysis might underestimate the true relations of cardio- and cerebrovascular risk.
We suggest using sequential analysis as a method for a robust measurement of IMT and cardiovascular risk.
Bots M, de Jong A, Grobbee DE: Left, right, near or far wall common carotid intima-media thickness measurements: associations with cardiovascular disease and lower extremity arterial atherosclerosis. J Clin Epidemiol. 1997, 50: 801-807. 10.1016/S0895-4356(97)00059-0
Craven TE, Ryu T, Espeland MA, Kahl F, McKinney W, Toole J, McMahan M, Thompson C, Heiss G, Crouse JR: Evaluation of the associations between carotid artery atherosclerosis and coronary artery stenosis: a case control study. Circulation. 1990, 82: 1230-1242.
Crouse JR, Tang R, Espeland MA, Terry J, Morgan W, Mercuri M: Associations of extracranial carotid atherosclerosis progression with coronary status and risk factors in patients with and without coronary artery disease. Circulation. 2002, 106: 2061-2066. 10.1161/01.CIR.0000033833.54884.34
O'Leary DH, Polak JF, Kronmal RA, Kittner S, Bond MG, Wolfson SK, Bommer W, Price T, Savage PJ: Distribution and correlates of sonographically detected carotid artery disease in the Cardiovascular Health Study. Stroke. 1992, 23: 1752-1760.
Chambless LE, Zhong MM, Arnett D, Folsom AR, Riley WA, Heiss G: Variability in B-mode ultrasound measurements in the atherosclerosis risk in communities (ARIC) study. Ultrasound Med Biol. 1996, 22: 545-554. 10.1016/0301-5629(96)00039-7
Burke GL, Evans GW, Riley WA, Sharrett AR, Howard G, Barnes RW, Rosamond W, Crow RS, Rautaharju PM, Heiss G: Arterial wall thickness is associated with prevalent cardiovascular disease in middle-aged adults. The Atherosclerosis Risk in Communities (ARIC) Study. Stroke. 1995, 26: 386-391.
Salonen R, Salonen JT: Progression of carotid atherosclerosis and its determinants: a population-based ultrasonography study. Atherosclerosis. 1990, 81: 33-40. 10.1016/0021-9150(90)90056-O
Bots M, van Swietjen M, Breteler M, de Jong A, Grobbee DE: Cerebral white matter lesions and atherosclerosis in the Rotterdam Study. Lancet. 1993, 341: 1232-1237. 10.1016/0140-6736(93)91144-B
Manolio T, Kronmal RA, Burke GL, Poirier V, O'Leary DH, Gardin M, Bryant MJ: Magnetic resonance abnormalities and cardiovascular disease in older adults. Stroke. 1994, 25: 318-327.
O'Leary DH, Polak JF, Kronmal RA, Manolio TA, Burke GL, Wolfson SK: Carotid-artery intima and media thickness as a risk factor for myocardial infarction and stroke in older adults. N Engl J Med. 1999, 340: 14-22. 10.1056/NEJM199901073400103
Hodis HN, Mack WJ, LaBree L, Selzer RH, Liu CR, Liu CH, Azen SP: The role of carotid arterial intima-media thickness in predicting clinical coronary events. Ann Intern Med. 1998, 128: 262-269.
Fried L, Kronmal RA, Newman A, Bild D, Mittelmark M, Polak JF, Gardin M: Risk factors for 5-year mortality in older adults. The Cardiovascular Health Study. J Am Med Assoc. 1998, 279: 585-592. 10.1001/jama.279.8.585. 10.1001/jama.279.8.585
Chambless LE, Heiss G, Sharatt A, Earp MW, Toole J: Prediction of ischemic stroke risk in the Atherosclerosis Riskin Communities Study. Am J Epidemiol. 2004, 160: 259-269. 10.1093/aje/kwh189
Terry J, Carr J, Tang R, Evans GW, Kouba W, Espeland MA, Mercuri M: Coronary artery calcium outperforms carotid artery intima-media thickness as a noninvasive index of prevalent coronary artery stenosis. Arterioscler Thromb Vasc Biol. 2005, 25: 1723-1728. 10.1161/01.ATV.0000173418.42264.19
Nagai Y, Matsumoto M, Metter EJ: The carotid artery as a noninvasive window for cardiovascular risk in apparently healthy individuals. Ultrasound Med Biol. 2002, 28: 1231-1238. 10.1016/S0301-5629(02)00578-1
Balkestein EJ, Staessen JA, Wang JG, van Der Heijden-Spek JJ, VanBortel M, Struijker-Boudier HA: Carotid and femoral artery stiffness in relation to three candidate genes in a white population. Hypertension. 2001, 38: 1190-1197. 10.1161/hy1101.095992
Kanters SDJM, Algra A, van Leeuwen MS, Banga JD: Reproducibility of In Vivo Carotid Intima-Media Thickness Measurements : A Review. Stroke. 1997, 28: 665-671.
Beux S, Carmassi F, Salvetti MV, Ghiadoni L, , Salvett A: Automatic evaluation of arterial diameter variation from vascular echographic images. Ultrasound Med Biol. 2001, 27: 1621-1629. 10.1016/S0301-5629(01)00450-1
Hayashi K, Handa H, Nagasawa S, Okumura A, Moritake K: Stiffness and elastic behavior of human intracranial and extracranial arteries. J Biomech. 1980, 13: 175-181. 10.1016/0021-9290(80)90191-8
Lehmann ED: Elastic properties of the aorta. Lancet. 1993, 342: 1417- 10.1016/0140-6736(93)92772-L
Sonesson B, Hansen F, Lanne T: Abdominal aortic aneurysm: a general defect in the vasculature with focal manifestations in the abdominal aorta? . J Vasc Surg. 1997, 26: 247-254. 10.1016/S0741-5214(97)70185-X
Murgo JP, Westerhof N, Giolma J, Altobelli S: Aortic input impedance in normal man: relationship to pressure waveforms. Circulation. 1980, 62: 105-116.
Murgo JP, Westerhof N, Giolma J, Altobelli S: Effects of exercise on aortic input impedance and pressure wave forms in normal humans. Circ Res. 1981, 48: 334-343.
Bots ML, Hoes AW, Koudstaal PJ, Hofman A, Grobbee DE: Common carotid intima-media thickness and risk of stroke and myocardial infarction: the Rotterdam Study. Circulation. 1997, 96: 1432-1437.
Selzer RH, Mack WJ, Lee PL, Kwong-Fu H, Hodis HN: Improved common carotid elasticity and intima-media thickness measurements from computer analysis of sequential ultrasound frames. Atherosclerosis. 2001, 154: 185-193. 10.1016/S0021-9150(00)00461-5
Schmidt-Trucksass A, Sandrock M, Cheng DC, Muller HM, Baumstark MW, Rauramaa R, Berg A, Huonker M: Quantitative measurement of carotid intima-media roughness--effect of age and manifest coronary artery disease. Atherosclerosis. 2003, 166: 57-65. 10.1016/S0021-9150(02)00245-9
Sandrock M, Cheng DC, Schmitz D, Schmidt-Trucksass A: Quantification of the wall inhomogeneity in B-mode sonographic images of the carotid artery. J Ultrasound Med. 2002, 21: 1395-1404.
Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB: Prediction of coronary heart disease using risk factor categories. Circulation. 1998, 97: 1837-1847.
Stary HC: Natural history and histological classification of atherosclerotic lesions: an update [comment] [see comments]. [Review] [5 refs]. Arterioscler Thromb Vasc Biol. 2000, 20: 1177-1178.
Stary HC: Lipid and macrophage accumulations in arteries of children and the development of atherosclerosis. [Review] [31 refs]. Am J Clin Nutr. 2000, 72: 1297S-1306S.
Stary HC: The development of calcium deposits in atherosclerotic lesions and their persistence after lipid regression. Am J Cardiol. 2001, 88: 16E-19E. 10.1016/S0002-9149(01)01713-1
This study is funded by the Deutsche Forschungs Gesellschaft (DFG) number DI-207-1 (German Research Foundation)
MS carried out the ultrasound examinations and drafted the manuscript, JH drafted the manuscript, JS developed the computer algorithm, DS carried out ultrasound examinations, AN read the manuscript, HB developed the algorithm and AS-T helped with the study design. All authors read and approved the final manuscript.
About this article
Cite this article
Sandrock, M., Hansel, J., Schulze, J. et al. Sequentially based analysis versus image based analysis of Intima Media Thickness in common carotid arteries studies - Do major IMT studies underestimate the true relations for cardio- and cerebrovascular risk?. Cardiovasc Ultrasound 6, 32 (2008). https://doi.org/10.1186/1476-7120-6-32
- Intima Medium Thickness
- Healthy Male Subject
- Image Base Analysis
- Sequential Base Analysis
- Cerebrovascular Risk Factor