The analysis of hydrogen stable isotopes (∂D) is a potentially powerful tool for studying animal ecology. Unlike other stable isotopes used in ecological research, however, we are less familiar with the physiological processes that influence the incorporation of hydrogen isotopes from dietary resources to animal tissues. Here we present the results of a controlled feeding experiment utilizing Japanese quail (Cortunix japonica) that was designed to: (1) estimate the relative contributions of diet to the ∂D signatures of blood plasma, red blood cells, intestine, liver, muscle and feathers; (2) investigate possible differences among these same tissues in diet to tissue discrimination; and (3) explore the differences in incorporation dynamics between deuterium (2H) and a well-studied isotope, 13C, for blood plasma solids and red blood cells. Tissues differed in both the relative contribution of diet to tissue ∂D and diet to tissue discrimination. The average residence time of both hydrogen and carbon was significantly lower in plasma solids than in red blood cells. The average residence time of hydrogen was significantly lower than that of carbon in plasma solids, but not in red blood cells. Although the average residence times of hydrogen and carbon were positively correlated, the correlation was weak. Hence the incorporation of carbon seems to be a poor predictor of that of hydrogen.

A thorough understanding of the processes that affect the incorporation of stable isotopes from resources to consumer tissues is fundamental to the reliable use of isotope analyses in animal ecology. This understanding is important for two reasons. First, because the rates at which consumer tissues incorporate the isotopic composition of resources differ (Phillips and Eldridge, 2005), knowledge of the dynamics of isotopic incorporation allows researchers to infer the temporal scope of resources used by an individual and to establish the contributions of intra- and inter-individual variation to the range of resources used by a population (Dalerum and Angerbjorn, 2005; Newsome et al., 2009). Second, knowledge of isotopic incorporation makes it possible to estimate the relative contributions of different resources to consumer tissue isotopic values, and to estimate discrimination factors between the isotopic composition of resources and consumer tissues (Wolf et al., 2011; Newsome et al., 2010) or between trophic levels (Post, 2002; Birchall et al., 2005). This need for a thorough understanding of the processes that influence the transfer of stable isotopes from the environment and resources to consumer tissues has motivated researchers to perform experiments and develop mathematical models to describe changes in the isotopic values of consumer tissues after a shift in dietary resources (Phillips and Gregg, 2001; Jardine et al., 2004; Carleton and Martínez del Rio, 2005; Martínez del Rio and Anderson-Sprecher, 2008). These investigations also examined the factors that determine variation in both incorporation rates, such as growth rate, temperature and protein turnover (Jardine et al., 2004; Carleton and Martínez del Rio, 2005; Carleton et al., 2008), and the discrimination between the isotopic values of resources and consumers (Caut et al., 2009).

When compared with carbon and nitrogen stable isotopes, the use of hydrogen stable isotopes in animal ecology is relatively nascent. Consequently, it is not surprising that our understanding of the processes that affect the incorporation of stable isotopes from resources to consumer tissues is not as complete for hydrogen as it is for carbon and nitrogen (reviewed by Martínez del Rio et al., 2009). This is not, however, to say that we know nothing about these mechanisms. Several studies have examined the relative contributions of resources to the hydrogen stable isotopes (∂2H) of consumer tissues (Estep and Dabrowski, 1980; Hobson et al., 1999; Ehleringer et al., 2008; Wang et al., 2009; Solomon et al., 2009; Nielson and Bowen, 2010) and the isotopic discrimination between resources and consumer tissues (Wolf et al., 2009). In addition, Podlesack et al. (Podlesack et al., 2008) conducted a study on the turnover of hydrogen isotopes in the body water, hair and tooth enamel of pack rats (Neotoma spp.). Although these studies have served to advance our knowledge on the relative contributions of resources to the ∂2H of consumer tissues and the incorporation rates and diet to tissue discrimination of 2H within individual tissues, we still lack information on how these processes may differ among tissues within an individual, among individuals within a species and between species. Without this information, we cannot yet confidently interpret the results of ∂2H analysis in a general comparative fashion. We performed a diet switch experiment using Japanese quail (Cortunix japonica) to examine the contributions of diet to the ∂2H values of quail body water, blood plasma, red blood cells, liver, intestine, muscle and feathers and estimated discrimination values for these tissues. We chose these tissues based on their prevalence in the methodological and mechanistic literature for carbon and nitrogen stable isotopes (Carelton et al., 2008; Caut et al., 2009). Previous attempts to determine the relative contributions of diet and drinking water to tissue ∂2H were calculated using a variety of mixing models that have yielded a wide range of results (Hobson et al., 1999; Ehleringer et al., 2008; Wang et al., 2009; Solomon et al., 2009; Nielson and Bowen, 2010). In an attempt to combat this complexity, Wolf et al. (Wolf et al., 2011) adopted a simple phenomenological strategy that, rather than attempting to determine the exact contributions of drinking water to tissue ∂2H, merely sought to characterize the relationship between drinking water ∂2H and the tissue ∂2H of house sparrows. Here we adopted this same strategy for Japanese quail.

A second function of our experiment was to compare incorporation rates for 2H and a well-studied isotope, 13C, in blood plasma and red blood cells. A thorough understanding of the differences in the dynamics of incorporation between these two isotopes will determine whether reliable comparisons can be made between the two isotopes. There are distinct differences in several aspects of the processes that may affect incorporation rates between hydrogen and carbon stable isotopes. For example, the carbon in the organic molecules that comprise animal tissues is entirely derived from diet (DeNiro and Epstein, 1978) whereas the hydrogen in these same molecules can come from several sources (Estep and Dabrowski, 1980; Hobson et al., 1999; Wang et al., 2009; Solomon et al., 2009; Nielson and Bowen, 2010; Wolf et al., 2011). Although the rate of hydrogen incorporation is not likely to be directly caused by that of carbon, it is reasonable to hypothesize that the processes that determine the magnitude of incorporation of carbon (e.g. protein turnover) will also have an effect on the incorporation of hydrogen.

