- Open Access
- Open Peer Review
Impact of contractile reserve on acute response to cardiac resynchronization therapy
© Moonen et al; licensee BioMed Central Ltd. 2008
- Received: 30 October 2008
- Accepted: 31 December 2008
- Published: 31 December 2008
Cardiac resynchronization therapy (CRT) provides benefit for congestive heart failure, but still 30% of patients failed to respond to such therapy. This lack of response may be due to the presence of significant amount of scar or fibrotic tissue at myocardial level. This study sought to investigate the potential impact of myocardial contractile reserve as assessed during exercise echocardiography on acute response following CRT implantation.
Fifty-one consecutive patients with heart failure (LV ejection fraction 27% ± 5%, 67% ischemic cardiomyopathy) underwent exercise Doppler echocardiography before CRT implantation to assess global contractile reserve (improvement in LV ejection fraction) and local contractile reserve in the region of the LV pacing lead (assessed by radial strain using speckle tracking analysis). Responders were defined by an increase in stroke volume ≥ 15% after CRT.
Compared with nonresponders, responders (25 patients) showed a greater exercise-induced increase in LV ejection fraction, a higher degree of mitral regurgitation and a significant extent of LV dyssynchrony. The presence of contractile reserve was directly related to the acute increase in stroke volume (r = 0.48, p < 0.001). Baseline myocardial deformation as well as contractile reserve in the LV pacing lead region was greater in responders during exercise than in nonresponders (p < 0.0001).
The present study showed that response to CRT largely depends not only on the extent of LV dyssynchrony and the severity of mitral regurgitation but also on the presence of contractile reserve.
- Leave Ventricular Ejection Fraction
- Mitral Regurgitation
- Cardiac Resynchronization Therapy
- Biventricular Pace
- Leave Ventricular Lead
The number of patients presenting with heart failure is increasing rapidly. Currently, cardiac resynchronization therapy (CRT) has emerged as an adjunctive treatment in patients with drug-refractory heart failure . The clinical benefit of CRT – evidenced by improvement in symptoms, quality of life, exercise capacity, and left ventricular (LV) systolic performance – has been largely demonstrated [2, 3]. Acutely, a CRT-increase in forward stroke volume by ≥ 15% has been shown to predict long-term clinical outcome in both ischemic and non ischemic cardiomyopathy . However, 20% to 30% of patients do not respond to CRT despite application of recommended selection criteria . Therefore, additional selection criteria are needed. Response to CRT largely depends on the extent of LV dyssynchrony, the severity of of mitral regurgitation and the possibility offers to the LV to recruit function [6–9]. CRT response could thus correlate with myocardial viability and inversely be associated with the extent of myocardial fibrosis [10, 11]. Furthermore, scar tissue or fibrosis – characterized by lack of contractile reserve – in the LV pacing lead region may prohibit response to CRT. Myocardial contractile reserve – important prognostic marker in heart failure – can be reliably assessed during exercise echocardiography . This study sought to investigate the potential impact of myocardial contractile reserve as assessed during exercise echocardiography on acute response following CRT implantation.
Consecutive patients (n = 51) with severe heart failure, scheduled for implantation of a permanent biventricular pacemaker, were included using established selection criteria for CRT: New York Heart Association class III or IV despite optimal medical therapy, LV ejection fraction ≤ 35%, QRS with left bundle branch block configuration or a duration >120 ms. The aetiologies of heart failure were idiopathic dilated cardiomyopathy in 17 patients and ischaemic heart disease in 34. All patients underwent exercise Doppler echocardiography before CRT implantation to assess global LV contractile reserve and local contractile reserve in the region of the LV pacing lead (assessed by longitudinal strain using speckle tracking analysis). No patients with ischaemic disease presented inducible ischaemia during test. The protocol was approved by the Human Ethical Committee of our University Hospital and all patients gave informed consent.
A symptom-limited graded bicycle exercise test was performed in a semi-supine position on a tilting exercise table. After an initial workload of 25 W maintained for 2 min, the workload was increased every 2 min by 25 W. Blood pressure and a 12-lead electrocardiogram were recorded every 2 min. Two-dimensional echocardiographic recordings were made throughout the test.
All echocardiographic parameters were obtained at rest and at peak exercise in the same cycling semi-supine position (Vivid 7 imaging device, GE Healthcare, United Kingdom) and were obtained in digital format and stored on optical disks for off-line analysis. Two-dimensional grayscale images (frame rates ≥ 70 s-1) and colour-tissue Doppler imaging (frame rates ≥ 115 s-1) were performed in the apical views using a narrow sector angle. Left ventricular end-diastolic and end-systolic volumes and ejection fraction were measured by the biapical Simpson disk method. Wall motion score analysis was applied to a 16-segment model of the LV with a semiquantitative scoring system (1 = normal, 2 = hypokinesia, 3 = akinesia, 4 = dyskinesia). Quantitation of MR was performed by the proximal isovelocity surface area method as previously described . Intra-LV dyssynchrony was referred to LV dispersion and was calculated as the difference between the longest and the shortest times to peak systolic velocities of 4 opposing basal walls (apical 2-and 4-chamber) . Global LV contractile reserve was expressed as the change in LV ejection fraction from rest to peak exercise. Regional contractile reserve in the pacing lead was assessed using speckle tracking analysis from LV long axis images (mid segments). After tracing the endocardial borders in the end-systolic frame, an automated tracking algorithm outlined the myocardial deformation in the dedicated LV segments . Peak systolic radial strain was only measured in the lateral, posterior and anterior regions, where the LV lead was positioned.
CRT implantation and setting
All patients received a biventricular pacing device for CRT with a right ventricular apical lead and LV pacing electrodes positioned through the coronary sinus in a LV epicardial vein. This coronary sinus lead was placed in a lateral position in 32 patients, in a postero-lateral position in 16 and in anterior position in 3. After a successful implant, echocardiography was used to optimize the atrio-ventricular delay in order to maximize LV filling time. Interventricular pacing interval was set to default value (V-V = 0 ms).
CRT response definition
Within the week after implantation, changes in forward stroke volume were assessed. Improvement was defined as an increase in stroke volume ≥ 15% after CRT  (see videos).
Data are expressed as mean ± 1 SD (STATISTICA version 6). Student's paired 2-tailed t test was used to compare measurements obtained at rest and during exercise. Categorical variables were compared with Fisher exact test. A value of p < 0.05 was considered significant. Categorical data are summarized as frequencies and percentages. Linear regression analysis was performed to evaluate the relation between the improvement in LV ejection fraction during exercise and the percentages of changes in forward stroke volume after CRT implantation. The optimal improvement in LV ejection fraction and in myocardial strain during exercise to predict response to CRT was determined by receiver operating characteristic curve analysis. To detect independent cofactors associated with response to CRT, a logistic multivariate analysis was performed. All significant variables were included in the multivariate model. Reproducibility of echocardiographic measurements has been previously published .
Responders versus non responders
Responders versus non responders
(n = 25, 49%)
(n = 26, 51%)
71 ± 8
69 ± 9
Male, n (%)
Ischemic cardiomyopathy, n (%)
QRS duration, ms
154 ± 24
169 ± 25
Diuretic, n (%)
β-Blockers, n (%)
ACEi, n (%)
AR Blockers, n (%)
Spironolactone, n (%)
LV-RV dyssynchrony, ms
54 ± 16
50 ± 17
LV dispersion, ms
109 ± 44
61 ± 35
Mitral effective regurgitant orifice, mm2
22 ± 11
15 ± 9
LV end-diastolic volume, ml
175 ± 34
193 ± 42
LV end-systolic volume, ml
128 ± 26
139 ± 35
LV ejection fraction, %
27.1 ± 5.0
27.7 ± 5.3
LV ejection fraction at exer., %
35 ± 5.7
31 ± 5.3
LV ejection fraction diff., %
8.4 ± 2.4
4.2 ± 2.6
Wall motion score index
2.2 ± 0.27
2.28 ± 0.35
Wall motion score index at exer
1.86 ± 0.34
2.15 ± 0.34
Wall motion score index diff
-0.3 ± 0.26
-0.13 ± 0.21
Strain target LV lead wall, %
15 ± 2.8
10.2 ± 6.3
Strain target LV lead wall at exer, %
19.3 ± 2.4
10.7 ± 8.5
Strain target LV lead wall diff., %
4.2 ± 1.6
0.44 ± 3.5
Area under the curves, sensitivity, specificity and optimal cutoff values of predictors for acute CRT response.
Data at inclusion
> 65 ms
Effective regurgitant orifice
> 15 mm2
Changes in LV ejection fraction at exercise
Changes in wall motion score index at exercise
Regional contractile reserve and response to CRT
The response to CRT is likely modulated by multiple factors. The present study confirms that response to CRT largely depends not only on the extent of LV dyssynchrony and the severity of mitral regurgitation but also on the possibility offers to the LV to recruit function, namely the presence of contractile reserve. The probability to acutely respond to CRT is highly likely in patients with moderate to severe mitral regurgitation, major LV dyssynchrony and significant contractile myocardial recruitment – global and in the LV lead target site – at exercise.
Role of LV dyssynchrony
Several studies have demonstrated that the major predictor of responsiveness to CRT is mechanical rather than electrical dyssynchrony [6, 13, 16, 17]. Among available techniques, tissue Doppler imaging has emerged as the most practical method for assessing LV dyssynchrony. As reported previously, we found that a significant mechanical delay between the basal segments of the LV, mainly the septum and the lateral wall, on tissue Doppler was highly predictive for response to CRT [6, 18]. Patients with a delay ≥ 65 ms showed a significant improvement in forward stroke volume under CRT. However, although this beneficial effect was rarely observed in patients with a delay < 65 ms, 34% (9 of 26) of the patients with severe LV dyssynchrony did not respond to CRT. This emphasizes that other factors play a role; LV dyssynchrony is necessary but not sufficient for CRT response. Indeed, delayed wall motion is mainly a marker of myocardial dysfunction and could be exhibited in viable and scar tissue according to loading conditions [19, 20].
Effect of mitral regurgitation
Response to CRT might be modulated by the presence of functional mitral regurgitation before implantation. Few authors have shown that patients with severe mitral regurgitation have less chance of improving with CRT . However, these studies included a limited number of patients. In the CARE-HF study, a randomized trial including a large number of patients, it was, conversely, shown that patients who were not improved were likely to have no significant mitral regurgitation as compared with responders . The results of the present study confirm and extend this observation. Indeed, we found that a mitral effective regurgitant orifice > 15 mm2 was associated with a higher likelihood of response to CRT. Patients with more severe mitral regurgitation have more LV asynchrony (98 ± 49 ms vs 67 ± 36 ms; p = 0.019). The interplay between LV dyssynchrony and mitral regurgitation is complex [4, 21–23]. In fact, the degree of mitral regurgitation is directly correlated to the asynchronous activation of papillary muscles and the decrease in mitral valve closing forces (mainly the consequence of global LV dyssynchrony). It could thus be expected that the beneficial effect of CRT would be greater in patients with significant mitral regurgitation.
Myocardial reserve and response to CRT
In patients with depressed LV function, the identification of contractile reserve during stress echocardiography has been shown to provide important prognostic information in heart failure patients [24, 25]. More specifically, in patients referred to CRT, few authors have reported that the presence of residual myocardial viability might modulate the response to CRT [11, 14, 26–30]. Contractile reserve is a specific marker of underlying myocardial viability, which can be reliably assessed during semi-supine exercise echocardiography . To the best of our knowledge, this study was the first to examine the value of exercise echocardiography for predicting response to CRT. Significant contractile myocardial recruitment – global and in the LV lead target site – characterized patients who presented the greater rise in forward stroke volume. An increase in LV ejection fraction by ≥ 6.7% and/or an increase in local strain by ≥ 2% were found to be predictive for response to CRT. Moreover, the extent of LV global contractile reserve was related to the percentage of acute changes in forward stroke volume. These data confirm and extend previous studies by demonstrating that a substantial amount of recruitable myocardium is needed to obtain improvement in LV function after CRT. It could be thus argued that in case of advanced myocardial remodeling process, fibrosis and loss of contractile material may severely altered myocardial conduction and contractile properties, which might in turn impede efficient biventricular pacing .
Some limitations should be acknowledged. Our data pertain only to patients with current criteria for CRT implantation. The population studied was not completely homogenous since it was composed of patients with myocardial dysfunction of ischaemic and non ischaemic origin. However, this represents our daily patients referred for CRT. Moreover, biventricular pacing is still a challenging therapy in both settings.
Heart failure patients referred to CRT have less chance of improving under therapy if they have no significant mitral regurgitation, no LV dyssynchrony and no contractile myocardial recruitment – global and in the LV lead target site – at exercise. It might thus be proposed that the assessment of these parameters would be routinely performed before submitting the patients to CRT.
- Vardas PE, Auricchio A, Blanc JJ, Daubert JC, Drexler H, Ector H, Gasparini M, Linde C, Morgado FB, Oto A, Sutton R, Trusz-Gluza M: Guidelines for cardiac pacing and cardiac resynchronization therapy: The Task Force for Cardiac Pacing and Cardiac Resynchronization Therapy of the European Society of Cardiology. Eur Heart J. 2007, 28: 2256-2295. 10.1093/eurheartj/ehm305.View ArticlePubMedGoogle Scholar
- Tournoux FB, Alabiad C, Fan D, Chen AA, Chaput M, Heist EK, Mela T, Mansour M, Reddy V, Ruskin JN, Picard MH, Singh JP: Echocardiographic measure of acute haemodynamic response after cardiac resynchronization therapy predicts long term clinical outcome. Eur Heart J. 2007, 28: 1143-1148. 10.1093/eurheartj/ehm050.View ArticlePubMedGoogle Scholar
- Cleland JG, Daubert JC, Erdmann E, Freemantle N, Gras D, Kappenberger L, Tavazzi L, Cardiac Resynchronization-Heart Failure (CARE-HF) Study Investigators: The effect of cardiac resynchronization on morbidity and mortality in heart failure. N Engl J Med. 2005, 352: 1539-1549. 10.1056/NEJMoa050496.View ArticlePubMedGoogle Scholar
- Lancellotti P, Mélon P, Sakalihasan N, Waleffe A, Dubois C, Bertholet M, Piérard LA: Effect of cardiac resynchronization therapy on functional mitral regurgitation in heart failure. Am J Cardiol. 2004, 94: 1462-5. 10.1016/j.amjcard.2004.07.154.View ArticlePubMedGoogle Scholar
- Abraham WT, Fisher WG, Smith AL, Delurgio DB, Leon AR, Loh E, Kocovic DZ, Packer M, Clavell AL, Hayes DL, Ellestad M, Trupp RJ, Underwood J, Pickering F, Truex C, McAtee P, Messenger J, MIRACLE Study Group. Multicenter InSync Randomized Clinical Evaluation: Cardiac resynchronization in chronic heart failure. N Engl J Med. 2002, 346: 1845-1853. 10.1056/NEJMoa013168.View ArticlePubMedGoogle Scholar
- Yu CM, Gorcsan J, Bleeker GB, Zhang Q, Schalij MJ, Suffoletto MS, Fung JW, Schwartzman D, Chan YS, Tanabe M, Bax JJ: Usefulness of tissue Doppler velocity and strain dyssynchrony for predicting left ventricular reverse remodeling response after cardiac resynchronization therapy. Am J Cardiol. 2007, 00 (8): 1263-1270. 10.1016/j.amjcard.2007.05.060.View ArticleGoogle Scholar
- Díaz-Infante E, Mont L, Leal J, García-Bolao I, Fernández-Lozano I, Hernández-Madrid A, Pérez-Castellano N, Sitges M, Pavón-Jiménez R, Barba J, Cavero MA, Moya JL, Pérez-Isla L, Brugada J, SCARS Investigators: Predictors of lack of response to resynchronization therapy. Am J Cardiol. 2005, 95: 1436-40. 10.1016/j.amjcard.2005.02.009.View ArticlePubMedGoogle Scholar
- Reuter S, Garrigue S, Barold SS, Jais P, Hocini M, Haissaguerre M, Clementy J: Comparison of characteristics in responders versus nonresponders with biventricular pacing for drug-resistant congestive heart failure. Am J Cardiol. 2002, 89: 346-50. 10.1016/S0002-9149(01)02240-8.View ArticlePubMedGoogle Scholar
- Richardson M, Freemantle N, Calvert MJ, Cleland JG, Tavazzi L, CARE-HF Study Steering Committee and Investigators: Predictors and treatment response with cardiac resynchronization therapy in patients with heart failure characterized by dyssynchrony: a pre-defined analysis from the CARE-HF trial. Eur Heart J. 2007, 28: 1827-34. 10.1093/eurheartj/ehm192.View ArticlePubMedGoogle Scholar
- Bleeker GB, Kaandorp TA, Lamb HJ, Boersma E, Steendijk P, de Roos A, Wall van der EE, Schalij MJ, Bax JJ: Effect of posterolateral scar tissue on clinical and echocardiographic improvement after cardiac resynchronization therapy. Circulation. 2006, 113: 969-976. 10.1161/CIRCULATIONAHA.105.543678.View ArticlePubMedGoogle Scholar
- Lim P, Bars C, Mitchell-Heggs L, Roiron C, Elbaz N, Hamdaoui B, Lellouche N, Dubois-Randé JL, Guéret P: Importance of contractile reserve for CRT. Europace. 2007, 9: 739-743. 10.1093/europace/eum117.View ArticlePubMedGoogle Scholar
- Lancellotti P, Hoffer EP, Piérard LA: Detection and clinical usefulness of a biphasic response during exercise echocardiography early after myocardial infarction. Am Coll Cardiol. 2003, 1 (7): 1142-7. 10.1016/S0735-1097(03)00049-4.View ArticleGoogle Scholar
- Bax JJ, Bleeker GB, Marwick TH, Molhoek SG, Boersma E, Steendijk P, Wall van der EE, Schalij MJ: Left ventricular dyssynchrony predicts response and prognosis after cardiac resynchronization therapy. J Am Coll Cardiol. 2004, 44: 1834-40. 10.1016/j.jacc.2004.08.016.View ArticlePubMedGoogle Scholar
- Ypenburg C, Sieders A, Bleeker GB, Holman ER, Wall van der EE, Schalij MJ, Bax JJ: Myocardial contractile reserve predicts improvement in left ventricular function after cardiac resynchronization therapy. Am Heart J. 2007, 154: 1160-1165. 10.1016/j.ahj.2007.07.035.View ArticlePubMedGoogle Scholar
- Lancellotti P, Cosyns B, Piérard LA: Dynamic left ventricular dyssynchrony contributes to B-type natriuretic peptide release during exercise in patients with systolic heart failure. Europace. 2008, 10: 496-501. 10.1093/europace/eun040.View ArticlePubMedGoogle Scholar
- Pitzalis MV, Iacoviello M, Romito R, Luzzi G, Anaclerio M, Forleo C: Ventricular asynchrony predicts a better outcome in patients with chronic heart failure receiving cardiac resynchronization therapy. J Am Coll Cardiol. 2005, 45: 65-69. 10.1016/j.jacc.2004.09.058.View ArticlePubMedGoogle Scholar
- Chung ES, Leon AR, Tavazzi L, Sun JP, Nihoyannopoulos P, Merlino J, Abraham WT, Ghio S, Leclercq C, Bax JJ, Yu CM, Gorcsan J, St John Sutton M, De Sutter J, Murillo J: Results of the Predictors of Response to CRT (PROSPECT) trial. Circulation. 2008, 117: 2608-16. 10.1161/CIRCULATIONAHA.107.743120.View ArticlePubMedGoogle Scholar
- Notabartolo D, Merlino JD, Smith AL, DeLurgio DB, Vera FV, Easley KA, Martin RP, León AR: Usefulness of the peak velocity difference by tissue Doppler imaging technique as an effective predictor of response to cardiac resynchronization therapy. Am J Cardiol. 2004, 94: 817-20. 10.1016/j.amjcard.2004.05.072.View ArticlePubMedGoogle Scholar
- Lancellotti P, Stainier PY, Lebois F, Piérard LA: Effect of dynamic left ventricular dyssynchrony on dynamic mitral regurgitation in patients with heart failure due to coronary artery disease. Am J Cardiol. 2005, 96: 1304-7. 10.1016/j.amjcard.2005.06.077.View ArticlePubMedGoogle Scholar
- Ryf S, Rutz AK, Boesiger P, Schwitter J: Is post-systolic shortening a reliable indicator of myocardial viability? An MR tagging and late-enhancement study. J Cardiovasc Magn Reson. 2006, 8: 445-51. 10.1080/10976640600599494.View ArticlePubMedGoogle Scholar
- Breithardt OA, Sinha AM, Schwammenthal E, Bidaoui N, Markus KU, Franke A, Stellbrink C: Acute Effects of Cardiac Resynchronization Therapy on Functional Mitral Regurgitation in Advanced Systolic Heart Failure. J Am Coll Cardiol. 2003, 41: 765-770. 10.1016/S0735-1097(02)02937-6.View ArticlePubMedGoogle Scholar
- Kanzaki H, Bazaz R, Schwartzman D, Dohi K, Sade LE, Gorcsan J: A Mechanism for Immediate Reduction in Mitral Regurgitation After Cardiac Resynchronization Therapy. J Am Coll Cardiol. 2004, 44: 1619-1625. 10.1016/j.jacc.2004.07.036.View ArticlePubMedGoogle Scholar
- Agricola E, Oppizzi M, Galderisi M, Pisani M, Meris A, Pappone C, Margonato A: Role of regional mechanical dyssynchrony as a determinant of functional mitral regurgitation in patients with left ventricular systolic dysfunction. Heart. 2006, 92: 1390-1395. 10.1136/hrt.2005.082115.View ArticlePubMedPubMed CentralGoogle Scholar
- Bax JJ, Wijns W, Cornel JH, Visser FC, Boersma E, Fioretti PM: Accuracy of currently available techniques for prediction of functional recovery after revascularization in patients with left ventricular dysfunction due to chronic coronary artery disease: comparison of pooled data. J Am Coll Cardiol. 1997, 30: 1451-60. 10.1016/S0735-1097(97)00352-5.View ArticlePubMedGoogle Scholar
- Ramahi TM, Longo MD, Cadariu AR, Rohlfs K, Slade M, Carolan S, Vallejo E, Wackers FJ: Dobutamine-induced augmentation of left ventricular ejection fraction predicts survival of heart failure patients with severe non-ischaemic cardiomyopathy. Eur Heart J. 2001, 22: 849-56. 10.1053/euhj.2001.2654.View ArticlePubMedGoogle Scholar
- Hummel JP, Lindner JR, Belcik JT, Ferguson JD, Mangrum JM, Bergin JD, Haines DE, Lake DE, DiMarco JP, Mounsey JP: Extent of myocardial viability predicts response to biventricular pacing in ischemic cardiomyopathy. Heart Rhythm. 2005, 2: 1211-1217. 10.1016/j.hrthm.2005.07.027.View ArticlePubMedGoogle Scholar
- Da Costa A, Thévenin J, Roche F, Faure E, Roméyer-Bouchard C, Messier M, Convert G, Barthélemy JC, Isaaz K: Prospective validation of stress echocardiography as an identifier of cardiac resynchronization therapy responders. eart Rhythm. 2006, 3: H406-413. 10.1016/j.hrthm.2005.12.017.View ArticleGoogle Scholar
- Ypenburg C, Schalij MJ, Bleeker GB, Steendijk P, Boersma E, Dibbets-Schneider P, Stokkel MP, Wall van der EE, Bax JJ: Extent of viability to predict response to cardiac resynchronization therapy in ischemic heart failure patients. J Nucl Med. 2006, 47: 1565-1570.PubMedGoogle Scholar
- Chalil S, Foley PW, Muyhaldeen SA, Patel KC, Yousef ZR, Smith RE, Frenneaux MP, Leyva F: Late gadolinium enhancement-cardiovascular magnetic resonance as a predictor of response to cardiac resynchronization therapy in patients with ischaemic cardiomyopathy. Europac. 2007, 9: 1031-1037. 10.1093/europace/eum133.View ArticleGoogle Scholar
- Ypenburg C, Schalij MJ, Bleeker GB, Steendijk P, Boersma E, Dibbets-Schneider P, Stokkel MP, Wall van der EE, Bax JJ: Impact of viability and scar tissue on response to cardiac resynchronization therapy in ischaemic heart failure patients. Eur Heart J. 2007, 28: 33-41. 10.1093/eurheartj/ehl379.View ArticlePubMedGoogle Scholar
- Kawara T, Derksen R, de Groot JR, Coronel R, Tasseron S, Linnenbank AC, Hauer RN, Kirkels H, Janse MJ, de Bakker JM: Activation delay after premature stimulation in chronically diseased human myocardium relates to the architecture of interstitial fibrosis. Circulation. 2001, 104: 3069-75. 10.1161/hc5001.100833.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.