Solving the woolly mammoth conundrum: amino acid ¹⁵N-enrichment suggests a distinct forage or habitat.

Understanding woolly mammoth ecology is key to understanding Pleistocene community dynamics and evaluating the roles of human hunting and climate change in late Quaternary megafaunal extinctions. Previous isotopic studies of mammoths' diet and physiology have been hampered by the 'mammoth conundrum': woolly mammoths have anomalously high collagen δ(15)N values, which are more similar to coeval carnivores than herbivores, and which could imply a distinct diet and (or) habitat, or a physiological adaptation. We analyzed individual amino acids from collagen of adult woolly mammoths and coeval species, and discovered greater  (15)N enrichment in source amino acids of woolly mammoths than in most other herbivores or carnivores. Woolly mammoths consumed an isotopically distinct food source, reflective of extreme aridity, dung fertilization, and (or) plant selection. This dietary signal suggests that woolly mammoths occupied a distinct habitat or forage niche relative to other Pleistocene herbivores.

Woolly mammoths (Mammuthus primigenius) were keystone herbivores in the Pleistocene mammoth steppe12. This megacontinental biome was inhabited by a now-extinct community of mammals, dominated by woolly mammoth, horse and bison. The mammoth steppe reached from north-western Canada, across the exposed Bering Isthmus, to Western Europe3. The ecological role of woolly mammoths within this ecosystem has been a subject of vigorous investigation345. Reconstructions of woolly mammoth behaviour and physiology have been largely based on morphology4. Isotopic studies of bulk tissues have provided independent tests of morphology-based hypotheses, as well as suggesting new ones36789. Compound-specific isotopic studies can provide a further level of understanding of ecosystem functioning within the mammoth steppe.

Bulk collagen nitrogen isotopic compositions (δ15NBulk) are commonly used in ecological studies to reveal the diet and trophic level of a species, as these values typically reflect the isotopic compositions of the plants at the base of the food web plus a 2–5‰ increase with each trophic level10. As a result, the δ15NBulk values of mammoth-steppe herbivore collagen are commonly ~+6‰ where the values of carnivores (~+9‰) are higher, reflective of this trophic enrichment3. The carnivore- rather than herbivore-like δ15NBulk values of woolly mammoth collagen (~+8‰3) are seemingly problematic and therefore require examination. The various hypotheses to explain this phenomenon (unique diet, niche feeding in a special habitat or distinct metabolic processes3678911121314) have different implications for our understanding of the now-vanished mammoth steppe ecosystem, woolly mammoth ecology, and related factors that contributed to extirpation of the woolly mammoth in this region.

Woolly mammoths may have consumed plants with higher δ15N values, such as graminoids and herbs rather than woody vegetation71213, as suggested by the morphology of their enamel plates4. However, an herbaceous diet alone is not sufficient to fully explain the woolly mammoth’s high δ15N values; some further form of habitat or plant selection is also required15. While modern Arctic graminoids and forbs from some sites have a δ15N range of –0.3 to +10‰, the average value of these species ranges from ~+1 to ~+4‰16, and still other studies have reported maximum δ15N values for modern sedges of +2‰17 and for modern herbs of +5.3‰18. The majority of these plants, therefore, are not sufficiently enriched in 15N to explain the woolly mammoth δ15NBulk values. Plants growing in drier habitats, however, have higher δ15N values than plants from a more mesic environment7, and woolly mammoths may have eaten plants experiencing water-stress38. Water stress can also cause 13C-enrichment of plants19. This enrichment, however, is unlikely to be directly observable in woolly mammoth collagen, because its carbon isotopic composition is likely dominated by the low δ13C values of fat reserves used to survive the winter7.

