|
| ![]() |
|
||
| Home Help Feedback Subscriptions Archive Search Table of Contents | ||||
First published online January 30, 2009
Journal of Experimental Biology 212, 471-482 (2009)
Published by The Company of Biologists 2009
doi: 10.1242/jeb.026377
Estimating energy expenditure of animals using the accelerometry technique: activity, inactivity and comparison with the heart-rate technique

1 School of Biological Sciences, University of Liverpool, Crown Street,
Liverpool, L69 7ZB, UK
2 School of Human and Life Sciences, Roehampton University, Holybourne Avenue,
London, SW15 4JD, UK
3 Institute of Environmental Sustainability, School of the Environment and
Society, University of Swansea, Singleton Park, Swansea, SA2 8PP, UK
4 Department of Zoology, La Trobe University, Bundoora, Melbourne, Victoria
3070, Australia
* Author for correspondence (e-mail: jonathan.green{at}liverpool.ac.uk)
Accepted 27 November 2008
| Summary |
|---|
|
|
|---|
Key words: accelerometry, heart rate, energetics, chicken, SDA, thermoregulation
| INTRODUCTION |
|---|
|
|
|---|
In addition to these two established methods, a newer method has been
developed, in which the overall dynamic body acceleration (ODBA) of an animal
is calculated and used as a calibrated proxy for the rate of energy
expenditure (Wilson et al.,
2006
). Using body motion as a proxy for energy expenditure is not
a new concept – studies in the medical literature date back to the 1960s
(e.g. Cavagna et al., 1963
).
However, its application to animal biology is a relatively new one,
facilitated by the ability of the latest generation of miniature data loggers
to record and store acceleration data in three axes at high frequencies
(Wilson et al., 2008
).
However, like all similar techniques, using ODBA has limitations. By
definition, this technique relies on recording ATP-fuelled movements of
animals and assumes that energy use increases with movement. Energy expended
in movement can represent the major portion of the energy budget of mammals
and birds (Weibel and Hoppeler,
2005
). Thus ODBA should provide a good estimate of energy
expenditure during periods of activity.
However, animals can spend hours not moving, yet they can be engaged in
energetically expensive physiological processes such as growth, production
(e.g. eggs or milk), digestion, thermoregulation or, conversely, energy-saving
processes such as torpor. For example, elephants spend less than 20% of their
time walking (Shannon et al.,
2008
), seabirds commonly spend over 50% of their day resting
(Falk et al., 2000
;
Grémillet et al., 1995
)
even when busy provisioning nestlings, and red deer rest for at least 45% of
the day, depending on the time of year
(Pépin et al., 2006
).
So while the relationship between ODBA and energy expenditure has been
investigated during activity for a number of species
(Fahlman et al., 2008
;
Halsey et al., 2008a
;
Halsey et al., 2008b
;
Halsey et al., 2008c
;
Wilson et al., 2006
), it is
well recognised that, during periods of inactivity, estimates of energy
expenditure might need to be derived from laboratory studies or allometric
estimates (Wilson et al.,
2008
).
In the present study, we establish whether body acceleration can be used to estimate the rate of energy expenditure across the full range of animal behaviours. As well as making the animals exercise (walk), we investigated whether an accelerometer might detect fine-scale movements associated with thermoregulation such as panting or shivering or, for example, whether specific dynamic action during the digestion of a meal increases energy expenditure without changing measures of body acceleration. We used the domestic chicken (Gallus gallus) as a model species as they are easily maintained in captivity and respond favourably to a variety of experimental manipulations. We also examined the relationship between heart rate and rate of energy expenditure in the same individuals across the same range of behaviours. This enabled the effectiveness of the accelerometry technique during different behaviours to be compared with that of a more established methodology.
Thus we set out to answer the following four questions:
| MATERIALS AND METHODS |
|---|
|
|
|---|
Respirometry
The rate of oxygen consumption
(
O2) is an
indirect measure of energy expenditure, assuming no anaerobic respiration, and
was determined using an open-circuit respirometry system (see
Green et al., 2008
). During
experiments, a bird was placed within a respirometer chamber
(470x480x433 mm high) made of clear plastic, which was suspended
above a variable-speed treadmill (Australian Treadmill Exerciser). Draught
excluders partially sealed the respirometer to the treadmill. Air was drawn
from the respirometer using a variable-speed pump (Recipritor) at
5 l
min–1. The flow rate was calculated before and after
experiments using a water-displacement technique. A sub-sample of air was
drawn from this main flow using a small air pump (Ametek R-1, Applied
Electrochemistry), passed through a drying column (Drierite, Hammond), and
analysed for the fractional content of O2 and CO2 by two
gas analysers (Ametek S-3A/I and ADI Instruments ML205). The outputs from the
gas analysers and a thermocouple located inside the respirometer were
collected at 100 Hz (Powerlab 4SP, ADInstruments) and displayed on a computer
using Chart software (ADInstruments). The rate of oxygen consumption was
determined from the rate of airflow from the respirometer and the difference
in the fractional oxygen concentration between ambient and out-flowing air.
Instantaneous corrections of the gas concentrations were calculated dry at
standard temperature (273 K) and pressure (101.3 kPa) using the method of
Frappell and colleagues (Frappell et al.,
1989
), assuming a first-order linear system (chamber volume=101 l;
flow=5 l min–1; tau=17 min, determined from a
semi-logarithmic plot of concentration against time following a perturbation,
R2=0.99).
O2 was
calculated with consideration of RQ-related errors
(Frappell et al., 1992
). Mean
O2 was
calculated every 5 min during experiments.
Heart rate and partial dynamic body acceleration
Heart rate (fH) was monitored using a customised
heart-rate transmitter (POLAR a3, Polar Electro Oy, Finland) with largest
dimensions 50x20x8 mm and a mass of 21 g. The transmitter was
attached to the upper back of the bird with paper tape (see
Wilson and Wilson, 1989
), and
custom-made brass electrodes were inserted under the skin. The transmitter
unit had a functional range of
1 m, and a receiver unit was placed on top
of the respirometer to ensure good reception. Outputs from the receiver unit
were collected by the Powerlab alongside respirometry data. Acceleration was
measured using the same data loggers as Wilson and colleagues
(Wilson et al., 2006
) and
Halsey and colleagues (Halsey et al.,
2008c
), also attached to the upper back with paper tape. The
loggers had largest dimensions of 42x36x13 mm, a mass of 24 g and
were set to record tri-axial acceleration (0–6 g) at 10
Hz with 22-bit resolution. This recording frequency is sufficiently high when
using measures of acceleration as a proxy for energy expenditure in this
species (Halsey et al.,
2008a
).
The absolute values (in g) for each axis were calculated
[abs(x), abs(y) and abs(z)], to represent the raw
accelerometry values. The x axis of the logger measured acceleration
laterally (coronal) – that is, wing to wing – and hence sway. The
y axis of the logger measured acceleration in the plane at right
angles to the x axis, approximately in the anteroposterior dimension
of the chicken (sagittal plane) – that is, head to tail – and
hence surge. The z axis measured acceleration in approximately the
dorsoventral dimension (transverse) and hence heave. From the downloaded
logger data, x, y and z, an approximation of absolute
g, resulting from only dynamic acceleration in that dimension,
was extracted (partial dynamic body acceleration, PDBAx,
PDBAy and PDBAz, respectively).
However, initial inspection of the data revealed unacceptable noise and
interference in the y-axis due to instrument failure, and ODBA was
not calculated. Instead, we summed PDBAx and
PDBAz, calculating partial dynamic body acceleration in
the x and z axes (PDBAxz).
PDBAxz is an accurate proxy for energy expenditure in
walking chickens (Halsey et al.,
2008a
).
Eating–digestion protocol
Eating and digesting were examined in a single protocol. Birds were
initially introduced to the respirometer and left until calm (<10 min). A
small dish containing 30 g of food pellets was then placed into the
respirometer. The birds immediately began eating and were left undisturbed for
10 min while their behaviour was recorded. The chickens had usually finished
eating by the end of this period, because either the dish was empty or they
were satiated. The dish was then removed, the amount of food consumed
calculated, the lights switched off, and the birds left undisturbed for
2
h to record specific dynamic action (SDA) during digestion. All experiments
took place at a constant 18°C.
Thermoregulation protocol
The rate of oxygen consumption and PDBAxz or
fH were measured in the dark at temperatures of: 1, 11,
18, 30 and 36°C. A chicken only experienced a single temperature on any
given day to avoid any complications of hysteresis in metabolism, heart rate
or body temperature. The temperature range was based on data from larger birds
(Meltzer, 1983
). In our
smaller chickens, we assumed that 1, 11 and 36°C would be outside the
thermoneutral zone (TNZ), whereas 18 and 30°C would be close to the lower
and upper critical temperatures, respectively. Birds were equipped with a
heart-rate transmitter or accelerometer and acclimatised to the experimental
temperature in a ventilated box for
1 h before being transferred to the
respirometer for recordings of
O2 and
PDBAxz or fH. The birds remained in
the respirometer for at least one hour, which allowed fH
and metabolism to settle to a steady level and resting rates to be
established.
Walking protocol
Birds, equipped with an accelerometer or heart-rate transmitter, were
introduced into the respirometer and allowed to rest for at least 15 min.
Birds were then run at randomly assigned speeds between 0.15 and 1.3 km
h–1 for between 8 and 16 min. The birds were allowed to rest
between speed adjustments. To encourage movement, food pellets were placed
outside the front of the respirometer; thus, some birds also exhibited pecking
behaviours. Indeed, other natural active behaviours were also undertaken,
particularly at the lower walking speeds, such as wing flapping and jumping.
Thus, the chickens exhibited a range of natural, active behaviours during the
trials. Again, experiments were conducted at 18°C. Some of these data on
dynamic body acceleration during walking have been reported previously
(Halsey et al., 2008a
;
Halsey et al., 2008c
).
Data analysis
Question 1 asked whether PDBAxz would change during
eating and digesting. First, 5 min periods were selected where the chicken was
eating continuously throughout and mean
O2,
PDBAxz and fH were calculated. Second,
during digestion, mean
O2,
PDBAxz and fH of the first 30 min
after eating had finished were selected for analysis as initial visual
inspection confirmed that the largest increase in metabolism was seen at this
time (e.g. Green et al.,
2006
). Data during both eating and digestion were compared with
the equivalent data obtained at the same temperature (18°C) during the
thermoregulation protocol using paired t-tests.
To investigate question 2, regarding whether changes in temperature would
induce a change in
O2,
PDBAxz and fH, the lowest five
recordings of
O2
were selected for each temperature for each bird, both during trials with the
heart-rate transmitter and those with the accelerometer. Two-way analysis of
variance with Tukey post-hoc comparisons were used to
investigate differences in
O2,
PDBAxz or fH between temperatures.
During walking, data were usually recorded for more than one 5 min period
at each speed for each animal. To determine the relationship between speed and
O2,
PDBAxz or fH, a mean was taken for
each animal at each speed. A grand mean and s.e.m. of these means was then
calculated for each speed, and weighted regression was used to look at
relationships between treadmill speed and each variable.
Subsequent analyses examined relationships between either
PDBAxz and
O2 or
fH and
O2 across all
behaviours in order to answer question 3. All 5 min periods of data from all
protocols were used in these analyses. These data were analysed using general
linear models (GLM) in which behaviour was a fixed factor, individual identity
was a random factor and PDBAxz or fH
was a covariate. A significant interaction between behaviour and
PDBAxz or fH indicated differences in
the slope of the relationship between PDBAxz or
fH and
O2 between the
different behaviours. Subsequent to this, in some analyses, behaviour was also
a random factor.
Predictive relationships were generated from these GLMs and used to answer
question 4 and compare the predictive power of relationships between
PDBAxz and
O2, or
fH and
O2. The standard
error of the estimates (s.e.e.) made using these equations were calculated
following the conventions established by Green and colleagues
(Green et al., 2001
). As Green
and colleagues (Green et al.,
2001
) outline, the s.e.e. and hence 95% prediction intervals of an
estimate made using a predictive equation will depend on the characteristics
of (1) the individuals used to derive the predictive relationship and (2) the
individuals for which further measurements of the independent variable were
made and for whom an estimate is being calculated (`sample' group). In
practice, this means that the number of individuals in the `sample' group and
number of measurements taken from those individuals will have a substantial
effect on the 95% prediction intervals surrounding an estimate
(Green et al., 2001
). When
dealing with predictive equations, the 95% confidence intervals associated
with an estimate are the 95% prediction intervals where a predictive equation
is used to make an estimate based on one further measurement of the
independent variable from one further individual. Consequently, the 95%
confidence intervals represent the maximum error range around the prediction.
Z-tests or proximate normal tests (for paired estimates) were used to
compare estimates made using predictive equations
(Green et al., 2001
).
|
|
| RESULTS |
|---|
|
|
|---|
|
|
O2 (paired
t-test; t(7)=2.98, P<0.05) when
compared with resting quietly at the same temperature (18°C)
(Fig. 2). Similarly, when
fH was recorded, both fH (paired
t-test; t(7)=7.27, P<0.01) and
O2 (paired
t-test; t(7)=5.28, P<0.01) were
significantly greater during eating. After eating, the magnitude of SDA during
digestion was evaluated (Fig.
3). When PDBAxz was recorded,
O2 was
significantly higher during digestion (paired t-test;
t(7)=3.24, P<0.05), but there was no
difference in PDBAxz (paired t-test;
t(7)=1.28, P=0.24). When fH
was recorded, both
O2 (paired
t-test; t(7)=3.78, P<0.01) and
fH (paired t-test; t(7)=3.66,
P<0.01) were significantly higher during digestion than while the
chickens rested at the same temperature.
Body acceleration during thermoregulation
When exposed to a range of different temperatures, the chickens showed
changes in metabolism, suggesting a TNZ with an upper critical temperature
around 30°C and a lower critical temperature as low as 20°C
(Fig. 4). When instrumented
with an accelerometer, the rate of oxygen consumption
(
O2) was at a
minimum at 30°C and was significantly elevated at 1°C and 11°C
(two-way ANOVA, with Tukey multiple comparisons;
F(4,28)=7.29, P<0.001). There were no
significant differences in PDBAxz between the different
temperatures (two-way ANOVA; F(4,28)=0.32,
P=0.86). When
O2 and heart
rate (fH) were recorded simultaneously, both
O2 and
fH were again at a minimum at 30°C, and both were
significantly elevated at 1°C and 11°C (two-way ANOVA, with Tukey
multiple comparisons; F(4,28)=44.3 and 22.7,
P<0.001).
Body acceleration during locomotion
The chickens exhibited a range of behaviours as well as walking while
undertaking the treadmill experiments. They also pecked the ground, jumped and
flapped their wings. During these trials, there was a positive, linear
relationship between
O2 and walking
speed in bantam chickens instrumented with an accelerometer (weighted
regression; P<0.05, R2=0.80)
(Fig. 5). This was matched by a
linear increase in PDBAxz with walking speed (weighted
regression; P<0.001, R2=0.96). Similarly when
O2 and
fH were recorded simultaneously, there was a linear
relationship between walking speed and both
O2 (weighted
regression; P<0.01, R2=0.90) and
fH (weighted regression; P<0.001,
R2=0.97) (Fig.
5).
|
Predicting the rate of oxygen consumption as a function of body acceleration
The behaviour-specific tests described above suggested that, although
PDBAxz might be an effective predictor of
O2 during active
behaviours such as locomotion or eating, it was likely to be less effective
while the chickens were inactive, either thermoregulating or digesting a meal.
To investigate this further,
O2 was plotted
as a function of PDBAxz for each animal (see
Fig. 6A for example). Visual
inspection confirmed this initial suspicion that PDBAxz
was likely to be an accurate predictor of
O2 during
walking but was unlikely to be so accurate while the chickens were inactive.
However, although there was a lot scatter in PDBAxz and
O2 during
inactivity, these data rarely overlapped with data recorded during activity.
Converting PDBAxz to a logarithmic scale suggested that,
with a single curvilinear model, it might be possible to predict
O2 from
PDBAxz across all behaviours
(Fig. 6C). An analysis of
covariance was conducted to investigate this further.
|
O2 did not
differ between behaviours. Therefore, a single function was used to predict
O2 from
PDBAxz across all behaviours
(Table 1). As both behaviour
and individual identity were significant factors in the model, two extra error
terms would need to be introduced when calculating s.e.e. values made using
this single model. Fig. 7
illustrates this effect and shows how using a single relationship might not be
the best approach. Calculation of 95% confidence intervals and 95% prediction
intervals shows that using this one-model approach will tend to underestimate
O2 during
walking. Even if
O2 were
predicted for a large sample of animals, the error of the estimate would be
considerable for estimates made during walking, despite the relatively close
relationship during this behaviour.
|
|
Further investigation was undertaken with analyses of covariance comparing
digestion with thermoregulation (`inactive' states), and walking with eating
(`active' states), to investigate the creation of two models. In the case of
thermoregulation and digestion, there was again no significant interaction
between behaviour and log(PDBAxz)
(F(1,588)=2.44, P=0.12) and a significant model
for inactivity was created (Table
2). While the model was relatively inaccurate
(R2=0.21), the relationship between
PDBAxz and
O2 during
inactivity was statistically significant. Similarly, there was no significant
interaction between behaviour and log(PDBAxz) when eating
and walking were compared (F(1,93)=0.22, P=0.64).
Furthermore, there was no significant effect of behaviour in this model
(F(1,108)=1.95, P=0.17), and so a model with only
one additional error term was created
(Table 3). Plotting predictions
made with this two-model approach (Fig.
8) shows that, although the 95% prediction intervals are still
relatively large during inactivity, they are substantially smaller during
activity than they are in the one-model approach. Furthermore, there is no
longer a systematic underestimation of
O2 during
walking, as in the one-model approach.
|
|
|
Predicting the rate of oxygen consumption as a function of heart rate
Heart rate appeared to be an effective predictor of
O2 during all
behaviours. To investigate this further,
O2 was plotted
as a function of fH during all behaviours for each animal
(see Fig. 6B for example).
Again, a logarithmic relationship appeared to best fit the data across all
behaviours. However, analysis of covariance indicated a significant
interaction between behaviour and log(fH)
(F(3,699)=5.23, P<0.001). Eliminating each of
the behaviours in turn and repeating this analysis revealed that the
relationship between fH and
O2 during
digestion was significantly different to the relationship between the other
three (Fig. 9). Thus, two
models were created: one during SDA and the other for all other behaviours. In
this two-model approach, there was a considerable overlap between the two
relationships (Fig. 9), and so
a one-model approach was also investigated. In this case, a single
relationship was constructed to enable prediction where behaviour was not
known, with behaviour included as an additional error term in the calculation
of the s.e.e. (Table 4;
Fig. 10). Despite adding this
potential for increased uncertainty, there was very little difference in the
95% prediction intervals when comparing the one-model and two-model approaches
(Figs 9 and
10).
|
|
|
| DISCUSSION |
|---|
|
|
|---|
O2 during
activity is not surprising. We reported this result after preliminary analysis
of some of the data collected in this study
(Halsey et al., 2008a
Question 1 – dynamic body acceleration during eating and digestion
Comparing PDBAxz and
O2 while eating
with the same measures during digestion neatly illustrates the potential
limitations in using PDBAxz to predict
O2 for any
particular behaviour. While eating, the chickens made repeated pecking
movements where the body pivoted from the hip, and the neck extended to reach
the food. This induced a 76% increase in
O2 and
commensurate 141% increase in PDBAxz
(Fig. 2), thus the difference
was easily detectable – an example of an increase in energy expenditure
being well detected by PDBAxz during active behaviours.
Furthermore, it should be possible to detect these non-locomotory active
behaviours from analysis of characteristics of the acceleration data
(Shepard et al., 2008
).
However, while the birds were digesting, no movement was seen by observers,
yet
O2 increased
by 42% from resting levels. Blaxter
(Blaxter, 1989
) notes that the
increase in metabolism during digestion is attributable predominantly to
upregulated biochemical processes, which one would not expect to be reflected
in increased body movement. Thus, although the increase in energy expenditure
was similar to that found in other birds
(Chappell et al., 1997
;
Green et al., 2006
;
Kaseloo and Lovvorn, 2003
;
Sedinger et al., 1992
), it was
accompanied by an increase in PDBAxz of only 10%, which
was not statistically significant.
Question 2 – dynamic body acceleration during thermoregulation
As with digestion, despite clear differences in energy expenditure
associated with thermoregulation, there was no significant effect of ambient
temperature on PDBAxz. In some respects, this is a
surprising finding, particularly with regard to ambient temperature below the
apparent lower critical temperature of the chickens. In the present study, we
did not determine the exact boundaries of the TNZ, but
Fig. 4 shows that it would be
in the range of 22–30°C, similar to that reported for larger broiler
chickens (Meltzer, 1983
).
Measurements made at 1°C and 11°C were clearly below the lower
critical temperature, as shown by the increase in energy expenditure compared
with that at 30°C. At 1°C and to a lesser extent at 11°C,
O2 was lower
while birds were equipped with an accelerometer than when they were equipped
with a heart-rate transmitter (Fig.
4). The experiments with the fH transmitter
were conducted first, and, even after repeating the measurements with the
accelerometer attached, we achieved the same result. Thus, we attributed this
slight difference to the chickens acclimatising to the cooler conditions in
their accommodation as the experiments progressed through the winter.
Although the potential for non-shivering thermogenesis has been
demonstrated (Toyomizu et al.,
2002
), the primary mechanism for regulatory thermogenesis in birds
is shivering (Dawson and Whittow,
2000
); thus, it is surprising that no increase in body movement
was registered by the accelerometer. Body tremors during shivering below the
LCT have been detected in small mammals and correlate well with ambient
temperature and metabolic rate (May,
2003
). Presumably, the magnitude of movement or tremors associated
with shivering in the chickens was either too small when compared with static
elements of acceleration or too low to be detected by the accelerometer.
Alternatively, extremely high-frequency movements might not have been detected
with our 10 Hz sampling protocol. In our experiments (and most likely
applications in birds), the accelerometer was fixed to the feathers with tape
(Wilson and Wilson, 1989
). It
is possible that this might have served as a buffer to the high-frequency
movements probably associated with shivering and that movements of the body
did not transmit through the feathers to the accelerometer. This could be
tested by using a harness attachment, but this is generally not recommended
for birds (Phillips et al.,
2003
). Finally, owing to instrument failure, we could not detect
small measures of acceleration in the y axis (surging acceleration)
and therefore could not calculate ODBA during our experiments. Although it
seems unlikely, it is possible that, during both thermoregulation (and even
digestion), noteworthy movement could have occurred exclusively in this plane
that we failed to detect.
Question 3 – dynamic body acceleration: comparing inactivity with walking
Both PDBAxz and
O2 increased
significantly with increasing walking speed
(Fig. 5). Coupled with the
significant changes in PDBAxz and
O2 while eating,
PDBAxz was evidently an excellent predictor of
O2 during
activity. The coefficient of determination (R2=0.77)
(Table 3) of this relationship
compared favourably with those from walking great cormorant (Phalacrocorax
carbo, R2=0.81) (Wilson et
al., 2006
), human (R2=0.91, 0.77)
(Halsey et al., 2008b
) and
coypu (Myocastor coypus, R2=0.91)
(Halsey et al., 2008c
).
During inactivity (resting during digestion and thermoregulation), we still
detected a relationship between PDBAxz and
O2 when all the
data were considered together. While the individual trials examining digestion
and thermoregulation did not detect these effects (see above), when all the
data were considered, there must have been small amounts of movement
correlated with both PDBAxz and
O2. The
coefficient of determination was much lower (R2=0.21)
(Table 2) during activity,
which can be attributed largely to intra-individual variation in
O2 at a given
value of PDBAxz and the small range of
PDBAxz during inactivity
(Fig. 6).
It is perhaps not surprising that PDBAxz predictions
are less precise for inactive behaviours than are fH
predictions. Although changes in heart rate can occur without commensurate
increases in metabolism (Blix et al.,
1974
), fH is still a more direct measure of
metabolism. The coefficient of determination of a combined relationship
between PDBAxz and
O2 was still
quite high (R2=0.69)
(Table 1), but this is
attributable to the range of PDBAxz during activity being
approximately 3–4 times greater than that during inactivity
(Fig. 6). Thus, the active data
have a disproportionate effect on the fit of the single relationship.
The rate of change of
O2 with
increasing PDBAxz is far greater during inactivity than
activity (Fig. 7). This
curvilinear relationship between PDBAxz and
O2 across all
activities and the large amount of intra-individual error coupled with a
truncated range of PDBAxz values during inactivity (see
Fig. 7) both have important
consequences that must be addressed if predictive relationships based on
measures of acceleration are to be used. First, small changes in measurements
of PDBAxz while inactive can lead to potentially large
changes in estimates of
O2; thus, it
will be important that calculations of PDBAxz or ODBA made
from inactive free-ranging animals are accurate and repeatable
(Wilson et al., 2006
). Second,
the relatively large s.e.e. and prediction intervals around estimates of
inactive
O2 will
mean that statistical comparison of two estimates of inactive
O2 are unlikely
to detect a significant difference (e.g.
Fig. 7). For example, a
difference in PDBAxz between resting and digesting could
not be detected (Fig. 3).
Predictions of
O2 made using
these PDBAxz estimates were not significantly different
either (proximate normal test, z=0.38, P=0.70). However,
this should not preclude the application of this method to inactive animals.
Applications of body acceleration currently use laboratory measurements and
time-budgets to estimate energy expenditure while free-ranging animals are
inactive (Wilson et al.,
2008
). Estimating energy expenditure while inactive from body
acceleration might generate large errors around prediction, but these errors
are at least quantifiable and depend on fewer assumptions. Certainly, both of
these methods are likely to be more accurate than using interspecific
allometric formulae as these can have very large errors of prediction, often
ignored or obscured by many orders of magnitude in mass data (e.g.
Nagy, 2005
). For example, an
analysis of data presented by McKechnie and colleagues
(McKechnie et al., 2006
) for
captive-bred animals (comparable to the domestic chickens in this study)
indicates that an estimate of basal metabolic rate for a body mass would have
a coefficient of variation (CV=100*s.e.e./Estimate) of 19.7%.
Furthermore, in many cases, using dynamic acceleration to estimate
O2 might be the
only practicable method available. A calibrated quantitative estimate is also
likely to be more useful than an un-calibrated qualitative measure where only
the magnitude of the proxy can be compared (e.g.
Weimerskirch et al.,
2001
).
In the current study, we concentrated on calibrations made over a
relatively short timescale (5 min). Preliminary work has shown that, if the
timescale over which averages for calibration relationships are constructed is
increased, s.e.e. will decrease (D. Thompson and C. Sparling, personal
communication). Furthermore, this will tend to reduce the number of data
points in the range of PDBAxz found during inactivity,
thus reducing the skew that these points currently exert on the one-model
relationship. However, predictions from a regression should be used only
speculatively to make estimates at a finer timescale than that used in
calibration (Green et al.,
2001
; Halsey et al.,
2008a
), which could limit the temporal resolution over which
O2 can be
estimated in the field. This might mean that it is not possible to calculate
activity-specific costs, but, for many applications where energy expenditure
is predicted on an hourly or daily basis (e.g.
Green et al., 2007
), this
might be the most appropriate approach.
Question 4 – comparing accelerometry with heart rate
Using heart rate to estimate energy expenditure in free-ranging animals is
a well-established technique (Butler et
al., 2004
). Like the use of body acceleration, it requires a
laboratory calibration in order to apply it in the field, and precautions
should be taken that the calibration procedure encompasses the full range of
behaviours undertaken by animals in the field (e.g.
Froget et al., 2002
;
Green et al., 2008
;
Ward et al., 2002
). In the
present study, it was easy to calibrate fH and
O2 in the same
experimental animals as those used in the calibration of
PDBAxz and
O2 and hence
compare these two proxies of energy expenditure. As in previous studies
(Froget et al., 2002
;
Green et al., 2008
),
concomitant changes in fH could be detected alongside
changes in
O2
during digestion (Fig. 3) and
thermoregulation (Fig. 4). The
relationship between fH and
O2 during
digestion was different to that during walking, whereas, in little penguins
(Eudyptula minor), there was no difference
(Green et al., 2008
). By
contrast, the relationship between fH and
O2 for
thermoregulation was the same as that during walking, whereas, in king
penguins (Aptenodytes patagonicus), it was different
(Froget et al., 2002
). These
two findings once again highlight the importance of species-specific
calibration procedures when applying the fH technique.
When comparing the two techniques directly, the coefficient of
determination of the one-model approach was higher for heart rate
(R2=0.84, Table
4 vs R2=0.69,
Table 1). To look at this in
more detail, we repeatedly simulated a day in the life of 10 chickens, where
O2 and s.e.e.
were estimated for 24 h using the one- and two-model approaches for both
PDBAxz and fH. In the simulation, the
amount of time spent active by the chicken was varied between 0 and 100%. We
assumed that, while inactive, PDBAxz and
fH were equivalent to the mean values for digestion
(Fig. 3) and that, while
active, PDBAxz and fH were equivalent
to the mean values while walking at 1 km h–1
(Fig. 5). We assumed that these
measurements had been taken 12 times each hour in each bird.
O2 was estimated
using the appropriate equations for both the one- and two-model approaches,
using both techniques, and s.e.e. and the coefficient of variation
(CV=100*s.e.e./Estimate) plotted as a function of the proportion of
time spent active (Fig.
11).
|
O2 with
approximately half the error of prediction than that for
PDBAxz, independent of the model selected. However, the
s.e.e. and CV of the PDBAxz estimate is not so low that it
could not be used in many applications. For example, CV never exceeds 14%,
which is less than the 19.7% calculated earlier for interspecific allometric
estimation of resting metabolic rate
(McKechnie et al., 2006
Considering different values of the percentage time active highlights the
importance of model selection once a technique has been selected. With
fH, the one-model approach provided a better estimate than
the two-model approach as the percentage time active increased over the
50% level. However, with PDBAxz, whatever the amount
of activity, the two-model approach provided an estimate with less error than
the one-model approach. As the percentage time active increased over the 50%
level, the s.e.e. and CV of two-model PDBAxz estimates
were very similar to the one-model fH estimates and better
than the two-model fH estimates. This shows once again
that, when animals are fairly active, using body acceleration can provide
estimates of
O2
with a similar accuracy to those made using heart rate
(Halsey et al., 2008b
;
Halsey et al., 2008c
).
Combining heart rate and accelerometry
The simulation shows the importance of selecting the best model available,
whichever proxy is being used. Crucially, this entails knowing at any
particular time what behaviour(s) the subject animals were undertaking. With
the heart-rate method, this can prove to be difficult unless the data logger
recording heart rate can store additional pertinent information and/or
observations are made of behaviour. Thus, in many applications using the
heart-rate method, a one-model approach has been recommended even if it means
a decrease in the precision of predictions (e.g.
Green et al., 2008
). However,
acceleration data can provide detailed information on behaviour as well as
being a proxy for energy expenditure. Different active behaviours produce
different characteristic traces in raw acceleration data
(Wilson et al., 2006
)
(Fig. 1), and thus analyses of
these data can reveal a minute-by-minute record of the behaviours undertaken
by animals (Gómez Laich et al.,
2008
; Tsuda et al.,
2006
; Yoda et al.,
1999
). Analyses thus far have tended to focus on active behaviours
(Gómez Laich et al.,
2008
), but future work should also concentrate on subtle changes
in acceleration traces that might be present during inactivity. For example,
during thermoregulation, animals will change their posture in order to
minimise or maximise heat loss (Dawson and
Whittow, 2000
), and this change in position should be detectable
by multi-axis accelerometers (Wilson et
al., 2008
; Yoda et al.,
1999
). Indeed, this is likely to explain the apparent change in
the baseline of the raw accelerometry traces while sleeping inside and outside
the TNZ (Fig. 1). A more crude
method would be simply to use dynamic body acceleration measures rather than
more-complex acceleration analyses. For instance, in the present study, a
threshold value of PDBAxz of 0.08 g could
be used to delineate activity and inactivity reliably
(Fig. 8), thus allowing
selection of the appropriate equation from a two-model approach. For example,
both Green and colleagues (Green et al.,
2008
) and Froget and colleagues
(Froget et al., 2002
) report
that the relationship between fH and
O2 during
thermoregulation or inactivity was different to that obtained during exercise
in penguins. However, in their studies, these authors did not suggest how this
difference in behaviour might be detected, because, as in the present study
(Fig. 9), the ranges of
fH in these behaviours overlapped substantially. Where
Green and colleagues (Green et al.,
2008
) recommended a one-model approach as a solution to this,
accelerometry could easily have detected whether birds were active or
inactive.
The inevitable conclusion is that the best way to estimate energy
expenditure accurately from free-ranging animals would be to record both heart
rate and three-axis acceleration simultaneously in the same animal. This
approach has been used in human studies
(Brage et al., 2005
) –
but with only one-axis acceleration. The combination of heart rate and body
acceleration would facilitate the choice of both technique and model to match
any circumstance encountered by study animals. Future studies should include
both fH and PDBAxz in predictive
models of
O2.
LIST OF ABBREVIATIONS
O2
| Acknowledgments |
|---|
| Footnotes |
|---|
Present address: School of Zoology, University of Tasmania, Hobart,
Tasmania 7001, Australia | References |
|---|
|
|
|---|
Blaxter, K. (1989). Energy Metabolism in Animals and Man. Cambridge: Cambridge University Press.
Blix, A. S., Stromme, S. B. and Ursin, H. (1974). Additional heart rate – an indicator of psychological activation. Aerosp. Med. 45,1219 -1222.[Medline]
Brage, S., Brage, N., Franks, P. W., Ekelund, U. and Wareham, N. J. (2005). Reliability and validity of the combined heart rate and movement sensor Actiheart. Eur. J. Clin. Nutr. 59,561 -570.[CrossRef][Medline]
Butler, P., Green, J., Boyd, I. and Speakman, J. (2004). Measuring metabolic rate in the field: the pros and cons of the doubly-labelled water and heart rate methods. Funct. Ecol. 18,168 -183.[CrossRef]
Cavagna, G. A., Saibene, F. and Margaria, R.
(1963). External work in walking. J. Appl.
Physiol. 18,1
-9.
Chappell, M. A., Bachman, G. C. and Hammond, K. A. (1997). The heat increment of feeding in house wren chicks: magnitude, duration and substitution for thermostatic costs. J. Comp. Physiol. B 167,313 -318.[CrossRef]
Dawson, W. R. and Whittow, G. C. (2000). Regulation of body temperature. In Sturkie's Avian Physiology (ed. G. C. Whittow), pp.343 -390. San Diego: Academic Press.
Fahlman, A., Wilson, R., Svärd, C., Rosen, D. A. S. and Trites, A. W. (2008). Activity and diving metabolism correlate in Steller sea lion Eumetopias jubatus. Aquat. Biol. 2,75 -84.[CrossRef]
Falk, K., Benvenuti, S., Dall'Antonia, L., Kampp, K. and Ribolini, A. (2000). Time allocation and foraging behaviour of chick-rearing Brünnichs Guillemots Uria lomvia in high-arctic Greenland. Ibis 142,82 -92.[CrossRef]
Frappell, P. B., Blevin, H. A. and Baudinette, R. V. (1989). Understanding respirometry chambers: what goes in must come out. J. Theor. Biol. 138,479 -494.[CrossRef][Medline]
Frappell, P. B., Lanthier, C., Baudinette, R. V. and Mortola, J. P. (1992). Metabolism and ventilation in acute hypoxia: a comparative analysis in a small mammalian species. Am. J. Physiol. 262,R1040 -R1046.[Medline]
Froget, G., Handrich, Y., Le Maho, Y., Rouanet, J.-L., Woakes,
A. J. and Butler, P. J. (2002). The heart rate/oxygen
consumption relationship during cold exposure of the king penguin: a
comparison with that during exercise. J. Exp. Biol.
205,2511
-2517.
Gómez Laich, A., Wilson, R. P., Quintana, F. and Shepard, E. L. C. (2008). Identification of imperial cormorant Phalacrocorax atriceps behaviour using accelerometers. Endanger. Species Res. (in press).
Green, J. A., Butler, P. J., Woakes, A. J., Boyd, I. L. and Holder, R. L. (2001). Heart rate and rate of oxygen consumption of exercising macaroni penguins. J. Exp. Biol. 204,673 -684.[Abstract]
Green, J. A., Frappell, P. B., Clark, T. D. and Butler, P. J. (2006). Physiological response to feeding in little penguins. Physiol. Biochem. Zool. 79,1088 -1097.[CrossRef][Medline]
Green, J. A., Boyd, I. L., Woakes, A. J., Green, C. J. and Butler, P. J. (2007). Feeding, fasting and foraging efficiency during chick-rearing in macaroni penguins. Mar. Ecol. Prog. Ser. 346,299 -312.[CrossRef]
Green, J. A., Frappell, P. B., Clark, T. D. and Butler, P. J. (2008). Predicting rate of oxygen consumption from heart rate while little penguins work, rest and play. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 150,222 -230.[CrossRef][Medline]
Grémillet, D. J. H., Schmid, D. and Culik, B. M. (1995). Energy requirements of breeding great cormorants Phalacrocorax carbo sinensis. Mar. Ecol. Prog. Ser. 121, 1-9.[CrossRef]
Halsey, L. G., Green, J. A., Wilson, R. P. and Frappell, P. B. (2008a). Accelerometry to estimate energy expenditure during activity: best practice with data loggers. Physiol. Biochem. Zool. (in press).
Halsey, L. G., Shepard, E. L. C., Hulston, C. J., Venables, M. C., White, C. R., Jeukendrup, A. E. and Wilson, R. P. (2008b). Acceleration versus heart rate for estimating energy expenditure and speed during locomotion in animals: Tests with an easy model species, Homo sapiens. Zoology 111,231 -241.[CrossRef][Medline]
Halsey, L. G., Shepard, E. L. C., Quintana, F., Gomez Laich, A., Green, J. A. and Wilson, R. P. (2008c). The relationship between oxygen consumption and body acceleration in a range of species. Comp. Biochem. Physiol. A Mol. Integr. Physiol. (in press).
Kaseloo, P. A. and Lovvorn, J. R. (2003). Heat increment of feeding and thermal substition in mallard ducks feeding voluntarily on grain. J. Comp. Physiol. B 173,207 -213.[Medline]
May, E. L. (2003). Application of piezoelectric sensor for measuring shivering in small marsupial. J. Therm. Biol. 28,469 -475.[CrossRef]
McKechnie, A. E., Freckleton, R. P. and Jetz, W.
(2006). Phenotypic plasticity in the scaling of avian basal
metabolic rate. Proc. Biol. Sci.
273,931
-937.
McNamara, J. M. and Houston, A. I. (1996). State-dependent life histories. Nature 380,215 -221.[CrossRef][Medline]
Meltzer, A. (1983). Thermoneutral zone and resting metabolic rate of broilers. Br. Poult. Sci. 24,471 -476.[CrossRef][Medline]
Nagy, K. A. (2005). Field metabolic rate and
body size. J. Exp. Biol.
208,1621
-1625.
Pépin, D., Renaud, P.-C., Dumont, B. and Decuq, F. (2006). Time budgets and 24-h temporal rest-activity patterns of captive red deer hinds. Appl. Anim. Behav. Sci. 101,339 -354.[CrossRef]
Phillips, R. A., Xavier, J. C. and Croxall, J. P. (2003). Effects of satellite transmitters on albatrosses and petrels. Auk 120,1082 -1090.[CrossRef]
Sedinger, J. S., White, R. G. and Hauer, W. E. (1992). Heat increment of feeding and partitioning of dietary energy in the yearling black brant. Can. J. Zool. 70,1047 -1051.[CrossRef]
Shannon, G., Page, B. R., MAckey, R. L., Duffy, K. J. and Slowtow, R. (2008). Activity budgets and sexual segregation in African elephants (Loxodonta africana). J. Mammal. 89,467 -476.[CrossRef]
Shepard, E. L. C., Wilson, R. P., Quintana, F., Gómez Laich, A., Liebsch, N., Albareda, D. A., Halsey, L. G., Gleiss, A., Morgan, D. T., Myers, A. E. et al. (2008). Identification of animal movement patters using tri-axial accelerometry. Endanger. Species Res. (in press).
Stephens, D. W. and Krebs, J. R. (1986). Foraging Theory. Princeton, NJ: Princeton University Press.
Toyomizu, M., Ueda, M., Sato, S., Seki, Y., Sato, K. and Akiba, Y. (2002). Cold-induced mitochondrial uncoupling and expression of chicken UCP and ANT mRNA in chicken skeletal muscle. FEBS Lett. 419,313 -318.
Tsuda, Y., Kawabe, R., Tanaka, H., Mitsunaga, Y., Hiraishi, T., Yamamoto, K. and Nashimoto, K. (2006). Monitoring the spawning behaviour of chum salmon with an acceleration data logger. Ecol. Freshw. Fish 15,264 -274.[CrossRef]
Ward, S., Bishop, C. M., Woakes, A. J. and Butler, P. J.
(2002). Heart rate and the rate of oxygen consumption of flying
and walking barnacle geese (Branta leucopsis) and bar-headed geese
(Anser indicus). J. Exp. Biol.
205,3347
-3356.
Weibel, E. R. and Hoppeler, H. (2005).
Exercise-induced maximal metabolic rate scales with muscle aerobic capacity.
J. Exp. Biol. 208,1635
-1644.
Weimerskirch, H., Martin, J., Clerquin, Y., Alexandre, P. and Jiraskova, S. (2001). Energy saving in flight formation. Nature 413,697 -698.[CrossRef][Medline]
Wilson, R. P. and Wilson, M.-P. (1989). Tape: a package attachment technique for penguins. Wildl. Soc. Bull. 17,77 -79.
Wilson, R. P., White, C. R., Quintana, F., Halsey, L. G., Liebsch, N., Martin, G. R. and Butler, P. J. (2006). Moving towards acceleration for estimates of activity-specific metabolic rate in free-living animals: the case of the cormorant. J. Anim. Ecol. 75,1081 -1090.[CrossRef][Medline]
Wilson, R. P., Shepard, E. L. C. and Liebsch, N. (2008). Prying into the intimate details of animal lives: use of a daily diary on animals. Endanger. Species Res. 4, 123-137.
Winship, A. J., Trites, A. W. and Rosen, D. A. S. (2002). A bioenergetic model for estimating the food requirements of Steller sea lions Eumetopias jubatus in Alaska USA. Mar. Ecol. Prog. Ser. 229,291 -312.[CrossRef]
Yoda, K., Sato, K., Niizuma, Y., Kurita, M., Bost, C.-A., Le Maho, Y. and Naito, Y. (1999). Precise monitoring of porpoising behaviour of Adélie penguins determined using acceleration data loggers. J. Exp. Biol. 202,3121 -3126.[Abstract]
![]()
CiteULike
Complore
Connotea
Del.icio.us
Digg
Reddit
Technorati
Twitter What's this?
This article has been cited by other articles:
![]() |
E. L. C. Shepard, R. P. Wilson, F. Quintana, A. Gomez Laich, and D. W. Forman Pushed for time or saving on fuel: fine-scale energy budgets shed light on currencies in a diving bird Proc R Soc B, September 7, 2009; 276(1670): 3149 - 3155. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Knight ENERGY EXPENDITURE CALCULATED FROM ACCELERATION J. Exp. Biol., February 15, 2009; 212(4): ii - ii. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||