Quail maintenance, experimental design and sample collection

Fifteen Japanese quail (Cortunix japonica Temminck and Schlegel 1849) [nine males and six females, initial mean body mass (Mb)=129 g, final mean Mb=121 g; paired t-test, t12=0.69, P=0.25] were kept on an initial diet of wheat (Triticum aestivum, ∂2H=–85‰, ∂13C=–26‰) and tap water (∂2H=–133‰) until their red blood cells and blood plasma achieved an asymptotic value that reflected that of the food and water (90 days). For the entire course of the experiment, the birds also received a vitamin and mineral supplement (Nekton S, Nekton GmbH, Pforzheim, Germany). Following this 90-day equilibration period, two female birds and three male birds were randomly assigned to each of three corn (Zea maize, ∂13C=–11‰) and distilled water (∂2H=–133‰) diets. The deuterium isotopic composition of corn was varied in each treatment. The levels of ∂2H in the corn diets were determined by three collections of corn harvested in three different geographical locations of differing latitude (Florida, USA: ∂2H=–22.3‰, Missouri, USA: ∂2H=–44.9‰ and Ontario, Canada: ∂2H=–54.0‰). Following the diet switch, blood samples were taken from each bird on days 0, 2, 4, 8, 16, 32, 64 and 90. The blood samples were centrifuged to separate blood cells from plasma. After 90 days, we plucked several primary feathers from each bird and allowed them to regrow. The newly grown feathers were sampled to examine the relationship of diet ∂2H and tissue ∂2H in feathers, a popular tissue in isotope-based avian movement studies. Following the feather sampling, the quail were euthanized and intestine, liver and muscle samples were collected to determine the relationship between dietary isotopic values and the isotopic values of a variety of tissues. We chose these tissues based on their prevalence in the isotope literature (Carleton et al., 2008; Caut et al., 2009).

Sample preparation and stable isotope analysis

All tissue samples except feathers were dried to a constant mass in an oven at 50°C and ground into a fine powder using a mortar and pestle. Feathers were rinsed with distilled water and air-dried. Using surgical scissors, we removed the feather barbs from the rachis and minced and homogenized the barbs. We did not include the rachis in the isotopic analysis. The powdered samples were then soaked in three 24-h treatments of petroleum ether to remove lipids. Following lipid removal, each powdered tissue sample was divided into two equal aliquots and a randomly selected sub-sample of each aliquot was loaded into a silver capsule (mean aliquot mass=0.14±0.03 mg, Costech Analytical Technologies Inc., Valencia, CA, USA) for deuterium analysis. To control for exchangeable hydrogen, encapsulated samples for deuterium analysis were equilibrated for at least 2 weeks prior to analysis in a 50°C oven under a constant stream of N2 alongside protein reference materials (chicken and turkey feathers) with known non-exchangeable ∂2H values (chicken feather: ∂2H=–92±2.0‰; turkey feather: ∂2H=–49±2.0‰). Despite the varied natures and large potential range of ∂2H values of the tissues we analyzed, we used these reference materials for all analyses because of a lack of availability of specific reference materials for other tissues. ∂2H values were determined using a Finnigan thermal conversion/elemental analyzer (TC/EA) (Thermo Fisher Scientific Inc., Waltham, MA, USA) coupled to a Thermo-Finnigan Delta Plus XL mass spectrometer (Thermo Fisher Scientific) at the Carnegie Institution of Washington (Washington, DC, USA). Samples were introduced into the TC/EA using a Costech Zero Blank autosampler (Costech Analytical Technologies). While in the autosampler, possible hydrogen exchange between ambient water vapor and samples was controlled for by flooding the sealed autosampler carousel with N2. In addition, all samples were equilibrated and analyzed using the same methods during the same time period. Consequently, we assume that the ∂2H composition of ambient water vapor was consistent throughout the analysis period and that any hydrogen exchange that occurred between ambient water vapor and tissues samples during the analysis process was also consistent. Analysis of two organic (keratin) and two inorganic (pump oil and mineral oil) reference materials showed that the precision of ∂2H analysis was ≤4‰. ∂13C values were determined using a Costech 4010 elemental analyzer (Costech Analytical Technologies) interfaced with a Thermo Finnigan DeltaPLUS XP Continuous Flow Stable Isotope Ratio Mass Spectrometer (Thermo Fisher Scientific) at the University of Wyoming Stable Isotope Facility (Laramie, WY, USA). Analysis of beef blood and acetanilide reference materials showed that the precision of ∂13C analysis was ≤0.1‰.

Statistical analysis

Ehleringer et al. (Ehleringer et al., 2008), Bowen et al. (Bowen et al., 2009) and Wolf et al. (Wolf et al., 2011) have used the slope of the linear relationship between the isotopic value of tissues and that of drinking water as a first-order approximation of the percent contribution of drinking water to tissue isotopic value. In this study, we used standard least-squares regression to characterize the relationship between the ∂2H of diet and that of individual tissues. Differences in the isotopic values among tissues were analyzed using repeated-measures (RM) ANOVA. RM ANOVAs were performed as univariate analyses using an RM MANOVA function. For each RM ANOVA, we completed a Mauchly’s test for sphericity. If this test was statistically significant, we used the Greenhouse–Geisser-corrected values for RM ANOVA. In cases where tissue values were missing, the entire bird was removed from the analysis. Post hoc means comparisons among means for all RM ANOVAs were carried out using Tukey’s honestly significant difference (HSD) tests. Differences in tissue ∂2H values between the sexes were determined using two tailed t-tests.