Several other factors could also have contributed to high δ15N values for woolly mammoth collagen. Woolly mammoths that had small ranges, or repeatedly travelled the same routes, could have deposited significant quantities of faeces in those areas12, causing 15N-enrichment in plants arising from this dung fertilization20. Partially decayed plant material can also have higher δ15N values than the original living plant21. Woolly mammoths may have removed snow and ice cover by trampling and (or) with their tusks5, allowing them to forage on winter-killed plants generally not utilized by other large herbivores that did not share the mammoth ecological niche. It has also been proposed that woolly mammoths had distinct metabolic processes, such as increased levels of nitrogen recycling associated with winter starvation2223 or poor quality food with low protein levels378911131424, or that woolly mammoths engaged in coprophagy12.

The nitrogen isotopic compositions of the individual amino acids in collagen, as opposed to bulk collagen, enable discrimination between 15N-enrichment occurring at the base of the food chain prior to consumption (source amino acids) versus that associated with metabolic processes (trophic amino acids) (Fig. 1). Phenylalanine (Phe) and glutamate (Glu) have been identified as characteristic of source and trophic amino acids, respectively2526. The δ15NPhe value reflects the isotopic composition of those amino acids in plants at the base of the food web, while the Δ15NGlu-Phe spacing (δ15NGlu – δ15NPhe) serves as a proxy for metabolic enrichment of 15N in the consumer’s body25.

Figure 1

Simplified pathway for nitrogen incorporation from soil to animal protein.

Arrows represent: uptake (dashed line), chemical transformations (solid line), and metabolic processes (solid curves). (a) Plant NO3– uptake. (b) NO3– converted to glutamine31. (c) NH4+ uptake. (d) NH4+ converted to glutamine by attachment to glutamate31. (e) Glutamine supplies nitrogen for synthesis of other amino acids. The associated shift in δ15N values depends on the specific amino acid, plant part and plant type31. (f) Consumption of amino acids by the animal. (g) Source amino acids are minimally involved in metabolic processes, undergoing small changes in δ15N values (e.g. increases in δ15NPhe values from diet to consumer tissue are commonly ≤ 2‰20). (h) Trophic amino acids are heavily involved in metabolic processes, undergoing enrichment in 15N (e.g. increases in δ15NGlu values from diet to consumer tissue are commonly 6–7‰20). Katherine Allan drew the images of grass and mammoth in Figure 1.

Results

Eight Pleistocene megafauna species were analyzed in this study. These included four herbivore species: woolly mammoth (Mammuthus primigenius), mastodon (Mammut americanum), horse (Equus sp.) and giant beaver (Castoroides ohioensis), and four carnivore species: brown bear (Ursus arctos), scimitar cat (Homotherium serum), wolf (Canis lupus), and short-faced bear (Arctodus simus) (see Supplementary Table S1). All samples were obtained from specimens collected near Old Crow, Yukon, Canada (latitude: 67°34’N; longitude: 139°48’W). A subset of these specimens was dated, including both herbivores and carnivores. Two horse specimens were dated to 18,370 and 27,180 14C years BP. The rest of the specimens yielded effectively non-finite radiocarbon dates ≥37,200 14C yr BP, and one specimen was dated by context to ~140,000 years BP (see Supplementary Table S1)2728. All collagen samples were considered well preserved based on their collagen yields, C/N ratios, and carbon and nitrogen contents29 (see Supplementary Table S2).

The δ15NBulk values for the Old Crow samples follow the pattern previously observed for Pleistocene megafauna3811; woolly mammoth collagen generally has δ15NBulk values similar to the carnivores and higher than the other herbivores, with some overlap with horses (Fig. 2a). The δ15NPhe values of woolly mammoth collagen, however, are higher than those of the carnivores and most of the other herbivores (Fig. 2b); horses with high δ15NBulk values for collagen show the most overlap with the δ15NPhe values of woolly mammoths. Woolly mammoth Δ15NGlu-Phe spacings overlap those of the other herbivores but are lower than the Δ15NGlu-Phe spacings of the carnivores, extending to negative values for most samples (Fig. 2c). Negative Δ15NGlu-Phe spacings have observed in terrestrial herbivores previously2530 and may be the result of relatively high δ15NPhe values in vascular plants31.

