- Research article
- Open Access
Gait analysis in clinically healthy sheep from three different age groups using a pressure-sensitive walkway
BMC Veterinary Researchvolume 8, Article number: 87 (2012)
Understanding normal gait requires allowing for variations in normal patterns by the sex, age, and species in question. Therefore, the aim of this study was to evaluate kinetic and temporospatial parameters in clinically healthy sheep from three different age groups with a pressure-sensing walkway. The sheep were judged to be healthy based on the results of complete physical and orthopaedic examinations and had no history of lameness. Twenty-one clinically healthy female Santa Ines sheep were divided into three groups: G1 – seven animals, aged from 8 to 12 months and weighing 19.5-33 kg; G2 - seven individuals, aged from 2 to 4 years and weighing 26.5-42 kg; and G3 - seven sheep, aged more than 5 years and weighing 37.3-45 kg. The animals were examined from two directions: first on the left side and then on the right side of the handler. The data from the first five valid trials in each direction were collected for each sheep and analysed using the designated software. A trial was considered valid if the sheep walked within the correct velocity (1.1-1.3 m/s) and acceleration (from −0.15 to 0.15 m/s2) ranges. The peak vertical force (PVF), vertical impulse (VI), gait cycle time, stance time, swing time, stride length, and the percentage body weight distribution among the four limbs were determined.
No significant differences were observed, in either the forelimbs or the hind limbs, between the left and right sides or between the two directions for any of the variables. No significant temporospatial differences were found among the groups. Significant PVF (%BW) differences were observed in the forelimbs (G1 > G3) and hind limbs (G1 > G3), and significant VI differences were observed in the forelimbs (G1 > G3).
Young healthy sheep differ from older sheep in the vertical forces they create when walking at the same velocity on a pressure-sensing walkway.
Gait analysis is usually used to directly aid patient treatment and to better understand locomotion [1–3]. Several methods of gait analysis are available, including visual analysis, which depends on the ability and perspicacity of the investigator, and specific analysis, which requires specialised equipment [2–4].
Force platforms and pressure distribution sensors, such as pressure-sensing walkways, may be used to measure forces, moments, and accelerations [2, 3]. Force platforms measure the orthogonal ground-reaction forces (vertical, mediolateral and craniocaudal) that result from limb positioning during locomotion, whereas pressure-sensing walkways measure vertical ground-reaction forces only [1–7]. Vertical forces have been commonly studied in animals due to their low variability and their magnitude relative to other orthogonal forces [1, 4, 6, 8, 9]. Craniocaudal and mediolateral forces may be affected by external factors and have smaller magnitudes [1, 9].
Different pressure-measuring devices have been used to perform gait analyses in animals [7, 10–14]. Some validation studies using dogs and horses have shown that certain kinetic measurements may differ when using a force platform rather than a pressure-sensing walkway [7, 10, 11], and highly accurate absolute force values may not be obtainable using a pressure plate . However, several pressure-sensing walkway systems have been proven to be useful in clinical settings and for repeated gait assessments [11–13].
Additionally, some authors have suggested that experiments using pressure-sensing walkways are less time consuming, as walkways capture multiple sequential steps and as the contralateral limbs may be evaluated in the same trial [7, 10, 11, 13, 15].
Understanding normal gait requires understanding that normal patterns may vary by sex, age , and the species in question. Sheep have been used as an in vivo experimental model for orthopaedic research studies due to their size, availability and ease of handling, and they can sometimes provide an alternative to dogs [16–18]; some aspects of gait analysis of the sheep require more comprehension, however.
For example, a previous study using a pressure-sensing walkway reported that female Suffolk-mix sheep (with body masses ranging from 69.3 to 103 kg) are not suitable for gait analysis assessments, largely due to their flight zone and flocking behaviour . In another study, however, female Merino-mix sheep (mean body mass 63 kg, range 46–90 kg) were trained to walk over a pressure-sensitive platform. The platform was used for nine weekly evaluations of a healing tibial defect .
Therefore, the aim of this study was to use a pressure-sensing walkway to evaluate kinetic and temporospatial parameters in clinically healthy female sheep from three different age groups. The hypothesis was that the older sheep would have altered kinetic parameters due to aging.
This study followed the guidelines for the care and use of laboratory animals and was approved by the Institutional Ethics Committee.
Twenty-one clinically healthy intact female Santa Ines sheep were divided into three groups: G1 – seven animals, aged from 8 to 12 months and weighing 19.5-33 kg (mean ± SD, 25.78 ± 4.75 kg); G2 - seven individuals, aged from 2 to 4 years and weighing 26.5-42 kg (mean ± SD, 31.52 ± 4.88 kg); and G3 – seven animals, aged more than 5 years and weighing 37.3-45 kg (mean ± SD, 42.21 ± 3.11 kg).
The sheep were judged to be healthy based on complete physical and orthopaedic examinations and had no history of lameness. Before the data collection, the sheep were trained in being led by a halter twice each day for a period of approximately three weeks. Subsequently, the sheep were trained to walk across the walkway twice each day for a period of one week. In addition, hoof trimming was performed one week before the exams. Food was used to motivate the sheep to walk in all of the training sessions. Each sheep was weighed on a single electronic scale immediately before the data collection. All of the sheep were handled by a single handler.
The kinetic and temporospatial gait parameters were measured using a 1.951 mm x 447 mm pressure-sensitive walkway (Walkway High Resolution HRV4; Tekscan, South Boston, Massachusetts, USA) containing 33,408 pressure-sensing cells. The sensors of the pressure-sensing walkway were equilibrated and calibrated using a phantom according to the manufacturer’s specifications. During each trial, the sheep were led in a straight line over the pressure-sensing walkway by a handler using a halter. The animals were handled from two different directions: first on the left side of the handler and then on the right side. Food was used to motivate the animals.
Forty trials (20 in each direction) were captured for each animal. The data from the first five valid trials in each direction were collected for each sheep and analysed using the designated software (Walkway 7.0; Tekscan). A trial was considered valid if the sheep walked within the correct velocity (1.1-1.3 m/s) and acceleration (from −0.15 to 0.15 m/s2) ranges without head movement or pulling on the halter and if all four limbs had contact with the surface of the walkway during each walk cycle.
The temporospatial gait cycle time (s), stance time (s) swing time (s) and stride length (m) parameters were determined for each limb. The stance time percentage was calculated as follows: (stance time/gait cycle time) x 100. The swing time percentage was calculated as (swing time/gait cycle time) x 100. The stride corresponded to the distance between two consecutive ground contacts of the same limb.
The peak vertical force (PVF) and vertical impulse (VI) kinetic parameters were also determined. Both were normalised to the sheep’s body weight and represented as percentages of body weight (%BW and %BW*s). The percentage body weight distribution among the four limbs was calculated by (PVF of the limb/total PVF of the four limbs) x 100.
Limb lengths and relative velocity
The limb lengths of the sheep were measure using retroreflective spherical markers (1.8 cm in diameter) that were tagged using a 3-camera kinematic system (Vicon MX-3+; Vicon, Oxford Metrics Group, Oxford, UK). The markers were placed on the skin by a single investigator using quick-drying glue. The forelimb markers were placed on the point of the cranial angle of the scapula, acromium of the scapulohumeral joint, lateral epicondyle of the humerus, styloid ulnar process, and distal lateral aspect of the III and IV metacarpus. On the hind limbs, the markers were placed on the greater trochanter of the femur, femorotibial joint between the lateral epicondyle of the femur and the fibular head, lateral malleolus and distal lateral aspect of the III and IV metatarsus. The length of each limb was determined by the sum of the distances between each pair of markers on that limb.
The relative velocities for the forelimbs and hind limbs were determined by the “Froude number”, which was defined as v2/gh, where v is the velocity, g is the gravitational acceleration and h is the height [12, 37].
The data were analysed using a linear repeated measurements model. For the temporospatial parameters and kinetic data, the animal side and the direction on the pressure-sensing walkway were considered to be the intra-sheep factors, while the group served as the inter-sheep factor. For the limb lengths, the animal’s side was considered to be the intra-sheep factor, and the group was the inter-sheep factor. The inter-sheep factor for body mass was the group. The sequential Bonferroni adjustment procedure was applied to our contrasts. The values were expressed as the mean ± standard deviation, and the coefficients of variation (CV) were calculated. A one-way ANOVA followed by Tukey's test was performed to compare the relative velocities of the forelimbs and hind limbs. Differences were considered significant at p < 0.05.
Pearson's correlation coefficients (r) were used to evaluate the linear relationships between the limb lengths and the PVF. The correlations were deemed significant at the 5% probability level.
After the training procedure, the sheep were able to walk appropriately on the pressure-sensing walkway.
No significant differences were found between the kinetic and temporospatial parameters of the left and right forelimbs or the left and right hind limbs when the group and the direction were not considered. None of the kinetic or temporospatial parameters, in either the forelimbs and or the hind limbs, differed significantly between the left and right sides or between directions when the group and the animal side were not considered.
The temporospatial values did not differ significantly among the groups in either the forelimbs or hind limbs (Tables 1 and 2). The forelimbs had significant differences in the PVF (G1 > G3) and VI (G1 > G3) (Table 3), while the PVF differed significantly (G1 > G3) in the hind limbs (Table 4).
The inter-group lengths differed in both the forelimbs (G1 < G3) and the hind limbs (G1 < G2 and G3) (Table 5). The velocity and the stride frequency did not differ significantly among the groups (Table 6), which was in contrast to a significant difference in body mass (G1 < G2 < G3) (Table 7). The relative velocity differed in both the forelimbs and hind limbs (G1 > G3). The average Froude numbers were 0.27 (G1), 0.25 (G2) and 0.23 (G3) for the forelimbs and 0.25 (G1), 0.24 (G2), and 0.23 (G3) for the hind limbs.
The duty factors for the complete study population ranged from 0.55 to 0.62, with a mean of 0.58 (SD 0.01), for the forelimbs, and from 0.58 to 0.66, with a mean of 0.61 (SD 0.02), for the hind limbs.
In our study, characteristics that are considered unfavourable for performing gait analysis in sheep, such as flight zone and flocking behaviour , were reduced by training, as has been previously reported , and with the same halter when inducing the animals to walk on a pressure-sensing walkway. Another important strategy was using food to motivate the sheep to traverse the pressure-sensing walkway in a straight line. Moreover, the body mass (mean 33.10 kg, range 19.5-45.0 kg, comparable to that of a large-sized dog) of the sheep used in our study facilitated handling.
Although sheep reach sexual maturity at approximately 7–12 months of age, depending on the breed, the closure of the physeal plates of the long bones may occur as late as 36 months . Thus, it is important to evaluate sheep of different ages, as in the present study, when evaluating skeletal growth.
Velocity measurements obtained by photoelectric cells and measurements from pressure-sensing walkway have been shown to be similar . For this reason, only the latter measurement method was used. The designated software calculated the velocity by dividing the distance between consecutive foot strikes by the time between them [7, 10, 21]. The velocity used, 1.1-1.3 m/s, was considered comfortable for the sheep. The mean velocity used in a study assessing healthy sheep walking on a pressure-sensing walkway was 1.06 m/s . In healthy dogs, the walking velocity on the same or similar type of pressure-sensing walkway has been reported to range from 0.5 to 1.14 m/s, depending on body sizes [10, 21–24]. For healthy cats, the mean velocity has been reported to range from 0.6 to 0.81 m/s [25–27].
The duty factor in this study was >0.5, with a mean of 0.58 for the forelimbs and a mean of 0.61 for the hind limbs. In a study of healthy Suffolk-mix sheep, the duty factors were 0.66 and 0.69 for the forelimbs and hind limbs, respectively . As has been reported previously , these values suggest a walking speed, as a duty factor >0.5 in the hind limbs is indicative of walking and a factor <0.5 is indicative of trotting or running .
Velocity and acceleration must be controlled because of their influence over the stance time, which is associated with the VI [2, 10, 29]. Studies of dogs have reported that velocity variations of 0.3 m/s may modify the ground reaction forces; as the velocity increased, the PVF increased and the VI decreased [4, 30]. In the present study, the velocity was similar among the groups and had a mean value of 1.17 m/s (SD 0.02). In another sheep study, the mean velocity was 1.06 m/s, but the velocity varied from 0.57 to 1.49 m/s in the forelimbs and from 0.57 to 1.76 m/s in the hind limbs .
The two directions in which the sheep walked on the pressure-sensing walkway did not show significant differences, suggesting that the side of the handler did not interfere with the values. The direction of locomotion relative to the camera did not influence any of the temporospatial or kinetic parameters evaluated in a prior sheep study .
In humans, the gait of children is different from that of adults. Characteristics such as the cycle time, stride length and velocity change with growth . In this study, no differences were detected in the temporospatial parameters among the groups, but the sheep moved at the same velocity. Stride length is diminished in healthy elderly humans [2, 31]. This finding was not observed in Group 3, probably because the animals were younger adults. Sheep life expectancy has been reported to be approximately 10 to 12 years .
In addition, a relationship between stride length and height has been found in humans . No stride length differences were found among the groups in our study, although the G1 animals had shorter limbs than those in other groups. The 8.2 cm and 6.7 cm disparities in the forelimbs and hind limbs, respectively, may not have been sufficient to cause detectable temporospatial differences. Additionally, a kinematic study of horses found that the duration of stance and stride in foals can be normalised using linear or dynamic scaling by the height at the withers and used to predict the values observed in adulthood. However, the height at the withers of the foals had increased by 30 cm when they reached adulthood .
In general, the stance and swing phases account for approximately 60% and 40%, respectively, of the gait cycle in healthy humans . In the present study, the distribution was 58.86% stance phase and 41.49% swing phase for the forelimbs versus 61.88% stance phase and 38.85% swing phase for the hind limbs. Another study of sheep reported a mean of 66.31% for the forelimbs and 68.89% for the hind limbs in stance phase and 33.69% for the forelimbs and 31.11% for the hind limbs in swing phase .
The walking PVI (%BW) is generally higher for the forelimbs than the hind limbs when measured using a pressure-sensitive walkway [10, 19, 22, 26, 34]. Although this finding is consistent with the results from all of the groups in our study (mean values of 47.35% for the forelimbs and 28.58% for hind limbs), the values showed some variations compared to previous studies of healthy sheep (means of 52.52% and 38.52% for the forelimbs and hind limbs, respectively) , healthy dogs (means of 58.11% and 42.30% for the forelimbs and hind limbs, respectively) , and healthy cats (means of 48.2% and 38.3% for the forelimbs and hind limbs, respectively) .
The differences in velocity and calibration  and the animal species in question may have contributed to these differences in reported values. Differences in the calibration protocol may influence the number of activated sensors and result in differing readings [14, 26, 35]. In the present study, the calibration was performed according to the manufacturer’s guidelines using a phantom (a short three-legged stool). Weight was added to the stool to match the weight of each animal, and hard rubber was attached to the bottom of the feet of the stool to mimic a sheep’s hoof. Few studies have specified which methods were used to calibrate their pressure-sensitive walkways [7, 11, 14, 26], which may explain the different measurements. In a study that evaluated vertical limb forces in dogs at a trotting velocity, for example, a 65 kg human subject placed one foot on the mat for calibration . In a kinetic evaluation of cats, the mat was calibrated by having a 50 kg human subject stand on it with both bare feet .
The body weight distribution among the limbs of healthy sheep and dogs during walking has been described as approximately 30% on each forelimb and 20% on each hind limb [19, 36]. Similar values were observed in our study, with 31.34% of the body weight placed on each forelimb and 18.79% placed on each hind limb.
In a study of healthy dogs walking on a force plate, the peak vertical forces were inversely correlated with physical size . A similar correlation was observed in the present study for the forelimbs (R = −0.57) and hind limbs (R = −0.44) Thus, the PFV values for the forelimbs and hind limbs were greater in Group 1 than in Group 3, with Group 3 being composed of larger sheep. Although body size of the Group 2 was similar to Group 3, no significant PFV differences were observed. Furthermore, the relative velocity was greater in Group 1 than in Group 3, which may have contributed to the differences in the PVF values.
However, the VI of the forelimbs was higher in the G1 animals than in the G2 animals. The velocity may influence the stance time and consequently the VI , as the impulse reflects the association between force and the time that the foot is on the ground . As a constant velocity was maintained in our study and as the stance time was statistically indistinguishable among the groups, the differences were associated with the different forces. These results are in contrast with those from another canine study in which the VI increased as the size of the dog increased , probably due to differences in velocity and individual dog size.
The vertical forces created by young healthy sheep walking on a pressure-sensing walkway at the same velocity differ from those of older sheep under similar circumstances. This finding may be a source of variation and should be controlled in locomotion research studies.
McLaughlin RM: Kinetic and kinematic gait analysis in dogs. Vet Clin North Am Small Anim Pract. 2001, 31: 193-201.
Whittle MW: Gait analysis: an introduction. Elsevier, Edinburgh; 2007:47-193.
Gillette RL, Angle TC: Recent development in canine locomotor analysis. A review. Vet J. 2008, 178: 165-176. 10.1016/j.tvjl.2008.01.009.
DeCamp CE: Kinetic and kinematic gait analysis and the assessment of lameness in the dog. Vet Clin North Am Small Anim Pract. 1997, 27: 825-840.
Anderson MA, Mann FA: Force plate analysis: a noninvasive tool for gait evaluation. Comp Cont Ed Pract Vet. 1994, 16: 857-864.
Budsberg SC, Thomas MW: Advanced diagnostic techniques. BSVA Manual of canine and feline musculoskeletal disorders. Edited by: Houlton JEF, Cook JL, Innes JF, Langley-Hobbs SJ. British Small Animal Veterinary Association, England; 2006:27-33.
Lascelles BD, Roe SC, Smith E, Reynolds L, Markham J, Marcellin-Little D, Bergh MS, Budsberg SC: Evaluation of a pressure walkway system for measurement of vertical limb forces in clinically normal dogs. Am J Vet Res. 2006, 67: 227-282.
Ishihara A, 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: 1805-1815. 10.2460/ajvr.2005.66.1805.
Fanchon L, Grandjean D: Accuracy of asymmetry indices of ground reaction forces for diagnosis of hind limb lameness in dogs. Am J Vet Res. 2007, 68: 1089-1094. 10.2460/ajvr.68.10.1089.
Besancon MF, Conzemius MG, Derrick TR, Ritter MJ: Comparison of vertical forces in normal greyhounds between force platform and pressure walkway measurement systems. Vet Comp Orthop Traumatol. 2003, 16: 153-157.
Oosterlinck M, Pille F, Huppes T, Gasthuys F, Back W: Comparison of pressure plate and force plate gait kinetics in sound Warmbloods at walk and trot. Vet J. 2010, 186: 347-351. 10.1016/j.tvjl.2009.08.024.
Oosterlinck M, Pille F, Back W, Dewulf J, Gasthuys F: A pressure plate study on fore and hindlimb loading and the association with hoof contact area in sound ponies at the walk and trot. Vet J. 2011, 190: 71-76. 10.1016/j.tvjl.2010.08.016.
Oosterlinck M, Pille F, Sonneveld DC, Oomen AM, Huppes T, Gasthuys F, Back W: Contribution of dynamic calibration to the measurement accuracy of a pressure plate system throughout the stance phase in sound horses. Vet J. 2012, in press. http://dx.doi.org/10.1016/j.tvjl.2012.01.029
Oosterlinck M, Bosmans T, Gasthuys F, Ingeborgh P, Van Ryssen B, Dewulf J, Pille F: Accuracy of pressure plate kinetic asymmetry indices and their correlation with visual gait assessment scores in lame and nonlame dogs. Am J Vet Res. 2011, 72: 820-825. 10.2460/ajvr.72.6.820.
Tekscan Boston. http://www.tekscan.com/medical/system-walkway.html.
Martini L, Fini M, Giavaresi G, Giardino R: Sheep model in orthopedic research: a literature review. Comp Med. 2001, 51: 292-299.
Pearce AI, Richards RG, Milz S, Schneider E, Pearce SG: Animal models for implant biomaterial research in bone: a review. Eur Cells Mater. 2007, 13: 1-10.
Potes JC, Reis JC, Silva FC, Relvas C, Cabrita AS, Simoes JA: The sheep as an animal model in orthopaedic research. Exp Pathol Health Sciences. 2008, 2: 29-32.
Kim J, Breur GJ: Temporospatial and kinetic characteristics of sheep walking on a pressure sensing walkway. Can Vet J. 2008, 72: 50-55.
Seebeck P, Thompson MS, Parwani A, Taylor WR, Schell H, Duda GN: Gait evaluation: A tool to monitor bone healing?. Clin Biomech. 2005, 20: 883-891. 10.1016/j.clinbiomech.2005.05.010.
Kim J, Kazmierczak KA, Breur GJ: Comparison of temporospatial and kinetic variables of walking in small and large dogs on a pressure-sensing walkway. Am J Vet Res. 2011, 72: 1171-1177. 10.2460/ajvr.72.9.1171.
Besancon MF, Conzemius MG, Evans RB, Ritter MJ: Distribution of vertical forces in the pads of greyhounds and Labrador retrievers during walking. Am J Vet Res. 2004, 65: 1497-1501. 10.2460/ajvr.2004.65.1497.
LeQuang T, Maitre P, Roger T, Viguier E: Is a pressure walkway system able to highlight a lameness in dog?. J Anim Vet Adv. 2009, 8: 1936-1944.
Light VA, Steiss JE, Montgomery RD, Rumph PF, Wright JC: Temporal-spatial gait analysis by use of a portable walkway system in healthy Labrador Retrievers at a walk. Am J Vet Res. 2010, 71: 997-1002. 10.2460/ajvr.71.9.997.
Romans CW, Gordon WJ, Robinson DA, Evans R, Conzemius MG: Effect of postoperative analgesic protocol on limb function following onychectomy in cats. J Am Vet Med Assoc. 2005, 227: 89-93. 10.2460/javma.2005.227.89.
Lascelles BD, Findley K, Correa M, Marcellin-Little D, Roe S: Kinetic evaluation of normal walking and jumping in cats, using a pressure sensitive walkway. Vet Rec. 2007, 160: 512-516. 10.1136/vr.160.15.512.
LeQuang T, Maitre P, Colin A, Viguier E: Evaluation spatial-temporal and pressure parameters of normal cats at a walk, using a pressure walkway. Proceedings of the Third International Conference on the Development of Biomedical Engineering in Vietnam. Edited by: Magjarevic R. 2010, Vietnam, 59-61.
Vilensky JA: Locomotor behavior and control in human and non-human primates: comparisons with cats and dogs. Neurosci Biobehav Rev. 1987, 11: 263-274. 10.1016/S0149-7634(87)80013-1.
McLaughlin RM, Roush JK: Effects of subject stance time and velocity on ground reaction forces in clinically normal Greyhounds at the trot. Am J Vet Res. 1994, 55: 1666-1671.
Riggs CM, DeCamp CE, Soutas-Litlle RW, Braden TD, Richter MA: Effects of subject velocity on force plate-measured ground reaction forces in healthy Greyhounds at the trot. Am J Vet Res. 1993, 54: 1523-1526.
Elble RJ, Thomas SS, Higgins C, Colliver J: Stride-dependent changes in gait of older people. J Neurol. 1991, 238: 1-5. 10.1007/BF00319700.
Khan BB, Iqbal A, Mustafa MI: Sheep and goat production. Department of Livestock Management: University of Agriculture Faisalabad 2003.
Back W: Predictive value of foal kinematics for the locomotor performance of adult horses. Res Vet Sci. 1995, 59: 64-69. 10.1016/0034-5288(95)90033-0.
Romans CW, Conzemius MG, Horstman CL, Gordon WJ, Evans RB: Use of pressure platform gait analysis in cats with and without bilateral onychectomy. Am J Vet Res. 2004, 65: 1276-1278. 10.2460/ajvr.2004.65.1276.
Brimacombe JM, Wilson DR, Hodgson AJ, Ho KCT, Anglin C: Effect of calibration method on Tekscan sensor accuracy. J Biomech Eng. 2009, 131: 034503-1-034503-4.
Budsberg SC, Verstraete MC, Soutas-Little RW: Force plate analysis of the walking gait in healthy dogs. Am J Vet Res. 1987, 48: 915-918.
Alexander RM, Jayes AS: A dynamic similarity hypothesis for the gaits of quadrupedal mammals. J Zool Lond. 1983, 201: 135-152.
The authors would like to thank the FAPESP (The State of São Paulo Research Foundation) and CNPq (National Council for Scientific and Technological Development) for financial support, the Coordination for the Improvement of Higher Level Personnel (CAPES) - PROCAD-NF No. 21/2009, and Dr. Guy Beauchamp of the Faculty of Veterinary Medicine, University of Montreal – Canada, for the statistical analysis. We would also like to thank the Universidade Estadual Paulista (UNESP) and the Universidade Federal Rural da Amazônia.
FSA and SCR conceived and designed the study; FAPA and RTC helped collected the data; CAH and AOE assisted in data collection methods and FOBM helped draft the manuscript; all authors read, contributed to and approved the final manuscript.