This article has Open Peer Review reports available.
Shear rate specific blood viscosity and shear stress of carotid artery duplex ultrasonography in patients with lacunar infarction
© Jeong and Rosenson; licensee BioMed Central Ltd. 2013
Received: 10 April 2012
Accepted: 9 March 2013
Published: 18 April 2013
This study describes a new method for determining site-specific vascular shear stress using dynamic measures of shear rate and blood viscosity (BV) in the carotid arteries, and examines characteristics of carotid arterial shear stress among patients with lacunar infarction.
Vascular shear stress measurements were conducted in 37 patients (17 lacunar infarction patients and 20 control subjects) using duplex ultrasonography. Vessel wall diameters and velocities were measured in each arterial segment at peak-systolic (PS) and end-diastolic (ED) phases, for calculation of PS/ED shear rates. PS/ED shear stresses [dyne/cm2] were determined with PS/ED shear rates and shear-rate dependent BV values. For comparison, both values of hematocrit-derived BV and BV measurements at 300 s-1 were used for calculation of shear stress.
All cardiovascular disease (CVD) risk factors including BV values were similar between the two groups. In both common carotid arteries, PS and ED shear stresses were significantly lower in the patients with lacunar infarction than in controls in multivariate models that included age, sex, and other major CVD risk factors. PS and ED shear stresses using the shear rate specific BV were 4.5% lower and 7.3% higher than those using the two other BVs, respectively.
Lacunar infarction was associated with reduced carotid arterial shear stress. The use of estimated BV for calculating carotid arterial shear stress provides more accurate assessment of the hemodynamic contribution of shear stress than previous models that have arbitrarily assigned a constant value to this dynamic flow property.
The vascular system can generally be viewed as a closed system, constraining a series of biochemical processes and mechanical stresses that are regulated by various homeostatic mechanisms . From this theoretical perspective, if mechanical stresses such as wall shear stresses and tensile stresses are not properly dispersed, then an artery may sustain injury. This mechanical injury ultimately results in atherothromboembolic disease and its clinical manifestations of stroke, myocardial infarction and claudication .
Blood is a fluid suspension of plasma and cells such as erythrocytes, leukocytes, and platelets, and it demonstrates non-Newtonian fluid mechanics. A non-Newtonian fluid has a variable, non-linear relation between blood viscosity and blood flow . Specifically, blood viscosity is higher at low shear rates () and is reduced as shear rate increases . This shear rate dependent aspect of blood viscosity (BV, μ) has presented challenges to the accurate calculation of vascular wall shear stress, the frictional force per unit area acting tangentially to the arterial wall. To date, nearly all previous studies on arterial wall shear stress have neglected the flow-dependency of blood  and improperly presented BV as a constant Newtonian fluid [6–8]. Although the non-Newtonian characteristics of blood have been assumed in some previous studies, BV over ranges of shear rates have not been measured directly, but the laws of viscosity were applied empirically [9, 10].
This study was designed to evaluate a new method for calculating vascular shear stresses along carotid arterial segments using duplex ultrasonography and dynamic BV measurements, and to examine whether the vascular shear stress of carotid artery was different between patients with lacunar infarction and control subjects using the new method (shear rate-specific BV and shear stress). This hypothesis is based on the observation that lacunar infarction is accompanied by endothelial dysfunction  and low cerebral arterial blood flow velocities . Finally, we compared the shear stresses of the new method with those of the conventional methods using other BV determinations such as hematocrit-derived BV or constant values at shear rate of 300 s-1.
Seventeen patients with lacunar infarction and 20 control subjects were examined for carotid duplex ultrasonography and measured for BV at the Department of Neurology, Chonbuk National University Hospital, Jeonju, South Korea. A diagnosis of lacunar infarction was made based upon acute onset of focal neurological deficit and a relevant lesion was seen on diffusion-weighted magnetic resonance imaging (MRI) of the brain . A lacunar infarct was diagnosed when the brain MRI identified a lesion less than 15 mm in diameter with characteristics of high signal intensity on diffusion-weighted imaging and low signal intensity on the apparent diffusion coefficient map. The location of the lacunar infarction included either the basal ganglia, corona radiata, thalamus, or brainstem, but not cortex . The control group was defined as patients with cardiovascular risk factors but no previous history of ischemic vascular diseases (e.g. stroke, myocardial ischemia, and lower extremity arterial disease). The institutional ethics committee at Chonbuk National University Hospital approved the present study. All subjects provided written informed consent.
Assessments and measurements
The medical history included data on prior CVD, type 2 diabetes mellitus, hypertension, dyslipidemia, and medication usage. Smoking status was recorded in pack-years. The patients were classified as hypertensive based on a persistent elevation of blood pressure (≥140/90 mmHg) or treatment with antihypertensive medications. Type 2 diabetes mellitus was defined by fasting blood glucose >7.0 mmol/L or previous use of glucose-lowering medications. Hypercholesterolemia was defined by a total cholesterol >6.2 mmol/L or low-density lipoprotein cholesterol >4.1 mmol/L, or current use of cholesterol-lowering medication. Plasma concentration of total homocysteine (tHcy) was measured by fluorescence polarization immunoassay (AxSYM, Abbott Laboratories, Abbott Park, IL).
Velocities were measured at both PS and ED with the sample volume reduced to the smallest possible size (1 mm) and placed in the center of flow, thereby enabling the maximum centerline velocity to be determined. The Doppler angle was generally maintained between 45 and 55 degrees. The sample volume box was placed in the mid-lumen parallel to the vessel wall in all of the studies, as recommended .
The degree of ICA stenosis at both gray-scale and Doppler ultrasonography was stratified into the following 6 categories: (1) normal (no stenosis), (2) less than 50% stenosis, (3) 50~69% stenosis, (4) ≥70% stenosis to near occlusion, (5) near occlusion, and (6) total occlusion, as recommended .
The patients with the following characteristics were excluded: 1) patients who showed ICA stenosis more than 50% in at least one side, 2) patients who have been performed for carotid artery revascularization including carotid endarterectomy or stenting. Finally, patients with lacunar infarction and control subjects, who had ICA of normal or less than 50% stenosis, were included in the present study.
where V and D are the maximum centerline velocity and local lumen diameter, respectively. The PS and ED shear rates were calculated with corresponding PS/ED velocities and diameters, which were designated as PS shear rate ( PS) and ED shear rate ( ED).
where Q is the volume flow rate of blood, d is the inside diameter of a capillary tube, △P is the pressure drop along the capillary tube length, μ is BV, and L is the length of the tube.
Shear rate () specific blood viscosity
Blood viscosity at shear rate of 300 s-1and of hematocrit-derived
where H meant hematocrit (%)/100.
Values of BV at 300 s-1 (μ 300 s-1) were chosen from the measurements using the computerized scanning capillary viscometer.
Descriptive data for the major characteristics were expressed as means ± standard deviations (SD) or percentage as appropriate. An independent t-test was used to determine the statistical differences in the continuous variables, whereas a chi-square test was used for categorical variables. For comparison of adjusted PS/ED shear stresses, analysis of covariance was performed separately in the carotid segments and expressed as adjusted means ± standard errors. The mean values of PS/ED shear stresses using 3 methods of BV determination were compared with paired t-test. All statistical analyses were conducted using PASW Statistics 18 (SPSS, Chicago, IL).
Characteristics of subjects
(n = 20)
(n = 17)
64.7 ± 8.0
68.3 ± 7.7
Smoker, ex- and current,%
Type 2 diabetes mellitus,%
Total cholesterol, mmol/L
4.8 ± 1.1
4.6 ± 0.9
HDL cholesterol, mmol/L
1.2 ± 0.3
1.0 ± 0.3
1.3 ± 0.6
2.4 ± 3.0
Plasma total Hcy, μmol/L
11.2 ± 4.0
12.3 ± 4.4
42.0 ± 5.0
41.7 ± 6.5
Hematocrit-derived BV, cP
4.0 ± 0.6
4.0 ± 0.8
BV, cP, 5 s-1
11.2 ± 2.4
11.4 ± 4.2
4.7 ± 0.8
5.0 ± 1.3
4.0 ± 0.7
4.3 ± 1.2
Carotid intima-media thickness along carotid arteries, and peak-systolic (PS) and end-diastolic (ED) shear rates in both common carotid arteries
(n = 20)
(n = 17)
Intima-media thickness, cm
0.9 ± 0.3
0.8 ± 0.2
1.0 ± 0.2
1.2 ± 0.2
1.0 ± 0.3
1.1 ± 0.3
0.9 ± 0.2
0.8 ± 0.2
1.1 ± 0.3
1.2 ± 0.2
1.0 ± 0.2
1.2 ± 0.3
Shear rate, s-1
Right CCA, PS
405.7 ± 211.5
288.4 ± 68.5
142.9 ± 98.8
97.0 ± 32.7
Left CCA, PS
503.8 ± 192.4
330.8 ± 109.9
201.2 ± 117.2
120.2 ± 55.9
Adjusted * (mean ± SE) peak-systolic and end-diastolic shear stresses (dyne/cm 2 ) of both common carotid arteries
(n = 20)
(n = 17)
Shear stress using shear-rate specific BV
Right CCA, PS
16.0 ± 1.0
12.2 ± 1.1
6.2 ± 0.4
4.7 ± 0.5
Left CCA, PS
19.5 ± 1.3
13.0 ± 1.3
9.1 ± 0.9
5.0 ± 1.0
Shear stress using BV at 300 s-1
Right CCA, PS
16.5 ± 1.2
12.0 ± 1.2
5.7 ± 0.5
3.9 ± 0.5
Left CCA, PS
21.0 ± 1.5
12.7 ± 1.7
8.9 ± 1.0
4.1 ± 1.0
Shear stress using Hct-derived BV
Right CCA, PS
16.6 ± 1.3
11.0 ± 1.4
5.8 ± 0.6
3.6 ± 0.6
Left CCA, PS
20.8 ± 1.6
11.9 ± 1.7
9.1 ± 1.2
3.8 ± 1.2
Comparisons of mean values (± SD, n=37) of PS and ED shear stresses and BVs according to applications of shear rate specific BV (A), BV at 300 s -1 (B), and hematocrit-derived BV (C)
Shear rate-specific BV (A)
BV at shear rate 300 s-1(B)
Hematocrit-derived BV (C)
A vs. B
A vs. C
B vs. C
Shear stress, dyne/cm2
14.0 ± 6.4
14.4 ± 7.2
14.4 ± 8.2
5.5 ± 3.0
4.9 ± 3.1
5.0 ± 3.6
16.2 ± 5.9
17.1 ± 6.8
17.1 ± 8.0
7.0 ± 3.7
6.6 ± 4.1
6.7 ± 4.9
Blood viscosity, cP
4.1 ± 1.0
4.2 ± 0.9
4.0 ± 0.7
4.8 ± 1.2
4.2 ± 0.9
4.0 ± 0.7
4.0 ± 1.0
4.2 ± 0.9
4.0 ± 0.7
4.6 ± 1.2
4.2 ± 0.9
4.0 ± 0.7
The present study introduces a novel method for calculating site-specific vascular shear at the CCA using dynamic measures of shear rate and BV. For the calculation, Casson model constants (k and τ y ) were derived from BV data obtained using a computerized scanning capillary viscometer. In this study, we demonstrate that shear stress is a dynamic parameter that varies widely along the carotid artery. With the novel method, vascular shear stress at the both CCAs was significantly lower in patients with lacunar infarction than the control subjects. Although CVD risk factors were not significantly different between the two groups, and even carotid IMT showed variable findings, vascular shear stresses at the CCAs showed discriminative features. Low shear stress along carotid artery was reported to be associated with ischemic stroke of large artery atherosclerosis, especially in the same side to the affected hemisphere . To our knowledge, this study is the first to compare vascular shear stress in the patients with lacunar infarction with control subjects, on a shear rate specific basis using the individual subject’s own BV profile.
The vascular shear stress of the CCAs, irrespective of the methods of shear stress calculation, showed similar patterns of difference between the patients with lacunar infarction and controls, as shown in Table 3. But, there were significant differences for the mean values of CCA shear stresses and blood viscosities according to the application methods of BV. Using the shear rate specific BV, the lowest values of BV were observed at the PS phase and the highest values were detected at the ED phase. As such, if blood is assumed to behave as a Newtonian fluid and a certain representative value is used, vascular shear stresses will be calculated as higher than their true values at the PS phases, and lower in the ED phases. In the present study, when shear rate specific BVs were used, PS and ED shear stresses were 4.5% higher and 7.3% lower than shear stresses using constant values of BV, as shown in Table 4. This discrepancy would become more pronounced in arterial segments with elevated pulsatility indexes . While the diagnostic methods described herein require further validation, this rationale underscores the non-Newtonian behavior of blood flow and makes a meritorious case for further research.
Our study results are consistent with the report that lacunar infarction is associated with endothelial dysfunction , which is controlled by vascular shear stress . Previously, we reported patients with lacunar infarction had lower blood flow velocities in cerebral arteries than patients with other ischemic stroke subtypes . The both low vascular shear stress in the CCAs and low blood flow velocities in cerebral arteries of patients with lacunar infarction may foster understanding of the pathophysiology of the disease.
Low carotid arterial shear stress in the patients with lacunar infarction who have no obvious carotid stenosis suggests increased resistance in small cerebral arteries  that may potentially increase arterial pressure in the more proximate carotid arteries. The transmission of increased pressure drop of cerebral arteries would cause an increase of carotid arterial pressure, a reduction of flow velocity, a compensatory increase in arterial diameter, and a decrease of vascular shear stress. When the increase of arterial pressure and the decrease of shear stress reinforce themselves through a feedback loop toward greater instability, a vicious cycle might cause lacunar infarction. Although thrombosis due to endothelial denudation under high shear could occur [24, 25], thrombosis in the milieu of low blood flow velocities and low shear stress may also increase the propensity for ischemic events, especially at vulnerable regions . Under low shear, increased transudation of blood or its components in small arteries of the deep cerebral tissue might also occur through the increased porosity of gap junctions, the morphological changes of endothelial cells , to cause lacunar (small-sized) hemorrhage or potentially lacunar infarction .
The present study has several limitations. First, we used a local centerline velocity to determine a local shear rate, which is expressed as 4 V/D. This equation is derived from a parabolic velocity profile for a fully developed laminar flow in a straight circular tube ; so, the assumptions in this model have limitations for estimation of shear stress in regions of complex flow even in the relatively straight CCAs. However, the use of the centerline velocity together with 4 V/D expression allows the estimation of local changes in the velocity gradient if there is no significant hemodynamic compromise . Further investigations are required to determine if centerline arterial shear stress is a parameter that has some critical threshold values for vascular diseases or pathologic conditions. Our study showed 76 carotid arteries in patients with lacunar infarction at CCA peak systolic shear stress of 12~13 dyne/cm2, which is slightly lower than the previous report of carotid arterial shear stress in the subjects with high risk vascular profiles  or with carotid plaque . Second, there were differences of shear rate and shear stress between the right and left CCAs. That might be largely due to the differences of cerebral arterial resistance, and in part, the proximal arterial geometry: the right CCA comes from a celiac trunk and might have more tortuosity than the left CCA which comes from aorta directly. Third, the present cross-sectional study could not address the causal relationship between low carotid artery shear stress and lacunar infarction.
We describe a new method to calculate dynamic vascular shear stress in human carotid arteries using flow rates obtained from ultrasonography and shear rate-specific BV. With the innovative method, the carotid arterial shear stress was significantly lower in patients with lacunar infarction than the control subjects. The prognostic contribution of decreased vascular shear stress to incident lacunar infarction requires further study.
This article was supported by Fund of Biomedical Research Institute of Chonbuk National University Hospital and the measurement of blood viscosity was supported by Medical Device Clinical Trial Center (MDCTC), Chonbuk National University Hospital.
- Cunningham KS, Gotlieb AI: The role of shear stress in the pathogenesis of atherosclerosis. Lab Invest. 2005, 85: 9-23.View ArticlePubMedGoogle Scholar
- Cecchi E, Giglioli C, Valente S, Lazzeri C, Gensini GF, Abbate R, Mannini L: Role of hemodynamic shear stress in cardiovascular disease. Atherosclerosis. 2010, 214: 249-256.View ArticlePubMedGoogle Scholar
- Soulis JV, Farmakis TM, Giannoglou GD, Chatzizisis YS, Giannakoulas GA, Parcharidis GE, Louridas GE: Molecular viscosity in the normal left coronary arterial tree. Is it related to atherosclerosis?. Angiology. 2006, 57: 33-40. 10.1177/000331970605700105.View ArticlePubMedGoogle Scholar
- Fahraeus R, Lindqvist T: The viscosity of blood in narrow capillary tubes. Am J Physiol. 1931, 96: 562-568.Google Scholar
- Cho YI, Kensey KR: Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows. Biorheology. 1991, 28: 241-262.PubMedGoogle Scholar
- Efstathopoulos EP, Patatoukas G, Pantos I, Benekos O, Katritsis D, Kelekis NL: Wall shear stress calculation in ascending aorta using phase contrast magnetic resonance imaging. Investigating effective ways to calculate it in clinical practice. Phys Med. 2008, 24: 175-181. 10.1016/j.ejmp.2008.01.004.View ArticlePubMedGoogle Scholar
- Box FMA, van der Grond J, de Craen AJM, Palm-Meinders IH, van der Geest RJ, Jukema JW, Reiber JHC, van Buchem MA, Blauw GJ, for the PSG: Pravastatin decreases wall shear stress and blood velocity in the internal carotid artery without affecting flow volume: results from the PROSPER MRI Study. Stroke. 2007, 38: 1374-1376. 10.1161/01.STR.0000260206.56774.aa.View ArticlePubMedGoogle Scholar
- Groen HC, Gijsen FJ, van der Lugt A, Ferguson MS, Hatsukami TS, van der Steen AF, Yuan C, Wentzel JJ: Plaque rupture in the carotid artery is localized at the high shear stress region: a case report. Stroke. 2007, 38: 2379-2381. 10.1161/STROKEAHA.107.484766.View ArticlePubMedGoogle Scholar
- Soulis JV, Giannoglou GD, Chatzizisis YS, Seralidou KV, Parcharidis GE, Louridas GE: Non-Newtonian models for molecular viscosity and wall shear stress in a 3D reconstructed human left coronary artery. Med Eng Phys. 2008, 30: 9-19. 10.1016/j.medengphy.2007.02.001.View ArticlePubMedGoogle Scholar
- Gijsen FJ, Wentzel JJ, Thury A, Mastik F, Schaar JA, Schuurbiers JC, Slager CJ, van der Giessen WJ, de Feyter PJ, van der Steen AF: Strain distribution over plaques in human coronary arteries relates to shear stress. Am J Physiol Heart Circ Physiol. 2008, 295: H1608-H1614. 10.1152/ajpheart.01081.2007.View ArticlePubMedGoogle Scholar
- Kim JS, Lee HS, Park HY, Kim SS, Kang HG, Kim NH, Park JS, Kim Y: Endothelial function in lacunar infarction: a comparison of lacunar infarction, cerebral atherosclerosis and control group. Cerebrovasc Dis. 2009, 28: 166-170. 10.1159/000226115.View ArticlePubMedGoogle Scholar
- Kim JT, Lee SH, Hur N, Jeong SK: Blood flow velocities of cerebral arteries in lacunar infarction and other ischemic strokes. J Neurol Sci. 2011, 308: 57-61. 10.1016/j.jns.2011.06.013.View ArticlePubMedGoogle Scholar
- Norrving B: Lacunar infarcts: no black holes in the brain are benign. Pract Neurol. 2008, 8: 222-228. 10.1136/jnnp.2008.153601.View ArticlePubMedGoogle Scholar
- Jiang Y, Kohara K, Hiwada K: Association between risk factors for atherosclerosis and mechanical forces in carotid artery. Stroke. 2000, 31: 2319-2324. 10.1161/01.STR.31.10.2319.View ArticlePubMedGoogle Scholar
- Tahmasebpour HR, Buckley AR, Cooperberg PL, Fix CH: Sonographic examination of the carotid arteries. Radiographics. 2005, 25: 1561-1575. 10.1148/rg.256045013.View ArticlePubMedGoogle Scholar
- Grant EG, Benson CB, Moneta GL, Alexandrov AV, Baker JD, Bluth EI, Carroll BA, Eliasziw M, Gocke J, Hertzberg BS: Carotid artery stenosis: gray-scale and Doppler US diagnosis–Society of Radiologists in Ultrasound Consensus Conference. Radiology. 2003, 229: 340-346. 10.1148/radiol.2292030516.View ArticlePubMedGoogle Scholar
- Kim S, Cho YI, Hogenauer WN, Kensey KR: A method of isolating surface tension and yield stress effects in a U-shaped scanning capillary-tube viscometer using a Casson model. J Non-Newtonian Fluid Mech. 2002, 103: 205-219. 10.1016/S0377-0257(02)00004-6.View ArticleGoogle Scholar
- Rosenson RS, Wolff D, Green D, Boss AH, Kensey KR: Aspirin. Aspirin does not alter native blood viscosity. J Thromb Haemost. 2004, 2: 340-341. 10.1111/j.1538-79333.2004.0615f.x.View ArticlePubMedGoogle Scholar
- Jeong SK, Cho YI, Duey M, Rosenson RS: Cardiovascular risks of anemia correction with erythrocyte stimulating agents: should blood viscosity be monitored for risk assessment?. Cardiovasc Drugs Ther. 2010, 24: 151-160. 10.1007/s10557-010-6239-7.View ArticlePubMedGoogle Scholar
- Guyton AC, Hall JE: Textbook of Medical Physiology. 2000, Philadelphia: W. B. Saunders, 10Google Scholar
- Carallo C, Lucca LF, Ciamei M, Tucci S, de Franceschi MS: Wall shear stress is lower in the carotid artery responsible for a unilateral ischemic stroke. Atherosclerosis. 2006, 185: 108-113. 10.1016/j.atherosclerosis.2005.05.019.View ArticlePubMedGoogle Scholar
- Lim MH, Cho YI, Jeong SK: Homocysteine and pulsatility index of cerebral arteries. Stroke. 2009, 40: 3216-3220. 10.1161/STROKEAHA.109.558403.View ArticlePubMedGoogle Scholar
- Kobayashi S, Okada K, Yamashita K: Incidence of silent lacunar lesion in normal adults and its relation to cerebral blood flow and risk factors. Stroke. 1991, 22: 1379-1383. 10.1161/01.STR.22.11.1379.View ArticlePubMedGoogle Scholar
- Fry DL: Acute vascular endothelial changes associated with Increased blood velocity gradients. Circ Res. 1968, 22: 165-197. 10.1161/01.RES.22.2.165.View ArticlePubMedGoogle Scholar
- Li ZY, Howarth SP, Tang T, Gillard JH: How critical is fibrous cap thickness to carotid plaque stability? A flow-plaque interaction model. Stroke. 2006, 37: 1195-1199. 10.1161/01.STR.0000217331.61083.3b.View ArticlePubMedGoogle Scholar
- Nesbitt WS, Westein E, Tovar-Lopez FJ, Tolouei E, Mitchell A, Fu J, Carberry J, Fouras A, Jackson SP: A shear gradient-dependent platelet aggregation mechanism drives thrombus formation. Nat Med. 2009, 15: 665-673. 10.1038/nm.1955.View ArticlePubMedGoogle Scholar
- Levesque MJ, Liepsch D, Moravec S, Nerem RM: Correlation of endothelial cell shape and wall shear stress in a stenosed dog aorta. Arterioscler Thromb Vasc Biol. 1986, 6: 220-229. 10.1161/01.ATV.6.2.220.View ArticleGoogle Scholar
- Labovitz DL, Boden-Albala B, Hauser WA, Sacco RL: Lacunar infarct or deep intracerebral hemorrhage. Neurology. 2007, 68: 606-608. 10.1212/01.wnl.0000254619.98089.43.View ArticlePubMedGoogle Scholar
- Sutera SP, Skalak R: The history of Poiseuille’s law. Annu Rev Fluid Mech. 1993, 25: 1-19. 10.1146/annurev.fl.25.010193.000245.View ArticleGoogle Scholar
- Irace C, Cortese C, Fiaschi E, Carallo C, Farinaro E, Gnasso A: Wall shear stress is associated with intima-media thickness and carotid atherosclerosis in subjects at low coronary heart disease risk. Stroke. 2004, 35: 464-468. 10.1161/01.STR.0000111597.34179.47.View ArticlePubMedGoogle Scholar
- Gnasso A, Irace C, Carallo C, De Franceschi MS, Motti C, Mattioli PL, Pujia A: In vivo association between low wall shear stress and plaque in subjects with asymmetrical carotid atherosclerosis. Stroke. 1997, 28: 993-998. 10.1161/01.STR.28.5.993.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2377/13/36/prepub
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.