The prognostic value of gray-white-matter ratio in cardiac arrest patients treated with hypothermia
- Michael Scheel†1, 4Email author,
- Christian Storm†3,
- Andre Gentsch1,
- Jens Nee3,
- Fridolin Luckenbach3,
- Christoph J Ploner2 and
- Christoph Leithner2
© Scheel et al.; licensee BioMed Central Ltd. 2013
Received: 21 January 2013
Accepted: 17 March 2013
Published: 8 April 2013
Mild therapeutic hypothermia alters the validity of a number of parameters currently used to predict neurological outcome after cardiac arrest and resuscitation. Thus, additional parameters are needed to increase certainty of early prognosis in these patients. A promising new approach is the determination of the gray-white-matter ratio (GWR) in cranial computed tomography (CCT) obtained early after resuscitation. It is not known how GWR relates to established outcome parameters such as neuron specific enolase (NSE) or somatosensory evoked potentials (SSEP).
Cardiac arrest patients (n = 98) treated with hypothermia were retrospectively analyzed with respect to the prognostic value of GWR, NSE and SSEP.
A GWR < 1.16 predicted poor outcome with 100% specificity and 38% sensitivity. In 62 patients NSE, SSEP and CCT were available. The sensitivity of poor outcome prediction by both NSE > 97 μg/L and bilateral absent SSEP was 43%. The sensitivity increased to 53% in a multi-parameter approach predicting poor outcome using at least two of the three parameters (GWR, NSE and SSEP).
Our results suggest a strong association of a low GWR with poor outcome following cardiac arrest. Determination of the GWR increases the sensitivity in a multi-parameter approach for prediction of poor outcome after cardiac arrest.
An increasing number of recent studies on outcome prognostication after cardiac arrest demonstrate that hypothermia treatment alters prognostic parameters that had been established in normothermic patients. For example, a level of neuron specific enolase (NSE) above 33 μg/L, obtained on day one to three after cardiac arrest, and motor reaction to painful stimuli not better than extension on day three, are associated with a substantial rate of false poor outcome predictions in these patients [1–8]. Individual patients have been reported who survived with good neurological outcome despite NSE levels of 97 μg/L , bilateral absent N20 of medianus somatosensory evoked potentials (SSEP) [4, 5], absent pupillary light responses or corneal reflexes [6, 7] and early status epilepticus . Before limitation of treatment is considered, a high level of certainty is needed and the combination of several prognostic parameters may increase the level of certainty for poor outcome prediction.
Brain imaging appears a self-evident choice for prediction of neurological outcome in patients with suspected hypoxic-ischemic encephalopathy. Severe cases are characterized by cerebral edema with decreased attenuation of gray matter in cranial computed tomography (CCT). A loss of contrast between gray and white matter has been shown to be associated with poor outcome [9, 10]. Recent studies have evaluated CT signs, such as cerebral gray-white matter ratio (GWR) as an early prognostic outcome parameter in hypoxic-ischemic encephalopathy [11–16]. GWR is obtained by measuring Hounsfield units (HU) in gray and white matter and calculating the ratio of both . In severe cases of cerebral edema, the ratio decreases to values below unity, also known as the ‘reversal sign’ .
There is currently no consensus on a distinct GWR cut-off value that may predict poor outcome with high specificity. A small number of studies on this subject suggest values between 1.15  and 1.22 , respectively. Furthermore, it is neither known whether hypothermia affects the prognostic significance of a reduced GWR, nor whether the sensitivity with which poor outcome is predicted by NSE and SSEP can be increased by measuring the GWR. We therefore retrospectively investigated the prognostic value of GWR in a cohort of cardiac arrest patients treated with hypothermia and compared GWR with the established outcome parameters NSE and SSEP.
Materials and methods
A total number of 353 consecutive patients after cardiac arrest were prospectively included into our database between 12/2005 and 10/2011. Since 2005, all patients have been treated with mild hypothermia in accordance with established guidelines (details of this treatment have been previously described ). The study protocol was approved by the local ethics committee on human research (Charité - Ethikkommission) and was conducted in accordance with the guidelines of the Declaration of Helsinki. Written informed consent to the use of routine clinical data is part of the standard contract between patients and the University Hospital Charité Berlin and was obtained from patients or their legal representatives. All data were collected from routine investigations in patients resuscitated from cardiac arrest. In our department, neurological assessment of these patients comprises several clinical, laboratory and electrophysiological investigations. Treatment decisions are always based on combined results from these examinations. Except for patients with confirmed brain death, treatment is continued for at least seven days in all patients .
In our database, we identified a subgroup of patients (n = 111) that had received cranial computed tomography (CCT) within the first seven days after cardiac arrest (median = 5 h; IQR = 2–24 hours). CCT was ordered independently of this study by the treating physician either to rule out a primary intracranial event leading to cardiac arrest (e.g. subarachnoid hemorrhage) or when clinical findings suggested intracranial complications (e.g. intracranial hemorrhage during anticoagulant therapy for myocardial infarction). Thirteen patients were excluded from further analysis for findings that would have biased measurement of Hounsfield units and GWR calculation (only contrast enhanced CCT available (3), hydrocephalus and shunt artifact (3), severe movement artifacts (2), intracerebral hemorrhage (3), old large ischemic lesion (1), massive calcification of the basal ganglia (1)).
Demographic data given as absolute numbers and percent or median and interquartile range (IQR)
Demographics of study population
Primary cause of arrest
Time to ROSC (min)
Total epinephrine dose (mg)
Length of ICU stay (days)
Time on ventilator (hours)
CCT acquisition, region of interest (ROI) measurements and GWR calculation
Only non-contrast enhanced CT scans were analyzed. CT scans were acquired on three different GE scanners (GE LightSpeed Pro 16 (n = 78), Lightspeed Ultra (n = 7), Lightspeed VCT (n = 13)-GE Healthcare, Little Chalfont, UK). All examinations followed a standard head CT protocol with a slice thickness of 5 mm.
For assessment of inter-raterreliability, ROIs were determined independently by two raters. ROIs were placed and GWR scores calculated according to criteria described previously . The ROIs consisted of circular shaped areas (0.1 cm2) placed bilaterally in caudate nucleus (CN), putamen (PU), thalamus (THL), posterior limb of internal capsule (PIC), forceps minor of the corpus callosum (CC), medial cortex (MC1) and medial white matter (MWM1) at the level of the centrum semiovale and high convexity area (MC2 and MWM2, respectively). Both raters were blind to clinical information and outcome of all patients as well as ROI placement and GWR scores of each other. GWR in the basal ganglia was calculated as GWR-BG = (CN + PU) / (CC + PIC). GWR at the cortical level was calculated as GWR-CE = (MC1 + MC2) / (MWM1 + MWM2). Average GWR was defined as the average of these two GWR-AV = (GWR-BG + GWR-CE)/2 .
Neuron specific enolase and somatosensory evoked potentials
NSE serum levels were determined in 84 out of 98 patients as described in an earlier study . We have previously found a cut-off value of 78.9 μg/L for poor outcome prediction with 100% specificity , and a cut-off value of 81.8 μg/L has recently been reported by a large prospective multicenter study . To our knowledge, the highest NSE cut-off value reported is 97 μg/L . We therefore applied this limit to the current study and have implemented it in our diagnostic pathway .
For 67 out of 98 patients results of SSEP were available. SSEPs were performed with stimulation of the median nerve at the wrist as described earlier . Results were classified as bilateral absent N20 or at least unilateral present. The majority of SSEPs were performed between 48 and 96 hours after cardiac arrest.
In 62 out of 98 patients, all three parameters (GWR, NSE and SSEP) were available. In this patient cohort, we determined sensitivity and specificity for poor outcome prediction for all three parameters separately as well as for the combination of SSEP and NSE (predicting poor outcome if both NSE > 97 μg/L and bilateral absent N20) and for the combination of SSEP, NSE and GWR (predicting poor outcome by at least two of three parameters, NSE > 97 μg/L, bilateral absent N20, GWR < 1.16).
We calculated the correlation coefficients for attenuation measurements determined by two different raters (MS and AG) for GWR measures to assess inter-rater reliability. Absolute attenuation values and GWR values were compared between patients with good outcome (CPC 1–2) and poor outcome (CPC 3–5) by means of a Mann-Whitney-U test. The central aim of the study was to determine whether a cut-off level of GWR exists that allows prediction of poor outcome with specificity close to 100%. Additionally, we conducted ROC analyses for the three different GWR measurements (GWR-BG, GWR-CE, GWR-AV) to illustrate the overall performance of GWR in discrimination of patients with good and poor outcome.
To evaluate a possible relationship of GWR with time between cardiac arrest and CCT, we divided all patients with poor outcome into three groups (CCT <6 h, 6–24 h, >24 h after cardiac arrest). GWR values of these three groups were compared by using a Kruskal-Wallis test.
Association of GWR and clinical outcome
Absolute attenuation in Hounsfield Units and GWR calculations
n = 37
n = 61
Single parameter approach-sensitivity and negative predictive values (at 100% specificity) with respect to poor outcome prediction
No. of patients
GWR AV < 1.16
GWR AV ≥ 1.16
NSE > 97 μg/L
NSE ≤ 97 μg/L
ROC analysis demonstrated an area under the curve (AUC) for GWR-AV = 0.81, GWR-BG = 0.77 and GWR-CE = 0.75 (Additional file 1: Figure S1). Therefore we used GWR-AV for all subsequent analysis and refer to it as GWR in the following.
GWR < 1.16 predicted poor outcome with a sensitivity of 38% and a specificity of 100%. Applying a GWR threshold of 1.2 would have misclassified 6 patients (CPC 1/2 = 2/4 patients).
Comparison of GWR, NSE and SSEP
Table 3 gives an overview on sensitivity and negative predictive value of GWR, NSE and SSEP for the whole study cohort. In our study, 62 Patients had received CCT, NSE and SSEP. In this subgroup, poor outcome was detected with a sensitivity of 53% by NSE > 97 μg/L, 55% by bilaterally absent SSEP and 35% by GWR < 1.16. No patient with bilateral absent SSEP, NSE > 97 μg/L or GWR < 1.16 survived with good outcome. NSE and SSEP detected poor outcome for only partly overlapping patient populations: 65% of patients (26/40) with poor outcome had absent SSEP and/or NSE > 97 μg/L and 70% of patients (28/40) with poor outcome had at least one of three parameters (SSEP, NSE, GWR) suggesting poor outcome. Hence, GWR only marginally increased sensitivity for poor outcome detection in a prediction approach that requires only one positive parameter.
Single and multi parameter approach-sensitivity and negative predictive values (at 100% specificity) for poor outcome prediction in the subgroup of 62 patients for which NSE, SSEP and CCT were available
No. of patients
GWR AV < 1.16
GWR AV ≥ 1.16
NSE > 97 μg/L
NSE ≤ 97 μg/L
SSEP absent and/or NSE > 97 μg/L
SSEP detected and NSE < 97 μg/L
SSEP absent and NSE > 97 μg/L
SSEP detected or NSE < 97 μg/L
2 or 3 abnormal (SSEP, NSE, GWR)
1 or 0 abnormal (SSEP, NSE, GWR)
One or more parameters pathologic
All parameters normal
Our study suggests that a gray-whitematter ratio (GWR) below 1.16 may be a specific parameter for prediction of poor outcome in patients treated with hypothermia. This corroborates previous findings that were obtained in normothermia and mixed patient cohorts (normothermia and hypothermia). None of 23 patients with a GWR < 1.16 survived with good outcome.
The majority of patients with GWR < 1.16 had at least one other parameter indicative of poor outcome (NSE > 97 μg/L and/or bilateral absent SSEP). Therefore, GWR did not substantially increase the sensitivity for poor outcome detection, if one accepts poor outcome prediction by a single parameter. However, because individual patients may survive despite poor outcome prediction by single parameters, it has been argued that poor outcome prediction should be based on a multi-parameter approach . We demonstrated in a subgroup of 62 patients who had received NSE, SSEP and CCT that a GWR < 1.16 increased the sensitivity of poor outcome detection (from 43% (17/40) to 53% (n = 21/40)) if at least two pathologic test results were required for poor outcome prediction. These results suggest that GWR is a useful additional parameter that may increase the level of certainty for poor outcome prediction in comatose patients after cardiac arrest.
GWR cut-off value and timing of CCT
Our cut-off value of 1.16 is similar to previously published values. Torbey et al. found a GWR < 1.18 to be 100% specific of poor outcome . The authors obtained GWR at the basal ganglia level in normothermic patients. None of the eleven patients with GWR < 1.18 survived. In the largest cohort so far by Metter et al., a GWR < 1.15 was found to be 100% specific for poor outcome . CCT was obtained within 24 hours and the majority of patients with GWR < 1.15 had been treated with hypothermia. Metter et al. concluded that a GWR < 1.2 had the best overall predictive value, but this threshold misclassified two patients. A of GWR < 1.2 in our patient cohort would have misclassified six patients as having poor outcome. Hence, our data suggest that a threshold of 1.2 can be associated with a substantial number of erroneous poor outcome predictions.
In the study by Metter et al. the number of patients who survived with little or no neurological deficit was small . This outcome distribution may have decreased the likelihood of detecting rare cases that survive with good neurological outcome despite a low GWR. In our study, the proportion of patients surviving with little or no neurological deficit was higher (38%, n = 37).
Inamasu et al. reported a high sensitivity and specificity of the ‘loss of boundary’ (LOB) sign and ‘sulcal effacement’ sign obtained prospectively very early after cardiac arrest . Only one of 52 patients with LOB sign and none of 20 patients with sulcal effacement sign survived with good outcome. These data corroborate the high specificity of early CT signs for poor outcome prediction.
Wu et al. investigated a semi-automated procedure, which used co-registration of CT to a brain atlas and determination of Hounsfield units for several brain regions . In a large cohort, the authors demonstrate that whole brain Hounsfield units obtained from early CCT (most within 24 hours) improved the predictive accuracy of a clinical examination. In this study, only 12% of patients had received hypothermia treatment, and therefore results should be interpreted with care for this patient group.
The majority of CCTs in our study were performed on admission, frequently within a few hours after cardiac arrest. Metter et al. did not find a significant association of attenuation with time since cardiac arrest within a narrow time frame of 24 hours . The distribution of GWR over time in our study, however, suggests that GWR may decrease over the days following cardiac arrest for patients with severe hypoxic-ischemic encephalopathy (Figure 2). Similarly, Wu et al. found a significant decrease in the density of the putamen for CT performed within 24–72 hours as compared to within the first 24 hours after cardiac arrest . No difference was found for other regions, but only few CTs (4%) were performed after day two, hence the power to detect a difference was low.
In our study, sensitivity of poor outcome detection by GWR < 1.16 was 23% for CCT performed within six hours, but increased to 60% for CCT performed between six hours and seven days. This finding, however, may be biased by patient selection and should be interpreted with care. Nonetheless our data argue for a prospective study investigating the prognostic value of GWR beyond six hours after cardiac arrest.
Since many CTs in our study were performed within the first few hours following admission, it is unlikely that hypothermia had substantially influenced the GWR in this subset of patients. However, as hypothermia changes the course of pathophysiological events leading to definite tissue damage in the days after cardiac arrest, it is still possible that hypothermia influences the relationship between early GWR and neurological outcome. Survival with good outcome may be possible in hypothermia patients despite a more severe early injury. Hence, re-evaluation of the prognostic value of GWR for hypothermia patients is mandatory, before GWR can be integrated into pathways for management of these patients.
Comparison with other outcome parameters
The lower sensitivity of GWR < 1.16 (35% as compared to 53% for NSE and 55% SSEP) and the fact that the vast majority of patients with GWR < 1.16 had either high NSE levels or bilateral absent SSEP may suggest that GWR < 1.16 identifies patients with more severe hypoxic-ischemic encephalopathy than NSE or SSEP. Most of the CCTs in our study were performed within a few hours, but SSEP and NSE were obtained around day three after cardiac arrest. The rather low sensitivity of GWR for the entire study cohort could therefore also be explained by the early time point. Indeed, our data suggest that the sensitivity of GWR increases substantially if CCT is performed > 24 h after admission.
If one accepts a prediction approach that uses a single parameter for poor outcome prediction, the combination of NSE and SSEP increases the sensitivity for detection of poor outcome. Apparently, a relevant number of patients with poor outcome are detected by NSE but not by SSEP and vice versa. If NSE and SSEP are performed, GWR < 1.16 did not significantly increase the number of patients for which poor outcome was detected. However, if one prefers a prediction approach that requires a second confirmatory parameter, GWR has an additional predictive value.
Limitations of our study
Due to the retrospective design of our study, we could only investigate the subset of patients who had received CCT (n = 111, i.e. 31% of all cardiac arrest patients). CCT was generally ordered by treating physicians if the cause of the cardiac arrest was deemed uncertain or if intracranial complications were suspected. Our results should therefore be extended to the entire cardiac arrest population with care and corroboration by a prospective study avoiding selection bias is desirable. The same caveat applies to the comparison of different prognostic parameters, which could only be performed in the subset of patients who had received NSE, SSEP and CCT.
The results of NSE, SEP and CCT were available to the treating physicians and influenced decisions on limitation of treatment. Therefore, we cannot exclude the possibility of a self-fulfilling prophecy. GWR was calculated for this study and not known to treating physicians. However, reduced GWR reflects brain edema and treatment decisions were influenced by CT if brain edema was reported by the radiologist. As a number of recent studies have indicated limited reliability of prognostic parameters in hypothermia patients, treatment withdrawal within the first days after cardiac arrest is restricted to patients with confirmed brain death at our institution. In the remainder, treatment is continued up to day seven . Hence, a major effect of a self-fulfilling prophecy from treatment withdrawal within the first days after cardiac arrest is unlikely in our study. Nonetheless, unexpected late recovery may occur in individual patients , and we cannot fully rule out that treatment withdrawal based on over-interpretation of prognostic parameters has prevented the detection of such exceptional cases.
Our study confirms the strong association of a low GWR with poor outcome. The optimal cut-off value of GWR was 1.16 in our study, which predicted poor outcome with 100% specificity and 38% sensitivity in our cohort of hypothermia treated patients. Most patients with GWR < 1.16 had absent SSEP or NSE > 97 μg/L. Hence, GWR did not substantially increase the sensitivity to predict poor outcome if these two parameters were available. However, GWR may be useful in a multi-parameter approach and may increase the certainty with which poor outcome can be predicted early after cardiac arrest in patients treated with hypothermia. Further prospective studies are needed to corroborate GWR thresholds and to evaluate the optimal timing of CCT for outcome prediction in cardiac arrest survivors.
M. Scheel is supported by the Friedrich C. Luft Clinical Scientist Pilot Program funded by Volkswagen Foundation and Charité Foundation.
- Al Thenayan E, Savard M, Sharpe M, Norton L, Young B: Predictors of poor neurologic outcome after induced mild hypothermia following cardiac arrest. Neurology. 2008, 71: 1535-1537. 10.1212/01.wnl.0000334205.81148.31.View ArticlePubMedGoogle Scholar
- Steffen IG, Hasper D, Ploner CJ, Schefold JC, Dietz E, Martens F, Nee J, Krueger A, Jörres A, Storm C: Mild therapeutic hypothermia alters neuron specific enolase as an outcome predictor after resuscitation: 97 prospective hypothermia patients compared to 133 historical non-hypothermia patients. Crit Care. 2010, 14: R69-10.1186/cc8975.PubMed CentralView ArticlePubMedGoogle Scholar
- Daubin C, Quentin C, Allouche S, Etard O, Gaillard C, Seguin A, Valette X, Parienti J-J, Prevost F, Ramakers M, Terzi N, Charbonneau P, Du Cheyron D: Serum neuron-specific enolase as predictor of outcome in comatose cardiac-arrest survivors: a prospective cohort study. BMC Cardiovasc Disord. 2011, 11: 48-10.1186/1471-2261-11-48.PubMed CentralView ArticlePubMedGoogle Scholar
- Leithner C, Ploner CJ, Hasper D, Storm C: Does hypothermia influence the predictive value of bilateral absent N20 after cardiac arrest?. Neurology. 2010, 74: 965-969. 10.1212/WNL.0b013e3181d5a631.View ArticlePubMedGoogle Scholar
- Bender A, Howell K, Frey M, Berlis A, Naumann M, Buheitel G: Bilateral loss of cortical SSEP responses is compatible with good outcome after cardiac arrest. J Neurol. 2012Google Scholar
- Bouwes A, Binnekade JM, Kuiper MA, Bosch FH, Zandstra DF, Toornvliet AC, Biemond HS, Kors BM, Koelman JHTM, Verbeek MM, Weinstein HC, Hijdra A, Horn J: Prognosis of coma after therapeutic hypothermia: a prospective cohort study. Ann Neurol. 2012, 71: 206-212. 10.1002/ana.22632.View ArticlePubMedGoogle Scholar
- Samaniego EA, Mlynash M, Caulfield AF, Eyngorn I, Wijman CAC: Sedation confounds outcome prediction in cardiac arrest survivors treated with hypothermia. Neurocrit Care. 2011, 15: 113-119. 10.1007/s12028-010-9412-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Rossetti AO, Oddo M, Liaudet L, Kaplan PW: Predictors of awakening from postanoxic status epilepticus after therapeutic hypothermia. Neurology. 2009, 72: 744-749. 10.1212/01.wnl.0000343006.60851.62.View ArticlePubMedGoogle Scholar
- Kjos BO, Brant-Zawadzki M, Young RG: Early CT findings of global central nervous system hypoperfusion. AJR Am J Roentgenol. 1983, 141: 1227-1232. 10.2214/ajr.141.6.1227.View ArticlePubMedGoogle Scholar
- Han BK, Towbin RB, De Courten-Myers G, McLaurin RL, Ball WS: Reversal sign on CT: effect of anoxic/ischemic cerebral injury in children. AJNR Am J Neuroradiol. 1989, 10: 1191-1198.PubMedGoogle Scholar
- Torbey MT, Selim M, Knorr J, Bigelow C, Recht L: Quantitative analysis of the loss of distinction between gray and white matter in comatose patients after cardiac arrest. Stroke. 2000, 31: 2163-2167. 10.1161/01.STR.31.9.2163.View ArticlePubMedGoogle Scholar
- Torbey MT, Geocadin R, Bhardwaj A: Brain arrest neurological outcome scale (BrANOS): predicting mortality and severe disability following cardiac arrest. Resuscitation. 2004, 63: 55-63. 10.1016/j.resuscitation.2004.03.021.View ArticlePubMedGoogle Scholar
- Choi SP, Park HK, Park KN, Kim YM, Ahn KJ, Choi KH, Lee WJ, Jeong SK: The density ratio of grey to white matter on computed tomography as an early predictor of vegetative state or death after cardiac arrest. Emerg Med J. 2008, 25: 666-669. 10.1136/emj.2007.053306.View ArticlePubMedGoogle Scholar
- Yanagawa Y, Un-no Y, Sakamoto T, Okada Y: Cerebral density on CT immediately after a successful resuscitation of cardiopulmonary arrest correlates with outcome. Resuscitation. 2005, 64: 97-101. 10.1016/j.resuscitation.2004.06.015.View ArticlePubMedGoogle Scholar
- Wu O, Batista LM, Lima FO, Vangel MG, Furie KL, Greer DM: Predicting clinical outcome in comatose cardiac arrest patients using early noncontrast computed tomography. Stroke. 2011, 42: 985-992. 10.1161/STROKEAHA.110.594879.PubMed CentralView ArticlePubMedGoogle Scholar
- Inamasu J, Miyatake S, Suzuki M, Nakatsukasa M, Tomioka H, Honda M, Kase K, Kobayashi K: Early CT signs in out-of-hospital cardiac arrest survivors: Temporal profile and prognostic significance. Resuscitation. 2010, 81: 534-538. 10.1016/j.resuscitation.2010.01.012.View ArticlePubMedGoogle Scholar
- Metter RB, Rittenberger JC, Guyette FX, Callaway CW: Association between a quantitative CT scan measure of brain edema and outcome after cardiac arrest. Resuscitation. 2011, 82: 1180-1185. 10.1016/j.resuscitation.2011.04.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Iyer RS, Thomas B: Teaching neuroimages: reversal sign on CT in a child with hypoxic-ischemic brain injury. Neurology. 2011, 76: e27-10.1212/WNL.0b013e31820af92b.View ArticlePubMedGoogle Scholar
- Leithner C, Storm C, Hasper D, Ploner CJ: Prognosis for Brain Function Following Cardiopulmonary Resuscitation and Therapeutic Hypothermia. Aktuelle Neurologie. 2012, 39: 145-154.View ArticleGoogle Scholar
- Jennett B, Bond M: Assessment of outcome after severe brain damage. Lancet. 1975, 1: 480-484.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.