To determine the dynamics of incorporation for hydrogen and carbon isotopes, we used the models described by Martínez del Rio and Anderson-Sprecher (Martínez del Rio and Anderson-Sprecher, 2008) to characterize the change in the isotopic composition of a tissue as a function of time (t):
formula
where X is the isotope in question (in our case, 2H or 13C), ∂X and ∂X0 are the initial and asymptotic values of the tissue, respectively, pi is the functional size of each pool in the model (in this case, pi=1 and τi is the average residence time of the isotope in the tissue in question. Following Carleton et al. (Carleton et al., 2008), we used small-sample Akaike’s information theoretic criterion (AICc) to determine whether the data were best described by one-compartment or two-compartment models. We chose the model with the lowest AICc value (Burnham and Anderson, 2002). The chosen model permitted estimation of the isotopic incorporation rates (in days), described as average residence time (τ, the reciprocal of the fractional incorporation rate), for carbon and hydrogen in blood plasma and red blood cells. Martínez del Rio and Anderson-Sprecher (Martínez del Rio and Anderson-Sprecher, 2008) provide the statistical details of these estimates. Differences in average residence times among tissues and isotopes were determined using paired t-tests. Two tailed t-tests were used to examine differences in average residence times between the sexes. All statistical analyses were performed using JMP 8.0 (SAS Institute, Cary, NC, USA) and had assumed alpha values of 0.05.
Fig. 1.

Relationships between the ∂2H of diet and that of seven tissue types in Japanese quail: liver, muscle, intestine, blood plasma, red blood cells and feathers. Regression lines are provided for significant relationships. The ∂2H values of muscle, intestine, blood plasma and red blood cells were significantly correlated with those of diet.

Fig. 1.

Relationships between the ∂2H of diet and that of seven tissue types in Japanese quail: liver, muscle, intestine, blood plasma, red blood cells and feathers. Regression lines are provided for significant relationships. The ∂2H values of muscle, intestine, blood plasma and red blood cells were significantly correlated with those of diet.

The ∂2H values of blood plasma, red blood cells, intestine and muscle showed significant linear correlations with the ∂2H values of the corn diets (blood plasma: t10=5.37, P=0.0003; red blood cells: t11=3.16, P=0.009; intestine: t11=3.75, P=0.003; muscle: t11=2.39, P=0.03). However, the ∂2H values of liver and feathers did not show significant linear relationships with those of diet (liver: t11=1.55, P=0.15; feathers: t6=0.34, P=0.74; Fig. 1). Two birds (one from treatment group 1 and one from treatment group 3) died before the end of the experimental course. Consequently, we did not include any data from these birds in our analyses. In addition, only one bird from treatment group 1 and two birds from treatment group 3 re-grew their feathers after they had been plucked. As a result, we had very small samples sizes for feathers for these two groups. All birds from treatment group 2 re-grew their feathers after they had been plucked. For blood plasma and red blood cells, we used the asymptote of the incorporation curve (∂X in Eqn 1) as the final tissue ∂2H value. For blood plasma, we were able to calculate these values for 12 of the 13 birds that finished the experimental course. For red blood cells, we were able to calculate these values for all 13 of the surviving birds. Using Ehleringer et al. (Ehleringer et al., 2008), Bowen et al. (Bowen et al., 2009) and Wolf et al.’s (Wolf et al., 2011) interpretation of the linear relationship between the isotopic values of tissues and those of drinking water, the estimated percent contributions (±s.e.m.) of diet to the ∂2H values of tissues were 74±15% for plasma, 48±17% for red blood cells, 79±21% for intestine and 30±12% for muscle (Fig. 1).

We observed significant differences in ∂2H values among blood plasma red blood cells, intestine, liver, muscle and feathers (Greenhouse–Geisser-corrected RM ANOVA, F1.97,15.75=44.76, P<0.0001; Table 1, Fig. 2). In addition, there was a significant interaction between tissue type and drinking water treatment group (Greenhouse–Geisser-corrected RM ANOVA, F3.94,15.75=5.58, P=0.0055). This interaction prevents comparisons across independent treatments. There was no significant effect of sex (Greenhouse–Geisser-corrected RM ANOVA, F1.97,15.75=0.88, P=0.43) on the tissue ∂2H values.

Table 1.

Mean (±s.d.) ∂2H values of Japanese quail tissues for the three diet treatments

Mean (±s.d.) ∂2H values of Japanese quail tissues for the three diet treatments
Mean (±s.d.) ∂2H values of Japanese quail tissues for the three diet treatments

For carbon, the one-compartment model had lower AICc values for both plasma and red blood cells. ΔAICc (ΔAICc=AICc,one–AICc,two) values ranged from –80.35 to –11.08. For hydrogen, the one-compartment model was also better supported by data for both tissues. However, for hydrogen, the one-compartment model was favored because the two-compartment model was unable to converge in most cases. In those cases where the two-compartment model was able to generate results, the one-compartment model had lower AICc values for both blood plasma and red blood cells. For hydrogen, ΔAICc values ranged from –40.66 to –21.20. Consequently, for the rest of our analyses, we used the one-compartment model. Using this one-compartment model, we were able to generate average residence times for carbon for all 13 birds that finished the experimental course for red blood cells and 12 birds for blood plasma. For hydrogen, we were able to generate average residence times for 12 birds for both red blood cells and blood plasma.

There were no significant differences in average residence times among corn diet treatments for either carbon or hydrogen isotopes in either tissue (ANOVA, carbon: red blood cells, F2,10=1.19, P=0.34, plasma, F2,9=1.38, P=0.29; hydrogen: red blood cells, F2,9=0.87, P=0.45, plasma, F2,9=2.70, P=0.12). In addition, there were no significant differences in average residence times between the sexes for either carbon or hydrogen isotopes in either tissue (carbon: red blood cells, t11=1.11, P=0.29, plasma, t10=0.71, P=0.49; hydrogen: red blood cells, t10=1.15, P=0.28, plasma, t10=1.59, P=0.14). Isotopic average residence times for hydrogen were lower than those for carbon in red blood cells (mean ± s.d. difference=7.72±12.69 days, paired t-test, t11=2.10, P=0.05; Fig. 3). Average residence times of hydrogen isotopes were also lower than those of carbon isotopes in blood plasma, but this difference was not statistically significant (mean ± s.d. difference=3.30±8.90 days, paired t-test, t11=1.30, P=0.22; Fig. 3). Average residence times for both carbon and hydrogen isotopes were significantly lower in blood plasma than in red blood cells (carbon: mean ± s.d. difference=17.18±4.88 days, paired t-test, t9=11.11, P<0.001; hydrogen: mean ± s.d. difference=10.16±10.37 days, paired t-test, t8=2.94, P=0.02; Fig. 4). The average residence times of carbon and hydrogen isotopes were weakly but significantly correlated (r2=0.19, N=24, P=0.03; Fig. 4).

Fig. 2.

2H values of different tissues for individual Japanese quail. Lines connect individual tissues from the same bird and are provided to obviate the overall pattern in the data. Feathers were more enriched in 2H than other tissues. This pattern was consistent among individual birds and diet treatments.

Fig. 2.

2H values of different tissues for individual Japanese quail. Lines connect individual tissues from the same bird and are provided to obviate the overall pattern in the data. Feathers were more enriched in 2H than other tissues. This pattern was consistent among individual birds and diet treatments.

Our results demonstrated a significant correlation between ∂2Htissue and ∂2Hdiet in four of six tissues. We also found differences in the ∂2H composition of different tissues. Finally, our results also showed a weak, but statistically significant, correlation between the mean residence times of hydrogen and carbon isotopes. The mean residence time of both carbon and hydrogen isotopes differed between plasma and red blood cells. In subsequent sections, we will examine possible explanations for these results, and the implications of these results for the use of ∂2H analysis in ecological studies.

The relationship between the ∂2H of diet and the ∂2H of quail tissues

Previous estimates of the contribution of diet to the ∂2H of tissues have yielded values ranging from 68 to 86% in a variety of avian tissues (Hobson et al., 1999; Wolf et al., 2011), 51–80% in aquatic invertebrates (Wang et al., 2009) and 88% in salmonid fish muscle (Solomon et al., 2009; Nielson and Bowen, 2010). We observed significant relationships between the ∂2H value of tissues and diet for four of six tissues: blood plasma, red blood cells, intestine and muscle (Fig. 3). The estimated percent contributions of diet to the hydrogen isotopic compositions of these tissues ranged from 30 to 86%. In addition, the correlation between the hydrogen isotopic value of these tissues and those of food varied. The coefficient of determination of these relationships ranged from 0.34 to 0.79 (Fig. 3). We observed these variations among tissues and individuals in the estimated contributions of diet to tissue ∂2H values despite the controlled nature of our experiment. This demonstrates that much of the variation observed by previous studies is the result of differences in the biochemical composition of tissues and the physiology and/or ecology of individuals. These results demonstrate the need for further exploration of the physiological mechanisms that influence the incorporation of hydrogen isotopes from resource to consumer tissues. Of particular interest is the influence of the ∂2H of diet on the ∂2H of individual amino acids using compound-specific analysis (Howland et al., 2003). Such an approach would provide insight into the relative contributions of dietary resources to the ∂2H values of individual amino acids and the influence of a tissue’s amino acid composition on its bulk ∂2H value.

Fig. 3.

Isotopic average residence times (τ) (±s.d.) for hydrogen and carbon in blood plasma and red blood cells of Japanese quail. Average residence times for hydrogen were lower than those for carbon in red blood cells. Average residence times of hydrogen isotopes were also lower than those of carbon isotopes in blood plasma, but this difference was not statistically significant. For comparison, we standardized the isotopic incorporation curves by setting ∂2H=1 and ∂2H0=0.

Fig. 3.

Isotopic average residence times (τ) (±s.d.) for hydrogen and carbon in blood plasma and red blood cells of Japanese quail. Average residence times for hydrogen were lower than those for carbon in red blood cells. Average residence times of hydrogen isotopes were also lower than those of carbon isotopes in blood plasma, but this difference was not statistically significant. For comparison, we standardized the isotopic incorporation curves by setting ∂2H=1 and ∂2H0=0.

An unexpected result of our experiment was the lack of a positive linear relationship between the ∂2H values of diet and tissue for liver and feathers. For feathers, this could have been the result of our relatively small sample size. However for liver, we suspect that our study’s failure to detect this relationship (Fig. 1) might be due to the relatively narrow range of ∂2H values employed on our diets and the lack of covariation between the ∂2H values of our diets and that of drinking water. The ∂2H values of the corn diets employed in our experiment only represent one-third of the variation seen in the ∂2H values of surface water (Waterisotopes.org, http://wateriso.eas.purdue.edu/waterisotopes/), and presumably producers, over the North American continent (<33‰, from –54 to –22‰). Although biologically relevant, this range may not be sufficient for establishing the relationship between the ∂2H value of diet and that of tissues, especially because the ∂2H value of drinking water was held constant across diet treatments in our experiments.

To assess the effect of both increasing the range of ∂2H values of diet treatments and the covariation between the ∂2H value of diet and drinking water, we performed a set of simulations. We varied the range of ∂2H values of diet from 30‰ (–52 to –22‰) to 60‰ (–82 to –22‰). We assumed that either the ∂2H values of drinking water were identical to those of diet, or that they were constant and independent of the value of diet (∂D=–135‰). Based on the estimated percent contribution of diet to intestinal tissue ∂2H observed in our quail, we assumed that food contributed 79% of the hydrogen in tissues. We chose to use the value for the estimated percent contribution of diet to intestinal tissue ∂2H because of its similarity to values found by previous studies (Hobson et al., 1999; Wang et al., 2009; Solomon et al., 2009; Nielson and Bowen, 2010; Wolf et al., 2011). To generate error, we assumed that the incorporation fractions of both diet and drinking water to the hydrogen content in tissues had a normally distributed error with a mean of 21% and a standard deviation of 29%. This mean value was based on the standard error of the slope of the relationship between the ∂2H values of diet and intestinal tissue found in our experiments. The standard deviation was calculated from the contributions of diet ∂2H to bird tissue ∂2H observed in our data and two additional sources (Wolf et al., 2011) (N.W., S.D.N., M.L.F. and C.M.d.R., submitted). Based on our data for intestinal tissue, we assumed that the discrimination factor from diet to tissues was normally distributed with a mean of –92‰ and a standard deviation of –16‰. The discrimination factor from drinking water to tissues was also normally distributed with a mean of –19‰ and a standard deviation of –17‰ (value for Japanese quail intestinal tissue from N.W., S.D.N., M.L.F. and C.M.d.R., submitted). To reflect the conditions of our experiment, we generated 15 points at regular intervals along the simulated range of ∂2Hdiet values. Accordingly, the model relating the ∂2H value of diet with that of tissues was:
formula
where εp is the error associated with the fractional contributions of diet and drinking water and ε is the error associated with diet to tissue discrimination. In these simulations, there was no significant relationship between the range of ∂2H values in diet treatments and the r2 value of the relationship between the ∂2H values of diet and tissue (r2=0.02, N=40, P=0.36, y=0.037+0.001x (Fig. 5). In addition, the r2 value for the simulations when the ∂2H value of water covaried with that of diet was not significantly different from when the ∂2H value of water was held constant (ANCOVA, F3,76=1.18, P>0.05; Fig. 5). Finally, most of the r2 values of the simulated relationships (r2=0.27) were lower than the P-values of the linear regression. These results demonstrate that the small number of tissues for which we found positive linear relationships between the ∂2H values of diet and tissue was not likely caused by the relatively narrow range of ∂2H values employed on our diets, nor was it caused by the lack of covariation between the ∂2H values of our diets and that of drinking water. Instead, the lack of significant relationships for several tissues was caused by the large intra-individual variation in both incorporation rates and discrimination factors exhibited by hydrogen.
Fig. 4.

Relationship between average residence times (τ) (±s.d.) of carbon and hydrogen isotopes in blood plasma and red blood cells of Japanese quail. Average residence times for both carbon and hydrogen isotopes were significantly lower in blood plasma than in red blood cells. The average residence times of carbon and hydrogen isotopes were weakly, but significantly, correlated. The dashed line represents a 1:1 line.

Fig. 4.

Relationship between average residence times (τ) (±s.d.) of carbon and hydrogen isotopes in blood plasma and red blood cells of Japanese quail. Average residence times for both carbon and hydrogen isotopes were significantly lower in blood plasma than in red blood cells. The average residence times of carbon and hydrogen isotopes were weakly, but significantly, correlated. The dashed line represents a 1:1 line.

Fig. 5.

Simulation to assess the effect of both increasing the range of ∂2H values of diet treatments and the covariation between the ∂2H value of diet and drinking water. (A) Drinking water ∂2H values covary with diet ∂2H values. (B) Drinking water ∂2H values are constant (–135‰). There was no significant relationship between the range of ∂2H values among diet treatments and the r2 value of the relationship between the ∂2H values of diet and tissue. In addition, the r2 values for the simulations when ∂2H values were held constant were not significantly different from those for the simulations when ∂2H values of water covaried with those of diet.

Fig. 5.

Simulation to assess the effect of both increasing the range of ∂2H values of diet treatments and the covariation between the ∂2H value of diet and drinking water. (A) Drinking water ∂2H values covary with diet ∂2H values. (B) Drinking water ∂2H values are constant (–135‰). There was no significant relationship between the range of ∂2H values among diet treatments and the r2 value of the relationship between the ∂2H values of diet and tissue. In addition, the r2 values for the simulations when ∂2H values were held constant were not significantly different from those for the simulations when ∂2H values of water covaried with those of diet.

This result has implications for the use ∂2H analysis for studying animal movements and migratory strategies. The most common use of hydrogen isotopes in the ecological literature compares the ∂2H value of an animal tissue (typically a bird feather) with the ∂2H value of precipitation in an ‘isoscape’ (West et al., 2006) to determine where the tissue was grown (Hobson, 2008; Wunder, 2010). Hobson et al. (Hobson et al., 1999) noted that large isotopic differences between drinking water and food at a single locality could lead to variation in the relationship between the isotopic content of an animal’s tissues and that of the precipitation of the area in which this animal lives. The results of our simulations showed that the variation in tissue ∂2H values was similar for scenarios when drinking water ∂D values covaried with diet ∂2H values and for scenarios when drinking water ∂D values were held constant. Our simulation demonstrates that for birds (e.g. Japanese quail) with high and variable drinking water rates, the reliability of ∂2H analysis for studying animal movements might not decrease when the ∂2H values of diet and drinking water differ.

Variation in ∂2H among tissues

In our experiment, feathers were more enriched in 2H than other tissues. In addition, the mean ∂2H values of all other tissues were significantly different (Table 1). This pattern of variation was consistent among individual birds from the same treatment (Fig. 2). In their drinking water switch experiment using house sparrows, Wolf et al. (Wolf et al., 2011) also found that feathers were significantly more enriched in 2H than other tissues. Although we see similar patterns in the discrimination values of hydrogen isotopes between house sparrows and Japanese quail, the variations in tissue-specific discrimination values seen among tissues in both species demonstrates the need to develop a mechanistic understanding of the processes that influence variations in hydrogen isotopic discrimination among tissues. Although we suspect that these processes might be similar to those that influence carbon and nitrogen [e.g. differences in nutrient routing, protein turnover and the macromolecular compositions of tissues (Martínez del Rio et al., 2009)], future ∂2H analysis of individual amino acids is necessary to determine both the general patterns in hydrogen isotopic discrimination among tissues and the processes that determine these patterns.

Similarities in the residence times of 13C and 2H

Similarities in the average residence times of carbon and hydrogen isotopes would allow for valuable comparisons between the two isotopes. Our results demonstrate a weak, yet statistically significant, correlation in the average residence times of carbon and hydrogen. However, the regression of the relationship of τ13C and τD only explained approximately 20% of the variation in deuterium retention (Fig. 4). In addition, the slope of the relationship of the residence times of the two isotopes was significantly lower than 1 (slope ± s.e.m.=0.51±0.23, t=2.08, P=0.04), and the intercept was not significantly different from 0 (intercept ± s.e.m.=2.27±4.36, t=0.62, P=0.53).

Differences in the average residence times of 13C and 2H

For red blood cells, the average residence times of hydrogen isotopes were significantly shorter than those of carbon isotopes. Average blood plasma residence times of hydrogen isotopes were also shorter than those of carbon isotopes, but this difference was not statistically significant (Fig. 3). The lack of a significant difference in τ2H and τ13C in blood plasma was likely the result of the large amount of variation in the retention time of hydrogen in plasma observed among individuals (Fig. 3). The coefficient of variation (CV) of residence time of hydrogen (CV=1.07) was approximately twice as high as that for carbon (CV=0.54) in blood plasma. This variation, and the shorter mean residence times observed for hydrogen isotopes, might be attributable to the greater number of potential sources for hydrogen in consumer tissues than for carbon in consumer tissues. The carbon in the molecules that comprise animal tissues is entirely derived from diet (DeNiro and Epstein, 1978), but hydrogen in animal tissues can come from either diet or body water (Hobson et al., 1999). In turn, the hydrogen in body water can come from two potential sources: pre-formed water (from food or drinking water) or metabolically produced water (Hobson et al., 1999; Wolf et al., 2011). Because the hydrogen in consumer tissues can be derived from multiple source pools, and each of these pools can have different dynamics (Ehleringer et al., 2008), the interaction of these pools can result in variations in the incorporation rates of hydrogen from each source. For example, in our experiment, birds were provided with water ad libitum. Variations in drinking water rates among individual birds could lead to variations in the incorporation of hydrogen derived from the diet into the body water pool (via metabolic water production). This would cause resultant variation in the incorporation of diet-derived hydrogen into consumer tissues. Similarly, in animals with low water intake requirements, the incorporation of diet-derived hydrogen into body water will be faster and, consequently, mean residence times will be shorter than those found in animals with higher drinking water requirements. In addition to the rate at which hydrogen from diet is incorporated into tissues, differences in drinking water intake will influence the fraction of the hydrogen found in consumer tissues that is derived from diet sources (Wolf et al., 2011).

Differences in average residence times of 13C and 2H between red blood cells and blood plasma

The average residence times of both carbon and hydrogen isotopes were significantly shorter in blood plasma than in red blood cells. This pattern agrees with the results of previous studies on the incorporation rates of carbon isotopes among tissues (Carleton et al., 2008). Isotopic ecologists have capitalized on the significant differences in 13C average residence times among tissues to study dietary ecology (Dalerum and Angerbjorn, 2005; Martínez del Rio et al., 2009). By conducting ∂13C analysis on multiple tissues from the same individual, ecologists can infer the timing of changes in resource use by an individual animal (Klassen et al., 2010). Although we observed significant differences in the average residence times of deuterium in blood plasma and red blood cells, our results prompt us to urge against using deuterium in a similar fashion. Estimating the timing of shifts in dietary resources using this method requires accurate estimates of the differences in isotopic incorporation rates between tissues (Klassen et al., 2010). In one-compartment incorporation models, such as those we used, average retention times are the reciprocal of fractional incorporation rates (Martínez del Rio and Anderson-Sprecher, 2008). The coefficient of variation of the difference between the fractional incorporation rate of blood plasma and red blood cells of hydrogen (CV=1.25) is almost double that of carbon (CV=0.71). Large variations in the difference between the fractional rates of incorporation of two tissues can lead to high error propagation and consequent large uncertainty in estimates of diet shift (Klassen et al., 2010). However, despite the variation we observed, we do see a significant difference in the average residence times of the two tissues. Consequently, we recommend further exploration of the average residence times of hydrogen isotopes in a variety of tissues before concluding that the average residence times of hydrogen isotopes in red blood cells and blood plasma are sufficiently variable so as to negate the reliable use of these tissues as independent indicators of hydrogen isotope incorporation in ecological studies.

Conclusions

Our results demonstrate the complexity of using ∂2H analysis for studying animal ecology. With our current state of knowledge on the physiological processes that influence the incorporation rates and diet to tissue discrimination of ∂2H, we face great difficulty in achieving accurate results using hydrogen isotope analysis. The results of our study show that variations in the ∂2H values among animal tissues, both among individuals and individual tissues in the same animal, is a natural phenomenon that is likely driven by physiological and ecological causes. This variation is high even in controlled experimental situations in which individuals are provided with food and water with identical and constant ∂2H values. In addition, our results demonstrate the potential differences in the incorporation of carbon and hydrogen isotopes. Given that hydrogen can be derived from two sources (diet and pre-formed water), it is likely that the hydrogen isotope composition of wild, free-ranging animals will be inherently more variable and unpredictable than for carbon isotopes. Hydrogen stable isotope analysis is a potentially powerful tool in animal ecology. Given this potential, we think that it is important that the limiting aspects of this tool are explored thoroughly. However, until more information is available on the reliability of the method, we recommend that future studies that aim to use ∂2H analysis in the field consider the variation demonstrated by our results in the design of their studies and the interpretation of their results.

     
  • 13C

    carbon isotope (6 protons, 7 neutrons) (moles)

  •  
  • 2H

    hydrogen isotope (1 proton, 1 neutron) (moles)

  •  
  • AICc

    Akaike’s information theoretic criterion for small samples

  •  
  • CV

    coefficient of variation

  •  
  • fA

    relative fractional contribution of source A

  •  
  • fB

    relative fractional contribution of source B

  •  
  • pi

    functional size of each pool in the incorporation model

  •  
  • 13C

    [(2 13C/1 13Csample2 13C/1 13Cstandard) 2 13C/1 13standard–1] × 103 (‰, Vienna Pee Dee belemnite)

  •  
  • 2H

    [(2H/1Hsample)–2H/1Hstandard) 2H/1Hstandard–1] × 103 (‰, standard mean ocean water)

  •  
  • ∂A

    isotopic composition of source A

  •  
  • ∂B

    isotopic composition of source B

  •  
  • ∂T

    isotopic composition of a consumer tissue

  •  
  • ∂X

    asymptotic value of an isotope (X) in a tissue

  •  
  • ∂X0

    initial value of an isotope (X) in a tissue

  •  
  • ε

    error associated with diet to tissue discrimination

  •  
  • εp

    error associated with the fractional contributions of diet and drinking water

  •  
  • τ

    average residence time of the isotope for a tissue

We thank Robert Carroll for his assistance with bird care and Roxanne Bowden for her assistance with ∂2H analysis.

FUNDING

N.W. was supported by the University of Wyoming Dr George E. Menken’s Memorial Scholarship. N.W., C.M.d.R. and S.D.N. were supported by the National Science Foundation (DIOS 0848028). M.L.F. was supported by the Carnegie Institution of Washington.

Birchall
J.
,
O’Connell
T. C.
,
Heaton
T. H. E.
,
Hedges
R. E. M.
(
2005
).
Hydrogen isotope ratios in animal body protein reflect trophic level
.
J. Anim. Ecol.
74
,
877
881
.
Bowen
G. J.
,
Ehleringer
J. R.
,
Chesson
L. A.
,
Thompson
A. H.
,
Podlesak
D. W.
,
Cerling
T.E.
(
2009
).
Dietary and physiological controls on the hydrogen and oxygen isotope ratios of hair from mid-20th century indigenous populations
.
Am. J. Phys. Anthrop.
139
,
494
504
.
Burnham
K. P.
,
Anderson
D. R.
(
2002
).
Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach
.
Springer
,
New York
.
Carleton
S. A.
,
Martínez del Rio
C.
(
2005
).
The effect of cold-induced increased metabolic rate on the rate of 13C and 15N incorporation in house sparrows (Passer domesticus)
.
Oecologia
144
,
226
232
.
Carleton
S. A.
,
Kelly
L.
,
Anderson-Sprecher
R.
,
Martínez del Rio
C.
(
2008
).
Should we use one-, or multi-compartment models to describe 13C incorporation into animal tissues?
Rapid Commun. Mass Spectrom.
22
,
3008
3014
.
Caut
S.
,
Angulo
E.
,
Courchamp
F.
(
2009
).
Variation in discrimination factors (Δ15N and Δ13C): the effect of diet isotopic values and applications for diet reconstruction
.
J. Appl. Ecol.
46
,
443
453
.
Dalerum
F.
,
Angerbjörn
A.
(
2005
).
Resolving temporal variation in vertebrate diets using naturally occurring stable isotopes
.
Oecologia
144
,
647
658
.
DeNiro
M. J.
,
Epstein
S.
(
1978
).
Influence of diet on the distribution of carbon isotopes in animals
.
Geochim. Cosmochim. Acta.
42
,
495
506
.
Ehleringer
J. R.
,
Bowen
G. J.
,
Chesson
L. A.
,
West
A. G.
,
Podlesak
D. W.
,
Cerling
T. E.
(
2008
).
Hydrogen and oxygen isotope ratios in human hair are related to geography
.
Proc. Natl. Acad. Sci. USA
105
,
2788
2793
.
Estep
M. F.
,
Dabrowski
H.
(
1980
).
Tracing food webs with stable hydrogen isotopes
.
Science
209
,
1537
1538
.
Hobson
K. A.
(
2008
).
Applying isotopic methods to tracking animal movements
. In
Tracking Animal Migrations With Stable Isotopes
(ed.
Hobson
K. A.
,
Wassenaar
L. I.
), pp.
45
78
.
London
:
Academic Press
.
Hobson
K. A.
,
Atwell
L.
,
Wassenaar
L. I.
(
1999
).
Influence of drinking water and diet on stable-hydrogen isotope ratios of animal tissues
.
Proc. Natl. Acad. Sci. USA
96
,
8003
8006
.
Howland
M. R.
,
Corr
L. T.
,
Young
S. M. M.
,
Jones
V.
,
Jim
S.
,
Van Der Merwe
N. J.
,
Mitchell
A. D.
,
Evershed
R. P.
(
2003
).
Expression of the dietary isotope signal in the compound-specific ∂13C values of pig bone lipids and amino acids
.
Int. J. Osteoarcheology
13
,
54
65
.
Jardine
T. D
,
McLatchy
D. L.
,
Fairchild
W. L.
,
Cukack
R. A.
,
Brown
S. B.
(
2004
).
Rapid carbon turnover during growth of Atlantic salmon (Salmo salar) in salt water, and evidence for reduced food consumption by growth stunts
.
Hydrobiologia
527
,
63
75
.
Klaasen
M.
,
Piersma
T.
,
Korthals
H.
,
Dekinga
A.
,
Dietz
M. W.
(
2010
).
Single point isotope measurements in blood cells and plasma to estimate the time since diet switches
.
Funct. Ecol.
24
,
796
804
.
Martínez del Rio
C.
,
Anderson-Sprecher
R.
(
2008
).
Beyond the reaction-progress variable: the meaning and significance of isotopic incorporation data
.
Oecologia
156
,
765
772
.
Martínez del Rio
C.
,
Wolf
N.
,
Carleton
S. A.
,
Gannes
L. Z.
(
2009
).
Isotopic ecology ten years after a call for more laboratory experiments
.
Biol. Rev.
84
,
91
111
.
Newsome
S. D.
,
Tinker
M. T.
,
Monson
D. H.
,
Oftedal
O.
,
Ralls
K.
,
Fogel
M. L.
,
Estes
J. A.
(
2009
).
Using stable isotopes to investigate individual diet specialization in California sea otters (Enhydra lutrisnereis)
.
Ecology
90
,
961
974
.
Newsome
S. D.
,
Clementz
M. R.
,
Koch
P. L.
(
2010
).
Using stable isotope biogeochemistry to study marine mammal ecology
.
Mar. Mamm. Sci.
26
,
509
572
.
Nielson
K. E.
,
Bowen
G. J.
(
2010
).
Hydrogen and oxygen in brine shrimp chitin reflect environmental water and dietary isotopic composition
.
Geochim. Cosmochim. Acta.
74
,
1812
1822
.
Phillips
D. L.
,
Eldridge
P. M.
(
2005
).
Estimating the timing of diet shifts using stable isotopes
.
Oecologia
147
,
195
203
.
Phillips
D. L.
,
Gregg
J. W.
(
2001
).
Uncertainty in source partitioning using stable isotopes
.
Oecologia
127
,
171
179
.
Podlesak
D. W.
,
Torregrossa
A. M.
,
Ehleringer
J. R.
,
Dearing
M. D.
,
Passey
B. H.
,
Cerling
T. E.
(
2008
).
Turnover of oxygen and hydrogen isotopes in the body water, CO2, hair, and enamel of a small mammal
.
Geochim. Cosmochim. Acta
72
,
19
35
.
Post
D. M.
(
2002
).
Using stable isotopes to estimate trophic position: models, methods, and assumptions
.
Ecology
83
,
703
718
.
Solomon
C. T.
,
Cole
J. J.
,
Doucett
R. R.
,
Pace
M. L.
,
Preston
N. D.
,
Smith
L. E.
,
Weidel
B. C.
(
2009
).
The influence of environmental water on the hydrogen stable isotope ratio in aquatic consumers
.
Oecologia
161
,
313
324
.
Wang
Y. V.
,
O’Brien
D. M.
,
Jenson
J.
,
Francis
D.
,
Wooller
M. J.
(
2009
).
The influence of diet and water on stable oxygen and hydrogen isotope composition of Chironomidae (Diptera) with paleoecological implications
.
Oecologia
160
,
225
233
.
West
J. B.
,
Bowen
G. J.
,
Cerling
T. E.
,
Ehleringer
J. R.
(
2006
).
Stable isotopes as one of nature’s ecological recorders
.
Trends Ecol. Evol.
21
,
408
414
.
Wolf
N.
,
Carleton
S. A.
,
Martinez del Rio
C.
(
2009
).
Ten years of experimental animal isotopic ecology
.
Funct. Ecol.
23
,
17
26
.
Wolf
N.
,
Bowen
G. J.
,
Martínez del Rio
C.
(
2011
).
The influence of drinking water on the ∂D and ∂18O values of house sparrow plasma, blood, and feathers
.
J. Exp. Biol.
214
,
98
103
.
Wunder
M. B.
(
2010
).
Using isoscapes to model probability surfaces for determining geographic origins
. In
Isoscapes: Understanding Movement, Pattern, and Process on Earth Through Isotope Mapping
(ed.
West
J. B.
,
Bowen
G. J.
,
Dawson
T. E.
,
Tu
K. P.
), pp.
251
272
.
New York
:
Springer
.