Isotopic partitioning by small mammals in the subnivium

Isotopic partitioning by small mammals in the subnivium Ivan Calandra 1 ,Ga € elle Labonne 2,3 , Olivier Mathieu 3 , Heikki Henttonen 4 , Jean L ev^ eque 3 , Marie-Jeanne Milloux 3 ,  Elodie Renvois e 5 , Sophie Montuire 2,3 & Nicolas Navarro 2,3 1 GEGENAA –EA 3795, Universit e de Reims Champagne-Ardenne, Reims, France 2 Laboratoire PALEVO,  Ecole Pratique des Hautes  Etudes, Dijon, France 3 Biog eosciences –UMR uB/CNRS 6282, Universit e Bourgogne Franche-Comt e, Dijon, France 4 Natural Resources Institute Finland, Vantaa, Finland 5 Evo-Devo Lab, Institute of Biotechnology, University of Helsinki, Helsinki, Finland Keywords Arvicolinae, carbon isotopes, foraging behavior, Lapland, nitrogen isotopes, seasonality. Correspondence Ivan Calandra, GEGENAA - EA 3795, Universit e de Reims Champagne-Ardenne, CREA - 2 esplanade Roland Garros, 51100 Reims, France. Tel: +33(0) 326 77 36 89; Fax: +33(0) 326 77 36 20; E-mail: [email protected] Funding Information This work was supported by the Universit e de Bourgogne grant no. 2011 BQR 087 to IC and a CNRS INSU INTERRVIE grant to NN. Received: 6 March 2015; Revised: 13 July 2015; Accepted: 25 July 2015 Ecology and Evolution 2015; 5(18): 4132–4140 doi: 10.1002/ece3.1653 Abstract In the Arctic, food limitation is one of the driving factors behind small mammal population fluctuations. Active throughout the year, voles and lemmings (arvicoline rodents) are central prey in arctic food webs. Snow cover, however, makes the estimation of their winter diet challenging. We analyzed the isotopic composition of ever-growing incisors from species of voles and lemmings in northern Finland trapped in the spring and autumn. We found that resources appear to be reasonably partitioned and largely congruent with phylogeny. Our results reveal that winter resource use can be inferred from the tooth isotopic composition of rodents sampled in the spring, when trapping can be conducted, and that resources appear to be partitioned via competition under the snow. Introduction Winter plays a key role in the population dynamics of small herbivorous mammals in the Arctic (Reid and Krebs 1996; Hansen et al. 1999a; Kausrud et al. 2008), but remains a poorly understood period of their annual cycle (Duchesne et al. 2011). Better knowledge of the resource allocations and foraging strategies used by these small mammals during winter is therefore important to improve our understanding of the ecology of these species in arctic food webs (Chappell 1980; Ims et al. 2013; Soininen et al. 2015). Arvicoline rodents (voles and lemmings) are key components of arctic ecosystems, as they constitute the sole prey of secondary consumers during winter (e.g., Ims et al. 2013). Aerial vegetative parts of plants seem to be the main constituents of the vole diet, while lemmings are known to be moss specialists (Stenseth and Ims 1993). It seems that each group of arvicoline shows specific feeding preferences; arvicoline diets would thus expected to be structured according to phylogeny, as has been shown for European rodents (Butet and Delettre 2011). Flexibility and seasonality in the diets of small arctic mammals have been described to some extent (Hansson 1971a; Hansson and Larsson 1978; Batzli and Henttonen 1990; Eskelinen 2002; Saetnan et al. 2009; Soininen et al. 2009, 2013a). In winter, as small arctic mammals forage in the subnivium (the interface between soil and snow; Pauli et al. 2013; Petty et al. 2015), this task has proven difficult because direct feeding observations and trapping are impossible (but see Bilodeau et al. 2013). Most dietary studies have therefore relied on gut or feces contents of animals trapped between June and September (but see Tast 1974; Soininen et al. 2015), which represent only 4132 ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. snapshots of an individual’s summer and autumn diets. So any method that can reconstruct the winter diet from animals trapped when the snow cover is gone would yield valuable information about the ecology of small arctic mammals. The isotopic composition of an animal’s tissue reflects the isotopic composition of the food ingested, shifted through fractionation by a given amount (i.e., discrimination factor), dependent on the animal’s metabolism (reviewed by Kohn and Cerling 2002; Cerling et al. 2010; Clementz 2012). In rodents, most studies in isotope ecology have relied on soft tissues or feces (e.g., Sare et al. 2005; Soininen et al. 2014), while stable isotope analyses on teeth have been restricted to reconstructing dietary ecologies in fossil taxa (e.g., Grimes et al. 2004; Hopley et al. 2006; Gehler et al. 2012; Ga zsiorowski et al. 2014). Unlike soft tissues, teeth do not decay and, in the case of rodent teeth, are generally conserved unaltered in owl pellets and can thus be used to reconstruct communities over time scales from days to years. Teeth grow progressively and retain the isotopic signal once formed, an interesting property for the study of short-term variations, such as seasonal differences (Dalerum and Angerbj€ orn 2005). The isotopic composition of tooth tissues corresponds to the diet during the maturation period (Kohn and Cerling 2002). In the case of arvicolines, the ever-growing incisors are completely renewed in 6– 8 weeks (Klevezal et al. 1990). Thus, the whole incisor should record the average diet of an individual over a period of 6–8 weeks before death, while the older half of the tooth should represent a signal 6–8 weeks old, but which does not include the last 3–4 weeks. By analyzing this specific part, it should be possible to quantify the April diet (corresponding to the late winter diet, as the snow cover is still present) of rodents trapped in June (when the snow has gone). While this assumption seems logical, there is much uncertainty about metabolic routing and how this affects the isotopic composition of a given tissue (Dalerum and Angerbj€ orn 2005). So, in order to test the feasibility and robustness of this approach, we selected a sample of voles from Finland with well-studied seasonal ecologies. In this study, stable carbon and nitrogen isotope analyses were carried out on the lower incisors of specimens from seven species of Finnish arvicolines, belonging to three different tribes, trapped at two separate seasons, in order to test the following hypotheses: (1) the diet of a given species will be closer to that of another species from the same tribe than to that of a species from a different tribe, and (2) it will be possible to infer the winter diet by analyzing the teeth of individuals trapped in spring. Answering these questions will further our understanding of seasonal ecologies in small arctic animals. Material and Methods Five Myodes glareolus females (“IBH” specimens in Table S1), raised in laboratory conditions, were used to estimate the discrimination factors in isotopic composition between food sources and tooth tissues. From weaning onwards (about 3 weeks after birth), the voles were fed exclusively with breeding diet pellets for rats/mice (number 1314 Fortified; Altromin Spezialfutter GmbH & Co. KG, Lippe, Germany). When sacrificed, all five females were healthy, and between four and 5 months old. The preweaning diet must have been completely overwritten by the postweaning pellet diet at the time of sacrifice, as incisors are completely renewed in 6–8 weeks (Klevezal et al. 1990). Sixty-two wild arvicolines (“UB” specimens in Table S1) were trapped at two sites in Finnish Lapland. Half of the individuals were trapped in spring (June, n=30), and half in autumn (September, n=32), most of them between 2010 and 2011. The first site, Pallasj€ arvi, is a boreal taiga zone. The second site, Kilpisj€arvi, in the north-westernmost part of Finland, is characterized at higher altitudes by alpine tundra, with subarctic mountain birch forests at lower altitudes, around the biological station. At the Kilpisj€ arvi site, only the lemmings were trapped in the tundra, while the voles could only be trapped in the forest habitat, which resembles the taiga of Pallasj€ arvi. Therefore, for subsequent analyses, arvicolines were not separated by trapping locality, but by species, tribe, and season. The specimens analyzed belong to seven arvicoline species (Microtus agrestis and M. oeconomus;Myodes glareolus, M. rufocanus, and M. rutilus;Lemmus lemmus [see Cover image]; and Myopus schisticolor) within three tribes (Arvicolini, Clethrionomyini, and Lemmini, respectively). The diets of these species are seasonally variable, species-specific and well known in the Arctic (Table 1). Samples of plants probably consumed by Finnish arvicolines were collected from Lapland (mostly from Kilpisj€ arvi; Table S2), in July 2012. The food pellets fed to the laboratory animals were also included in the analyses. We also included reviews by Ben-David et al. (2001) and Drucker et al. (2010, 2012) and reviewed the literature on mosses (see Fig. 1A). The shift from one trophic level to the next is classically thought to be approximately +3& in d 15 N (e.g., Ben-David and Flaherty 2012). The heads of the wild specimens and of the laboratory voles were prepared following the protocol in Appendix S1. As teeth sometimes contain inorganic carbonate, acidification tests were performed. We found that the inorganic carbonate content of the teeth was low enough to have no effect on the carbon and nitrogen isotopic compositions (Appendix S1, Fig. S1). Both lower incisors were extracted from the mandibles. The enamel and dentine of ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 4133 I. Calandra et al. Isotopic Partitioning in Small Mammals the oldest halves of the incisors (which had formed at least 3–4 weeks before trapping; Klevezal et al. 1990) were manually ground into a fine homogeneous powder, with a mortar and pestle. The resulting powder was then dried overnight at 60°C. All plant samples were frozen (75°C) for at least 24 h, then lyophilized (<1.5 mbar, 132°C) for at least another 24 h, and then ground to a fine powder. All the plant and pellet samples were dried overnight at 60°C, the day before the preparation of the tin capsules. Three tin capsules per individual (arvicolines, plants and pellets) were analyzed for elemental carbon and nitrogen proportions and for stable carbon (d 13 C) and nitrogen (d 15 N) isotopic compositions. Isotopic compositions are expressed relative to the conventional standards, Vienna Pee Dee Belemnite (VPDB) for carbon (Craig 1953) and atmospheric N 2 for nitrogen (Mariotti 1983). To calibrate the elemental and isotopic analyses, a glutamic acid standard (L-USGS 40: C =40.8%, N=9.5%, d 13 C=26.389 0.042&,d 15 N=4.5 0.1&) was used. The precision of analysis was 0.15&(d 13 C) and 0.2&(d 15 N). The values for the three tin capsules per individual were averaged before analysis. The variation between these three samples did not exceed the precision of analysis. Interspecific and seasonal dietary variations were tested. An unsupervised model-based classification (Gaussian finite mixture model, fitted by an expectation–maximization algorithm; Fraley and Raftery 2002; Fraley et al. 2012) was run on the wild voles, to test for a grouping based on the isotopic composition of the teeth. The optimal model was chosen according to the highest Bayesian information criterion. The magnitude of the agreement between this classification and the known heterogeneities existing in the sample (taxonomy and season) was then calculated using Cohen’s Kappa statistic (Cohen 1960). A two-way ANOVA with type II sums of squares was performed on each stable isotope ratio. A square root transformation on the d 15 N values, √(d 15 N+1), was required prior to statistical analysis. The factors were tribe, species nested within tribe, and season; all interactions between these factors were also tested. As the sample size was in some cases quite small, the species factor Table 1. Dietary data from the literature, for the arvicolines studied. Tribe Species Spring/summer diet Autumn/winter diet References Arvicolini Microtus agrestis Grass/sedge (ca. 50%, up to 80%), forbs (ca. 40%), shrubs (ca. 10%); complemented by invertebrates, bark, berries, fungi Grass/sedge (ca. 70%, up to 90%), forbs (ca. 15%, up to 65%); complemented by bark, invertebrates, berries, fungi Hansson (1971a,b), Stenseth et al. (1977), Hansson and Larsson (1978), Saetnan et al. (2009), Butet and Delettre (2011) Microtus oeconomus Grass shoots (ca. 50%, up to 80%), forbs (ca. 20%, up to 65%), horsetail (up to 20%), shrubs (5–10%) Underground rhizomes and shoots of grass/sedge (up to 95%); complemented by forbs and shrubs Tast (1974), Batzli and Henttonen (1990), Soininen et al. (2009, 2013a) Clethrionomyini Myodes glareolus Forbs (40–50%), invertebrates (30–40%), grass/sedge (ca. 15%, up to 30%), berries (ca. 10%, up to 35%) Forbs (ca. 20%, up to 45%), grass/sedge (ca. 20%), lichen (ca. 20%, up to 35%), fungi (ca. 10%, up to 55%), shrubs (ca. 10%, up to 45%), berries (ca. 10%, up to 25%) Hansson (1969, 1971a), Hansson and Larsson (1978), Sulkava (1978, in Viro and Niethammer 1982), Butet and Delettre (2011) Complemented by fungi Complemented by invertebrates Myodes rufocanus Shrubs (ca. 50%, especially Vaccinium shoots), forbs (ca. 25%), grass/sedge (ca. 5%), horsetail (ca. 5%) Shrubs (ca. 60%, especially Vaccinium), grass (ca. 15%), forbs (ca. 10%); complemented by Betula bark, seeds/berries Hansson and Larsson (1978), Henttonen and Viitala (1982), Henttonen et al. (1992), Soininen et al. (2009, 2013a) Myodes rutilus Fungi (30–65%), fruits/seeds (10–15%), invertebrates (5–20%), lichen (ca. 10%); complemented by Vaccinium shoots Fungi (ca. 60%), lichen (ca. 25%), fruits/seeds (ca. 10%); complemented by Vaccinium shoots Grodzinski (1971, in Hansson 1985), Henttonen and Peiponen (1982), Bangs (1984) Lemmini Lemmus lemmus Moss (ca. 60%, up to 90%), grass/sedge (ca. 20%, up to 80%), dicots (ca. 10%, up to 50%) Moss (ca. 80%, up to 100%), grass/sedge (ca. 10%) Koshkina (1961, in Batzli 1993), Tast (1991), Saetnan et al. (2009), Soininen et al. (2013b) Myopus schisticolor Moss (ca. 90%); complemented by leaves Moss (ca. 90%, up to 100%); complemented by leaves Bondrup-Nielsen (1993), Eskelinen (2002) 4134 ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. Isotopic Partitioning in Small Mammals I. Calandra et al. was not tested alone. The significance of interesting linear contrasts between groups was assessed using pairwise t-tests with a Welch correction for heteroscedasticity and a Holm correction on the P-values. These contrasts were as follows: the mean differences between any two tribes, the mean differences between any two species of the same tribe, and the seasonal differences within a tribe. Even though mixing models in general, and Bayesian mixing models in particular, are often used (Wolf et al. 2009; Newsome et al. 2012; Phillips 2012), it was not possible to apply them to this dataset (see Appendix S2 for details). The open-source software R 3.1.0 (R Development Core Team 2014) was used with the following packages: car (Fox and Weisberg 2011), doBy (Højsgaard et al. 2013), irr (Gamer et al. 2012), mclust (Fraley et al. 2012), R.utils (Bengtsson 2014), RSvgDevice (Luciani et al. 2014), and xlsx (Dragulescu 2013). Results The five laboratory specimens show very little deviation from the mean, lower than any of the other populations (Table S3). The discrimination factors (the difference between the isotopic composition of the teeth and that of the food sources) are D 13 C tooth-food =+2.40 0.13&and D 15 N tooth-food =+6.01 0.16&(Fig. S1). These values fall within the ranges of discrimination factors in hairs and bone collagen published for voles (see Appendix S2). The carbon and nitrogen compositions of plants in northern Finland are given in Table S2. Our plant samples from Lapland have isotopic compositions that fall within the ranges of published data from northern America (Fig. 1A). The unsupervised classification (cluster analysis) grouped the isotopic data into four clusters (Fig. 2A), based on taxonomic affinity and season, although there is some overlap (Fig. 2B). The agreement between clustering and taxonomy plus season is significant (Cohen’s Kappa: j=0.537, P<0.001), with the fourth cluster probably lowering the jvalue. The clustering into tribes (plus season) is also supported by the two-way ANOVA, where the tribe effect is significant (d 13 C: F=155.26, P<0.001 and d 15 N: F=54.41, P<0.001). The six contrasts between tribes show that each tribe differs from the other two (Table S4). There are also significant differences between species within each tribe (effect of species nested within tribe for d 13 C: F=21.12, P<0.001 and for d 15 N: F=8.61, P<0.001; Fig. 1B–D, Tables S3–S4). Both the seasonal effect (d 13 C: F=24.89, P<0.001 and d 15 N: F=56.55, P<0.001) and the interaction between tribe and season (d 13 C: F=4.07, P=0.023 and d 15 N: F=8.01, P<0.001) are significant. The interaction between species (nested within tribe) and season is also significant (d 13 C: F=3.47, P=0.014 and d 15 N: F=7.28, P<0.001). The large seasonal differences observed between animals trapped in spring and their autumn Trees/shrubs Lichen Forbs Gram Fungi Mosses δ13C (‰) δ15N (‰) 10 0 4 –4 –24–28–32 –8 –28 –26 –24 –22 0 2 4 6 8 δ13C (‰) δ15N (‰) L. lemmus M. schisticolor Lemmini Autumn Spring Plants (review) δ13C (‰) –28 –26 –24 –22 0 2 4 6 8 Arvicolini M. agrestis M. oeconomus δ15N (‰) δ13C (‰) –28 –26 –24 –22 0 2 4 6 8 Clethrionomyini M. glareolus M. rufocanus M. rutilus δ15N (‰) (A) (B) (C) (D) Figure 1. (A) Review of carbon (d 13 C) and nitrogen (d 15 N) isotopic compositions (in &)of plant types. The ranges for mosses are based on Nadelhoffer et al. (1996), Brooks et al. (1997), McLeman (2006), and Loisel et al. (2009). Ranges for other plant types are summarized from Ben-David et al. (2001) and Drucker et al. (2010, 2012). gram =graminoids. (B–D) Carbon and nitrogen isotopic compositions of the teeth studied. Black symbols for specimens trapped in spring; white symbols for autumn. (B) Lemmini: Lemmus lemmus (squares) and Myopus schisticolor (diamonds). (C) Clethrionomyini: Myodes glareolus (upright triangles), M. rufocanus (inverted triangles) and M. rutilus (stars). (D) Arvicolini: Microtus agrestis (pentagons) and M. oeconomus (octagons). Note that (A) and (B–D) are not drawn in the same isotopic space. ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 4135 I. Calandra et al. Isotopic Partitioning in Small Mammals counterparts (Fig. 1B–D, Table S3) are driven solely by Arvicolini on d 13 C(+1.6&from autumn to spring), and by Clethrionomyini on d 15 N(3.8&from autumn to spring) (Table S4). Testing for these seasonal contrasts within each tribe will require additional samples, so interspecific seasonal differences are only described qualitatively here. There are marked increases (1–3&)ind 13 C in spring in species of both Arvicolini and Clethrionomyini (Fig. 1C–D). Additionally, a very large decrease in d 15 N (approximately 5&) is observed in spring in Myodes glareolus and M. rutilus (Fig. 1C). Discussion Our isotopic data reveal that each tribe has a specific diet, sometimes modulated by seasonality, thus extending the hypothesis that diet is structured according to phylogeny (Butet and Delettre 2011) to include Arctic arvicolines. Lemmini have very specialized diets based on mosses, where d 13 C values span a very large range, but are generally the lowest of all plant types (Nadelhoffer et al. 1996; Brooks et al. 1997; McLeman 2006). The more negative position of Lemmini in the isotopic space is consistent with this diet. The higher proportion of graminoids, fungi, lichen, and invertebrates in the diets of Arvicolini and Clethrionomyini than in that of Lemmini explains the higher d 13 C and d 15 N values in these two tribes. Clethrionomyini have the highest d 13 C values, consistent with their substantial consumption of fungi and lichens. Within Lemmini, the higher values for both isotopes (d 13 C and d 15 N) in L. lemmus can be explained by the fact that it is less dependent on mosses than M. schisticolor, consuming more graminoids to complement its diet. Within Arvicolini, the higher d 15 N values of M. oeconomus compared to M. agrestis seem to contradict the dietary data available in the literature, as M. oeconomus consumes only green plant parts (generally low d 15 N, but higher values for horsetail, for example), while M. agrestis may include some invertebrates and fungi (high d 15 N) in its diet. Within Clethrionomyini, M. rutilus has higher d 13 C values than M. rufocanus, consistent with its greater consumption of fungi and lichens (high d 13 C). High d 15 N values were expected for M. glareolus and M. rutilus, based on their consumption of fungi and invertebrates. The broad range of d 15 N values observed in these species could result from two complementary factors: (1) the low d 15 N values of lichen partially compensate for the high values of fungi and invertebrates and (2) the proportions of lichens, fungi, and invertebrates are seasonally very variable, and therefore, average values are somewhat artificial. The isotopic composition of ever-growing rodent incisors can be used to assess the isotopic composition of the diet, with a delay of one to 2 months. We found no sign of isotopic seasonality in Lemmini, consistent with their constant reliance on mosses and with observed food storage behavior (Eskelinen 2002). In contrast, seasonality was recorded in both Arvicolini and Clethrionomyini. Arvicolini rely more on graminoids, at the expense of forbs and shrubs, in winter. The winter diet should therefore have a higher carbon isotopic composition. Individuals trapped in spring do indeed have higher d 13 C values 2468 – 660 – 640 – 620 – 600 –28 –26 –24 –22 0 2 4 6 8 EII VVI VII EEE EEI EEV VEI VEV EVI VVV Number of components BIC (A) (B) δ15N (‰) δ13C (‰) Groups Arvicolini Clethrionomyini Lemmini Figure 2. Unsupervised model-based classification. (A) Bayesian information criterion (BIC) relative to the number of components (i.e., clusters) tested by each model. The arrow highlights the best model (highest BIC): spherical variance–covariance matrices and varying volume (model VII, white triangles) with 4 components. (B) Clustering of each arvicoline tooth sample according to the unsupervised classification (symbols), with taxonomic affinities identified by gray shading. The samples are shown in the carbon (d 13 C) and nitrogen (d 15 N) isotopic space of arvicoline incisors. 4136 ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. Isotopic Partitioning in Small Mammals I. Calandra et al. than autumn specimens, thus reflecting their winter diet. Microtus agrestis and M. oeconomus are known to show interference competition in summer, favoring the latter (Henttonen et al. 1977). Coexistence is thought to be linked to relaxed competition in winter, when M.oeconomus shows a higher level of individual activity than M. agrestis (Hoset and Steen 2007). Stable isotopes also suggest resource partitioning during winter, as M.oeconomus seems to have higher d 15 N values. Alternatively, the greater activity of M. oeconomus in winter, despite potential food shortage, may result in nutritional stress. This may induce impaired nitrogen balance, leading to higher d 15 N values unrelated to any change in food supply (Fuller et al. 2005; Petzke et al. 2010). As the lighter isotope is preferentially excreted and not replaced by external sources when the organism is fasting, protein synthesis during tooth renewal should lead to an increase in d 15 N (Vanderklift and Ponsard 2003; Mart ınez del Rio et al. 2009; Lee et al. 2012). Myodes rufocanus is close to Arvicolini, in the isotopic space, but without any perceptible seasonal shift. Its winter signature is very similar to that of M. agrestis, suggesting that some resource competition occurs. This observation is in agreement with the over-winter population dynamics of M. rufocanus, related to the size of the vole community (Hansen et al. 1999b). The winter diets of M. glareolus and M. rutilus include fewer or no invertebrates (high d 15 N), but more lichens (low d 15 N, high d 13 C), most likely because of food availability, and should thus result in lower d 15 N but higher d 13 C values. These species indeed show a marked decrease in d 15 N and an increase in d 13 C between animals trapped in autumn and those trapped in spring. Interestingly, resource use between these two closely related species, which produce viable F1 individuals (Grant 1974) and present mtDNA introgression (Boraty nski et al. 2014), seems closer in winter than in summer, suggesting probable resource competition under the snow. Stable isotope analysis suggests strong resource partitioning between these species and M. rufocanus. This is in agreement with the suggested absence of competition between M. rufocanus and M. rutilus (Mal’kova and Yakimenko 2007). The present study is, to our knowledge, the only ecological analysis of carbon and nitrogen isotopes from rodent teeth (see Ga zsiorowski et al. 2014 for palaeoenvironmental reconstructions). It is also the only study on isotopic ecology based on the teeth of small arctic mammals. We show that food resources are partitioned among arvicoline tribes and species. By analyzing only the apical part of ever-growing incisors, it is possible to access the late winter diet of arvicolines trapped in June, when the snow cover is gone, thus bringing to light resource use of this key guild in the Arctic food web. This methodology could be used on other rodent species, and in other contexts where observations and trapping are difficult during certain periods. This approach also opens up the possibility of using stable isotope analysis on pellets to track the ecology of both rodents and birds of prey. Acknowledgments We thank J. Guigue, L. Bruneau, and N. Guichard (Universit e Bourgogne Franche-Comt e) for their help with the experimental procedures. R. Virtanen (Botanical Museum, University of Oulu, Finland) and E. Luokkanen (Lapland Regional Environment Centre, Finland) kindly identified our plant samples. Carmela Chateau-Smith (Universit e Bourgogne Franche-Comt e) corrected the English manuscript. Laboratory voles were raised at the University of Helsinki, under conditions approved by the National Animal Experiment Board (license numbers ESLH-200705823/Ym-23). This work is a contribution of the Biog eosciences and EPHE laboratories. Conflict of Interest The authors declare no conflict of interest. References Bangs, E. E. 1984. Summer food habits of voles, Clethrionomys rutilus and Microtus pennsylvanicus, on the Kenai Peninsula, Alaska. Can. Field Nat. 98:489–492. Batzli, G. O. 1993. Food selection by lemmings. Pp. 281–301 in N. C. Stenseth and R. A. Ims, eds. The biology of lemmings. Academic Press, San Diego. Batzli, G. O., and H. Henttonen. 1990. Demography and resource use by Microtine rodents near Toolik Lake, Alaska, USA. Arct. Alp. Res. 22:51–64. Ben-David, M., and E. A. Flaherty. 2012. Stable isotopes in mammalian research: a beginner’s guide. J. Mammal. 93:312–328. Ben-David, M., E. Shochat, and L. G. Adams. 2001. Utility of stable isotope analysis in studying foraging ecology of herbivores: examples from moose and caribou. Alces 37:421–434. Bengtsson, H. (2014) R.utils: Various programming utilities. R package version 1.29.8. Bilodeau, F., A. J. Kenney, B. S. Gilbert, E. Hofer, G. Gauthier, D. G. Reid, et al. 2013. Evaluation of a technique to trap lemmings under the snow. Arctic 66:32–36. Bondrup-Nielsen, S. 1993. Food preference and diet of the wood lemming (Myopus schisticolor). Pp. 303–309 in N. C. Stenseth and R. A. Ims, eds. The biology of lemmings. Academic Press, San Diego. Boraty nski, Z., J. Melo-Ferreira, P. C. Alves, S. Berto, E. Koskela, O. T. Pentik€ ainen, et al. 2014. Molecular and ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 4137 I. Calandra et al. Isotopic Partitioning in Small Mammals ecological signs of mitochondrial adaptation: consequences for introgression? Heredity 113:277–286. Brooks, J. R., L. B. Flanagan, N. Buchmann, and J. R. Ehleringer. 1997. Carbon isotope composition of boreal plants: functional grouping of life forms. Oecologia 110:301–311. Butet, A., and Y. Delettre. 2011. Diet differentiation between European arvicoline and murine rodents. Acta Theriol. 56:297–304. Cerling, T. E., J. M. Harris, M. G. Leakey, B. H. Passey, and N. E. Levin. 2010. Stable Carbon and Oxygen Isotopes in East African Mammals: Modern and Fossil. Pp. 949–960 in L. Werdelin and W. J. Sanders, eds. Cenozoic mammals of Africa. University of California Press, Berkeley. Chappell, M. A. 1980. Thermal energetics and thermoregulatory costs of small arctic mammals. J. Mammal. 61:278–291. Clementz, M. T. 2012. New insight from old bones: stable isotope analysis of fossil mammals. J. Mammal. 93:368–380. Cohen, J. 1960. A coefficient of agreement for nominal scales. Educ. Psychol. Measur. 20:37–46. Craig, H. 1953. The geochemistry of the stable carbon isotopes. Geochim. Cosmochim. Acta 3:53–92. Dalerum, F., and A. Angerbj€ orn. 2005. Resolving temporal variation in vertebrate diets using naturally occurring stable isotopes. Oecologia 144:647–658. Dragulescu, A. A. 2013. xlsx: Read, write, format Excel 2007 and Excel 97/2000/XP/2003 files. R package version 0.5.5. Drucker, D. G., K. A. Hobson, J.-P. Ouellet, and R. Courtois. 2010. Influence of forage preferences and habitat use on 13 C and 15 N abundance in wild caribou (Rangifer tarandus caribou) and moose (Alces alces) from Canada. Isot. Environ. Health Stud. 46:107–121. Drucker, D. G., K. A. Hobson, S. C. M€ unzel, and A. Pike-Tay. 2012. Intra-individual variation in stable carbon (d 13 C) and nitrogen (d 15 N) isotopes in mandibles of modern caribou of Qamanirjuaq (Rangifer tarandus groenlandicus) and Banks Island (Rangifer tarandus pearyi): implications for tracing seasonal and temporal changes in diet. Int. J. Osteoarchaeol. 22:494–504. Duchesne, D., G. Gauthier, and D. Berteaux. 2011. Habitat selection, reproduction and predation of wintering lemmings in the Arctic. Oecologia 167:967–980. Eskelinen, O. 2002. Diet of the wood lemming Myopus schisticolor. Ann. Zool. Fenn. 39:49–57. Fox, J., and S. Weisberg 2011. An {R} Companion to Applied Regression. R package version 0.5.5. Fraley, C., and A. E. Raftery. 2002. Model-based clustering, discriminant analysis and density estimation. J. Am. Stat. Assoc. 97:611–631. Fraley, C., A. E. Raftery, T. B. Murphy, and L. Scrucca. 2012. Mclust Version 4 for R: Normal Mixture Modeling for Model-Based Clustering, Classification, and Density Estimation. Technical Report, 597. Fuller, B. T., J. L. Fuller, N. E. Sage, D. A. Harris, T. C. O’Connell, and R. E. M. Hedges. 2005. Nitrogen balance and d 15 N: why you’re not what you eat during nutritional stress. Rapid Commun. Mass Spectrom. 19:2497–2506. Gamer, M., J. Lemon, I. Fellows, and P. Singh. 2012. Various coefficients of interrater reliability and agreement. R package version 0.84. Ga zsiorowski, M., H. Hercman, and P. Socha. 2014. Isotopic analysis (C, N) and species composition of rodent assemblage as a tool for reconstruction of climate and environment evolution during Late Quaternary: a case study from Bi snik Cave (Cze zstochowa Upland, Poland). Quatern. Int. 339–340:139–147. Gehler, A., T. T€ utken, and A. Pack. 2012. Oxygen and carbon isotope variations in a modern rodent community – implications for palaeoenvironmental reconstructions. PLoS One 7:e49531. Grant, P. R. 1974. Reproductive compatibility of voles from separate continents (Mammalia: Clethrionomys). J. Zool. 174:245–254. Grimes, S. T., M. E. Collinson, J. J. Hooker, D. P. Mattey, N. V. Grassineau, and D. Lowry. 2004. Distinguishing the diets of coexisting fossil theridomyid and glirid rodents using carbon isotopes. Palaeogeogr. Palaeoclimatol. Palaeoecol. 208:103–119. Hansen, T. F., N. C. Stenseth, and H. Henttonen. 1999a. Multiannual vole cycles and population regulation during long winters: an analysis of seasonal density dependence. Am. Nat. 154:129–139. Hansen, T. F., N. C. Stenseth, H. Henttonen, and J. Tast. 1999b. Interspecific and intraspecific competition as causes of direct and delayed density dependence in a fluctuating vole population. Proc. Natl Acad. Sci. USA 96:986–991. Hansson, L. 1969. Spring populations of small mammals in central Swedish Lapland in 1964-68. Oikos 20:431–450. Hansson, L. 1971a. Small rodent food, feeding and population dynamics. A comparison between granivorous and herbivorous species in Scandinavia. Oikos 22:183–198. Hansson, L. 1971b. Habitat, food and population dynamics of the field vole Microtus agrestis (L.) in south Sweden. Viltrevy 8:267–373. Hansson, L. 1985. Clethrionomys food: generic, specific and regional characteristics. Ann. Zool. Fenn. 22:315–318. Hansson, L., and T.-B. Larsson. 1978. Vole diet on experimentally managed reforestation areas in Northern Sweden. Holarctic Ecol. 1:16–26. Henttonen, H., and V. A. Peiponen. 1982. Clethrionomys rutilus (Pallas, 1779) –Polarr€ otelmaus. Pp. 165–176 in J. Niethammer and F. Krapp, eds. Handbuch der S€ augetiere Europas. Akademische Verlagsgesellschaft, Wiesbaden. Henttonen, H., and J. Viitala. 1982. Clethrionomys rufocanus (Sundevall, 1846) –Graur€ otelmaus. Pp. 147–164 in J. Niethammer and F. Krapp, eds. Handbuch der S€ augetiere Europas. Akademische Verlagsgesellschaft, Wiesbaden. 4138 ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. Isotopic Partitioning in Small Mammals I. Calandra et al. Henttonen, H., A. Kaikusalo, J. Tast, and J. Viitala. 1977. Interspecific competition between small rodents in subarctic and boreal ecosystems. Oikos 29:581–590. Henttonen, H., L. Hansson, and T. Saitoh. 1992. Rodent dynamics and community structure: Clethrionomys rufocanus in northem Fennoscandia and Hokkaido. Ann. Zool. Fenn. 29:1–6. Højsgaard, S., U. Halekoh, J. Robison-Cox, K. Wright, and A. A. Leidi. 2013. Doby: doBy –Groupwise summary statistics, LSmeans, general linear contrasts, various utilities. R package version 4.5-10. Hopley, P. J., A. G. Latham, and J. D. Marshall. 2006. Palaeoenvironments and palaeodiets of mid-Pliocene micromammals from Makapansgat Limeworks, South Africa: a stable isotope and dental microwear approach. Palaeogeogr. Palaeoclimatol. Palaeoecol. 233:235–251. Hoset, K. S., and H. Steen. 2007. Relaxed competition during winter may explain the coexistence of two sympatric Microtus species. Ann. Zool. Fenn. 44:415–424. Ims, R. A., J.-A. Henden, A. V. Thingnes, and S. T. Killengreen. 2013. Indirect food web interactions mediated by predator, € A ırodent dynamics: relative roles of lemmings and voles. Biol. Lett. 9:20130802. Kausrud, K. L., A. Mysterud, H. Steen, J. O. Vik, E. Ostbye, B. Cazelles, et al. 2008. Linking climate change to lemming cycles. Nature 456:93–97. Klevezal, G. A., M. Pucek, and L. I. Sukhovskaja. 1990. Incisor growth in voles. Acta Theriol. 35:331–344. Kohn, M. J., and T. E. Cerling. 2002. Stable Isotope Compositions of Biological Apatite. Pp. 455–488 in M. L. Kohn, J. Rakovan and J. M. Hughes, eds. Phosphates – geochemical, geobiological and materials importance. Mineralogical Society of America and Geochemical Society, Chantilly and Saint Louis, USA. Lee, T. N., C. L. Buck, B. M. Barnes, and D. M. O’Brien. 2012. A test of alternative models for increased tissue nitrogen isotope ratios during fasting in hibernating arctic ground squirrels. J. Exp. Biol. 215:3354–3361. Loisel, J., M. Garneau, and J.-F. H elie. 2009. Modern Sphagnum d 13 C signatures follow a surface moisture gradient in two boreal peat bogs, James Bay lowlands, Quebec. J. Quat. Sci. 24:209–214. Luciani, T. J., M. Decorde, and V. Lise 2014. RSvgDevice: An R SVG graphics device. R package version 0.6.4.4. Mal’kova, M. G., and V. V. Yakimenko. 2007. Spatial population structure in forest voles of the genus Clethrionomys in the southern taiga of the Middle Irtysh region. Russian J. Ecol., 38:190–197. Mariotti, A. 1983. Atmospheric nitrogen is a reliable standard for natural 15N abundance measurements. Nature 303:685– 687. Mart ınez del Rio, C., N. Wolf, S. A. Carleton, and L. Z. Gannes. 2009. Isotopic ecology ten years after a call for more laboratory experiments. Biol. Rev. 84:91–111. McLeman, C.I.A. 2006. Determining the relationships between forage use, climate and nutritional status of barren ground caribou, Rangifer tarandus groenlandicus, on Southampton Island, Nunavut, using stable isotopes analysis of d13C and d15N. Msc thesis, University of Waterloo, Ontario, Canada. Nadelhoffer, K., G. Shaver, B. Fry, A. Giblin, L. Johnson, and R. McKane. 1996. 15N natural abundances and N use by tundra plants. Oecologia 107:386–394. Newsome, S. D., J. D. Yeakel, P. V. Wheatley, and M. T. Tinker. 2012. Tools for quantifying isotopic niche space and dietary variation at the individual and population level. J. Mammal. 93:329–341. Pauli, J. N., B. Zuckerberg, J. P. Whiteman, and W. Porter. 2013. The subnivium: a deteriorating seasonal refugium. Front. Ecol. Environ. 11:260–267. Petty, S. K., B. Zuckerberg, and J. N. Pauli. 2015. Winter Conditions and Land Cover Structure the Subnivium, A Seasonal Refuge beneath the Snow. PLoS One 10:e0127613. Petzke, K. J., B. T. Fuller, and C. C. Metges. 2010. Advances in natural stable isotope ratio analysis of human hair to determine nutritional and metabolic status. Curr. Opin. Clin. Nutr. Metab. Care 13:532–540. Phillips, D. L. 2012. Converting isotope values to diet composition: the use of mixing models. J. Mammal. 93:342– 352. R Development Core Team. 2014. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Reid, D. G., and C. J. Krebs. 1996. Limitations to collared lemming population growth in winter. Can. J. Zool. 74:1284–1291. Robovsky, J., V.  Ri c ankov a, and J. Zrzavy. 2008. Phylogeny of Arvicolinae (Mammalia, Cricetidae): utility of morphological and molecular data sets in a recently radiating clade. Zoolog. Scr. 37:571–590. Saetnan, E. R., J. O. Gjershaug, and G. O. Batzli. 2009. Habitat use and diet composition of Norwegian lemmings and field voles in central Norway. J. Mammal. 90:183–188. Sare, D. T. J., J. S. Millar, and F. J. Longstaffe. 2005. Tracing dietary protein in red-backed voles (Clethrionomys gapperi) using stable isotopes of nitrogen and carbon. Can. J. Zool. 83:717–725. Soininen, E. M., A. Valentini, E. Coissac, C. Miquel, L. Gielly, C. Brochmann, et al. 2009. Analysing diet of small herbivores: the efficiency of DNA barcoding coupled with high-throughput pyrosequencing for deciphering the composition of complex plant mixtures. Front. Zool. 6:16. Soininen, E. M., V. T. Ravolainen, K. A. Br  athen, N. G. Yoccoz, L. Gielly, and R. A. Ims. 2013a. Arctic small rodents have diverse diets and flexible food selection. PLoS One 8: e68128. Soininen, E. M., L. Zinger, L. Gielly, E. Bellemain, K. A. Br  athen, C. Brochmann, et al. 2013b. Shedding new light on ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 4139 I. Calandra et al. Isotopic Partitioning in Small Mammals the diet of Norwegian lemmings: DNA metabarcoding of stomach content. Polar Biol. 36:1069–1076. Soininen, E. M., D. Ehrich, N. Lecomte, N. G. Yoccoz, A. Tarroux, D. Berteaux, et al. 2014. Sources of variation in small rodent trophic niche: new insights from DNA metabarcoding and stable isotope analysis. Isot. Environ. Health Stud. 50:361–381. Soininen, E. M., G. Gauthier, F. Bilodeau, D. Berteaux, L. Gielly, P. Taberlet, et al. 2015. Highly overlapping winter diet in two sympatric lemming species revealed by DNA Metabarcoding. PLoS ONE 10:e0115335. Stenseth, N. C., and R. A. Ims. 1993. Food selection, individual growth and reproduction –an introduction. Pp. 263–280 in N. C. Stenseth and R. A. Ims, eds. The biology of lemmings. Academie Press, San Diego, CA. Stenseth, N. C., L. Hansson, and A. Myllym€ aki. 1977. Food selection of the field vole Microtus agrestis. Oikos 29:511–524. Tast, J. 1974. The food and feeding habits of the root vole, Microtus oeconomus, in Finnish Lapland. Aquilo Ser. Zool. 15:25–32. Tast, J. 1991. Will the Norwegian Lemming become endangered if climate becomes warmer? Arct. Alp. Res. 23:53–60. Vanderklift, M., and S. Ponsard. 2003. Sources of variation in consumer-diet d 15 N enrichment: a meta-analysis. Oecologia 136:169–182. Viro, P., and J. Niethammer. 1982. Clethrionomys glareolus (Schreber, 1780) –R€ otelmaus. Pp. 128–146 in J. Niethammer and F. Krapp, eds. Handbuch der S€ augetiere Europas. Akademische Verlagsgesellschaft, Wiesbaden. Wolf, N., S. A. Carleton, and C. Mart ınez del Rio. 2009. Ten years of experimental animal isotopic ecology. Funct. Ecol. 23:17–26. Supporting Information Additional Supporting Information may be found in the online version of this article: Figure S1. Carbon (d 13 C) vs. nitrogen (d 15 N) isotopic compositions (in &) of the food pellet (grey triangle), and of laboratory vole teeth before (white dots) and after (black dots) decarbonation with HCl. Dotted lines connect samples before and after HCl treatment. Table S1. Information about the specimens studied: tribe (following Robovsky et al. 2008); species; inventory number (ID); trapping locality, season and year; carbon (d 13 C) and nitrogen (d 15 N) isotopic compositions (in &), and carbon ([C]) and nitrogen ([N]) elemental compositions (in %) of the teeth studied. Table S2. Information on the food samples studied: species; sampling locality; part of plant analysed; carbon (d 13 C) and nitrogen (d 15 N) isotopic compositions (in &), carbon ([C]) and nitrogen ([N]) elemental compositions (in %), and carbon to nitrogen ratio (C:N =[C]/ [N]). Table S3. Mean and standard deviation (SD) of the carbon (d 13 C) and nitrogen (d 15 N) isotopic compositions (in &) of the arvicoline teeth studied. Table S4. Results of the pairwise t-tests (df: degree of freedom, P:Pvalue, t:t-statistics). Appendix S1. Details of the method: preparation of arvicoline heads and pre-treatment. Appendix S2. Bayesian mixing models and discrimination factors. 4140 ª2015 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. Isotopic Partitioning in Small Mammals I. Calandra et al.