Figure 2

Nitrogen isotopic compositions of Pleistocene Old Crow megafauna.

(a) Bulk collagen nitrogen isotopic compositions (δ15NBulk). Results for woolly mammoths are displayed in blue, other herbivores in purple, and carnivores in red. (b) Phenyalanine (source) amino acid nitrogen isotopic compositions (δ15NPhe) of each species. (c) Difference between the nitrogen isotopic composition of glutamate and phenylalanine (Δ15NGlu-Phe) for each species.

Discussion

The high δ15NPhe values of the woolly mammoth imply that they selectively consumed plants more enriched in 15N than forage consumed by most of the other herbivores. The fact that the δ15NPhe values of woolly mammoths are higher than those of carnivores suggests that the latter consumed herbivores subsisting on less 15N-rich forage than consumed by woolly mammoths. In short, the carnivores did not consume significant quantities of woolly mammoth. Partial overlap between horse and woolly mammoth δ15NBulk and δ15NPhe values likely indicates that horses exploited a similar niche to the woolly mammoth in some cases. The low Δ15NGlu-Phe spacings of woolly mammoths indicate that their δ15NBulk values arise from the higher δ15N values of the plants they consumed, and not from a specialized metabolic process.

It seems that woolly mammoths occupied a specialized dietary or habitat niche. A dietary niche implies that woolly mammoths selected particular herbaceous plants or consumed large quantities of decayed plants in winter, while a habitat niche suggests that woolly mammoths occupied more arid habitats, or lived in distinct ranges where they left considerable quantities of dung that fertilized the plants growing there. While some Old Crow horses appear also to have exploited such a niche, it was not generally shared by other mammoth steppe megafauna. The Old Crow samples likely represent various time intervals through the late Pleistocene and potentially varied climate regimes. The fact that the relative differences in average δ15NBulk values among herbivore species are generally consistent across the mammoth steppe3 suggests that most herbivore species ate the same forage types regardless of climatic differences or time period. This implies that mammoth steppe herbivores targeted specific forage types.

Two significant conclusions arise from these observations. First, the woolly mammoth occupied a distinct niche from other contemporaneous herbivores. This unique habitat or forage existed across the entirety of the mammoth steppe, although its size may have varied with changing climate across geographic or temporal zones. Other evidence of the woolly mammoth’s dependence on a specialized niche may be provided by the retraction of woolly mammoth populations into small, isolated refugia during the last interglacial warm period (MIS 5e, 130-116 kyr BP), and the subsequent re-expansion upon return to glacial conditions32. An investigation of the isotopic compositions of woolly mammoths from a variety of time periods and sites across the mammoth steppe could reveal the extent of adaptability of the woolly mammoth to disruptions in its niche, such as may have occurred with the onset of climatic shifts at the end of the Pleistocene33.

Second, the horse δ15NBulk and δ15NPhe values overlap those of the woolly mammoth and the other herbivores. This implies that horses fed from a wide diversity of habitats or forage types, including the woolly mammoths’ niche. Such behaviour would suggest that the mammoth steppe ecosystem supported herbivores with a variety of ecological adaptations, and that even in the Pleistocene Arctic, resources were sufficiently abundant to support both specialist and generalist strategies.

Methods

Radiocarbon dating

Radiocarbon dates for six woolly mammoths discussed here have been published previously2728. A further subset of samples was dated for this study; these included two other herbivores and two carnivores (see Supplementary Table S1). Collagen was extracted, combusted, graphitized and radiocarbon dated at the University of Arizona Accelerator Mass Spectrometry (AMS) Laboratory. All dates are presented as uncalibrated radiocarbon years before present (mean ±1SD).

Bulk collagen nitrogen isotope (δ 15NBulk) analysis

