Defining the mid-diastolic imaging period for cardiac CT – lessons from tissue Doppler echocardiography
© Otton et al.; licensee BioMed Central Ltd. 2013
Received: 20 August 2012
Accepted: 25 January 2013
Published: 1 February 2013
Aggressive dose reduction strategies for cardiac CT require the prospective selection of limited cardiac phases. At lower heart rates, the period of mid-diastole is typically selected for image acquisition. We aimed to identify the effect of heart rate on the optimal CT acquisition phase within the period of mid-diastole.
We utilized high temporal resolution tissue Doppler to precisely measure coronary motion within diastole. Tissue-Doppler waveforms of the myocardium corresponding to the location of the circumflex artery (100 patients) and mid-right coronary arteries (50 patients) and the duration and timing of coronary motion were measured. Using regression analysis an equation was derived for the timing of the period of minimal coronary motion within the RR interval. In a validation set of 50 clinical cardiac CT examinations, we assessed coronary motion artifact and the effect of using a mid-diastolic imaging target that was adjusted according to heart rate vs a fixed 75% phase target.
Tissue Doppler analysis shows the period of minimal cardiac motion suitable for CT imaging decreases almost linearly as the RR interval decreases, becoming extinguished at an average heart rate of 91 bpm for the circumflex (LCX) and 78 bpm for the right coronary artery (RCA). The optimal imaging phase has a strong linear relationship with RR duration (R2 = 0.92 LCX, 0.89 RCA). The optimal phase predicted by regression analysis of the tissue-Doppler waveforms increases from 74% at a heart rate of 55 bpm to 77% at 75 bpm. In the clinical CT validation set, the optimal CT acquisition phase similarly occurred later with increasing heart rate. When the selected cardiac phase was adjusted according to heart rate the result was closer to the optimal phase than using a fixed 75% phase. While this effect was statistically significant (p < 0.01 RCA/LCx), the mean effect of heart-rate adjustment was minor relative to typical beat-to-beat variability and available precision of clinical phase selection.
High temporal resolution imaging of coronary motion can be used to predict the optimal acquisition phase in cardiac CT. The optimal phase for cardiac CT imaging within mid-diastole increases with increasing heart rate although the magnitude of change is small.
KeywordsCardiac CT Image quality Heart rate Tissue Doppler Echocardiography
Improvements in cardiac computed tomography technology have enabled fast and reliable imaging of the coronary arteries. Nevertheless, temporal resolution remains a limiting factor and the timing of image acquisition is critical for optimal imaging quality. While various empiric principles have been adopted to guide coronary image acquisition, to date, studies of timing strategies have generally focused on preferable timing bands of CT acquisition [1–4] during systole or diastole according to heart rate. Mid-diastole is typically selected at lower heart rates. Precise guidance as to the optimal phase within this period has been limited due to the temporal resolution of coronary tracking techniques available.
Echocardiography, and in particular tissue Doppler, is particularly suited due to its high temporal resolution, accuracy and ability to interrogate the myocardial motion at the arterial segments most vulnerable to imaging artifact . Tissue Doppler provides motion assessment of the coronary arteries with a temporal resolution of less than 5 ms  as opposed to the much poorer 30-50 ms resolution of invasive angiography [7, 8], 10-20 ms resolution for MRI  or the 70-160 ms temporal resolution of CT imaging itself . Tissue Doppler may therefore be used to distinguish between coronary motion present at single-percentage phase intervals within the diastolic period.
Optimal timing of cardiac CT is important firstly to minimize or eliminate motion artifact, and secondly to guide imaging parameters with regards to aggressive dose minimization strategies, which frequently employ limited cardiac phases [11, 12]. It is well recognized that the fraction of the cardiac cycle occupied by diastole decreases with increasing heart rate and comes proportionally later in the cardiac cycle. It would therefore be expected that the optimal window for cardiac imaging may vary with heart rate and that small adjustments to the acquired phase may result in improved image quality.
We sought to investigate the motion of the coronary arteries throughout the cardiac cycle in order to provide both a qualitative and mathematical description of the period of optimal diastolic cardiac image acquisition in order to guide coronary imaging by low-dose, limited-phase cardiac CT. We compared the resultant analysis and prediction equations to motion artifact outcomes within a clinical cardiac CT cohort.
Echocardiographic measurements were obtained from anonymised routine clinical data and participants within the CT cohort gave written consent for the use of clinical data for research purposes. The study was performed in accordance with the principles of the Helsinki declaration and was approved by the St Vincent’s Hospital Human Ethics Committee (approval number 10/044).
In order to record physiological motion of the coronary arteries, consecutive clinical echocardiograms performed at our institution were retrospectively reviewed from June 2010. Echocardiograms were performed using either a Philips ie33 or Sonos 7500 echocardiography system by an experienced technician. One hundred echocardiograms with clear tissue Doppler traces of the left ventricular lateral wall and fifty with right ventricular lateral wall tissue Doppler were selected from a sequential cohort of 324 patients. Patients were referred from both inpatient and outpatient settings for standard indications including the assessment of ischemic heart disease, heart failure and valvular pathology. Exclusion criteria were poor image quality and the presence of a cardiac rhythm other than sinus rhythm.
The tracked volume was selected to encompass the lateral tricuspid annulus corresponding to the course of the mid right coronary artery and the lateral mitral valve annulus, corresponding to the mid-distal position of the circumflex artery.
Utilizing the recorded ECG trace simultaneous to echocardiographic image acquisition, the time interval from the mid-point of the QRS trace to the start and end of the isovolumic relaxation time (IVRT), center of the E’ and A’ waves and start and end of the period minimal cardiac motion (PMCM) was measured. The period of minimal cardiac motion was defined as the diastolic phase between the E’ and A’ during which time the spectrum of the tissue Doppler velocity trace indicated tissue velocity of less than 2 cm/s. The duration of the QRS interval, current and preceding RR interval and ejection fraction were also noted.
Clinical CT validation study
In a validation cohort of 50 consecutive cardiac CT studies we prospectively assessed motion artifact within mid-diastole. The prospective cohort of 50 consecutive clinical cardiac CT patients included all scans performed for standard clinical indications, and included a broad range of heart rates. Patients were excluded from the validation cohort if the underlying rhythm was not sinus rhythm. Patients within the derivation echocardiography set were not included within the CT validation cohort.
All CT scans were obtained utilizing a 320 detector-row CT (Toshiba Aquilion One) with z-axis coverage adjusted to ensure complete cardiac coverage within a single gantry rotation and image acquisition of the entire coronary tree at identical cardiac phase. The x-ray tube current (mA), voltage (kV) and contrast volume (Ultravist370) were adjusted according to body size and standard clinical protocols. A FC03 filter kernel was used for reconstruction, without the use of iterative noise reduction. The resultant CT images were displayed at 0.5 mm  isotropic resolution. Where multiple cardiac cycles were acquired (for example due to very high heart rates or multiple volume acquisition), only the first RR interval was used for timing analysis.
The acquired cardiac phase was defined as the percentage of the RR interval corresponding to the midpoint of the half-scan reconstruction. All scans were acquired with a minimum of 170 ms of widened data acquisition, enabling reconstruction of a range of 70-80% of cardiac phase or greater, during which the mid-diastolic period is expected to occur at all physiological heart rates. In patients with a baseline heart rate of greater than 65 beats per minute or greater, a 30-80% range of cardiac phases were acquired. In all scans it was anticipated that the period of mid-diastole would be fully captured.
Tissue Doppler analysis
Utilizing separate ordinary least squared linear regression models within STATA statistical software v10, StataCorp, USA we assessed the relationship between the RR interval and the dependent variables: time to the center of the IVRT, time to the E’ wave, time to the A’ wave and time to center of the period of minimal cardiac motion. The relationship and regression residuals were inspected visually. In addition, non-linearity was assessed using the Ramsay RESET test for omitted variables. Using simple algebraic substitution (noting Heart Rate = 60 seconds/RR interval) the regression equation relating the RR interval and the time to the center of the PMCM was used to form a prediction equation relating heart rate to the ideal target of cardiac imaging. This heart rate adjusted ideal phase equation was subsequently tested in the CT validation cohort.
As the duration of the PMCM, or imaging window, is limited by a lower bound of zero, truncated linear regression was used to generate an equation relating the RR interval and PMCM. This equation was used to forecast the PMCM based on the RR interval. The standard deviation of the forecast was calculated and used to create 90% confidence intervals of the forecast, which was then graphically displayed. The forecast equation of the PMCM was then displayed against heart rate within heart rates clinically relevant to diastolic imaging. This allowed the relationship between the CT acquisition time and the probability that the acquisition time falls within the PMCM to be graphically demonstrated.
CT validation cohort
The relationship between heart rate and ideal imaging target (the mid-point of the artifact free CT image time period) was assessed using ordinary least square regression. The null hypothesis, that there was no relationship between heart rate and ideal imaging target, was tested using a two-sided T test of the regression slope coefficient.
A fixed 75% phase
The heart-rate predicted ideal imaging phase
was measured. The error of each method of phase selection was compared using a paired t-test. The two-sided threshold for significance of the T-test was set at p < 0.05.
Timing of tissue Doppler features and the cardiac imaging window
The center of the right coronary IVRT showed only a modest correlation with heart rate (R2 = 0.38) irrespective of QRS duration or ejection fraction. The center of the IVRT occurred within the range (95% confidence interval) 32–63% of the RR duration, occurring later as a percentage of the RR interval as the heart rate increased
Correlation of the time to E’ was strongly related to the RR duration (LCX R2 = 0.60; RCA R2 = 0.52 p < 0.001 for both). The A’ was very strongly related to the RR duration (LCX R2 = 0.96; RCA R2 = 0.92. P < 0.001 for both). (Figure 4).
The coefficient of determination (R2) of the relationship between the ideal imaging target and the RR interval was 0.92 for the circumflex and 0.86 for the right coronary artery. The regression equation for the optimal target of CT acquisition was found to be:
Optimal CT trigger time(ms) Circumflex = RR interval x 0.67 + 80 ms.
Optimal CT trigger time(ms) Right Coronary Artery = RR interval x 0.68 + 70 ms.
Heart rate adjusted phase targets for both the right and circumflex coronary arteries
Optimal diastolic target (% phase)
Doppler study of the duration of the cardiac imaging window in mid-diastole
As expected, tissue Doppler recording of both the right and left atrio-ventricular groove demonstrated greatest motion during systole. The mean duration of the IVRT, or end-systolic imaging window, was 55 ms and only a weak positive relation with the RR interval was noted (R2 = 0.28, p < 0.001).
Validation patient demographic data
Body Mass Index
In 38% of patients the stated reason for referral for CT was assessment due to cardiovascular risk factors or equivocal prior investigations. In 34% of patients the stated reason for referral for cardiac CT was of chest pain, and 14% were referred for prior arrhythmia. No referral reason was recorded for 14% of patients. The median radiation dose of all scans was 6.2 mSv.
The ideal cardiac phase for both the right and left coronary arteries increased with heart rate. The regression coefficient for the circumflex was 0.17 (95% CI 0.090-0.25, T-test P < 0.01) and 0.15 (95% CI 0.038-0.27, T-test P < 0.05) for right coronary artery indicating a 1.7% or 1.5% increase in the ideal imaging target in mid-diastole for every 10 beat increase in heart rate.
Analysis of time between the ideal cardiac imaging target, 75% phase target and the heart rate adjusted target (as per Table 1) demonstrated that the heart rate adjusted target was statistically closer to the ideal target than the fixed 75% target for both the right coronary artery and left anterior descending/circumflex arteries, (p < 0.01 for both). The absolute benefit of the heart rate adjusted strategy was 0.5% (SD 0.14%) for the circumflex and 0.4% (SD 0.13%)for the RCA. This partially reflects minor differences between the two prediction models at the most common heart rate 55–65 bpm. At heart rates greater than 65 the heart rate adjusted strategy was a mean of 1% for the RCA and 0.6% for the circumflex, closer to the ideal imaging phase.
This study demonstrates the role of careful examination of the physiology of coronary motion in understanding the optimal timing of cardiac CT within diastole and the impact of image acquisition speed on motion artifact. The study also demonstrates the difficulties of cardiac imaging timing given the inherent physiological variability in coronary motion. While we primarily assessed the effects of coronary motion on cardiac CT, the research may also have a role in other gated imaging modalities such as cardiovascular MRI.
The temporal resolution of half-scan reconstruction of current generation single source cardiac CT lies between 140 and 175 ms. While current dual-source imaging reduces the temporal resolution to 70 ms-83 ms [1, 13], in its lowest dose, high-pitch exemplification, the time taken to image the entire axial cardiac volume remains in the order of 280 ms. Retrospective gating of cardiac CT, and increasing the window of acquisition can go some way towards reducing the effects of cardiac motion at higher heart rates, but at the cost of greatly increased radiation dose. Retrospective gated cardiac CT may deliver 4 or 5 times the ionizing radiation of prospective CT  and even when prospective cardiac CT protocols are used, for every extra 100 ms of the cardiac cycle acquired, radiation dose is increased by up to 45% . On the other hand, very low dose prospectively gated cardiac CT acquired during minimal phase acquisition enables only one or a very narrow band of phases to be analyzed, compromising diagnosis in the case of motion artifact. Accurate timing of cardiac CT is therefore vital for both radiation dose reduction and prevention of imaging artifacts.
High temporal resolution imaging of coronary artery motion reveals interesting insights about coronary motion, which may be useful for effective cardiac imaging. Previous measurements of coronary artery motion using electron beam CT [16, 17] or coronary angiography  are inadequate for the precise definition of the cardiac rest period within mid-diastole. This is because the low temporal resolution of these techniques can only provide information as to which broad band of phase acquisition is most appropriate according to heart rate. The issue as to whether end-systole or diastole is superior at a given heart rate is clinically important and has been previously studied [3, 4, 18–20]. Likewise the interaction of gantry rotation time, heart rate and image artifact requires careful attention in order to optimize image quality [21, 22]. The current research adds to these prior works by further elucidating the nature, effects and variability of coronary motion during the imaging phase of mid-diastole.
The equations and tables provided are intended as a guide for gated cardiac imaging across all CT types. They indicate the heart rates at which diastolic imaging becomes viable, the window available for imaging and the chance of motion artifact for given heart rates. In the validation cohort, a strategy of varying the targeted phase according to heart rate was statistically superior to a fixed percentage acquisition in terms of proximity to the ideal imaging phase, although the extent of benefit was small. These tools should prove useful for the achievement of high quality cardiac CT imaging. Nevertheless, there are several limitations to the formulas and conclusions provided and they should not be applied without consideration.
Firstly, it should be noted that even small differences in absolute timing make a large difference in the optimal phase selection. Variation can occur due to triggering parameters, QRS measurement and biological variability. Beat to beat variation of more than 50 ms can be expected at most heart rates and no system of CT triggering can ever account for the unpredictability associated with ectopic beats or atrial fibrillation. The wide variation in the ideal imaging target within the CT validation cohort may have reflected ectopy during the period of heart-rate assessment, although as we were only able to capture the two heartbeats immediately prior to imaging we were unable to assess the impact of ectopy in this study. While measurement of time intervals on tissue Doppler has generally good reliability , the use of a single measurement in our study should also be considered a limitation and may have added to variability of tissue Doppler measurements.
Secondly, the definition of motion artifact used in this study was conservative and it is likely that a degree of motion artifact may be tolerated without affecting clinical interpretability. The effect of motion artifact on diagnosis is inherently subjective and the precise heart-rate boundaries of CT interpretability may vary from the predictions of this study.
Lastly, the difference between a fixed 75% phase trigger and a heart rate adjusted strategy was small, and both strategies benefit from a degree of extra phase acquisition at most heart rates and gantry speeds. CT coronary angiography performed as a part of this study was designed to ensure coverage of mid-diastole and analysis of systolic or end-systolic phases could not be performed. Adjustment of the precise imaging target within mid-diastole should be a consideration secondary to the more fundamental issue of whether end-systolic or diastolic imaging is required.
While the benefit of heart rate adjustment for coronary CT phase target selection within any individual patient is small, over a large cohort where dose minimization through single-phase acquisition is mandated, even minor adjustments in diastolic timing may result in a reduced chance of image motion artifact. Knowledge of physiological coronary motion reinforces the importance of heart rate control, acquisition time and phase selection in cardiac imaging.
Left circumflex artery
Right coronary artery
Period of minimal cardiac motion
Isovolumic relaxation time.
Grants and funding: Database infrastructure for the study was provided via a grant from the St Vincent’s Clinic Foundation, Sydney Australia, a charitable medical research foundation. http://www.clinic.stvincents.com.au/clinic-foundation.
James Otton would like to acknowledge scholarship support from the National Health and Medical Research of Council of Australia, and grant support from the St Vincent’s Clinic Foundation of Australia
Chung-Yao Yu is the current Swire Fellow in Advanced Cardiac Imaging and would like to acknowledge the Swire company for their generous financial support.
- Achenbach S, Manolopoulos M, Schuhback A, Ropers D, Rixe J, Schneider C, Krombach GA, Uder M, Hamm C, Daniel WG, Lell M: Influence of heart rate and phase of the cardiac cycle on the occurrence of motion artifact in dual-source CT angiography of the coronary arteries. J Cardiovasc Comput Tomogr. 2012, 6: 91-98. 10.1016/j.jcct.2011.11.006.View ArticlePubMedGoogle Scholar
- Leschka S, Husmann L, Desbiolles LM, Gaemperli O, Schepis T, Koepfli P, Boehm T, Marincek B, Kaufmann PA, Alkadhi H: Optimal image reconstruction intervals for non-invasive coronary angiography with 64-slice CT. Eur Radiol. 2006, 16: 1964-1972. 10.1007/s00330-006-0262-x.View ArticlePubMedGoogle Scholar
- Seifarth H, Wienbeck S, Pusken M, Juergens KU, Maintz D, Vahlhaus C, Heindel W, Fischbach R: Optimal systolic and diastolic reconstruction windows for coronary CT angiography using dual-source CT. AJR Am J Roentgenol. 2007, 189: 1317-1323. 10.2214/AJR.07.2711.View ArticlePubMedGoogle Scholar
- Herzog C, Abolmaali N, Balzer JO, Baunach S, Ackermann H, Dogan S, Britten MB, Vogl TJ: Heart-rate-adapted image reconstruction in multidetector-row cardiac CT: influence of physiological and technical prerequisite on image quality. Eur Radiol. 2002, 12: 2670-2678.PubMedGoogle Scholar
- Zamorano J, Wallbridge DR, Ge J, Drozd J, Nesser J, Erbel R: Non-invasive assessment of cardiac physiology by tissue Doppler echocardiography. A comparison with invasive haemodynamics. Eur Heart J. 1997, 18: 330-339. 10.1093/oxfordjournals.eurheartj.a015236.View ArticlePubMedGoogle Scholar
- Van de Veire NR, De Sutter J, Bax JJ, Roelandt JR: Technological advances in tissue Doppler imaging echocardiography. Heart. 2008, 94: 1065-1074. 10.1136/hrt.2007.120758.View ArticlePubMedGoogle Scholar
- Wang Y, Vidan E, Bergman GW: Cardiac motion of coronary arteries: variability in the rest period and implications for coronary MR angiography. Radiology. 1999, 213: 751-758.View ArticlePubMedGoogle Scholar
- Johnson KR, Patel SJ, Whigham A, Hakim A, Pettigrew RI, Oshinski JN: Three-dimensional, time-resolved motion of the coronary arteries. J Cardiovasc Magn Reson: official journal of the Society for Cardiovascular Magnetic Resonance. 2004, 6: 663-673. 10.1081/JCMR-120038086.View ArticleGoogle Scholar
- Tangcharoen T, Bell A, Hegde S, Hussain T, Beerbaum P, Schaeffter T, Razavi R, Botnar RM, Greil GF: Detection of coronary artery anomalies in infants and young children with congenital heart disease by using MR imaging. Radiology. 2011, 259: 240-247. 10.1148/radiol.10100828.View ArticlePubMedGoogle Scholar
- Flohr TG, Raupach R, Bruder H: Cardiac CT: how much can temporal resolution, spatial resolution, and volume coverage be improved?. J Cardiovasc Comput Tomogr. 2009, 3: 143-152. 10.1016/j.jcct.2009.04.004.View ArticlePubMedGoogle Scholar
- Husmann L, Valenta I, Gaemperli O, Adda O, Treyer V, Wyss CA, Veit-Haibach P, Tatsugami F, von Schulthess GK, Kaufmann PA: Feasibility of low-dose coronary CT angiography: first experience with prospective ECG-gating. Eur Heart J. 2008, 29: 191-197.View ArticlePubMedGoogle Scholar
- Achenbach S, Marwan M, Ropers D, Schepis T, Pflederer T, Anders K, Kuettner A, Daniel WG, Uder M, Lell MM: Coronary computed tomography angiography with a consistent dose below 1 mSv using prospectively electrocardiogram-triggered high-pitch spiral acquisition. Eur Heart J. 2010, 31: 340-346. 10.1093/eurheartj/ehp470.View ArticlePubMedGoogle Scholar
- McCollough CH, Schmidt B, Yu L, Primak A, Ulzheimer S, Bruder H, Flohr TG: Measurement of temporal resolution in dual source CT. Medical physics. 2008, 35: 764-768. 10.1118/1.2826559.View ArticlePubMedPubMed CentralGoogle Scholar
- Maruyama T, Takada M, Hasuike T, Yoshikawa A, Namimatsu E, Yoshizumi T: Radiation dose reduction and coronary assessability of prospective electrocardiogram-gated computed tomography coronary angiography: comparison with retrospective electrocardiogram-gated helical scan. J Am Coll Cardiol. 2008, 52: 1450-1455. 10.1016/j.jacc.2008.07.048.View ArticlePubMedGoogle Scholar
- Labounty TM, Leipsic J, Min JK, Heilbron B, Mancini GB, Lin FY, Earls JP: Effect of padding duration on radiation dose and image interpretation in prospectively ECG-triggered coronary CT angiography. AJR Am J Roentgenol. 2010, 194: 933-937. 10.2214/AJR.09.3371.View ArticlePubMedGoogle Scholar
- Achenbach S, Ropers D, Holle J, Muschiol G, Daniel WG, Moshage W: In-plane coronary arterial motion velocity: measurement with electron-beam CT. Radiology. 2000, 216: 457-463.View ArticlePubMedGoogle Scholar
- Mao S, Lu B, Oudiz RJ, Bakhsheshi H, Liu SC, Budoff MJ: Coronary artery motion in electron beam tomography. J Comput Assist Tomogr. 2000, 24: 253-258. 10.1097/00004728-200003000-00012.View ArticlePubMedGoogle Scholar
- Mok GS, Yang CC, Chen LK, Lu KM, Law WY, Wu TH: Optimal systolic and diastolic image reconstruction windows for coronary 256-slice CT angiography. Acad Radiol. 2010, 17: 1386-1393. 10.1016/j.acra.2010.06.011.View ArticlePubMedGoogle Scholar
- Sun G, Li M, Li L, Li GY, Zhang H, Peng ZH: Optimal systolic and diastolic reconstruction windows for coronary CT angiography using 320-detector rows dynamic volume CT. Clin Radiol. 2011, 66: 614-620. 10.1016/j.crad.2011.02.007.View ArticlePubMedGoogle Scholar
- Husmann L, Leschka S, Desbiolles L, Schepis T, Gaemperli O, Seifert B, Cattin P, Frauenfelder T, Flohr TG, Marincek B, Kaufmann PA, Alkadhi H: Coronary artery motion and cardiac phases: dependency on heart rate – implications for CT image reconstruction. Radiology. 2007, 245: 567-576. 10.1148/radiol.2451061791.View ArticlePubMedGoogle Scholar
- Greuter MJ, Dorgelo J, Tukker WG, Oudkerk M: Study on motion artifacts in coronary arteries with an anthropomorphic moving heart phantom on an ECG-gated multidetector computed tomography unit. Eur Radiol. 2005, 15: 995-1007. 10.1007/s00330-004-2602-z.View ArticlePubMedGoogle Scholar
- Greuter MJ, Flohr T, van Ooijen PM, Oudkerk M: A model for temporal resolution of multidetector computed tomography of coronary arteries in relation to rotation time, heart rate and reconstruction algorithm. Eur Radiol. 2007, 17: 784-812. 10.1007/s00330-006-0228-z.View ArticlePubMedGoogle Scholar
- Fraser AG, Payne N, Madler CF, Janerot-Sjoberg B, Lind B, Grocott-Mason RM, Ionescu AA, Florescu N, Wilkenshoff U, Lancellotti P, Wutte M, Brodin LA: Feasibility and reproducibility of off-line tissue Doppler measurement of regional myocardial function during dobutamine stress echocardiography. Eur J Echocardiogr. 2003, 4: 43-53. 10.1053/euje.4.1.43.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/13/5/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.