Validation of the UNESP-Botucatu unidimensional composite pain scale for assessing postoperative pain in cattle
© de Oliveira et al.; licensee BioMed Central Ltd. 2014
Received: 4 February 2014
Accepted: 15 August 2014
Published: 6 September 2014
The recognition and measurement of pain in cattle are important in determining the necessity for and efficacy of analgesic intervention. The aim of this study was to record behaviour and determine the validity and reliability of an instrument to assess acute pain in 40 cattle subjected to orchiectomy after sedation with xylazine and local anaesthesia. The animals were filmed before and after orchiectomy to record behaviour. The pain scale was based on previous studies, on a pilot study and on analysis of the camera footage. Three blinded observers and a local observer assessed the edited films obtained during the preoperative and postoperative periods, before and after rescue analgesia and 24 hours after surgery. Re-evaluation was performed one month after the first analysis. Criterion validity (agreement) and item-total correlation using Spearman's coefficient were employed to refine the scale. Based on factor analysis, a unidimensional scale was adopted.
The internal consistency of the data was excellent after refinement (Cronbach;s α coefficient = 0.866). There was a high correlation (p < 0.001) between the proposed scale and the visual analogue, simple descriptive and numerical rating scales. The construct validity and responsiveness were confirmed by the increase and decrease in pain scores after surgery and rescue analgesia, respectively (p < 0.001). Inter- and intra-observer reliability ranged from moderate to very good. The optimal cut-off point for rescue analgesia was > 4, and analysis of the area under the curve (AUC = 0.963) showed excellent discriminatory ability.
The UNESP-Botucatu unidimensional pain scale for assessing acute postoperative pain in cattle is a valid, reliable and responsive instrument with excellent internal consistency and discriminatory ability. The cut-off point for rescue analgesia provides an additional tool for guiding analgesic therapy.
KeywordsCattle Castration Orchiectomy Reliability Responsiveness Validity
The assessment of pain in animals is challenging due to their lack of verbal expression . This challenge is intensified in cattle because, as they are prey in their natural state, they may avoid expressing pain to limit vulnerability . Cattle are routinely subjected to surgical procedures related to management and production, such as dehorning and orchiectomy, usually without adequate analgesia -. In surveys of veterinarians concerning the use of analgesics in cattle practice, lack of knowledge in recognising pain , the belief that farm animals feel less pain than smaller animals , economic reasons , and the lack of valid and reliable instruments to assess pain have been cited as the main reasons why analgesics are not used more frequently .
The assessment of pain in cattle is important in determining the need for analgesic intervention, in evaluating the effectiveness of treatment and in comparing the effects of various analgesics. Difficulty in the assessment of pain is not as serious a problem in other domestic species, as there are validated scales in the literature for the assessment of acute pain in dogs -, cats - and horses -. In cattle, a scoring system for gait has been validated and found to be reliable and sensitive for identifying cows with severe hoof lesions . However, there are no validated scales for the assessment of acute postoperative pain in cattle.
To develop an accurate tool to assess pain, it is necessary that the scale show validity, responsiveness and reliability ,-. The primary tool used to assess pain in animals is species-specific behaviour . Although pain signals in ruminants may not be easy to recognise, changes associated with pain can often be seen in the animals' appearance, posture, gait, appetite, weight, interaction with other animals and with the environment and in the frequency of movement and vocalisation. In addition, animals in pain may protect, lick or bite the wound area ,,-. When experiencing pain after orchiectomy, cattle may stomp the ground with their feet, ease their quarters, directing attention to the lesion site, stand still with their limbs extended backward or apart or lie down with their hind limbs extended ,. Such behaviour is absent or reduced when the animal receives adequate analgesia .
The scales most often used for measuring postoperative pain are ordinal in nature and can be classified as either unidimensional or multidimensional . Unidimensional scales, such as the simple descriptive (SDS), numerical rating (NRS) and visual analogue (VAS) scales ,, only measure pain intensity , whereas multidimensional or composite assessment scales also take into account the sensory and affective qualities of pain ,.
In developing an instrument to assess pain, it is important to determine the minimum score that indicates a need for analgesic intervention. For this purpose, ROC (Receiver Operating Characteristic) curve analysis may be used. This methodology was first used in veterinary medicine in a pain assessment scale for cats ,.
The primary objective of this study was to validate a scale to assess acute pain in cattle subjected to orchiectomy. To that end, the following specific objectives were defined: 1) develop a record of pain-related behaviour; 2) correlate the number of steps taken, time spent lying down and number of lying bouts with the degree of postoperative pain; 3) refine the proposed scale; 4) evaluate inter- and intra-observer reliability, validity and responsiveness; and 5) define a cut-off point for analgesic intervention.
The study was approved by the Ethics Committee for Animal Use of FMVZ-UNESP-Botucatu under protocol number 147/2011. Forty Nellore cattle two to three years of age weighing 365 ± 51 kg were used. The animals were considered to be healthy based on clinical and laboratory assessment (complete blood count and biochemical tests).
Prior to the assessment of normal behaviour, the animals were maintained in groups of 3 to 4 for 24 hours in a pasture paddock, with hay and feed placed in troughs and with water ad libitum, to permit them to adapt to the new environment. At this point, pedometers (Pedometer Plus®)a were placed on the right forelimbs of the animals to provide a means of measuring the number of steps, time spent lying down and number of lying bouts for 24 hours prior to (baseline - D0) and 24 hours after (D1) surgery . The behaviour of the animals was then filmed intermittently over a 24-hour period to determine the normal behaviour of each animal. Following the 24-hour period of baseline behavioural assessment, an experienced surgeon performed an orchiectomy on each animal using the open technique. Each animal received 0.025 mg/kg xylazine (Anasedan®)b intramuscularly (IM) followed after 10 minutes by the injection of 10 mL of 1% lidocaine without vasoconstrictor (Xylestesin®)c in each spermatic cord and 5 mL in the incision line. Ten minutes thereafter, the orchiectomy was performed. Immediately after xylazine administration, benzathine penicillin (Pentabiótico®)d was administered IM at a dose of 30,000 IU/kg.
After the procedure, the animals were kept in the same enclosure in which they had been housed prior to surgery. Rescue analgesia was applied to all animals after the four-hour assessment at the end of the surgical procedure, with 3 mg/kg ketoprofen (Ketojet®)e and 0.2 mg/kg morphine (Dimorf®)c, both administered intravenously (IV) in separate syringes. For the rescue analgesia, the animals were placed in a restraining chute, and the analgesic was administered in the marginal ear vein using a butterfly cannula (19G) after restraint of the head. The analgesic doses were selected based on the literature and on information provided by the manufacturer ,.
One observer evaluated and filmed the animals' behaviour for 25 minutes in each of the following periods: 24 and 16 hours before and 1, 2, 4, 5, 8 and 24 hours after orchiectomy. The animals were filmed using two cameras positioned near the observation paddock. To minimise potential effects of the observer on the animals' behaviour, the observer and cameras were positioned behind a black plastic screen. There were two openings in the screen through which filming and observation could be performed.
The pain scoring scale was developed based on previous studies ,,, on the results of the pilot study and on the analysis of footage taken during the experiment. During the experiment, categories or levels were incorporated, modified or excluded from the original scale. Behaviour was assessed before surgery (M1), at the anticipated time of greatest pain, between 1 and 4 hours after orchiectomy (M2) (according to the pilot study, this was the period during which the animals expressed the most intense pain-related behaviour), one hour after administration of the rescue analgesic (M3) to investigate the efficacy of analgesia and 24 hours after the surgical procedure (M4), totalling 66 hours of video assessment. Based on analysis of the videos obtained at each time period, the percentage duration of specific behaviours was recorded, including time spent eating, ruminating, drinking, walking, standing and in recumbency. Also noted were changes in locomotion and in standing and recumbency posture, interactions with the environment and with other animals and position of the head. This analysis was used to identify behaviours related to pain.
To assess content validity, the scale was sent to three evaluators with expertise in cattle behaviour. These evaluators analysed and scored the behaviour by degree of importance according to the following scale: -1 = irrelevant item; 0 = do not know; 1 = relevant item. The total item correlation was evaluated, and items that achieved a score of ≥ 0.5 were accepted .
The proposed instrument yielded a variable score scale composed of behavioural categories. The variables were ordinal in nature and exhibited three descriptive levels, to which a score (a numerical value) was assigned. In the scoring, zero reflected normality and one or two represented changes related to pain, with a maximum score of 16 points.
The same three evaluators who performed the content validity analysis analysed the edited footage to validate the scale. For this purpose, they received a hard drive containing four films, each approximately three minutes long, for each of the 40 cattle evaluated, corresponding to the time points previously described. The chronological order of the videos was randomised so that the evaluators were blinded with respect to the assessed time points, and the descriptions of the animals' behaviour did not include the pain scores. The researcher responsible for the study, who was considered the local observer, also analysed the videos so that the agreement between the blinded observers and the local observer could be compared to determine inter-observer reliability. One month after the first evaluation, the blinded observers evaluated the videos again with the order of the cattle and videos changed to establish intra-observer reliability.
After watching each film, the observers specified, based on their clinical experience, whether rescue analgesia should be performed and provided sequentially determined pain scores using the Visual Analogue (VAS), Numerical Rating (NRS), Simple Descriptive (SDS) scales and the proposed scale. The data regarding the application of rescue analgesia were used to determine the minimum score related to the need for rescue analgesia.
The Shapiro-Wilk normality test was used to compare the percentage of time during which the animals engaged in specific behaviours (states) during each of the time periods observed in the footage. Because all variables were nonparametric, Friedman's test was used. The paired t-test was used to compare the number of steps taken, the time spent lying down (minutes) and the number of lying bouts during the pre (D0) and postoperative (D1) periods. Differences were considered significant when p < 0.05.
The criterion validity was evaluated based on the agreement between the scores determined by the blinded observers and those determined by the local observer. First, the percentage of absolute agreement for each scale item was determined considering only M2. The percentage of absolute agreement was considered satisfactory when it was ≥ 60%. The weighted kappa coefficient was then calculated with a 95% confidence interval (CI)  for each scale item, considering all of the assessment times as a group (MA = M1, M2, M3 and M4). The kappa coefficient results were interpreted according to Altman's classification : 0.81 to 1.0, very good; 0.61 to 0.8, good; 0.41 to 0.6, moderate; 0.21 to 0.4, fair; and < 0.2, poor.
Factor analysis was used to define the number of factors (dimensions or domains) determined by different variables to establish the dimensionality of the scale . Exploratory factor analysis was performed based on principal component analysis, and factors were identified based on the Kaiser criterion, which recommends retaining all components with eigenvalues > 1 . The factor structure was determined by attributing each item with a factor loading and communality > 0.5 to a factor.
The item-total correlation was evaluated using Spearman's non-parametric correlation coefficient between each item and the sum of all scale items. This correlation coefficient was used to evaluate the relevance of each item to the instrument and to identify items that contributed strongly to the total scale score. Items with a correlation coefficient < 0.4 were rejected .
The internal consistency of the scale after refinement was evaluated by calculating Cronbach's alpha coefficient . Internal consistency with a value of > 0.70 was considered adequate . The concurrent (criterion) validity was evaluated by comparing the scores obtained using the scale with the scores determined using the VAS, NRS and SDS. Spearman's correlation coefficient was calculated for each blinded observer and for the local observer as well as for the blinded observers as a group.
To determine the inter-observer reliabilities with or without the local observer and the intra-observer reliability for each scale item, the intraclass correlation coefficient (ICC) with 95% CI was used for MA and for M2 and M4 grouped together . The two-factor model with absolute agreement criterion was selected, and the values obtained were interpreted using Altman's classification, described earlier.
The construct validity was established based on hypothesis-testing methodology. The first hypothesis was that if the scale actually measures pain, the scores after surgery should be higher than the preoperative scores (M1 versus M2). The second and third hypotheses were that the scores should decrease after the administration of analgesics and over time (M2 versus M3 and M2 versus M4, respectively). The scores were expressed as medians, and the Wilcoxon test was used for the analysis of significance (p < 0.05) . This analysis evaluated the responsiveness of the scale.
To determine a minimum score related to the need for intervention or rescue analgesia, ROC (Receiver Operating Characteristic) curve analysis was conducted to provide a graphical representation of the relationship between "true positives" (sensitivity) and "false positives" (1-specificity). The area under the curve (AUC) was also determined; this value indicates the discriminatory ability of the test , with AUC values above 0.9 representing high accuracy .
Mean values (±SD) of percentage duration of the behaviour of 40 cattle during the perioperative period
55.3 ± 33.4A
6.2 ± 14.4B
39.0 ± 32.3A
34.1 ± 36.8A
2.8 ± 19.9
5.1 ± 14.3
0.4 ± 97.5
0.07 ± 0.42
0.25 ± 0.8
0.1 ± 0.3
8.8 ± 6.7A
0.3 ± 0.7B
10.0 ± 17.1A
6.7 ± 8.2A
4.8 ± 4.8A
0.4 ± 1.8B
0.5 ± 1.9B
31.7 ± 28.8B
25.5 ± 26.7B
43.7 ± 32.2A
50.5 ± 33.5A
16.6 ± 20.7A
1.7 ± 11.1B
34.4 ± 24.0A
4.7 ± 20.8B
1.8 ± 7.6B
Abnormal lying down ‡
7.5 ± 15.3
Lying down with head resting on/close to the ground
4.6 ± 6.9A
1.0 ± 1.9
0.07 ± 0.2
0.25 ± 1.2
1.2 ± 4.3
Head below the line of the spinal column
2.5 ± 6.9
0.05 ± 0.3
0.04 ± 0.3
Regarding the frequency of the occurrence of specific behaviours, arching the back and extending the neck cranially were observed only at M2 in 4 and 13, respectively, of the 40 animals studied. Kicking, wagging the tail abruptly, looking at and licking the wound were observed more frequently at M2 (21, 7, 14 and 7 of 40 animals, respectively) than at M3 (6, 2, 3 and 0 of 40 animals, respectively) and M4 (7, 0, 1, and 6 of 40 animals, respectively) and were not observed at M1.
No difference between D0 and D1 was observed in the number of steps (D0 - 5401 ± 1142; D1 - 4702 ± 1737; p = 0.08) or in time spent lying down (D0 - 486 ± 265 minutes; D1 - 509 ± 213; p = 0.64); however, the animals lay down more often in the period 24 hours after surgery (D0 - 17 ± 4; D1 - 20 ± 4.9; p = 0.0007).
All of the scale items proposed yielded scores higher than 0.5 regarding the item-total correlation and were thus accepted.
Refinement of the proposed scale
Percentage of absolute agreement and agreement by the weighted kappa reliability coefficient - criterion validity
The absolute agreement between the scores assigned by the blinded observers and those assigned by the local observer at M2 was considered unsatisfactory for the items standing posture and head position (<60%) based on the analysis of all evaluators. Only one evaluator did not find satisfactory agreement for the items locomotion, interactive behaviour and miscellaneous behaviours, and these items were therefore retained in the scale.
Agreement between the local evaluator and blinded observers for each item on the scale
0.60 (0.43 - 0.75)
0.26 (0.14 - 0.38)
0.40 (0.26 - 0.54)
0.51 (0.37 - 0.66)
0.11 (0.03 - 0.19)
0.23 (0.12 - 0.34)
0.56 (0.41 - 0.71)
0.47 (0.34 - 0.61)
0.61 (0.41 - 0.75)
0.65 (0.52 - 0.77)
0.67 (0.58 - 0.75)
0.77 (0.69 - 0.86)
0.67 (0.55 - 0.79)
0.56 (0.44 - 0.68)
0.74 (0.63 - 0.85)
0.86 (0.80 - 0.93)
0.76 (0.67 - 0.86)
0.81 (0.72 - 0.89)
Attention to surgical wound
0.68 (0.55 - 0.82)
0.67 (0.51 - 0.81)
0.85 (0.75 - 0.95)
0.80 (0.72 - 0.81)
0.62 (0.50 - 0.74)
0.84 (0.78 - 0.92)
Factor analysis (construct validity)
Exploratory factor analysis of the acute postoperative pain scale in cattle
Attention to surgical wound
% Cumulative variance
Correlation coefficient of item score with total score
Spearman's correlation coefficient between item score and total score
p = 0.000
p = 0.000
p = 0.000
p = 0.000
Attention to surgical wound
p = 0.000
p = 0.000
Behaviour included in the miscellaneous behaviour item after refining the acute postoperative pain scale in cattle
Percentage of absolute agreement between the local and the blinded observers
Hind limbs extended caudally
Head below the line of spinal column
Licks surgical wound
Validation of the UNESP-Botucatu unidimensional pain scale for assessing postoperative pain in cattle
Evaluation of the internal consistency of the scale
UNESP-Botucatu unidimensional pain scale for acute postoperative pain assessment in cattle
▪ (0) Walking with no obviously abnormal gait.
▪ (1) Walking with restriction, may be with hunched back and/or short steps.
▪ (2) Reluctant to stand up, standing up with difficulty or not walking.
▪ (0) Active; attention to tactile and/or visual and/or audible environmental stimuli; when near other animals, can interact with and/or accompany the group.
▪ (1) Apathetic: may remain close to other animals, but interacts little when stimulated.
▪ (2) Apathetic: may be isolated or may not accompany the other animals; does not react to tactile, visual and/or audible environmental stimuli.
▪ (0) Moves normally.
▪ (1) Restless, moves more than normal or lies down and stands up with frequency.
▪ (2) Moves less frequently in the pasture or only when stimulated.
▪ (0) Normorexia and/or rumination.
▪ (1) Hyporexia.
▪ (2) Anorexia.
▪ Wagging the tail abruptly and repeatedly.
▪ Licking the surgical wound.
▪ Moves and arches the back when in standing posture.
▪ Kicking/foot stamping.
▪ Hind limbs extended caudally when in standing posture.
▪ Head below the line of spinal column.
▪ Lying down in ventral recumbency with full or partial extension of one or both hind limbs.
▪ Lying down with the head on/close to the ground.
▪ Extends the neck and body forward when lying in ventral recumbency.
(0) All of the above described behaviours are absent.
(1) Presence of 1 of the behaviours described above.
(2) Presence of 2 or more of the behaviours described above.
Additional file 1: Video 1.: Normal behaviours and pain behaviours in cattle. (MP4 17 MB)
The Cronbach's α coefficient of the scale after refinement was 0.866, indicating that the instrument has excellent internal consistency and lending weight to the feasibility of using the total score to interpret the results obtained.
Concurrent validity (criterion validation)
Correlation between the UNESP-Botucatu unidimensional pain scale and the VAS, NRS and SDS
Agreement between blinded observers, including or not including the local observer, for each scale item
Local and blinded observers
M2 and M4
M2 and M4
0,50 (0,34 - 0,62)
0,44 (0,23 - 0,60)
0,52 (0,42 - 0,62)
0,45 (0,31 - 0,59)
0,68 (0,56 - 0,75)
0,63 (0,50 - 0,74)
0,69 (0,62 - 0,75)
0,63 (0,52 - 0,73)
0,57 (0,45 - 0,67)
0,54 (0,38 - 0,67)
0,61 (0,52 - 0,69)
0,56 (0,44 - 0,67)
0,80 (0,72 - 0,85)
0,76 (0,91 - 0,84)
0,80 (0,75 - 0,85)
0,77 (0,68 - 0,84)
0,66 (0,58 - 0,73)
0,72 (0,63 - 0,80)
0,71 (0,65 - 0,76)
0,75 (0,68 - 0,82)
Intra-observer reliability for each scale item and observer
0.67 (0.57 - 0.75)
0.62 (0.52 - 0.71)
0.69 (0.60 - 0.77)
0.81 (0.75 - 0.86)
0.74 (0.66 - 0.80)
0.77 (0.69 - 0.82)
0.80 (0.74 - 0.85)
0.63 (0.52 - 0.71)
0.69 (0.60 - 0.77)
0.96 (0.95 - 0.97)
0.82 (0.77 - 0.87)
0.81 (0.75 - 0.85)
0.82 (0.77 - 0.87)
0.61 (0.50 - 0.70)
0.78 (0.71 - 0.83)
Total scores determined by local and blinded observers
1 (0-5) †
0 (0-7) †
3.5 (1-8) †
0 (0-6) †
Determination of the cut-off point - ROC curve
Determination of optimum cut-off based on ROC curve analysis
Sensitivity (95% CI)
Specificity (95% CI)
100.00 (98.3 - 100.0)
0.00 (0.0 - 0.5)
100.00 (98.3 - 100.0)
51.14 (47.5 - 54.8)
99.54 (97.5 - 100.0)
65.55 (62.0 - 69.0)
98.16 (95.3 - 99.5)
75.24 (72.0 - 78.3)
96.77 (93.5 - 98.7)
81.29 (78.3 - 84.0)
95.85 (92.3 - 98.1)
87.35 (84.7 - 89.7)
89.86 (85.1 - 93.5)
91.66 (89.4 -- 93.5)
76.04 (69.8 - 81.6)
95.56 (93.8 - 96.9)
61.29 (54.5 - 67.8)
97.04 (95.6 - 98.1)
35.48 (29.1 - 42.2)
99.06 (98.1 - 99.6)
17.51 (12.7 - 23.2)
99.87 (99.3 - 100.0)
0.00 (0.0 - 1.7)
100.00 (99.5 - 100.0)
Because animals are unable to report their pain as humans can , the recognition of pain in animals requires the ability to understand the behaviour of the target species, the behavioural changes typically observed in animals in pain and the specific changes that occur in each animal's behaviour in response to pain. In this context, the video records obtained in this study served as an initial survey of items that might be appropriate for use in building and subsequently validating a scale for the assessment of pain in cattle.
The use of video recording for the validation of scales and for behavioural assessment is a common tool -, that permits the simultaneous analysis of an animal by multiple evaluators to be performed as often as necessary. In this study, following evaluation of the videos by the researcher, the films were edited according to the behaviour observed at different points in time. After reviewing the films, changes not covered in the initial scale were identified and items deemed irrelevant were excluded. The behaviour of animals in pain (M2) showed a reduction in eating and moving around, and when animals in pain did move around, they did so with restrictions and/or short steps and/or hunched backs. In addition, animals in pain spent more time lying down with their heads on or near the ground. When in the standing position, these animals assumed an abnormal posture, e.g., hunched and rigid and/or with the hind limbs extended caudally. Arched-back movements were also observed more frequently in animals in pain, along with cranial extensions of the neck while lying down, kicking, wagging the tail abruptly and looking at and licking the surgical wound. Given their relationship to pain, these behaviours were incorporated into the scale.
Some of the behaviours observed in this study have been described previously in cattle subjected to orchiectomy. These behaviours include remaining idle for longer periods, assuming an abnormal standing posture , and exhibiting gait changes involving shorter, more cautious steps ,. With respect to time spent lying down, the results reported in the literature vary according to whether xylazine and analgesics were used. It would seem that animals spend more time in the standing position when xylazine is not used ,,, and that they spend more time lying down and less time moving around when it is used ,. In the present study, although xylazine might contribute to lying down after the surgery, the fact that administration of additional analgesia resulted in less rather than more lying behaviour suggests that this behaviour is due to pain rather than to the sedative effect of xylazine.
The use of xylazine in the present study was necessary because the Nellore breed is skittish and difficult to control. No animal showed recumbency after orchiectomy when leaving the restraining chute, once again demonstrating that the sedation was mild. Low doses of xylazine (0.015 to 0.025 mg/kg IV or IM) generally promote sedation without recumbency in ruminants . In conclusion, decreased activity in cattle may be a good indicator of pain.
The reduced time spent eating observed in our study is consistent with other findings in the literature that report reductions in grazing time ,, eating frequency ,, and, in the case of calves, suckling time . The benefit of rescue analgesia with respect to this behaviour was also evident from the fact that the time spent feeding increased after rescue analgesia was performed.
Kicking and abrupt wagging of the tail were observed more frequently at M2, as described previously ,,. These events may occur after the local anaesthetic effect has lost its effectiveness  but may also be related to the presence of flies, which might represent one limitation of the study. Although it is impossible to completely eliminate flies from the environment, care was taken to reduce the number of flies present by using fly repellent, and the study was conducted in the winter when there is a low incidence of insects. One possible indication that flies had little effect on the kicking and abrupt tail wagging observed in this study was that kicking the abdomen was not observed at baseline; furthermore, the behaviour of wagging the tail abruptly and repeatedly is very characteristic and differs from the motion the animal makes to ward off flies. Thus, it would appear that these two behaviours are also related to pain in cattle .
Although it might be expected that the number of steps taken by the animals would be reduced after orchiectomy, no significant difference was observed in the number of steps recorded before and after orchiectomy. This finding differs from results previously reported in the literature . The difference may be explained by the fact that the animals in the previous study did not receive analgesia , in contrast to this study, in which there was only a short (4-hour) span during which analgesia was not provided. The short period of pain experienced by the animals in the present study was most likely insufficient to influence the data obtained during the 24-hour evaluation.
The larger number of lying bouts observed after orchiectomy may be related to the restlessness and discomfort of the animals in the period prior to the application of the rescue analgesic. A similar phenomenon was observed in previous studies in cattle ,. Lying-down behaviour is also evaluated on the pain scales commonly used for dogs ,, and cats ,, showing that although it is important to develop species-specific tools to assess pain, some pain behaviours are common among species.
Because the pedometer is not expensive and is relatively easy to handle, it can be a useful tool in the assessment of pain in cattle, especially when data analysis is carried out over relatively short periods.
Methods other than pain scales have also been used to investigate pain following castration in cattle. These methods include the assessment of physiological and neuroendocrine changes, such as serum cortisol concentration, and infrared thermography . Facial expression of pain and kinematic and force platform gait analysis have been used in mice and in horses and dogs, respectively -. However, these methods have not yet been validated in cattle, and they either do not provide information in real time or require special equipment that is not currently available and/or is impractical under field conditions.
Validity and reliability are the key attributes of a scale that can be used to identify and quantify pain in animals. Reliability demonstrates the ability of the scale to reproduce the results regardless of the evaluator and at different times by the same evaluator . In this study, the assessment of content validity was performed using the judgment of experts in the field who analysed the representativeness of each item in relation to the scale as a whole . This methodology, which is well accepted ,,, refers to the scope and adequacy with which the instrument reflects the phenomenon of interest, in this case, pain .
Criterion validity tests the effectiveness of a scale's measurement by comparing results obtained using that scale to results obtained using a previously validated method . Criterion validity can be predictive when evaluating the criterion after testing and concurrent when evaluating the instrument and the criterion simultaneously . In tests of criterion validity, the correlation between the scale and another instrument, ideally the gold standard ,, is evaluated.
Considering that, to our knowledge, no gold standard instrument has been developed to evaluate pain in cattle and that correlation of the total scores obtained using our proposed scale with the scores determined by VAS may be questionable, an alternative method was used to investigate criterion validity in this study. The method involved comparing the agreement between pain scores assigned by blinded evaluators and a "gold standard" evaluator, in this case, the local evaluator. This method has been used with instruments designed for use in cats  and in young children .
Although the VAS, SDS and NRS may not show inter-observer reliability when tested on animals, they are nonetheless widely used to validate veterinary pain scales ,,, because the gold standards of verbal expression and self-assessment evaluation are not available in animals. Although inter-observer reliability may not be adequate when using VAS , intra-observer agreement or reliability is consistent over time  and may be a good option for measuring and comparing pain assessed by the same trained observer over time, as was done in this study .
The same methodology used for criterion validity in studies in cats , was used to refine the scale by comparing the pain scores determined by blinded observers with those determined by the local evaluator. Subsequently, criterion validity should be evaluated by correlating the results obtained using the proposed scale and another instrument considered the gold standard (concurrent validity) . Given the absence in the literature of validated scales for pain assessment in cattle, the pain scores on the scale proposed in this study were compared with the scores obtained using three other classical scales used in animals, the VAS, the NRS and the SDS. There was a high correlation between the results obtained using the four scales. Although these scales have not been validated in animals, this approach has been widely used to evaluate pain scales in veterinary medicine ,,,.
Using factor analysis, it is possible to determine the dimensionality of the scale , i.e., the number of factors (dimensions or domains) represented by different variables . Because the scale in question generated only one factor, it was considered unidimensional, in contrast to the scales validated for cats, which were considered multidimensional based on this analysis ,. Factor analysis is commonly used to develop an instrument and to relate a large number of variables such that the items that define specific parts of the construct are grouped together .
Despite the low reliability observed for the items standing posture and head position and the low correlation of the item attention to the surgical wound with the total scale score in this study, it was deemed important to retain a behaviour for each item.
The behaviours rigid hind limbs, hunched back and head below the line of spinal column may not have been clearly visible in the videos, and this may have resulted in the observed poor correlation between the blinded observers and the local observer on the items standing posture and head position. Conversely, the behaviour head below the line of spinal column obtained satisfactory agreement when considering only M2. Regarding the item attention to the surgical wound, the description of the behaviour looking at the surgical wound may not have been wholly appropriate because it produced different results when assessed by the local evaluator and the blinded observers. A description such as moves the snout in the direction of the surgical wound might have clarified observation of this behaviour; the description looking at the surgical wound was subjective because it could denote looking at the abdomen and/or to the side for another reason, a fact that may have confused the observers.
According to the Cronbach's α value, the scale employed in this work has excellent internal consistency ,. Internal consistency ensures that the scores of the items comprising the scale can be summed to produce a total score related to the overall assessment of pain intensity .
The moderate to good inter-observer agreement found in this work demonstrates the consistency of the results obtained by different evaluators and the ability of the instrument to produce consistent results . The lower level of agreement for the item locomotion may be due to the short video analysis time and the animals' way of walking, which may have hindered the definition of the category. Thus, the results of both inter- and intra-observer reliability tests demonstrated good repeatability and stability of the scale.
The analysis grouping of M2 and M4 was important to confirm the reliability of the scale because these points represent the two most challenging times for pain assessment. A similar approach was used to validate a postoperative acute pain scale in cats ,, but in that case, only M2 was considered separately. In our study, M4 was included because it also represents a challenging time, given the reduction in analgesic effect and the manifestation of pain-related behaviour that typically occurs after 24 hours.
Construct validity examines whether a given instrument detects predictable changes in the construct . It can be evaluated by the well-known group method. This method determines whether the instrument detects differences between groups and is based on testing the hypothesis that time and intervention, both surgical and analgesic, should alter the pain scores . The observed differences between pain scores at the time of greatest pain (M2) and the scores at other time points confirm the construct validity used in this work by verifying the reduction in pain scores in response to analgesia and over time . This method has also been used to validate scales in veterinary medicine , and attests to the responsiveness of the scale using a similar approach.
ROC curve analysis was used to determine the minimum score required for analgesic intervention , as was previously performed for a pain assessment scale in cats ,. The determination of scores that suggest a need for the use of analgesics assists the professional's clinical decision, affirms the effectiveness of analgesic treatment  and helps avoid unnecessary suffering in animals. Based on the balanced sensitivity and specificity criteria observed in this study, an optimum cut-off of > 4 was identified, i.e., additional analgesia is recommended when the pain score is ≥ 5 (0-10 point scale). It should be emphasised that according to clinical evaluation, additional analgesia must be performed if deemed necessary even if the score is lower than the cut-off point.
The high AUC observed (0.963) in this study indicates that the scale has excellent discriminatory ability and high accuracy, i.e., the instrument can correctly classify subjects with or without pain ,. Similar results were observed in the validation of a pain scale in cats ,.
A possible limitation of this study is the absence of a control or uncastrated group of animals. The inclusion of a control group was considered when the study was designed, and a pilot study was performed to address this point. Subsequently, the authors decided to use only one castrated group and a larger number of animals based on the rationale that the animals' behaviour during the time period immediately prior to surgery could be considered a control because, at this point, the animals had already adapted to the environment and no management changes were performed during the study that could influence the results. This methodology has previously been used in cats ,, dogs , and horses . The results observed here, which show significant changes in pain scores before surgery, after surgery and after analgesia, support the validity of the construct as well as the responsiveness of the scale. Additional support for the idea that the pre-surgical time period provided an appropriate control comes from the fact that the observers were blinded to the test moments and the order of the videos was randomised to avoid any bias. Although, in a very few cases, it was possible to observe the region of the testicles in the videos, it was not possible to determine from the video footage whether the animal had already been castrated. Another consideration is that it would be difficult to compare a different, uncastrated control group of animals with a group of castrated animals because the response to pain varies according to each individual.
The results of this study allow us to state that the UNESP-Botucatu unidimensional pain scale for assessing acute postoperative pain is valid and reliable. However, clinical tests with different analgesics and surgical protocols are recommended to assess the scale's clinical applicability.
It is concluded that, following the refinement of the originally proposed scale, the UNESP-Botucatu unidimensional pain scale for assessing acute postoperative pain in cattle is a valid, reliable and responsive instrument with excellent internal consistency and discriminatory ability. The cut-off point for rescue analgesia provides an additional tool for guiding analgesic therapy.
aSAE Afimilk, Israel.
bVetbrands, Paulínia, São Paulo, Brazil.
cCristália, Itapira, São Paulo, Brazil.
dFort Dodge, Campinas, São Paulo, SP.
eAgener, Embu Guaçu, São Paulo, SP.
The financial support of the São Paulo Research Foundation (FAPESP) is acknowledged (Process 2010/08967-0 and 2011/02472-1).
- Anil SS, Anil L, Deen J: Challenges of pain assessment in domestic animals. J Am Vet Med Assoc. 2002, 220 (3): 313-319. 10.2460/javma.2002.220.313.View ArticlePubMedGoogle Scholar
- Watts S, Clarke K: A survey of bovine practitioners' attitudes to pain and analgesia in cattle. Cattle Pract. 2000, 10: 361-362.Google Scholar
- Whay H, Huxley J: Pain relief in cattle: a practitioner's perspective. Cattle Pract. 2005, 13: 81-85.Google Scholar
- Hewson CJ, Dohoo IR, Lemke KA, Barkema HW: Factors affecting Canadian veterinarians' use of analgesics when dehorning beef and dairy calves. Can Vet J. 2007, 48 (11): 1129-1136.PubMed CentralPubMedGoogle Scholar
- Fajt VR, Wagner SA, Norby B: Analgesic drug administration and attitudes about analgesia in cattle among bovine practitioners in the United States. J Am Vet Med Assoc. 2011, 238 (6): 755-767. 10.2460/javma.238.6.755.View ArticlePubMedGoogle Scholar
- Lorena SERS, Luna SPL, Lascelles BDX, Corrente JE: Attitude of Brazilian veterinarians in the recognition and treatment of pain in horses and cattle. Vet Anaest Analg. 2013, 40 (4): 410-418. 10.1111/vaa.12025.View ArticleGoogle Scholar
- Raekallio M, Heinonen K, Kuussaari J, Vainio O: Pain alleviation in animals: attitudes and practices of Finnish veterinarians. Vet J. 2003, 165: 131-135. 10.1016/S1090-0233(02)00186-7.View ArticlePubMedGoogle Scholar
- Anil L, Anil SS, Deen J: Pain detection and amelioration in animals on the farm: issues and options. J Appl Anim Welfare Sci. 2005, 8 (4): 261-278. 10.1207/s15327604jaws0804_3.View ArticleGoogle Scholar
- Hewson CJ, Dohoo IR, Lemke KA, Barkema HW: Canadian veterinarians' use of analgesics in cattle, pigs, and horses in 2004 and 2005. Can Vet J. 2007, 48 (2): 155-164.PubMed CentralPubMedGoogle Scholar
- Firth AM, Haldane SL: Development of a scale to evaluate postoperative pain in dogs. J Am Vet Med Assoc. 1999, 214 (5): 651-659.PubMedGoogle Scholar
- Holton L, Reid J, Scott EM, Pawson P, Nolan A: Development of a behaviour-based scale to measure acute pain in dogs. Vet Rec. 2001, 148 (17): 525-531. 10.1136/vr.148.17.525.View ArticlePubMedGoogle Scholar
- Morton CM, Reid J, Scott EM, Holton LL, Nolan AM: Application of a scaling model to establish and validate an interval-level pain scale for assessment of acute pain in dogs. Am J Vet Res. 2005, 66 (12): 2154-2166. 10.2460/ajvr.2005.66.2154.View ArticlePubMedGoogle Scholar
- Murrell JC, Psatha EP, Scott EM, Reid J, Hellebrekers LJ: Application of a modified form of the Glasgow pain scale in a veterinary teaching centre in the Netherlands. Vet Rec. 2008, 162 (13): 405-410. 10.1136/vr.162.13.403.View ArticleGoogle Scholar
- Brondani JT, Luna SPL, Minto BW, Santos BPR, Beier SL, Matsubara LM, Padovani CR: Validity and responsiveness of a multidimensional composite scale to assess postoperative pain in cats. Arq Bras Med Vet Zoo. 2012, 64 (6): 1529-1538. 10.1590/S0102-09352012000600019.View ArticleGoogle Scholar
- Brondani JT, Luna SPL, Minto BW, Santos BPR, Beier SL, Matsubara LM, Padovani CR: Reliability and cut-off point related to the analgesic intervention of a multidimensional composite scale to assess postoperative pain in cats. Arq Bras Med Vet Zoo. 2013, 65 (1): 153-162. 10.1590/S0102-09352013000100024.View ArticleGoogle Scholar
- Brondani JT, Mama KR, Luna SPL, Wright BD, Niyom S, Ambrosio J, Vogel PR, Padovani CR: Validation of the English version of the UNESP-Botucatu multidimensional composite pain scale for assessing postoperative pain in cats. BMC Vet Res. 2013, 9: 1-15. 10.1186/1746-6148-9-143.View ArticleGoogle Scholar
- Bussiéres G, Jacques C, Lainay O, Beauchamp G, Leblond A, Cadore JL, Desmaizieres LM, Cuvelliez SG, Troncy E: Development of a composite orthopaedic pain scale in horses. Res Vet Sci. 2008, 85 (2): 294-306. 10.1016/j.rvsc.2007.10.011.View ArticlePubMedGoogle Scholar
- Graubner C, Gerber V, Doherr M, Spadavecchia C: Clinical application and reliability of a post-abdominal surgery pain assessment scale (PASPAS) in horses. Vet J. 2011, 188 (2): 178-183. 10.1016/j.tvjl.2010.04.029.View ArticlePubMedGoogle Scholar
- Sutton GA, Dahan R, Turner D, Paltiel O: A behaviour-based pain scale for horses with acute colic: scale construction. Vet J. 2013, 196: 394-401. 10.1016/j.tvjl.2012.10.008.View ArticlePubMedGoogle Scholar
- Sutton GA, Paltiel O, Soffer M, Turner D: Validation of two behaviour-based pain scales for horses with acute colic. Vet J. 2013, 197 (3): 646-650. 10.1016/j.tvjl.2013.04.007.View ArticlePubMedGoogle Scholar
- Flower F, Weary D: Effect of hoof pathologies on subjective assessments of dairy cow gait. J Dairy Sci. 2006, 89: 139-146. 10.3168/jds.S0022-0302(06)72077-X.View ArticlePubMedGoogle Scholar
- Crellin D, Sullivan TP, Babl E, O'Sullivan R, Hutchinson A: Analysis of the validation of existing behavioral pain and distress scales for use in the procedural setting. Pediat Anest. 2007, 17: 720-733. 10.1111/j.1460-9592.2007.02218.x.View ArticleGoogle Scholar
- Figueiredo KMOB, Lima KC, Guerra RO: Instrumentos de avaliação do equilíbrio corporal em idosos. [Instruments for the assessment of physical balance in the elderly]. Rev Bras Cineantropom Desempenho Hum. 2007, 9 (4): 408-413.Google Scholar
- Alexandre NMC, Coluci MZO: Validade de conteúdo nos processos de construção e adaptação de instrumentos de medidas. [Content validity in the development and adaptation processes of measurement instruments]. Ciência & Saúde Coletiva. 2011, 16 (7): 3061-3068. 10.1590/S1413-81232011000800006.View ArticleGoogle Scholar
- Flecknell P: Analgesia from a veterinary perspective. Brit J Anaest. 2008, 101 (1): 121-124. 10.1093/bja/aen087.View ArticleGoogle Scholar
- Robertson IS, Kent JE, Molony V: Effect of different methods of castration on behavior and plasma cortisol in calves of three ages. Res Vet Sci. 1994, 56 (1): 8-17. 10.1016/0034-5288(94)90189-9.View ArticlePubMedGoogle Scholar
- Molony V, Kent JE, Robertson IS: Assessment of acute and chronic pain after different methods of castration of calves. Appl Anim Behav Sci. 1995, 46: 33-48. 10.1016/0168-1591(95)00635-4.View ArticleGoogle Scholar
- Molony V, Kent JE: Assessment of acute pain in farm animals using behavioral and physiological measurements. J Anim Sci. 1997, 75 (1): 266-272.PubMedGoogle Scholar
- Fischer A, Knight T, Cosgrove G, Death A, Anderson C, Duganzich D, Mathews L: Effects of surgical or banding castration on stress responses and behaviour of bulls. Aust Vet J. 2001, 79 (4): 279-284. 10.1111/j.1751-0813.2001.tb11981.x.View ArticleGoogle Scholar
- Ting STL, Earley B, Hughes JML, Crowe MA: Effect of ketoprofen, lidocaine local anesthesia, and combined xylazine and lidocaine caudal epidural anesthesia during castration of beef cattle on stress responses, immunity, growth, and behavior. J Anim Sci. 2003, 81 (5): 1281-1293.PubMedGoogle Scholar
- Ting STL, Earley B, Crowe MA: Effect of repeated ketoprofen administration during surgical castration of bulls on cortisol, immunological function, feed intake, growth, and behavior. J Anim Sci. 2003, 81 (5): 1253-1264.PubMedGoogle Scholar
- Currah JM, Hendrick SH, Stookey JM: The behavioral assessment and alleviation of pain associated with castration in beef calves treated with flunixin meglumine and caudal lidocaine epidural anesthesia with epinephrine. Can Vet J. 2009, 50 (4): 375-382.PubMed CentralPubMedGoogle Scholar
- González LA, Schwartzkopf-Genswein KS, Caulkett NA, Janzen E, McAllister TA, Fierheller E, Schaefer AL, Haley DB, Stookey JM, Hendrick S: Pain mitigation after band castration of beef calves and its effects on performance, behavior, Escherichia coli, and salivary cortisol. J Anim Sci. 2010, 88 (2): 802-810. 10.2527/jas.2008-1752.View ArticlePubMedGoogle Scholar
- Millman ST: Behavioral responses of cattle to pain and implications for diagnosis, management, and animal welfare. Vet Clin N Am-Food A. 2013, 29 (1): 47-58. 10.1016/j.cvfa.2012.11.007.View ArticleGoogle Scholar
- Pereira LV, Sousa FAEF: Estimação em categorias dos descritores da dor pós-operatória. [Estimation in categories of postoperative pain descriptors]. Rev latino-am enfermagem. 1998, 6 (4): 41-48.Google Scholar
- Von Baeyer CL, Spagrud LJ: Systematic review of observational (behavioral) measures of pain for children and adolescents aged 3 to 18 years. Pain. 2007, 127 (1-2): 140-150. 10.1016/j.pain.2006.08.014.View ArticlePubMedGoogle Scholar
- Holton LL, Scott EM, Nolan AM, Reid J, Welsh E, Flaherty D: Comparison of three methods used for assessment of pain in dogs. J Am Vet Med Assoc. 1998, 212 (1): 61-66.PubMedGoogle Scholar
- Riebold TW: Ruminants. Veterinary Anesthesia and Analgesia. Edited by: Tranquilli WJ, Thurmon JC, Grimm KA. Blackwell Publishing, Oxford; 2007:731-746.Google Scholar
- Suraseranivongse S, Santawat U, Kraiprasit K, Petcharatana S, Prakkamodom S, Muntraporn N: Cross-validation of a composite pain scale for preschool children within 24 hours of surgery. Brit J Anaest. 2001, 87 (3): 400-405. 10.1093/bja/87.3.400.View ArticleGoogle Scholar
- Cohen J: Weighted kappa: nominal scale agreement with provision for scaled disagreement or partial credit. Psychol Bull. 1968, 70: 213-220. 10.1037/h0026256.View ArticlePubMedGoogle Scholar
- Altman D: Some common problems in medical research. Practical Statistics for Medical Research. Chapman & Hall, London; 1991:404-408.Google Scholar
- Brondani JT, Loureiro Luna SP, Padovani CR: Refinement and initial validation of a multidimensional composite scale for use in assessing acute postoperative pain in cats. Am J Vet Res. 2011, 72 (2): 174-183. 10.2460/ajvr.72.2.174.View ArticlePubMedGoogle Scholar
- Kaiser H: The varimax criterion for analytic rotation in factor analysis. Psychometrika. 1958, 23: 187-200. 10.1007/BF02289233.View ArticleGoogle Scholar
- Cronbach L: Coefficient alpha and the internal structure of tests. Psychometrika. 1951, 16: 297-333. 10.1007/BF02310555.View ArticleGoogle Scholar
- Jensen M: Questionnaire validation: a brief guide for readers of the research literature. Clin J Pain. 2003, 19: 345-352. 10.1097/00002508-200311000-00002.View ArticlePubMedGoogle Scholar
- Bartko J: Measurement and reliability: statistical thinking considerations. Schizophrenia Bull. 1991, 17 (3): 483-489. 10.1093/schbul/17.3.483.View ArticleGoogle Scholar
- Deyo R, Diehr P, Patrick D: Reproducibility and responsiveness of health status measures. Control Clin Trials. 1991, 12: 142S-158S. 10.1016/S0197-2456(05)80019-4.View ArticlePubMedGoogle Scholar
- Streiner D, Cairney J: What's under the ROC? An introduction to receiver operating characteristics curves. Can J Psychiatry. 2007, 52 (2): 121-128.PubMedGoogle Scholar
- Webster H, Morin D, Jarrell V, Shipley C, Brown L, Green A, Wallace R, Constable P: Effects of local anesthesia and flunixin meglumine on the acute cortisol response, behavior, and performance of young dairy calves undergoing surgical castration. J Dairy Sci. 2013, 96 (10): 6285-6300. 10.3168/jds.2012-6238.View ArticlePubMedGoogle Scholar
- White BJ, Coetzee JF, Renter DG, Babcock AH, Thomson DU, Andresen D: Evaluation of two-dimensional accelerometers to monitor behavior of beef calves after castration. Am J Vet Res. 2008, 69 (8): 1005-1012. 10.2460/ajvr.69.8.1005.View ArticlePubMedGoogle Scholar
- Pauly C, White BJ, Coetzee JF, Robert B, Baldridge S, Renter DG: Evaluation of analgesic protocol effect on calf behavior after concurrent castration and dehorning. Int J Appl Res Vet Med. 2012, 10 (1): 54-61.Google Scholar
- Sutherland M, Ballou M, Davis B, Brooks T: Effect of castration and dehorning singularly or combined on the behavior and physiology of Holstein calves. J Anim Sci. 2013, 91: 935-942. 10.2527/jas.2012-5190.View ArticlePubMedGoogle Scholar
- Coetzee J: Assessment and management of pain associated with castration in cattle. Vet Clin N Am-Food A. 2013, 29: 75-101. 10.1016/j.cvfa.2012.11.002.View ArticleGoogle Scholar
- Langford DJ, Bailey AL, Chanda ML, Clarke SE, Drummond TE, Echols S, Glick S, Ingrao J, Klassen-Ross T, LaCroix-Fralish ML, Matsumiva L, Sorge RE, Sotocinal SG, Tabaka JM, Wong D, van den Maagdenberg AMJM, Ferrari MD, Craig KD, Mogil JS: Coding facial expressions of pain in the laboratory mouse. Nat Methods. 2010, 7: 447-449. 10.1038/nmeth.1455.View ArticlePubMedGoogle Scholar
- Ishihara I, Bertone AL, Rajala-Schultz PJ: Association between subjective lameness grade and kinetic gait parameters in horses with experimentally induced forelimb lameness. Am J Vet Res. 2005, 66 (10): 1805-1815. 10.2460/ajvr.2005.66.1805.View ArticlePubMedGoogle Scholar
- Evans R, Horstman C, Conzemius M: Accuracy and optimization of force platform gait analysis in Labradors with cranial cruciate disease evaluated at a walking gait. Vet Surg. 2005, 34 (5): 445-449. 10.1111/j.1532-950X.2005.00067.x.View ArticlePubMedGoogle Scholar
- Pasquali L: Psicometria: [Psychometrics]. Rev Esc Enferm USP. 2009, 43: 992-999. 10.1590/S0080-62342009000500002.View ArticleGoogle Scholar
- Gauvain-Piquard A, Rodary C, Rezvani A, Serbouti S: The development of the DEGRR: a scale to assess pain in young children with cancer. Eur J Pain. 1999, 3: 165-176. 10.1053/eujp.1999.0118.View ArticlePubMedGoogle Scholar
- Reid J, Nolan A: A comparison of the postoperative analgesic and sedative efects of flunixin and papaveretum in the dog. J Small Anim Pract. 1991, 32: 603-608. 10.1111/j.1748-5827.1991.tb00900.x.View ArticleGoogle Scholar
- DeVon HA, Block ME, Moyle-Wright P, Ernst DM, Hayden SJ, Lazzara DJ, Savoy SM, Kostas-Polston E: A psychometric toolbox for testing validity and reliability. J Nurs Scholarship. 2007, 39 (2): 155-164. 10.1111/j.1547-5069.2007.00161.x.View 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.