Collagen for δ15NBulk analysis was extracted at the Laboratory for Stable Isotope Science (LSIS) following previously published methods28, or was previously extracted and analyzed for another study2728 (see Supplementary Table S2). The δ15NBulk values were obtained using a Costech Elemental Combustion System (ECS 4010) attached to a ThermoFisher Delta Plus XL IRMS or to a ThermoFisher Delta V Plus IRMS. The δ15NBulk values were measured over three analytical sessions. In the first two analytical sessions, the δ15N values were calibrated to AIR using USGS40 (L-glumatic acid; accepted value –4.52‰34) and IAEA-N2 (ammonium sulfate; accepted value +20.3‰35), while the third analytical session substituted USGS41 (L-glumatic acid; accepted value +47.57‰34) for IAEA-N2. The same standards were used to create calibration curves for determining carbon and nitrogen contents (wt%) of each sample, from which C/N ratios were calculated. Keratin (MP Biomedicals Inc., Cat No. 90211, Lot No. 9966H) was used as an internal standard in each analytical session. For a total of 18 keratin measurements over the three analytical sessions, average values (mean ±1 SD) were δ15N = +6.4 ± 0.2‰ (accepted value = +6.4‰), and C/N = 3.6 ± 0.4 (accepted value = 3.7). The standard deviation of a sample analyzed as an instrumental duplicate was δ15NBulk = ±0.0‰, and C/N = ±0.1 (1 SD). The standard deviations (1 SD) for method duplicates of δ15NBulk values ranged from 0.0 to 0.2‰, and for C/N ratios, from 0.0 to 0.3. All samples were considered to be well preserved based on their extraction yields, C/N ratios, and carbon and nitrogen contents2936. Eight samples had high carbon and/or nitrogen contents, but this anomaly likely arises from a weighing error as they were well preserved by other measures.

Amino acid nitrogen isotope (δ 15NAmino Acid) Analysis

Using collagen first extracted for δ15NBulk measurements, amino acids were hydrolysed, derivatized to their N-acetyl-methyl ester derivative, and their individual δ15NAmino Acid values measured using an Aligent 6890N-ThermoFisher Gas Chromatograph-Combustion 3-ThermoFisher Delta Plus XL IRMS. An Agilent Technologies VF-23MS column was used in the GC. We followed published methods2537 with only slight modifications: (i) the quantity of collagen hydrolysed was increased from 2 to 6 mg, and the quantity derivatized was increased from 0.25 to 1.5 mg; and (ii) the initial GC column temperature was set at 60 °C instead of 40 °C, and its final temperature of 250 °C was held for 15 min instead of 20 min. A representative chromatogram is shown in Supplementary Fig. S1. Three reference gas pulses were introduced into the IRMS at the beginning of each analytical session and one pulse was introduced at the end of each session. The isotopic composition of the reference gas was calibrated using four amino acid standards. Three of these standards, alanine, leucine and phenylalanine, were purchased as their NACME derivative from Sigma Aldrich. The fourth, proline, was purchased as an amino acid and derivatized in-house. The nitrogen isotopic compositions of the derivatized standards were established by multiple measurements performed in the same manner as used for isotopic analysis of bulk collagen, and calibration to AIR using international standards. The amino acid reference standards were injected every three to five runs. All samples were analyzed a minimum of three times, and the average variation was ±0.7‰ (1 SD) for phenyalanine, and ±0.7‰ (1 SD) for glutamate, with a range of 0.0–3.9‰. An internal standard, norleucine, was also analyzed. Its nitrogen isotopic compositions were offset from expected values by an average of +1.3‰.

Additional Information

How to cite this article: Schwartz-Narbonne, R. et al. Solving the woolly mammoth conundrum: amino acid 15N-enrichment suggests a distinct forage or habitat. Sci. Rep. 5, 09791; doi: 10.1038/srep09791 (2015).