Feasibility of perfusion cardiovascular magnetic resonance in paediatric patients
© Buechel et al; licensee BioMed Central Ltd. 2009
Received: 18 June 2009
Accepted: 30 November 2009
Published: 30 November 2009
As coronary artery disease may also occur during childhood in some specific conditions, we sought to assess the feasibility and accuracy of perfusion cardiovascular magnetic resonance (CMR) in paediatric patients.
Methods and results
First-pass perfusion CMR studies were performed under pharmacological stress with adenosine and by using a hybrid echo-planar pulse sequence with slice-selective saturation recovery preparation. Fifty-six perfusion CMR examinations were performed in 47 patients. The median age was 12 years (1 month-18 years), and weight 42.8 kg (2.6-82 kg). General anaesthesia was required in 18 patients. Mean examination time was 67 ± 19 min. Diagnostic image quality was obtained in 54/56 examinations. In 23 cases the acquisition parameters were adapted to patient's size. Perfusion CMR was abnormal in 16 examinations. The perfusion defects affected the territory of the left anterior descending coronary artery in 11, of the right coronary artery in 3, and of the circumflex coronary artery in 2 cases. Compared to coronary angiography, perfusion CMR showed a sensitivity of 87% (CI 52-97%) and a specificity of 95% (CI 79-99%).
In children, perfusion CMR is feasible and accurate. In very young children (less than 1 year old), diagnostic image quality may be limited.
Coronary artery disease (CAD) in paediatric patients may occur after surgery for congenital heart disease involving the coronary arteries, after surgery for congenital coronary anomalies or in patients with an inflammatory disease affecting the mid- and small-size arteries such as Kawasaki disease or Takayasu arteritis [1–3]. Therefore assessment of myocardial ischemia is important not only in the adult population, but increasingly in the paediatric population as well, in which early detection and therapy of CAD may help preventing irreversible myocardial dysfunction .
Exercise-electrocardiography, stress echocardiography, quantitative x-ray coronary angiography (QCA), single-photon emission computed tomography (SPECT) and positron emission tomography (PET) have been so far the available techniques to evaluate myocardial perfusion. QCA is still the reference standard for detecting or quantifying CAD. Radionuclide perfusion imaging suffers from some attenuation artefacts and suboptimal spatial resolution; both techniques necessitate exposure to ionising radiation. As children may need repeated evaluation during long-term follow-up, a non-invasive and radiation-free technique for evaluating the myocardial perfusion is advantageous. Perfusion cardiovascular magnetic resonance (perfusion CMR) using the first-pass kinetics of a gadolinium bolus for detecting myocardial ischemia has shown promising results in clinical trials in adult patients with known or suspected CAD [5–7]. In children perfusion CMR may present some technical challenges, due to higher heart rate (HR), shorter circulation time and high spatial resolution required. On the top of all these challenges patient's cooperation may be sometimes suboptimal, and younger children need to be examined in general anaesthesia.
Aims of this study were to assess the feasibility of performing perfusion CMR in paediatric patients and to determine the diagnostic performance of CMR compared to QCA.
Cardiac diagnoses/surgery and interventions on the coronary arteries
N = 56
Cardiac diagnoses and procedures
Arterial switch operation
Interventions on the coronary arteries
Transfer of the coronary arteries during surgery for CHD
Direct coronary artery surgery
Coronary artery stenting
Coronary surgery included aortocoronary bypass to the left coronary artery in 5 patients, a venous patch augmentation of the left coronary artery in 4, an aortocoronary bypass to the right coronary artery (RCA) in 2 and reduction of a coronary aneurysm with patch augmentation of a nearby stenosis of the RCA in one. Two patients underwent catheter-guided stenting of a coronary stenosis following Kawasaki disease.
All CMR examinations were performed using a 1.5 T scanner (Signa MR/I EchoSpeed, General Electric Medical Systems, Milwaukee, Wisconsin;USA) and a phased-array cardiac coil. The images were acquired during breath holding in all patients, by asking the patients to stop breathing in end-inspiration (n = 38) or by stopping mechanical ventilation in the same position (n = 18). Patient's HR and ECG tracing were continuously monitored; blood pressure was measured at 2 minutes intervals.
The protocol included cine steady-state free precession sequences in a vertical and horizontal long-axis plane, as well as a stack of short-axis acquisitions covering the entire left ventricle for evaluation of left ventricular volumes and function. The steady-state free precession images were acquired at rest with following parameters: TE minimum, flip angle 45°, bandwidth 125 kHz, matrix 224 × 224, number of excitations 1, field of view 26-45 cm, slice thickness 6-8 mm with 1-3 mm gap depending on the body size, retrospective gating with 20 images reconstructed per cardiac cycle.
Myocardial perfusion was assessed, in accordance to our standard protocol, under pharmacological stress with intravenous adenosine in a dose of 0.14 mg/kg/min during 3 minutes. The myocardial signal increase during first-pass was evaluated during a GADOLINIUM bolus injection (Gadopentetate dimeglumine, Magnevist®, Schering, Berlin Germany) in a dose of 0.1 mmol/kg injected in a cubital vein through a large bore cannula at flow rate of 5 ml/sec and followed by a saline flush of 25 ml, flow 5 ml/sec. Five to seven slices were acquired every second heart beat in a short axis plane (slice thickness 8 mm, gap 4 mm) using an hybrid echo-planar pulse sequence with selective notch pulse preparation [5, 6] with following parameters: repetition time 6.5 ms, echo time 1.5 ms, echo train length 4, flip angle 25°, bandwidth ± 125 kHz, matrix 128 × 128, rectangular field of view of 26 × 20 cm to 38 × 28 cm (depending on patient's size), number of excitations 1, acquisition every second RR interval. The images were acquired during a breath-hold of 60 cardiac beats (where patients were instructed to hold their breath as long as possible) or during suspended mechanical ventilation in children requiring general anaesthesia.
Late gadolinium enhancement imaging (LGE) was performed 10-15 minutes after injection (Gadopentetate dimeglumine, Magnevist®, Schering, Berlin Germany; total dose 0.2 mmol/kg) by acquiring images in long- and short-axis planes with a inversion-recovery-prepared gradient-recalled echo pulse sequence and with an adjusted inversion time (180-210 ms) in order to null viable myocardium. Acquisition parameters: echo time minimum full, matrix 256 × 192, number of excitations 2, flip angle 20°, bandwidth ± 31.25 KHz, 24 views per segment, rectangular field of view 26 × 20 cm to 38 × 28 cm depending on patient's size, slice thickness 8 mm, trigger delay depending on HR.
Image analysis was performed on a commercially available workstation (Advantage Windows 4.0, General Electric). The ventricular function was calculated as previously described with the MASS analysis software package (Magnetic Resonance Analytical Software System Version 4.0, MEDIS medical imaging systems, Leiden, Netherlands) . First-pass and LGE images were reviewed as cine loops and stand-images respectively, by two experienced cardiac imagers (JS, EVB), blinded to the results of quantitative coronary angiography. Perfusion abnormalities were assessed qualitatively by visual comparison of the contrast enhancement in different myocardial regions and the final diagnosis was reached by consensus. A significant perfusion deficit was present if gadolinium wash-in was delayed in the subendocardial layer or transmurally, inducing a signal reduction in the myocardium that persisted throughout the entire first-pass of gadolinium, i.e. typically for 5-6 seconds. If a signal reduction was seen in the basal part of the interventricular septum only (membranous portion), this was not considered to represent a perfusion deficit. Small and transient signal reductions not extending over more than one slice were not considered as perfusion deficit.
X-ray coronary angiography
X-ray coronary angiographies were performed during cardiac catheterisation with a biplane Integris BH 5000 cardiovascular system (Philips Medical System, Leiden, The Netherlands). Selective coronary angiography was performed by hand injection of an ionized contrast agent (Iopamiro®, Sintetica SA, Mendrisio, Switzerland). In small children or in cases in whom anatomy precluded selective cannulation of the coronary arteries, contrast was injected into the aortic root at a dose of 1-1.5 ml/kg and with the maximal flow velocity allowed by the size of the catheter (1-2 ml/sec) by using a power injector. The images were acquired with a frame rate of 25/sec in the standard projections for visualisation of the right and the left coronary artery as previously described .
Quantitative coronary angiography (QCA). X-ray coronary angiographies were analysed using Inturis Cardio View software, Release 1.2 (Philips Medical Systems, Leiden, The Netherlands). Before measuring the diameter of the coronary arteries, calibration was obtained by comparing the diameter of the catheter shown on the angiographic images with the catheter diameter in French supplied by the manufacturer. In order to increase the accuracy of the measurements the external diameter of the catheter was measured when empty of contrast. Stenosis was defined as a narrowing of a coronary artery with a reduction of the vessel diameter of ≥ 50%.
An experienced cardiologist (CB) blinded to the results of perfusion CMR analyzed the QCAs.
Numerical data are reported as mean values and standard deviations or as median values and ranges where appropriate. Accuracy was evaluated by calculating sensitivity and specificity (and their 95% confidence interval) of perfusion CMR vs. QCA. A p value < 0.05 was considered as statistically significant. The Research Ethics Board of our institution approved the study. Parents gave their informed consent for the perfusion examination.
Fifty-six perfusion CMR examinations were performed in 47 consecutive patients. Perfusion CMR was successfully performed and diagnostic image quality was obtained in 54 of 56 examinations. The image quality of two examinations was insufficient, due to lack of compliance in a 10 years-old patient and to the very small body size (2.7 kg) in a 3 weeks-old newborn. The overall mean examination (scanning) time was 67 ± 19 min. Baseline mean heart rate (HR) was 78 ± 19 bpm, systolic blood pressure (BP) 102 ± 17 mmHg, diastolic BP 58 ± 13 mmHg. Peak HR under hyperaemia was 100 ± 19 bpm, systolic BP 99 ± 22 mmHg, diastolic BP 53 ± 15 mmHg.
Perfusion CMR was performed in general anaesthesia in 18 patients. In 23 examinations the acquisition parameters were modified according to patient's size. Slice thickness was reduced from 8 to 7 mm in 16, the gap between slices from 4 to 3 mm in 19 studies and to 2 mm in 4. A field of view smaller than 30 × 23 cm was used in 19 studies, all patients younger than 10 years. The injection velocity of the CM bolus was adapted to the size of the cannula and reduced from 5 to 4 ml/sec in 4 cases, to 3 ml/sec in 5 and to 2 ml/sec in the youngest child.
LGE was performed in addition to the first-pass acquisitions in 13 cases. Indication for performing LGE acquisitions consisted of a patient's clinical history suggesting potential occurrence of myocardial scar, i.e. a positive history for myocardial infarction, myocarditis or cardiomyopathy. In one case the image quality of LGE was insufficient due to arrhythmias.
Minor technical problems occurred during 3 examinations and consisted of triggering difficulties, software problems, and initial IV line dysfunction in one patient each. The technical difficulties were transient and the examination could be successfully concluded in all 3 cases.
All examinations were performed without any significant medical complications. In one child under general anaesthesia atropine was given to treat bradycardia (not related to adenosine), one patient showed a self-limiting mild arterial hypotension, and a girl with a total occlusion of left anterior descending coronary artery (LAD) experienced chest pain during adenosine infusion, but remained haemodynamically stable.
Left ventricular volumes and function
Data about left ventricular (LV) volume and function were available for 53 examinations. Mean LV end-diastolic volume (LVEDV) was 79.9 ± 24.6 ml/m2 and mean ejection fraction (EF) was 54.5 ± 9.8%. A dilated LV (LVEDV > 100 ml/m2)  was found in 4 patients: two of them had a combined aortic valve disease, one a LV subvalvular outflow tract obstruction, and another one was performing high-endurance sports.
Reduced function with EF < 45% was measured in 9 patients: 7 of these presented with abnormal coronary arteries, including Kawasaki disease in 5, Bland-White Garland syndrome in one, and LAD narrowing after arterial switch operation in one. Five of these patients had previous surgery of the coronary arteries and 3 had a history of myocardial infarction.
Quantitative X-ray coronary angiography
QCA was performed in 31 patients. The median time interval between perfusion CMR and QCA was 13.5 ± 19 months. Nineteen angiographies were performed within a time interval of less than 3 months (median 2.4 months). Normal coronary arteries were found in 20 patients and a coronary artery stenosis in 8, consisting of a LAD stenosis in 5, RCA stenosis in 2, and stenosis of the LCX in one case. Aneurysmatic enlargement of the coronary arteries was described in 3 children after Kawasaki disease.
Correlation of the findings between perfusion-CMR and coronary angiography
Quantitative coronary x-ray angiography
Compared to QCA perfusion CMR showed a sensitivity of 87% (CI 52%-97%), and a specificity of 95% (CI 79%-99%). Agreement between the perfusion CMR results and the intraoperative findings was 100%.
The clinical impact of perfusion CMR as a non-invasive technique for the assessment of ischemic heart disease is rapidly increasing . In children CAD is much less frequent than in adults and mainly related to congenital heart defects and their surgical repair, or to inflammatory disease affecting the coronary arteries. Nevertheless, particularly young patients may benefit from the advantages of perfusion CMR, such as the lack of exposure to ionising radiation , and its non-invasiveness overcoming the potential complications of cardiac catheterisation .
To our knowledge, this is the largest study reporting the use of perfusion CMR in paediatric patients. Our data demonstrate a feasibility of 96%, and a high accuracy, with a sensitivity of 87% and a specificity of 95%. These results are similar to those reported in adults in single-centre [5, 7] and multi-centre trials [6, 13]. The largest multicentre multivendor perfusion CMR study published recently has shown that perfusion CMR has a better diagnostic performance than SPECT, when results from QCA were used as the reference standard . In the paediatric population a comparison with nuclear perfusion techniques is problematic. PET and SPECT are not well-established diagnostic tools, since ethical reasons preclude the establishment of normal database in infants and children . Moreover, selected congenital conditions such as cyanosis, with possible hypoxia-related effects on the myocardium may preclude a correct interpretation of the perfusion metabolism pattern shown by radionuclide techniques .
The paediatric experience reported so far in assessing ischemic disease with perfusion CMR is scarce [16, 17]. One study demonstrated the feasibility of performing CMR coronary angiography and late enhancement imaging in children . Prakash et al. reported their preliminary experience showing good agreement between perfusion CMR and QCA or SPECT in a limited number of patients and without providing data on sensitivity and specificity. The authors concluded that further validation of the technique with additional studies is needed .
Several larger perfusion CMR trials established key components of perfusion pulse sequences [5, 14]. However, paediatric patients require a particular technical approach with some modifications of pulse sequences parameters. Children have the characteristics of a small body size, higher HR, and, at least in younger subjects, limited compliance. In younger patients examined in general anaesthesia, mechanical ventilation with the option of suspend ventilation during image acquisition is crucial for obtaining a sufficient study quality. Correspondingly, we studied all patients younger than 8 years in general anaesthesia and during breath-holding. This recommendation seems to be valid not only for perfusion CMR, but also for CMR coronary angiography and LGE imaging in children .
We utilized a gadolinium dose of 0.1 ml/kg, which has been demonstrated to give the best results in adults, also in all our paediatric patients. In contrast, the gadolinium injection rate needed to be reduced and adapted to the size of the intravenous cannula in 10 patients. Since Gebker et al. showed that myocardial enhancement and upslope are largely independent from the injection rate of gadolinium bolus as long as the injection speed is not below 3 ml/sec , we targeted to keep injection speed above this limit. Correspondingly in the only case, in which we had to use a slower injection rate of 2 ml/sec, the images resulted not diagnostic.
High HR in children, caused by both, an elevated baseline HR and a more sensitive response to adenosine, which we could consistently observe in our study, require an efficient and robust ECG-triggering system. The relatively new vector cardiac gating systems usually yield to a good triggering, even in higher heart rate, without need for any additional adjustment.
In spite of all these challenges, by carefully tailoring the examination conditions in each individual case, we were able to obtain a diagnostic imaging quality in 96% of the examinations. Perfusion CMR was misleading in two cases: in 3 weeks-old newborn weighting 2.7 kg in whom the rate of gadolinium injection needed to be reduced to 2 ml/sec, and in another child unable to be compliant to hold his breath.
Limitations of the study
This is a retrospective study concealing some limitations. Firstly, our data present a rather wide range of time intervals between perfusion CMR and QCA. The reason for this is that some QCA data used for comparison had been obtained before perfusion CMR was available as a diagnostic technique in our institution. In these patients catheterisation was not repeated if no clinical suspicion of ischemic heart disease was present at perfusion CMR. Other patients were selected for QCA only, when perfusion CMR findings were abnormal. This more clinically oriented approach may have introduced a selection bias in the patient group undergoing QCA, and may have influenced accuracy. This kind of potentially serious limitation is often encountered in paediatric research, as in children ethical reasons preclude the liberal use of invasive examinations or anaesthesia for research purposes.
Secondly, as this study was performed in children, we sought to keep the examinations as short as possible and the patients compliant. Therefore, LGE was not performed consistently in every case, but only if indicated for clinical reasons or if perfusion CMR findings suggested the presence of ischemia and/or scar. Nevertheless, in spite of having applied this single perfusion CMR pulse sequence for evaluation of ischemia, we found a high sensitivity and specificity for the detection of CAD. However, in clinical practice with the suspicion of myocarditis or a storage disease, LGE should be performed as the preferred CMR approach .
Finally, as the CMR scans were performed for clinical reasons, we endeavoured to keep the sequence robust and performed only careful changes in the acquisition parameter for adapting them to the patient's size. Testing different "paediatric protocols" in a standardized way was not the aim of our study, but may be needed in the future for further optimisation and validation of the method in the paediatric population
In the paediatric population, perfusion CMR is feasible and accurate. Age and body size are still limiting factors, and very young children, e.g. < 1 year, are considered suboptimal candidates for this non-invasive technique. However, in older children with suspected coronary alterations perfusion CMR can be recommended as a highly valuable first test as well as an important adjunct to invasive x-ray coronary angiography.
- Hauser M, Bengel FM, Kuhn A, Sauer U, Zylla S, Braun SL, Nekolla SG, Oberhoffer R, Lange R, Schwaiger M, Hess J: Myocardial Blood Flow and Flow Reserve After Coronary Reimplantation in Patients After Arterial Switch and Ross Operation. Circulation. 2001, 103: 1875-1880.View ArticlePubMedGoogle Scholar
- Vogel M, Smallhorn JF, Gilday D, Benson LN, Ash J, Williams WG, Freedom RM: Assessment of Myocardial Perfusion in Patients After the Arterial Switch Operation. J Nucl Med. 1991, 32: 237-241.PubMedGoogle Scholar
- Mavrogeni S, Papadopoulos G, Douskou M, Kaklis S, Seimenis I, Baras P, Nikolaidou P, Bakoula C, Karanasios E, Manginas A, Cokkinos DV: Magnetic resonance angiography isequivalent to X-Ray coronary angiography for the evaluation of coronary arteries in kawasaki disease. Journal of the American College of Cardiology. 2004, 43: 649-652. 10.1016/j.jacc.2003.08.052.View ArticlePubMedGoogle Scholar
- Mavroudis CBC, Muster AJ, Pahl E, Sanders JH, Zales VR, Gevitz M: Expanding indications for pediatric coronary artery bypass. J Thorac Cardiovasc Surg. 1996, 111: 181-189. 10.1016/S0022-5223(96)70415-4.View ArticlePubMedGoogle Scholar
- Schwitter J, Nanz D, Kneifel S, Bertschinger K, Buchi M, Knusel PR, Marincek B, Luscher TF, von Schulthess GK: Assessment of Myocardial Perfusion in Coronary Artery Disease by Magnetic Resonance: A Comparison With Positron Emission Tomography and Coronary Angiography. Circulation. 2001, 103: 2230-2235.View ArticlePubMedGoogle Scholar
- Giang TH, Nanz D, Coulden R, Friedrich M, Graves M, Al-Saadi N, Luscher TF, von Schulthess GK, Schwitter J: Detection of coronary artery disease by magnetic resonance myocardial perfusion imaging with various contrast medium doses: first european multi-centre experience. Eur Heart J. 2004, 25: 1657-1665. 10.1016/j.ehj.2004.06.037.View ArticlePubMedGoogle Scholar
- Nagel E, Klein C, Paetsch I, Hettwer S, Schnackenburg B, Wegscheider K, Fleck E: Magnetic Resonance Perfusion Measurements for the Noninvasive Detection of Coronary Artery Disease. Circulation. 2003, 108: 432-437. 10.1161/01.CIR.0000080915.35024.A9.View ArticlePubMedGoogle Scholar
- Maceira AMPS, Khan M, Pennell DJ: Normalized left ventricular systolic and diastolic function by steady state free precession cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2006, 8: 417-426. 10.1080/10976640600572889.View ArticlePubMedGoogle Scholar
- Takahashi MSR, Wishner SH, Ritchie GW, Francis PS: Selective coronary arteriography in infants and children. Circulation. 1983, 68: 1021-1028.View ArticlePubMedGoogle Scholar
- Sakuma H: Magnetic resonance imaging for ischemic heart disease. J Magn Reson Imaging. 2007, 26: 3-13. 10.1002/jmri.20976.View ArticlePubMedGoogle Scholar
- Bacher K, Bogaert E, Lapere R, De Wolf D, Thierens H: Patient-Specific Dose and Radiation Risk Estimation in Pediatric Cardiac Catheterization. Circulation. 2005, 111: 83-89. 10.1161/01.CIR.0000151098.52656.3A.View ArticlePubMedGoogle Scholar
- Vitiello R, McCrindle BW, Nykanen D, Freedom RM, Benson LN: Complications associated with pediatric cardiac catheterization. Journal of the American College of Cardiology. 1998, 32: 1433-1440. 10.1016/S0735-1097(98)00396-9.View ArticlePubMedGoogle Scholar
- Yates RWM, Marsden PK, Badawi RD, Cronin BF, Anderson DR, Tynan MJ, Maisey MN, Baker EJ: Evaluation of Myocardial Perfusion Using Positron Emission Tomography in Infants Following a Neonatal Arterial Switch Operation. Pediatric Cardiology. 2000, 21: 111-118. 10.1007/s002469910015.View ArticlePubMedGoogle Scholar
- Schwitter J, Wacker CM, van Rossum AC, Lombardi M, Al-Saadi N, Ahlstrom H, Dill T, Larsson HBW, Flamm SD, Marquardt M, Johansson L: MR-IMPACT: comparison of perfusion-cardiac magnetic resonance with single-photon emission computed tomography for the detection of coronary artery disease in a multicentre, multivendor, randomized trial. Eur Heart J. 2008, 29: 480-489. 10.1093/eurheartj/ehm617.View ArticlePubMedGoogle Scholar
- Hernandez-Pampaloni M, Allada V, Fishbein MC, Schelbert HR: Myocardial perfusion and viability by positron emission tomography in infants and children with coronary abnormalities: correlation with echocardiography, coronary angiography, and histopathology. Journal of the American College of Cardiology. 2003, 41: 618-626. 10.1016/S0735-1097(02)02867-X.View ArticlePubMedGoogle Scholar
- Prakash A, Powell AJ, Krishnamurthy R, Geva T: Magnetic resonance imaging evaluation of myocardial perfusion and viability in congenital and acquired pediatric heart disease. Am J Cardiol. 2004, 93: 657-661. 10.1016/j.amjcard.2003.11.045.View ArticlePubMedGoogle Scholar
- Taylor AM, Dymarkowski S, Hamaekers P, Razavi R, Gewillig M, Mertens L, Bogaert J: MR Coronary Angiography and Late-Enhancement Myocardial MR in Children Who Underwent Arterial Switch Surgery for Transposition of Great Arteries. Radiology. 2005, 234: 542-547. 10.1148/radiol.2342032059.View ArticlePubMedGoogle Scholar
- Gebker RPI, Neuss M, Schnackenburg B, Bornstedt A, Jahnke C, Gooma O, Fleck E, Nagel E: Determinants of myocardial response in CMR perfusion imaging using Gd-BOPTA (Multihance). J Cardiovasc Magn Reson. 2005, 7: 565-572. 10.1081/JCMR-200060647.View ArticlePubMedGoogle Scholar
- Kellman P, Arai AE: Imaging Sequences for First Pass Perfusion. A Review. J Cardiovasc Magn Reson. 2007, 9: 525-537. 10.1080/10976640601187604.View ArticlePubMedGoogle Scholar
- Mahrholdt H, Goedecke C, Wagner A, Meinhardt G, Athanasiadis A, Vogelsberg H, Fritz P, Klingel K, Kandolf R, Sechtem U: Cardiovascular Magnetic Resonance Assessment of Human Myocarditis: A Comparison to Histology and Molecular Pathology. Circulation. 2004, 109: 1250-1258. 10.1161/01.CIR.0000118493.13323.81.View ArticlePubMedGoogle Scholar
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.