Quantification of functional hand grip using electromyography and inertial sensor-derived accelerations: clinical implications
© Martin-Martin and Cuesta-Vargas; licensee BioMed Central Ltd. 2014
Received: 1 September 2014
Accepted: 4 December 2014
Published: 11 December 2014
Assessing hand injury is of great interest given the level of involvement of the hand with the environment. Knowing different assessment systems and their limitations generates new perspectives. The integration of digital systems (accelerometry and electromyography) as a tool to supplement functional assessment allows the clinician to know more about the motor component and its relation to movement. Therefore, the purpose of this study was the kinematic and electromyography analysis during functional hand movements.
Ten subjects carried out six functional movements (terminal pinch, termino-lateral pinch, tripod pinch, power grip, extension grip and ball grip). Muscle activity (hand and forearm) was measured in real time using electromyograms, acquired with the Mega ME 6000, whilst acceleration was measured using the AcceleGlove.
Electrical activity and acceleration variables were recorded simultaneously during the carrying out of the functional movements. The acceleration outcome variables were the modular vectors of each finger of the hand and the palm. In the electromyography, the main variables were normalized by the mean and by the maximum muscle activity of the thenar region, hypothenar, first interosseous dorsal, wrist flexors, carpal flexors and wrist extensors.
Knowing muscle behavior allows the clinician to take a more direct approach in the treatment. Based on the results, the tripod grip shows greater kinetic activity and the middle finger is the most relevant in this regard. Ball grip involves most muscle activity, with the thenar region playing a fundamental role in hand activity.
Relating muscle activation, movements, individual load and displacement offers the possibility to proceed with rehabilitation by individual component.
The hand is one of the fundamental elements of evolution for human beings and their interaction with the environment . One element of great importance in the development of activities of daily living (ADLs) is the maintenance of sensory and motor skills (required for object manipulation) . Changes in the functional capabilities of the hand could have a direct impact on the development of ADLs . Therefore, the correct assessment of these skills is very important for planning rehabilitation processes . One of the most complete tests of function is the Southampton Hand Assessment Procedure (SHAP) . This is based on carrying out six fundamental tasks, namely power, ball and extension grip, and terminal, termino-lateral and tripod pinch. They are all performed in three modalities (light objects, heavy objects and interaction). The main outcome variable is the action time in each of the tasks .
However, the evaluation of the functional capabilities of any musculoskeletal structure – and hence the hand – should be more complete. It requires the use of instruments that are reliable and sensitive, allowing the analysis of muscle activity during contraction . Surface electromyography allows the analysis of muscle activation through electrodes placed on the skin , a method used in previous pertinent studies [8–11], and the analysis of the function of the hand in healthy subjects based on the electrical activity of the muscles (electromyography) could provide reference values [8–11]. Furthermore, the use of inertial sensors (accelerometers) provides an overview of normal movement patterns  and inertial sensors are instruments that have shown high reliability and validity for the kinematic analysis of human movement  and own functional hand gestures [13, 14].
The analysis of the hand function based on electromyography and accelerometers with a gesture of reference has previously been undertaken [13, 14] to describe the maximum electromyography variables and the variation of acceleration in a task. Also, a combined electromyography and accelerometry system has been used by other authors to recognize movement . However, no study has been found that performs an analysis of the function of the hand in which the electromyography signals and inertial sensors are integrated in the course of the six functional tasks in the SHAP. Linking both elements (electromyography and accelerometers) with the various functional tests of the hand could be the point of origin for a full functional assessment.
Therefore the objective of this study is to conduct a descriptive analysis of the six functional tasks of SHAP, integrating electromyographic and kinematic signals collected during the execution of each gesture. A secondary objective of this study is to compare the values of muscle activation and functional kinematics during the tasks included in the SHAP. Our hypothesis is that the kinematic and muscle activation records will be different for each of the functional gestures that compose the SHAP.
Material and method
A descriptive study was carried out based on muscle activity and acceleration variables.
This study was conducted at the Faculty of Health Sciences, University of Málaga. Data were collected during the months of January to March, 2014.
The sample used to conduct this study comprised 10 healthy subjects (five men and five women). The inclusion criteria were: aged between 18 and 35 years, right-handed, no skin disorders, the absence of any type of motor impairment in the upper extremity and the cognitive ability to understand the orders given. Exclusion criteria were left-handed (the glove used is adapted to the right hand, see material section), skin and upper extremity motor disorders, fractures and/or surgery of the right arm in the 12 months prior to the study, cognitive impairment that would prevent understanding of the tasks and implanted electrical medical devices (neurostimulator, pacemaker).
The participants were randomly selected and they had to be able to read, understand and sign the informed consent provided by the researcher. The informed consent showed the purpose of the study and the procedure to be followed, together with the questionnaires to be completed.
The study has the ethical support of the Ethics and Research Committee of the institution, in line with the Helsinki declaration . In addition, the participant data have been treated in accordance with the Organic Law of Protection of Personal Data.
Sample descriptive tools
The Upper Limb Functionality Index (ULFI) shows reliability values of 0.96 in the tests-retests, with a confidence interval of 95%, and the internal consistency based on Cronbach’s alpha is 0.89 in the Spanish version . The psychometric properties of the QuickDASH questionnaire show a Cronbach’s alpha of 0.94 and a test-retest reliability of 0.95. The questionnaire consists of three parts: general (11 items), work (four items) and sports/music (four items) . In addition, the dynamometer used was of a palmar pressure hydraulic type with a Jamar “Sammons Preston Rolyan” . The palmar pressure force was obtained in kg/cm2 in two different positions: with the elbow in flexion and in extension. The Jamar hand dynamometer was adjusted based on the metacarpal measurements.
Accelerometry was recorded with the Acceleglove  device, together with Acceleglove Visualizer recording software. The lycra AcceleGlove is fitted with six inertial accelerometer-type sensors: one on the back of the second phalanx of each finger and a sixth on the back of the hand. The sampling frequency of the device was 120 Hz. Each accelerometer records three spatial axes (X, Y, Z) from the AcceleGlove within a range of ±1.5 g. The unit of measurement of the device is “g”, (1 g = 9.8 m/s2), the value of universal gravitation. According to the manufacturer, the relationship between the axes of the glove is as follows: if the hand is horizontal, the Z axis is the gravity vector, which is perpendicular to the surface of the earth; the X (flexion and extension movements) and Y (adduction and abduction movements) axes are perpendicular to each other and to Z (ratio of movements). The hardware provides continuous acceleration values for each of the elements of the hand (fingers and palm). The reference unit for time measure is the time stamp of Unix 1 January 1970 .
The time in seconds was calculated by subtracting from the first time record and dividing by the sample rate. The magnitudes of the accelerations were calculated based on the vector module of each of the elements, . The time unit was transformed into seconds.
The electromyography variables were registered using the MEGA ME6000 MT-M6T8-0-10 , measured in microvolts (μV), with the capture and data processing software Megawin 3.0.1, manufactured by Mega Electronics Ltd. (Kuopio, Finland). The sampling frequency of the device was 1000 Hz. Data capture was performed with ECG Lessa electrodes in two sizes : child-size for the hand and adult-size for the forearm region. The electrodes were positioned and arranged in accordance with Perotto et al.  and the European Society of Electromyography (SENIAM) . The skin was prepared to minimize the impedance experienced in electromyographic recording according to criteria established by SENIAM . In addition, for this purpose, we excluded subjects whose BMI was equal to or greater than 35 kg/m 2 from the study.
Description of the sample
Dynamo Max Ext (kg/cm2)
Dynamo Max Flex (kg/cm2)
QuickDash Sport (score)
Subsequently, the AcceleGlove (outcomes recording device) and the electromyography electrodes were positioned. The muscles for which electrical activity was recorded were: the hypothenar muscles (opponens digiti minimi), the thenar muscles (flexor pollicis brevis), the first dorsal interosseous muscle, the flexors of the wrist (palmaris longus), the flexor carpi ulnaris muscle and the extensors of the wrist (extensor carpi radialis). To avoid a crosstalk effect in the electromyographic recording, the guidelines established by SENIAM  and Merletti  were followed. In addition, an exhaustive study of palpatory anatomy was conducted with each participant to improve electromyographic recording.
Tests of maximal voluntary muscle contraction (MVC) were performed  on each selected muscle. This test was performed three times for each muscle. The maximum recorded value of each muscle was considered the maximum value of activation and used as a reference for the normalization of muscle activation during the functional task.
Participants placed the middle finger of the right hand on a reference mark located 20 cm from the interaction object. Participants received four acoustic alerts with five seconds of rest between them. The first alert was for the participants to remain at rest. The second alert was to move the hand to the interaction object and hold it. The third alert was to release the object and return to the starting point. The fourth alert was the end of the test. All participants performed each task three times. The second repetition of each functional movement was used for data analysis. The process was carried out with each of the functional movements. Electromyographic and kinematic data were recorded during movement and the rest time. Both devices were synchronized in real time and the data recorded simultaneously using DV Trigger (synchronization hardware), manufactured by Mega Electronics Ltd. (Kuopio, Finland). This element allows the recording process to be started simultaneously in two devices.
The outcome variables analysed in this study were of two different types: kinematic and electromyographic. The kinematic variables considered in this study were average acceleration of each finger and the palm of the hand (values for each of the functional gestures) and average acceleration of the gesture (values for each of the gestures made). Accelerometer values were measured in “g”, based on module vector of each sensor.
The electromyographic variables analysed in this study were mean muscle activation (the average muscle activation recorded for each muscle during the execution of the six functional gestures) and mean gesture activation (the average of all recorded muscle measurements analysed during the functional activation task). EMG values are presented in absolute values (uV) and normalized values (%). The normalized value represents the actual load or activation level required by a muscle with reference to the MVC value recorded. These variables provide an overview of each gesture.
The quantitative outcome variables based on mean acceleration of the resultant vector and the surface electromyography values relative to MVC and absolute terms (microvolts) were extracted. Data were extracted by a blinded researcher with more than 10 years of experience in kinematic and electromyographic data analysis.
A descriptive analysis (mean and standard deviation) was carried out independently on the selected anthropometric data, questionnaires and dynamometric values. A descriptive analysis was then performed with the outcome variables. The software used was SPSS 15.0  for Windows.
This study has 10 participants (five men and five women), with a mean age of 26.8 (±3.67) years old. The dynamometry tests performed did not reflect major changes in relation to the position of the arm; accordingly, no significant differences were observed in action between the two positions (elbow flexed at 90° and extended elbow). Furthermore, a description of the functional capacity of the upper limbs was undertaken; the questionnaires used were ULFI and QuickDash, with mean values of 0.85 (±1.80) and 4.31 (±10.30) respectively. The remaining descriptive variables for the sample can be observed in Table 1.
Mean values of acceleration by vector and task (mean ± sd), in [g]
Mean values of muscle activity (mean ± sd) in [μV]
Normalized mean values of maximum EMG activation (mean ± sd) in [%MVC]
Addressing the primary objective of this study, the joint use of electrical activity and acceleration recording systems is a feasible way to analyse hand function, based on the SHAP tasks. This allows the comparison of the electromyographic and kinematic values during task performance. Both elements could be used to identify which functional tasks have greater acceleration based on the values of the hand segments when carrying out the tasks. Recording electrical activity by means of electromyography provides knowledge of muscle action in each movement.
On the basis of the results obtained, we demonstrate that muscle activation and acceleration of the segments of the hand change according to the task being performed. Therefore, in line with the hypothesis in this study, kinematic and muscle activation records will be different for each of the functional gestures that compose the SHAP.
The study reflects different patterns of movement and activation for each of the gestures. These patterns of muscle activation and movement (acceleration) can be used for functional recovery of the hand. The analysis of data derived from accelerometry and electromyography will allow the clinician to take a more direct approach in treatment, with understanding of which muscles should be recovered first, based on the level of activation and the elements with greater mobility.
The mean acceleration values (Table 2) indicate that the tripod movement requires the greatest hand movement. This may be due to the placement which must be performed by the first three fingers of the hand in order to make the grip and the compensatory flexion carried out by the others. The middle finger is the hand element with the greatest acceleration. The highest mean values produced in the vector of the middle finger (1.28 g) should be treated with caution as being indicative of functionality (2 degrees of freedom). Individually, the thumb has a greater range of movement (3 degrees of freedom), despite registering lower acceleration values (1.15). The tripod grip pinch (1.28 g) has the highest average acceleration value, followed by the extension grip (1.27 g).
From a functional perspective, rehabilitation treatment should not be performed directly on the movement with greater mobility (in this study, the tripod grip) as this may hinder functional achievement. Subjects who are incapable of reaching the first set of goals may decide to drop out. In consequence, work on movements must be progressive and, based on the results of this study, start with the terminal pinch and ball grip, movements which produce less variation in the acceleration.
There is a direct relationship between EMG absolute values (Table 3) and the normalized results (Table 4), as greater electrical activation implies greater individual load. However, there is a second relationship: force generation (and hence EMG values) is not the same for each muscle. In consequence, not all muscles can achieve the same absolute values; if they do, specific muscles may be subject to direct overload, as shown in Table 3. The higher muscle load values, absolute activity values (Table 3) and normalized values (Table 4) coincide in showing the importance of these elements. As mentioned above, it is not always the case that the muscles with greater electrical activity (Table 4) present the highest normalized values (Table 3). A single absolute activation value may present different normalized values depending on the muscle in question.
The muscle load record (Table 4) in healthy subjects shows how muscles behave with normal motor control and without alterations to activation. The muscles with the highest degree of activation in the series of six functional movements are the short flexor of the thumb (thenar region), followed by the flexor carpi ulnaris (Table 4); the ball grip requires the greatest muscle activity, whereas the lateral grip requires the least. This information can be used to guide the functional recovery process from least to greatest muscle activity.
Kinematic and electromyography uses
Electromyography and accelerometry variables were analysed in a previous study [13, 14] that aimed to generate normal movement patterns by examining the changes in the different phases of a task. In addition, as in this study, electromyography (forearm) and accelerometry (forearm and wrist) have been used jointly by other authors for movement recognition. The results obtained could be used to estimate acceleration movements by means of electromyography values . This shows the overall behaviour of both methods and their effectiveness for the analysis of movement, which can be used for prosthetic control. Gazzoni et al.  provide an analysis of the change in the activation and localization of muscle activity through electromyography according to the hand function and its angle. Furthermore, Ngeo et al.  use a simultaneous accelerometry and electromyography system for detecting muscle activation patterns. These patterns may be predictors of acceleration segments made by hand. The use of electromyography for the recognition of movement patterns has also been employed by Birdwell et al. . In this case, the electromyographic signals provide the reference for virtual hand movements.
The forearm muscles selected for this study (flexor carpi ulnaris and palmaris longus) have been analysed by other authors to examine power grip at different intensities . The results (Table 3 of this study) show that the flexor muscle with the greatest activity in mean values in power grip is the ulnar flexor carpi (28.83 μV ± 11.03), just as determined in the study carried out by Oskouei et al. .
Oskouei et al.  agree on the importance of the extensor muscles, although these are not assessed in their study; this aspect is analysed by Hoozemans et al.  as a predictor of strength (calibration procedure using fluctuating bursts of grip; the most effective predictions were between 27N and 47N), establishing the existence of a high relation between the extensor muscles, power grip and ability to counteract the wrist flexion. This factor is related to the importance of the wrist extensors, extensor carpi radialis (36.62 μV ± 14.43), as the forearm muscle with the most activity (Table 3).
A variable to take into account in rehabilitation treatment is muscle fatigue in repetitive tasks (three-tip pinch) analysed in twelve hand muscles . In this regard, the intrinsic muscles of the thumb (abductor pollicis brevis and flexor pollicis brevis) show greater activity. In terms of activation, this is in line with the results of Table 4, which shows the thenar region (flexor pollicis brevis) as the most active in the hand in the three-tip grip (6.09%MVC ± 1.86).
Moreover, muscle stress changes according to the position adopted by the forearm to perform the different tasks, with the position of least fatigue for the forearm being neutral . The neutral position is used in this study as a starting point in rest position and at the end of the task, with an approximate duration of two seconds in each period. However, the muscular load may be reduced through use of support surfaces. Onyebeke et al.  discuss the use of these elements of support while performing tasks of long duration with a mouse on a computer (similar to the pleota tennis grip), finding that these reduce muscle stress. The way in which the electrodes are arranged on the skin and their initial position may affect the electromyography results . In this study, as in the results of Takala et al. , the electrodes were applied in the neutral position, after identifying the muscles by means of maximum voluntary contractions and in line with the indicated bibliography.
The angular opening of the wrist is one of the elements which most affects activation in the extensors digitorum, extensor carpi ulnaris and flexor carpi radialis . These values are lower at high speed and the results may be affected; grip angle should be considered as a variable in studies related to the tasks.
In consequence, based on the results obtained, the functional rehabilitation process should begin with the least demanding movements in terms of mobility (terminal pinch) and electrical load (lateral pinch). Individually, rehabilitation may start with the muscle with the greatest functional load, such as the short flexor of the thumb and the flexor carpi ulnaris. Based on mobility (accelerations), rehabilitation would start with the middle finger and ring finger. However, starting rehabilitation with the thumb should be considered, given its high functionality and integration in all hand movements.
- Young RW: Evolution of the human hand: the role of throwing and clubbing. J Anat 2003, 202: 165–174. 10.1046/j.1469-7580.2003.00144.xView ArticleGoogle Scholar
- Hill EL, Khanem F: The development of hand preference in children: the effect of task demands and links with manual dexterity. Brain Cogn 2009, 71: 99–107. 10.1016/j.bandc.2009.04.006View ArticleGoogle Scholar
- Rosberg H-E, Carlsson KS, Cederlund RI, Ramel E, Dahlin LB: Costs and outcome for serious hand and arm injuries during the first year after trauma - a prospective study. BMC Public Health 2013, 13: 501. 10.1186/1471-2458-13-501View ArticleGoogle Scholar
- Shelton FD, Volpe BT, Reding M: Motor impairment as a predictor of functional recovery and guide to rehabilitation treatment after stroke. Neurorehabil Neural Repair 2001, 15: 229–237. 10.1177/154596830101500311View ArticleGoogle Scholar
- Light CM, Chappell PH, Kyberd PJ: Establishing a standardized clinical assessment tool of pathologic and prosthetic hand function: normative data, reliability, and validity. Arch Phys Med Rehabil 2002, 83: 776–783. 10.1053/apmr.2002.32737View ArticleGoogle Scholar
- Drost G, Stegeman DF, van Engelen BGM, Zwarts MJ: Clinical applications of high-density surface EMG: a systematic review. J Electromyogr Kinesiol 2006, 16: 586–602. 10.1016/j.jelekin.2006.09.005View ArticleGoogle Scholar
- Merletti R, Aventaggiato M, Botter A, Holobar A, Marateb H, Vieira TMM: Advances in surface EMG: recent progress in detection and processing techniques. Crit Rev Biomed Eng 2010, 38: 305–345. 10.1615/CritRevBiomedEng.v38.i4.10View ArticleGoogle Scholar
- Muceli S, Jiang N, Farina D: Multichannel surface EMG based estimation of bilateral hand kinematics during movements at multiple degrees of freedom. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Conf 2010, 2010: 6066–6069.Google Scholar
- Hashemi Oskouei A, Paulin MG, Carman AB: Intra-session and inter-day reliability of forearm surface EMG during varying hand grip forces. J Electromyogr Kinesiol 2012,23(1):216–222.View ArticleGoogle Scholar
- Hoozemans MJM, van Dieën JH: Prediction of handgrip forces using surface EMG of forearm muscles. J Electromyogr Kinesiol Off J Int Soc Electrophysiol Kinesiol 2005, 15: 358–366. 10.1016/j.jelekin.2004.09.001View ArticleGoogle Scholar
- Takala E-P, Toivonen R: Placement of forearm surface EMG electrodes in the assessment of hand loading in manual tasks. Ergonomics 2013, 56: 1159–1166. 10.1080/00140139.2013.799235View ArticleGoogle Scholar
- Cuesta-Vargas A, Galán-Mercant A, Williams J: The use of inertial sensors system for human motion analysis. Phys Ther Rev 2010, 15: 462–473. 10.1179/1743288X11Y.0000000006View ArticleGoogle Scholar
- Martin-Martin J, Cuesta-Vargas AI: A biomechanical study of spherical grip. SpringerPlus 2013, 2: 591. 10.1186/2193-1801-2-591View ArticleGoogle Scholar
- Martin-Martin J, Cuesta-Vargas AI: A kinematic and electromyographic study of grip force. Eur J Physiother 2014, 16: 183–191.Google Scholar
- Li Y, Chen X, Zhang X, Wang K, Wang Z: A sign-component-based framework for chinese sign language recognition using accelerometer and sEMG data. IEEE Trans Biomed Eng 2012,59(10):2695–704.View ArticleGoogle Scholar
- WMA declaration of Helsinki - ethical principles for medical research involving human subjects. http://www.wma.net/es/30publications/10policies/b3/
- Cuesta-Vargas AI, Gabel PC: Cross-cultural adaptation, reliability and validity of the Spanish version of the upper limb functional index. Health Qual Life Outcomes 2013, 11: 126. 10.1186/1477-7525-11-126View ArticleGoogle Scholar
- Hervás MT, Navarro Collado MJ, Peiró S, Rodrigo Pérez JL, López Matéu P, Martínez Tello I: Spanish version of the DASH questionnaire. Cross-cultural adaptation, reliability, validity and responsiveness. Med Clínica 2006, 127: 441–447.View ArticleGoogle Scholar
- Patterson medical - strength dynamometer jamar hand. https://www.pattersonmedical.co.uk/app.aspx?cmd=get_item&id=1174
- AcceleGlove. http://www.anthrotronix.com/
- LLC B: Time measurement systems: Iso 8601, metric time, unix time, 12-hour clock, 24-hour clock, decimal time, thai six-hour clock. Memphis (Tennessee), USA: General Books LLC; 2010.Google Scholar
- Mega Electronics Ltd pioneers in biosignal monitoring technology. http://www.megaemg.com/
- Registration papers lessa - products. http://www.lessap.com/productos.php
- Perotto A, Delagi EF: Anatomical guide for the electromyographer: the limbs and trunk. Springfield: Charles C Thomas Publisher; 2005.Google Scholar
- Hermens HJ, Freriks B, Merletti R: SENIAM: European Recommendations for Surface Electromyography. Volume 8. Enschede, the Netherlands: Roessingh Research and Development; 1999.Google Scholar
- Merletti R: The electrode-skin interface and optimal detection of bioelectric signals. Physiol Meas 2010., 31: 3 p preceding S157Google Scholar
- PT JMM: Daniels and Worthingham’s muscle testing: techniques of manual examination, 8e (Daniel’s & Worthington’s Muscle Testing. 8th edition. California: Saunders; 2007.Google Scholar
- IBM SPSS statistics. http://www-01.ibm.com/software/analytics/spss/products/statistics/
- Gazzoni M, Celadon N, Mastrapasqua D, Paleari M, Margaria V, Ariano P: Quantifying forearm muscle activity during wrist and finger movements by means of multi-channel electromyography. PLoS One 2014, 9: e109943. 10.1371/journal.pone.0109943View ArticleGoogle Scholar
- Ngeo JG, Tamei T, Shibata T: Continuous and simultaneous estimation of finger kinematics using inputs from an EMG-to-muscle activation model. J Neuroeng Rehabil 2014, 11: 122. 10.1186/1743-0003-11-122View ArticleGoogle Scholar
- Birdwell JA, Hargrove LJ, Weir RF, Kuiken TA: Extrinsic finger and thumb muscles command a virtual hand to allow individual finger and grasp control. IEEE Trans Biomed Eng 2014. Epub ahead of printGoogle Scholar
- Danna-Dos Santos A, Poston B, Jesunathadas M, Bobich LR, Hamm TM, Santello M: Influence of fatigue on hand muscle coordination and EMG-EMG coherence during three-digit grasping. J Neurophysiol 2010, 104: 3576–3587. 10.1152/jn.00583.2010View ArticleGoogle Scholar
- Qin J, Chen H, Dennerlein JT: Wrist posture affects hand and forearm muscle stress during tapping. Appl Ergon 2013, 44: 969–976. 10.1016/j.apergo.2013.03.013View ArticleGoogle Scholar
- Onyebeke LC, Young JG, Trudeau MB, Dennerlein JT: Effects of forearm and palm supports on the upper extremity during computer mouse use. Appl Ergon 2014, 45: 564–570. 10.1016/j.apergo.2013.07.016View ArticleGoogle Scholar
- Dai W, Zhao X, Wang Z, Yang L: Electromyographical study on muscle fatigue in repetitive forearm tasks. J Huazhong Univ Sci Technol Med Sci 2007, 27: 358–361. 10.1007/s11596-007-0403-8View ArticleGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.