Proceedings of the Japan Academy, Series B
Online ISSN : 1349-2896
Print ISSN : 0386-2208
ISSN-L : 0386-2208
Review
A new era of isotope ecology: Nitrogen isotope ratio of amino acids as an approach for unraveling modern and ancient food web
Naohiko OHKOUCHI
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2023 Volume 99 Issue 5 Pages 131-154

Details
Abstract

Food web research is rapidly expanding through study of natural fractional abundance of 15N in individual amino acids. This paper overviews the principles of this isotope approach, and from my perspective, reanalyzes applications, and further extends the discussion. It applies kinetic isotope effects that enriches 15N in certain amino acids associated with the metabolic processes, which was clearly demonstrated by observations of both natural ecosystem and laboratory experiments. In trophic processes ‘trophic amino acids’ such as glutamic acid that significantly enrich 15N, whereas ‘source amino acids’ such as phenylalanine and methionine show little 15N enrichment. Through various applications conducted over the years, the principles of the method have shown to operate well and disentangle complex food webs and relevant problems. Applications include food chain length estimate, nitrogen resource assessment, tracking fish migration, and reconstruction of paleodiet. With this approach, considerations of a wide range of classical issues have been reinvigorated, while in the same time, new challenging frontiers are emerging.

Introduction

Attempts to detect small fluctuations of the nitrogen isotope ratio (15N/14N) in natural amino acids date back to as early as 1930s, when applied isotope chemistry was still in its infancy.1) The pioneer study of Schenheimer and Rittenberg at Columbia University was followed by decades of efforts to isolate the individual amino acids from living organisms in order to determine their nitrogen isotope ratio.2)4) To precisely detect the very small isotopic variations mostly ranging from 0.36 to 0.39% as fractional abundance of 15N, it is a prerequisite that the amino acids are isolated from the samples with high degree of purity without modifying the isotope ratios. In 1980s to early 1990s a few laboratories overcame the technical difficulties, and outstanding progress was made that provided insights into microbial processes and the diet of fossil animals.5)11) Around 1990, a continuous-flow analytical system, ‘gas chromatography/combustion/isotope-ratio mass spectrometry (GC/C/IRMS)’ was newly developed and led to accelerate studies. Even though the isotopic compositions of all 20 amino acids cannot be determined due to the chromatographic limitations,12),13) merits of this instrument are significant; both high-throughput and high sensitivity. Not surprisingly, this method soon became mainstream.14)

Meanwhile, natural-level stable nitrogen isotopic composition of whole body of living organisms contributed significantly as a natural tracer to investigate the food chains in a particular ecosystem. The phenomenon that 15N of animal body is enriched through the food chain was first recognized in late 1960s by E. Wada then at the Tokyo University of Education.15) Subsequent studies revealed that the magnitude of 15N enrichment averages 3–4‰ per trophic step.16)19) With the knowledge that natural-abundance of carbon isotopic ratio (i.e., 13C/12C) that exhibits smaller elevation within food chains, carbon-nitrogen dual isotopic diagrams provide a useful tool for understanding and visualizing the trophic relationship among organisms in the ecosystem.20)23) Numerous applications of this approach have been made over the last several decades, and the field has come to be known as ‘isotope ecology’ or ‘stable isotope ecology’.23),24) Already in the relatively early phase of the research on the isotope ecology, the applications further expanded to other related fields such as paleoanthropology and archaeology.8),11),20),25),26) Considering the nitrogen in living organisms is derived mainly from amino acids, it would have been inevitable for isotope biogeochemists to focus on the 15N abundance of amino acids to further refine this approach.

As a successful integration of the above two research streams, McClelland and Montoya at Georgia Institute of Technology reported earlier this century that the 15N abundance of some amino acids like phenylalanine from marine zooplankton do not vary much with the trophic process, whereas other amino acids substantially enrich 15N relative to phytoplankton as their diet.27) The paper stimulated several groups including my group at JAMSTEC, where the use of nitrogen isotopic composition of amino acids began to be evaluated around 2005 by Y. Chikaraishi.28) It was followed by the proposition of an equation relating the nitrogen isotopic compositions of some amino acids and ‘trophic position’ based on the controlled feeding experiments of aquatic organisms.29),30) We viewed such formulation as a crucial step on the way toward various applications of this novel biogeochemical tool. At about the same time, teams at least from University of California Santa Cruz, University of Hawaii, and Bristol University shared this research.31)33) This approach has spread and received attentions from a wide range of research fields such as ecology, fisheries, anthropology, and others as a potential tool for solving long-standing issues in their fields.34) It should be noted that in the same time, other research groups applied the same method to various topics.35)39)

Below, I describe the principles of this approach, and reanalyze applications related to the isotope ecology during the last two decades, and further extend the discussion. It is not my intention to summarize extensive literature in this growing field, but rather to aim to present a useful introduction of this topic. I also note that I arbitrarily select and focus on some specific topics with which I have been deeply involved, rather than covering the full spectrum of isotope ecology applications of this approach. The readers may find other recent reviews also useful.34),39)

Principles

Amino acid metabolism and 15N.

In animals, amino acids are absorbed primarily as short peptides in the gut. After transport to the individual cell, the peptides are broken down into individual amino acids, joining the free amino acid pool (Fig. 1). The free amino acids in the pool are utilized for the synthesis of various proteins, but after a certain period of time, most proteins are broken down again to amino acids and returned to the pool. A portion of amino acids in the pool are constantly metabolized and undergo deamination (i.e., removal of amino group, -NH2).40),41)

Fig. 1.

A schematic view of major (but not all) metabolic processes of amino acids in an animal cell highlighted in nitrogen dynamics. Amino group of many amino acids are transferred to glutamic acid catalyzed by transaminase, and then the amino group of the glutamic acid is cleaved by glutamate dehydrogenase (GDH) and released as ammonia. The ammonia is excreted from the organisms as ammonia or urea. Note that both transaminases and GDH can catalyze the reverse reactions.

The deamination catalyzed by transaminase (or aminotransferase) is a common metabolic process for major amino acids to form the corresponding α-keto carboxylic acids (e.g., aspartic acid and α-oxaloacetate) as a half-reaction (Fig. 1). In contrast, the released amino group is transferred to α-ketoglutaric acid in the other half reaction to form glutamic acid that play a central role in amino acid metabolism. Some mechanistic details of this enzymatic reaction42) are known, and several studies have demonstrated that the entire reaction catalyzed by transaminase produces kinetically 15N-depleted glutamic acid.4),43) The glutamic acid, one of the major amino acids composing the pool, is oxidatively deaminated to form ammonia and α-ketoglutaric acid catalyzed by glutamate dehydrogenase (GDH). We still know little of the exact magnitude of isotopic fractionation associated with the GDH deamination,44) but the overall processes kinetically produce somewhat 15N-depleted ammonia and leave 15N-enriched amino acids in the pool. The ammonia is excreted from the organisms as it is in most aquatic organisms like fish, whereas in amphibians and mammals it is excreted after transformation to urea. Loss of 15N-depleted nitrogen leaves organisms and amino acids enriched in 15N.

A simplified model of trophic isotopic fractionation.

Figure 2 shows a simplified dynamic model of the above amino acid metabolism, especially considering nitrogen isotope ratio. The size of the amino acid pool in vivo is considerably smaller than the protein reservoir, with typical example commonly found in metabolic systems.45),46) The model assumes that amino acids brought into the pool from both the diet and protein reservoirs are mixed instantaneously. If zero-order kinetics is applied, the differential equations for the amino acid pool in the metabolic system are as follows:   

\begin{equation} \frac{\mathrm{d}[\text{AAs}]}{\mathrm{dt}} = f_{0} - f_{1}^{+} + f_{1}^{-} - f_{2}, \end{equation} [1]
  
\begin{align} \frac{\mathrm{d}[\text{AAs}]\delta_{d}}{\mathrm{dt}} &= f_{0}(\delta_{\text{d}} - \varepsilon_{0}) - f_{1}^{+}(\delta_{\text{AA}} - \varepsilon_{1}^{+}) \\ &\quad + f_{1}^{-}(\delta_{\text{p}} - \varepsilon_{1}^{-}) - f_{2}(\delta_{\text{AA}} - \varepsilon_{2}) \end{align} [2]
where f terms are fluxes, and (δ − ε) terms are isotopic values for fluxes. When considering a mature animal that the nitrogen mass balance and isotope balance are both achieved, above equations are described, respectively, as follows:   
\begin{equation} 0 = f_{0} - f_{1}^{+} + f_{1}^{-} - f_{2}, \end{equation} [3]
  
\begin{align} 0& = f_{0}(\delta_{\text{d}} - \varepsilon_{0}) - f_{1}^{+}(\delta_{\text{AA}} - \varepsilon_{1}^{+}) \\ &\quad + f_{1}^{-}(\delta_{\text{p}} - \varepsilon_{1}^{-}) - f_{2}(\delta_{\text{AA}} - \varepsilon_{2}). \end{align} [4]
In a mature animal, it is reasonably assumed that the nitrogen inflow (f0) and outflow (f2) are in equilibrium (i.e., zero growth rate). In this case, from Eq. [1], protein synthesis flux ($f_{1}^{ + }$) and degradation flux ($f_{1}^{ - }$) are also balanced as well.   
\begin{equation} f_{0} = f_{2}\ \textit{and}\ f_{1}^{+} = f_{1}^{-}(= f_{1}). \end{equation} [5]
Such a balance can be widely observed in healthy adult animals or human beings.41),47) The isotopic fractionations associated with the peptide hydrolysis catalyzed by both protease and peptidase to individual amino acids (ε0 and $\varepsilon_{1}^{ - }$) and the synthesis of protein ($\varepsilon_{1}^{ + }$) are assumed to be zero for simplification. To my knowledge, there is no literature that estimates at least $\varepsilon_{1}^{ + }$. In contrast, as some experimental studies have indicated, the amino acid metabolic step including transamination and deamination described above accompanies the significant isotopic fractionation (ε2). The isotope balance Eq. [4] can be arranged to:   
\begin{equation} \frac{f_{1}}{f_{0}} = \frac{\delta_{\text{d}} - \delta_{\text{AA}} + \varepsilon_{2}}{\delta_{\text{AA}} - \delta_{\text{p}}}. \end{equation} [6]
In the cases of investigated healthy higher animals including human, $\frac{f_{1}}{f_{0}}$ ratio for the maintenance of proteins generally falls in a range 2 to 3.41),45),48)50) A considerable amount of amino acids released by protein catabolism is re-utilized for the synthesis of new proteins. Because in most animals the size of protein pool is 2 to 3 orders of magnitude larger than that of the amino acid temporal pool, δp can be regarded as the bulk isotope ratio of the whole organism. As numerous observations have indicated that the bulk organic matter of the predator is somewhat enriched in 15N relative to that of the diet (here the magnitude is referred to ‘TDF’ or trophic discrimination factor), we apply δp = δd + TDF to Eq. [6]. Then, the magnitude of 15N enrichment along the trophic process is given by:   
\begin{equation} \mathit{TDF} = \varepsilon_{2} + \left( 1 + \frac{f_{1}}{f_{0}} \right)(\delta_{\text{p}} - \delta_{\text{AA}}). \end{equation} [7]
Because it takes some time for the incorporated amino acids to be metabolized, the snapshot results should be treated in caution.51)55) If δp − δAA is sufficiently small, Eq. [7] implies:   
\begin{equation} \mathit{TDF} \approx \varepsilon_{2}. \end{equation} [8]
The Eq. [8] provides a theoretical background that when the growth rate is zero, TDF primarily reflects the magnitude of isotopic fractionation associated with the amino acid metabolism such as deamination. However, the Eq. [7] also suggests that in the fasting state (f0 → 0), for example, TDF value fluctuates, which may also explain a part of variation of TDF value observed in natural environment.19),20),56),57) The relationship between TDF and ε2 needs to be carefully discussed.

Fig. 2.

The simple model for intake and metabolism of nitrogen by animals. ‘AA’ denotes amino acid pool in vivo. Nitrogen isotopic ratios of diet (d), temporal pool of amino acids (AA), protein reservoir (p), and excrement (ex) are denoted by δd, δAA, δp, and δex, respectively. Flows of nitrogen denoted by f and isotope effects by ε. Isotopic ratios of nitrogen transferred the indicated pathways are given by δd − ε0, $\delta_{\text{p}} - \varepsilon_{1}^{ - }$, $\delta_{\text{AA}} - \varepsilon_{1}^{ + }$, and δAA − ε2.

Unfortunately, only a limited number of studies have reported the isotope effects associated with the metabolism of individual amino acids, making it difficult to rigorously verify the above considerations. Nevertheless, several studies have demonstrated that the transamination such as the synthesis of glutamic acid from aspartic acid results in an average isotopic fractionation of ∼6‰.7),43) Although we do not know if this is representative value for isotopic fractionation of various transaminations and is a kinetic isotopic fractionation that does not necessarily have to match ε2 in the model, the TDF value strongly reflects the ε2 value, given that the amino acid pool is constantly being supplied by feeding and protein degradation. In practice, the TDF values range 2 to 5‰ (3.4‰ on average) in the various food webs and, $( 1 + \frac{f_{1}}{f_{0}} )$p − δAA) term in Eq. [7] is expected to vary negatively to some extent, probably −4‰ to −1‰. If the amino acid distribution of the prey is not significantly different from the predator (which warrants the constancy of ε2), the δex value should have a slightly (∼1 to 2‰) higher than δp, and both the $\frac{f_{1}}{f_{0}}$ and δex terms regulate the TDF value.

The same scheme can be applied to explain the isotopic ratios of individual amino acids, the subject of this paper. It has been known that amino acid composition among various organisms (even among organisms at different trophic levels) are relatively similar (Fig. 3). Thus, the $\frac{f_{1}}{f_{0}}$ ratios of individual amino acids may not differ significantly as a whole. However, when we focus on individual organs, even within a single organism, as in the eye lens discussed below, there is diversity in amino acid composition and metabolic rate (i.e., varying $\frac{f_{1}}{f_{0}}$ ratio). Therefore, the TDF values for individual amino acids potentially differ among organs. In any case, the magnitude of 15N-enrichment of individual amino acids in the protein pools is not constant, governed primarily by ε2 value and $( 1 + \frac{f_{1}}{f_{0}} )$p − δAA) term unique to each individual amino acid.

Fig. 3.

Hydrolyzable amino acid compositions as weight % of four types of aquatic organisms; a) natural phytoplankton (n = 4),166) b) natural zooplankton (n = 4),166) c) cultured bacteria (n = 3),166) and d) Atlantic salmon (n = 2).167) Glx and Asx indicate the sums of glutamic acid and glutamine, and aspartic acid and asparagine, respectively. n.d. denotes ‘no data’.

With another model, the dynamics of nitrogen isotope ratios in living organisms were carefully examined not only in mature stage (i.e., steady state), but also during growth and fasting stages (in section 4.7 of Ref. 23). Based on this study, it was demonstrated that δp increases during fasting (f0 < f2) and decreases during growth (f0 > f2). The magnitude of change in δp is also shown to be sensitive to the balance between f0 and f2. The fact that the magnitude of change in δp in ecological studies is consistent means that animal metablism in the food webs is regulated in highly reproducible manner by ingestion and excretion. Predation may control growth and reproductive patterns across animal taxa, forcing metabolic dynamics into a mode consistent with that seen in nitrogen isotope ratio.

Estimation of trophic position.

In the photoautotrophic cell as a start of the food chain, synthesis and degradation of 20 proteinaceous α-amino acids are balanced so that metabolism can take place normally. Under such a balance, it is one of the beauties of Nature or reflection of ‘isotopically ordered world’58),59) that 15N is distributed among amino acids according to certain rules, which was practically confirmed by many studies. For photoautotrophs in the aquatic environment (i.e., primary producers like algae and cyanobacteria), for example, glutamic acid is 3 to 4‰ enriched in 15N relative to phenylalanine.27),29),60),61) A similar regularity, albeit in different offset is observed in terrestrial photoautotrophs, in which glutamic acid is on average 8 to 9‰ depleted in 15N relative to the phenylalanine.62)64) Although the details of how such isotopic differences are to be explained are not yet known, these isotopic rules are principally independent of species or biota.

In terms of 15N enrichment magnitude along the trophic transfer, proteinaceous amino acids are practically classified as either ‘source’ or ‘trophic’ amino acids. The former group shows little 15N enrichment in trophic processes, whereas the latter substantially enriches 15N in the same process. Generally, phenylalanine and methionine and others such as lysine, tyrosine, and threonine are grouped into the source amino acids, and other most amino acids like glutamic acid, aspartic acid, alanine, valine, leucine, isoleucine and proline are grouped into the trophic amino acids. This grouping is, clearly, not determined simply by essential or nonessential designations that primarily reflect carbon relationships, but instead nitrogen metabolic pathways.65)

Much discussion ensued based on various experimental results, for accurately estimating the trophic position. The idea is that the offset in nitrogen isotope ratios between trophic and source amino acids is strongly related to the trophic position (TP), a measure of the position it occupies in a food web. This has been generalized as Fig. 4 and the following equation29),34);   

\begin{equation} \mathit{TP} = \frac{\delta^{15}\text{N}_{\text{trophic}} - \delta^{15}\text{N}_{\text{source}} - \beta}{\mathit{TDF}_{\text{trophic}} - \mathit{TDF}_{\text{source}}} + 1 \end{equation} [9]
where TDF and β represents, respectively the trophic discrimination factor and the δ15N offset between trophic and source amino acids from photoautotroph at the base of the food web described above. The weighted average of the TDFs of the 20 amino acids can be viewed as TDF as explained above. Trophic position is a quasi-continuous variable that can take decimal values, the basis of which is 1.0 for the primary producer and 2.0 for the primary consumers. The nitrogen isotopic ratio is expressed as conventional δ notation relative to the atmospheric N2 or AIR standard (15N/14N ≃ 0.003677)66): δ15N ≡ [(Rsample/RAIR) − 1] × 103, where R is 15N/14N ratio. Units for δ are parts per thousand or permil (‰).

Fig. 4.

Schematic diagram showing the relationship between the δ15N offset of trophic and source amino acids and trophic position (TP). Typically, the trophic amino acid used in TP calculations is glutamic acid (Glu) and source amino acid is phenylalanine (Phe). Trophic discrimination factor (TDF) and β are shown graphically. Especially in the aquatic system, there are two scenarios have been proposed, and extreme cases are shown here. Scenario I indicates the case that TDF of trophic amino acids is constant through the food chain (TDFGlu = 8.0‰, TDFPhe = 0.4‰),29) whereas scenario II indicates that TDF of trophic amino acids decreases in high TP range (TDFGlu = 6.4‰).39)

Inasmuch as nitrogen isotopic variation is produced as a result of kinetic isotope effects, TDF and β should be, in theory, expressed as isotope fractionation or ratio of isotopic rate constants (i.e., ε) rather than the ‘difference of δ15N values’. However, numerically these two terms are very close, and thus in most cases there will not make any meaningful difference in the interpretation at this stage.

Thus far, δ15N values of glutamic acid and phenylalanine have been empirically recognized to be most valuable for assessing TP of organisms. Based on both field observations and diet-controlled culture experiments, it was proposed that the TDFGlu is as large as 7 to 8‰ on average in particular for TP 1 to 2 and TP 2 to 3 (Fig. 4, Scenario I). Subsequent some studies confirmed the TDFs of trophic amino acids nearly constant between aquatic and terrestrial predation processes.67)69) In contrast, the other studies have suggested that the universality of the TDFs of trophic amino acids do not hold, but are apparently variable at least in aquatic food chains with some factors such as diet quality and mode of nitrogen excretion.39),59),70)75) Some studies have suggested that TDFGlu in high TP organisms may be lower than the above values (7–8‰).39),60),70)72) For example, (TDFGluTDFPhe) value was estimated to be ∼5.7‰ based on a comparison of δ15N and gut contents in the teleosts across a wide range of TP (Fig. 4, Scenario II).60) Although the solution to this problem, which is difficult to know, is not straightforward, it is not surprising that TDF potentially fluctuates, since cellular amino acids form a strongly regulated dynamic system that are constantly synthesized and metabolized. Even among organisms at different trophic levels, amino acid compositions are relatively similar (Fig. 3), but minor compositional changes in AAs may be related to the TDF variability.76),77) In practice, however, if the TDFGlu is as small as ∼4‰ or lower as suggested by some studies,39),69),70) organisms with extremely high TP such as TP = 6 or higher need to be assumed to explain the results presented in the following sections. In a ‘healthy population’ within a single species, the balance of amino acid metabolism would fall within a certain level, and thus the TDF values for the amino acids would also fall in a certain narrow range. If the metabolism fluctuates due to external or internal factors, the TDF is also expected to fluctuate.

In contrast to glutamic acid, experimentally determined TDF of phenylalanine as a source amino acid is close to zero. Phenylalanine is hydroxylated to form tyrosine as the first step in its metabolism, which does not accompany the cleavage of the carbon-nitrogen covalent bond.28),77) Reasonably, TDFPhe is quite small. Coupled with the fact that the phenylalanine is an essential amino acid for all animals, δ15NPhe values even in higher-order predators are expected to be close to those of autotrophs such as algae or plants, which form the basis of the ecosystem. In the following discussion, to simplify the discussion, I assume the TDF values of glutamic acid and phenylalanine are 8.0‰ and 0.4‰, respectively, if not mentioned.

The essence of above equation can be graphically visualized in the δ15NPhe15NGlu diagram (Fig. 5). Organisms belonging to a single trophic chain can be, in theory, plotted on a line roughly parallel to the vertical axis in this diagram, even if the TDF of glutamic acid varies along the food chain. In contrast, organisms with equivalent TP are plotted on a line parallel to the broken lines drawn from the lower left to the upper right in the diagram. Alike δ13C-δ15N diagram that has been applied classically, the horizontal axis of the δ15NPhe15NGlu diagram indicates the source factor, whereas the vertical axis indicates the trophic factor. However, compared to the δ13C-δ15N diagram, an essential difference exists in this amino acid isotope approach. This recently developed approach is able to separate two major δ15N-controlling factors that originate from completely distinctive processes; environmental δ15N-baseline and trophic effect.78),79)

Fig. 5.

δ15NPhe15NGlu diagram for interpreting the food chains of aquatic and terrestrial food webs. Aquatic food chains are shown in blue and a terrestrial food chain in green. The diagonal dashed lines indicate trophoclines, or iso-TP lines. In this diagram, TP = (δ15NGlu − δ15NPhe − 3.4)/7.6 + 1 is applied to the aquatic food chain, and TP = (δ15NGlu − δ15NPhe − 8.4)/7.6 + 1 to the terrestrial food chain.29),62) Three hypothetical food chains are shown in the diagram; Aquatic food chains 1 and 2 start at (δ15NPhe, δ15NGlu) = (4.0‰, 7.4‰) and (−3.0‰, 0.4‰), respectively, whereas terrestrial food chain starts at (δ15NPhe, δ15NGlu) = (11.0‰, 2.6‰). Hypothetical organism D feeds on TP 2 organism A belonging to aquatic food chain 1, and TP 3 organism B belonging to aquatic food chain 2, with a ratio of 1:1. Hypothetical organism E feeds on TP 2 organism A belonging to aquatic food chain 1, and TP 2 organism C belonging to terrestrial food chain, with a ratio of 7:3.

A few caveats of the approach.

A few caveats should be noted in this approach. First and foremost, ‘glutamic acid’ referred throughout this paper is a mixture of glutamic acid and glutamine, because glutamine is chemically deamidated to form glutamic acid in the process of amino acid extraction, or hydrolysis of peptide bonds. Currently, no efficient protocol is known to liberate glutamine from the peptides. These two amino acids are closely related metabolically, so that the nitrogen isotope ratios of the amino groups bonded to the α-carbons of glutamic acid and glutamine would be expected to have close values.

As a similar issue, what we call the nitrogen isotopic ratio of “amino acids” except glycine, is the weighted average of those of amino acid enantiomers.80)83) In the natural environment, d-amino acids are not only produced by racemization of l-amino acids, but are also abundantly synthesized by prokaryotes as peptidoglycan, a cell wall component. For alanine and glutamic acid, it has been demonstrated that d-enantiomers are up to a few permil depleted in 15N relative to the corresponding l-enantiomers in both cultured and natural organisms.84)86) If such case studies are widely true for other amino acids and in various environments, then the contribution of d-amino acids from microbes would lower the δ15N values of the amino acids to some extent. Lack of evidence hampers further discussion, but for the precise interpretation of isotopic records that strongly involved in microbial processes, this problem must be carefully evaluated in the future.

Secondly, the nitrogen isotopic composition of organic tissue reflects only that of the diet only absorbed and not all the contents found in the gut of organisms. Even though organic substances with high carbon-to-nitrogen ratio such as lipids and carbohydrates are digested and absorbed, they lack nitrogen and must have little effect on the TP defined above. The TP estimated through the nitrogen isotope ratio should not be strictly equivalent to the TP viewed in terms of carbon (or energy). Flows of carbon and nitrogen are decoupled not only in environmental processes but also in metabolic processes. This fact is recognized, for example, by considering the assimilations of carbon and nitrogen by autotrophs through separate mechanisms; For carbon, ambient CO2 is fixed through enzymatic reaction catalyzed by rubisco, whereas for nitrogen, inorganic nitrogen such as NO3 and ammonia is fixed through enzymatic reaction catalyzed by either glutamine synthase or GDH. Furthermore, N2 fixation and denitrification are highly important processes for the nitrogen cycle, but are not directly related to the carbon cycle. To be exact, the TP discussed here should be defined as ‘nitrogen-based TP’, and ‘carbon-based TP’ requires separate evidence and considerations. For example, in the case of humans, the main source of nitrogen is from land animals such as cattle and fish, so the nitrogen-based TP is 3 or higher. On the other hand, carbon is mainly consumed from grains such as rice and wheat, so the carbon-based TP is about 2.

Thirdly, care should be taken when selecting samples. It has been long known that individual organs and tissues even from a single organism have distinctive bulk isotopic signatures.26),54),87),88) Such isotopic variabilities among organs and tissues even in a single organism may be ascribed mainly to the difference in metabolic rate (i.e., $\frac{f_{1}}{f_{0}}$ ratio in Eq. [7]) among organs89); muscle is constantly metabolized (i.e., produced and destroyed constantly), whereas organs such as bones, the eye lens and brain are in principal mostly accumulating amino acids and are metabolized much slowly. As described in the preceding section, even in the case of amino acids, differences in metabolic rates among organs and tissues result in differences in isotopic ratios.

Grazing food chain in the aquatic environment

A survey of the aquatic food web, especially in a confined space such as a lake, provides a good example to start with. The lake and its watershed form materially semi-enclosed system, and most organic matter produced there is consumed in the system. In other words, the organisms that inhabit that closed space are materially connected with prey-predator relationship. We have been studying three different types of lakes, Lake Superior, Lake Baikal, and Lake Biwa. Partial results were shown in elsewhere54) and studies of Lake Biwa and Lake Suwa are still on going.

Lake Superior.

Figure 6 shows the amino acid isotopic results from Lake Superior, one of the Great Lakes in northern North America, studies in collaboration with a research group at the University of Minnesota.90) We examined 8 animal species out of ∼80 species inhabiting this oligotrophic lake, which provided intriguing aspects as well as compared with the knowledge based on classical methods such as gut content analysis. Mean TP of Daphnia is ∼2.0, indicating this small zooplankton to be a typical herbivorous and feed only on photoautotrophs. The TP of small-sized copepod Diporeia is ∼2.4 on average, suggesting omnivorous with nearly half of its diet may depend on primary producers. The result of opossum shrimp Mysis (TP = 2.6 to 2.9) indicate somewhat more carnivorous than Diporeia, consistent with the previous observations.91)93) The copepod Limnocalanus macrurus registered at carnivorous TP 3.1 on average. In this way, herbivores, omnivores, and carnivores are clearly discriminated, and the degree of their carnivory would be estimated with some precision. The accuracy of the dietary evaluation is, in theory, dependent on the stability of TDF of glutamic acid and β value, but in application, it may also be dependent on the microbial activity, which will be discussed later.

Fig. 6.

Food web analysis of Lake Superior based on the nitrogen isotopic ratios of glutamic acid and phenylalanine. a) δ15NPhe15NGlu diagram with trophoclines, and b) trophic position (TP) for 8 species collected from three sites (Bap, WM, and Ont) in Lake Superior. Of three sites, Bap and Ont located coastal, and WM offshore; Species analyzed are; zooplankton, Daphnia, Diporeia, opossum shrimp Mysis, copepod (Limnocalanus macrurus), lake herring (Coregonus artedi), kiyi (Coregonus kiyi), deep water sculpin (Myoxocephalus thompsonii), and lake trout (Salvelinus namaycush). The two axes for estimating TP are shown at the top and bottom of diagram b. The top axis corresponds to Scenario I shown in Fig. 4, with TDFGlu set to 8.0‰,29) and the bottom axis corresponds to Scenario II, with TDFGlu set to 6.4‰.39) In both cases, TDFPhe was set to 0.4‰. The original data was reported and discussed in elsewhere.90)

The estimated TPs of two small planktivorous coregonids, lake herring (Coregonus artedi) and kiyi (Coregonus kiyi) are both slightly higher than 3.0, suggesting them potential resource competitors. However, they appear to prefer different diets and are thought to be well segregated. The TP of benthic deep water sculpin (Myoxocephalus thompsonii) and lake trout (Salvelinus namaycush) are 3.8 and 4.0 on average, respectively, nearly one unit above the coregonids. The lake trout has been thought to be a top predator of the lake, whose stomach content analysis have suggested heavy reliance on other fish species, particularly deep water sculpin and coregonids.93),94) The amino acid isotopic evidences support the interpretation of bulk isotopic results from Lake Ontario.92) Collectively, these observations revealed that the food chain in Lake Superior basically has four trophic steps. Accurate estimation of the TP would be helpful for understanding the food chain in detail, which potentially track persistent organic pollutants (POPs) and heavy metals intensively discussed in Great Lakes.95)

δ15NPhe values of above 8 species widely vary over 7‰ (−5‰ to 2‰), especially at coastal Bap and Ont sites (Fig. 6). It would reflect the seasonal changes in δ15N at the base of the food web associated with water dynamics, in addition to the existence of discrete ecosystems in coastal and offshore areas. Moreover, the contribution of allochthonous resources to the coastal areas may have an impact on the δ15NPhe. Theoretically, such variabilities of δ15NPhe can be particularly observed in low-TP organisms.96),97) The food chain converges toward higher-order, long-lived predators, leading that the δ15NPhe values of the predators also tend to focus in a certain range.

Lake Baikal.

An intriguing exercise is to compare the food web of Lake Superior (Fig. 6) with that of Lake Baikal (Fig. 7). Lake Baikal is the largest and deepest (∼1,642 m at maximum depth) freshwater lake on Earth and located in a similar latitudinal range. In this lake, many endemic species inhabit both coastal and offshore surface water, forming the characteristic ecosystem. In the offshore water, the ecological diversity is remarkably low with a limited number of species.98) The TP of mesozooplankton Episula baikalensis even collected in different years is 2.0 on average, suggesting them to depend entirely on photoautotrophs like Daphnia in Lake Superior. In contrast, microzooplankton gammarid Macrohectopus branikii occupies TP of ∼2.3 on average, suggesting omnivorous with about two thirds of its diet depend on primary producers. Benthic sculpin (Batrachocottus multiradiatus) registered clearly omnivorous at TP ∼2.7, which is a strong contrast with the deep water (DW) sculpin in Lake Superior whose TP obviously shows it carnivore. Commercially important pelagic whitefish omul (Coregonus autumnalis migratorius) has a mean TP of 3.6, which is consistent with behavior observations that its staple diet is the gammarid amphipod.98) Planktic sculpin (Comephorus dybowski) are higher order predators than omul, with an average TP of 4.3. Top predator of the lake is Baikal seal (Pusa sibirica) with TP 5.1 on average. In offshore ecosystem of Lake Baikal, the amino acid isotopic results showed that the TP of the top predator is about one unit longer than that estimated by δ15Nbulk (∼4.0),99) which might be the result of more variable TDF of bulk tissue.

Fig. 7.

Food web analysis of Lake Baikal based on the nitrogen isotopic ratios of glutamic acid and phenylalanine. a) δ15NPhe15NGlu diagram with trophoclines, and b) trophic position (TP) for 7 species collected from offshore in Lake Baikal; phytoplankton (Melosira), zooplankton (Episula baikalensis), amphipod (Macrohectopus branikii), benthic sculpin (Batrachocottus multiradiatus), omul fish (Coregonus autumnalis migratorius), planktic sculpin (Comephorus dybowski), and seal (Pusa sibirica). The two axes for estimating TP are shown at the top and bottom of diagram b. The top axis corresponds to Scenario I shown in Fig. 4, with TDFGlu set to 8.0‰,29) and the bottom axis corresponds to Scenario II, with TDFGlu set to 6.4‰.39) In both cases, TDFPhe was set to 0.4‰.

The food chain length (FCL) that corresponds to the TP of the top predator, is a central characteristic of ecological communities. Based on the amino acid isotopic results under the assumption of constant TDF through the food chain, the FCL are ∼4.0 for Lake Superior, ∼5.1 for Lake Baikal, ∼3.9 for Lake Biwa (Ishikawa, N. et al., unpublished results), and ∼3.0 for Lake Suwa (Urai, A. et al., unpublished results). In Lake Superior, the result appears to be consistent or somewhat shorter than those of previous studies.100),101) If the TDF of glutamic acid were smaller than that assumed here (i.e., 8.0‰), the estimated FCL would be systematically longer than those described above. Any FCL estimates have specific assumptions inherent to a particular approach and should be carefully assessed. Nevertheless, the estimated FCL in the largest Lake Baikal (water volume: 23,600 km3) was significantly longer than the other three lakes, whereas it is by far the lowest in the smallest Lake Suwa (water volume: 0.063 km3), which does not contradict the hypothesis that ecosystem size determines FCL in lakes.101) The inclusion of freshwater mammal seal in the Lake Baikal ecosystem is interpreted to push up the FCL nearly one unit from planktic sculpin.

Given all species comprising the ecosystem have the same opportunity to be preyed on their predators, ecosystem generates a triangle with the top-predator at the apex in the δ15NPhe15NGlu diagram. It is theoretically anticipated that the δ15NPhe of the high-TP organisms converges to a narrow range corresponding to averaged δ15NPhe of biomass. The overall distribution of δ15NPhe is −0.6 ± 1.9‰ (1σ, n = 30) for Lake Superior and 2.6 ± 1.9‰ (1σ, n = 53) for Lake Baikal. The ∼3‰ difference is attributed to the 15N-enriched nitrogen pool in Lake Baikal relative to Lake Superior due to geochemical processes such as ammonia evaporation or denitrification in the watershed.102) Together with the mean δ15NPhe value of the top predator, the ecosystem triangle in the δ15NPhe15NGlu diagram can be strictly defined by three components, including the FCL and the variation in δ15NPhe of primary producers.

Feeding habit variabilities among closely related species.

How do dietary changes relate to other biological properties such as habitat, behavior, and even evolution? It often happens that similar organisms occupy different trophic niche, as seen in examples of deep-inhabiting sculpin in Lake Superior and Lake Baikal. It might be ascribed simply to ecological reasons such as the local extinction of predators or lack of food associated with the environmental change. Although sufficient evidence has not yet been accumulated, it appears that closely related species from different ecosystems often occupy distinct trophic niches. Here we compare 5 species of righteye flounder (albeit at different genus levels) that inhabit the soft muddy areas of the seafloor with similar behavioral patterns, collected at water depth of 220–501 m off the Pacific coast of northeast Japan.95),103) Trophic positions of these species are roughly consistent among samples in the single species, but vary strikingly from 3.1 and 4.1 among species (Fig. 8). Of these, the lowest TP Rikuzen sole, Dexistes rikuzenius, may feed on detritus sinking from the overlying water column, whereas the highest TP flathead flounder, Hippoglossoides dubius, is clearly carnivorous. Such TP variations among righteye flounder species do not appear to be related to the size of specimens within a single species. Despite similar habitats and behavior, these righteye flounders feed on markedly different diet. In nature, feeding habit is not conservative and may be strongly regulated by several factors, including availability of the prey and competition with other species. Differences in species mean size might be related to food quality issues. Also, the dietary change is a potential driving force of sympatric speciation. Accumulation of the results would provide an insight into such an intriguing topic.

Fig. 8.

Summary of trophic position of 5 species of righteye flounder along the total length of the specimens; They are Dexistes rikuzenus (rikuzen sole), Clidoderma asperrimum (roughscale sole), Hippoglossoides dubius (flathead flounder), Glyptocephalus stelleri (blackfin flounder), and Microstomus achne (slime flounder). All these specimens were collected in the water depth range of 220–501 m from off Pacific coast of northeast Japan.95)

Ecological consequences of the Anthropocene.

As partly discussed above, the δ15NPhe value is informative for reconstructing the δ15N of the nitrogen pool in the system. It provides clues for understanding environmental issues such as the ecological consequence of eutrophication. In eutrophic lakes around the world, university laboratories and local museums often preserve diverse collections of formalin-fixed biological archives. Such biological archives potentially provide valuable samples for reconstructing the recent changes in lake environment and understanding its ecological consequence. The fixation by formalin (the aqueous solution of formaldehyde) has been one the most common fixatives used for the preservation of organisms. In the aqueous solution formaldehyde forms methanediol or methylene glycol, which is empirically known to react with the residues of some proteinaceous amino acids, including serine, threonine, glutamine, lysine, arginine, cysteine, and tyrosine to form reactive hydroxymethyl groups.104) The chemical process neither adds nor removes nitrogen from the protein (but adds carbon). Indeed, it was experimentally proven that the biological specimens fixed by formalin encode reliable nitrogen isotopic information in individual amino acids.105),106)

This approach was applied to formalin-fixed specimens of gobiid fish (Isaza, Gymnogobius isaza, year 2nd), which were collected from Lake Biwa throughout the 20th century and stored in Center for Ecological Research, Kyoto University (Fig. 9). Based on the δ15Nbulk records of Isaza specimens as well as sediment collected from the lake, a gradual 15N enriching trend had been recognized, especially during the latter half of the 20th century when human activities seriously perturbed the lake environment.107),108) For this temporal 15N enriching trend, two conflicting interpretations had been proposed; change in trophic niche (i.e., TP elevation) of Isaza along with the eutrophication, and elevation of δ15N of nitrogen pool. Given the elasticity of ecosystems, even in a single species TP would have a certain fluctuation depending on the environment. The results clearly demonstrated that the temporal change in the nitrogen isotopic composition was mostly attributable to that of nitrogen pool in the system.105) Through the 20th century, the estimated TP of Isaza remained nearly constant (∼3.2). In contrast, δ15NPhe at the base of the food web (i.e., δ15Nbase = δ15NPhe − 0.4 × (TP − 1)) systematically increased from ∼5‰ in early 20th century to ∼9‰ in the latest 20th century (Fig. 9). Such a ∼4‰ elevation of δ15Nbase during the 20th century has been explained by the reflection of enhanced denitrification along with the eutrophication of the lake.109),110) Despite the deterioration of water quality and the almost simultaneous invasions of exotic species such as black bass and bluegill, trophic function of Isaza has remained largely unchanged. More formalin-fixed specimens of various species are currently being investigated.

Fig. 9.

A case study in Lake Biwa on the effects of lake eutrophication in the 20th century on the feeding habits of fish and the reconstruction of nitrogen cycles in the lake. a) Concentrations of nitrate (black) and phosphate (red) observed in the hypolimnetic water in the north basin of Lake Biwa,168) b) estimated trophic position of formalin-fixed gobiid fish Isaza (Gymnogobius isaza, Tanaka), and c) estimated δ15Nbase, or δ15NPhe value at the base of the food web (i.e., δ15Nbase = δ15NPhe − 0.4 × (TP − 1)).

Numerous lakes, ponds, and rivers in the world that have experienced or are presently experiencing the eutrophication like Lake Biwa. Studies of the formalin-fixed specimens archived around the world will provide unique but crucial information for discussing the ecological consequences of the Anthropocene.

Nitrogen resource apportionment between aquatic and terrestrial organic matters.

Besides glutamic acid and phenylalanine, the δ15N value of methionine (δ15NMet), another source amino acid has been investigated its utility for the ecosystem study. As with phenylalanine, the metabolic processes of methionine start not with deamination, but with the methyl group removal to form homocysteine. Consequently, it is reasonably expected that the catabolism of methionine does not change its nitrogen isotope ratio. Taking advantage of the stability of the nitrogen isotope ratios of source amino acids along the food chain (∼0.4‰/TP for phenylalanine and ∼0.5‰/TP for methionine),29) an attempt was made to estimate the respective contributions of terrestrial and aquatic resources of nitrogen based on the δ15N offset between two source amino acids, phenylalanine and methionine (Fig. 10).111) By subtracting the nitrogen isotope ratios of the two source amino acids, the variation in the δ15N baseline can be canceled. The utility and accuracy of this estimate strongly depends on the robustness of the δ15N relationship between these two amino acids, although the relationship needs marginal corrections by TP. Our current evidence indicated as large as ∼12‰ isotopic offset (= δ15NMet − δ15NPhe) between aquatic and terrestrial ecosystems; the former is ∼−5‰, whereas the latter ∼−17‰ (Fig. 10). This could be successfully used to accurately assess TP in animals dependent on both terrestrial and aquatic resources by correcting for β values.112),113)

Fig. 10.

Nitrogen resource apportionment (aquatic vs. terrestrial) based on nitrogen isotopic difference between phenylalanine and methionine (i.e., δ15NMet − δ15NPhe): From top to bottom, aquatic organisms, terrestrial organisms, and stream organisms that feed on both.111) At the bottom, the scale of terrestrial contribution (%) is shown.

The results of applying this approach to 8 species of animals including fish and insect larvae from a headwater stream food web are also shown in the bottom of Fig. 10, which are not contradicting with current knowledge on feeding behavior of these species.111) Of these, larvie of mayfly (Heptageniidae spp.) exhibits strong reliance (∼90%) on aquatic nitrogen resources. Nitrogen for the three fish species of goby (Rhindogobius flumineus), trout (Oncorhynchus masou ishikawae), and minnow (Rhynnchocypris oxycephalus jouyi) in the stream was found to be obtained roughly half from both the terrestrial and aquatic ecosystems. By focusing on protein resources, it allows to delve deeply into a classical topic for biogeochemistry, the interaction between distinctive terrestrial and aquatic ecosystems.114),115) Research fields like paleoanthropology would benefit from the uniqueness of this approach. For example, the consumption of aquatic resources by inland hunter-gatherer communities has been thought to be strongly related to environmental adaptation, but it is challenging to assess accurately with classical approaches. This novel approach, if further confirmed, has the potential to contribute to even cultural issues.

Microbial processes and endosymbiosis

The outcomes of recent biogeochemistry research left little doubt about the fundamental role of microbes in the environment.116),117) It is obviously necessary to know how such microbial processes affect the nitrogen isotope ratios of amino acids presented here. Interpretation of TP estimated by this approach often ignores microbial processes and may be somewhat less quantitative. When applied to natural ecosystems where degradation dominates, such as the deep sea, the effects of microbial activity must be carefully taken into account. In food web studies, the ‘detritus food chain’ has often been overlooked in comparison with ‘grazing food chain’, but these two are functionally interacted and not neatly separated.

Based on the laboratory culture experiments of microbial species including bacteria, archaea, and eukarya (fungi) in synthetic media aerobically,118) the TP of microbes grown with free amino acid mixture (i.e., casamino acids) elevates one unit from their ‘diet’ (Fig. 11). In contrast, for the microbes grown with ammonium as the sole nitrogen source, δ15N of amino acid pattern including δ15NPhe15NGlu relationship is very similar to that of algae with TP of 1.118),119) For heterotrophic microbes including bacteria and fungi in the terrestrial environment, the amino acid approach was also demonstrated to be effective.120) These observations may reflect the fact that amino acid metabolism in microorganisms is essentially the same as in other organisms, including photoautotrophs. An ecological perspective is that the apparent TP of detritivores that consume microbe-rich organic matter would be not simply one unit higher that of the macroscopic ‘prey’, but even higher. The TP values found in deep-sea organisms may represent an example of microbial links in the natural food chain. The observed TP of deep sea (200–550 m water depth) fish off Japan coast starts at TP ∼3.2.95) On the other hand, if a mixture of phytoplankton (TP 1) and zooplankton (TP 2) biomasses directly sinks to the deep ocean and is fed by the deep-sea fish, their TP should theoretically start at ∼2.5. The difference between these observational and theoretical values is explained potentially by the microbial links in the water column.

Fig. 11.

The nitrogen isotopic difference between glutamic acid and phenylalanine (i.e., δ15NGlu − δ15NPhe) from various microbes cultured in the laboratory.118) Culture experiments were conducted by both grown with ammonia (closed symbols) and with casamino acids, a mixture of amino acids (open symbols). Because δ15NGlu and δ15NPhe values of the casamino acids are 10.2‰ and 5.8‰, respectively, the (δ15NGlu − δ15NPhe) value of microbes is expected to be 12.0‰ (dotted line) when they assimilate casamino acids. Broken line indicates the expected value (δ15NGlu − δ15NPhe = 3.4‰) when the microbes grow using free NH3. Cultured microbes are selected from major three domains of organisms. Bacteraia: Escherichia coli, Archaea: Methanothermobacter thermautotrophicus, Sulfolobus tokodaii, Halobacterium salinarum, and Eukarya: Saccharomyces cerevisiae. Details of cultured conditions were described in elsewhere.118)

An interesting application of this knowledge is host-symbiont nutritional relationship in endosymbiosis observed widely in natural environment. So far, the nutritional relationship between microbial symbionts and hosts such as foraminifera, sea slugs, and deep-sea mussel have been investigated (Tame, A. et al., unpublished results).121)124) In all these cases, in favorable environments, microbes grown with inorganic nitrogen (i.e., TP = 1) as symbionts supply nutrition (i.e., amino acids) they produced as photosynthates to the hosts. The host is supported by the amino acids produced by the symbionts as well as by its own food harvesting, and its trophic position takes values between 1 and 2. In contrast, in unfavorable environment such as shortage of food, the hosts directly feed on the symbionts, and TP of the hosts substantially increases up to 2. These studies clearly demonstrated that the hosts switch the nitrogen source between being supplied by the symbiont and consuming the symbiont, dependent on the environment. Amino acid isotopic approach provides insight into the strategies of the hosts for survival in a fluctuating environment.

Tracking migration routes of marine fish

Since many marine fishes that are important as fishery resources are migratory, knowing the details of their migration is essential for fisheries management and marine conservation. Isoscape, which maps isotopes, is a relatively recent approach for potentially contributing to this issue.125)128) A utility of δ15NPhe that reflects the base of the food chain is its application to marine biological samples through isoscape for tracking the migration of fish in time and space.129) This approach is clearly demonstrated to be helpful for reconstructing the long-distance migration route of chum salmon (Oncorhynchus keta) in the northern North Pacific.130) The δ15NPhe as well as δ15Nbulk recorded in the growth ring of vertebral centra131),132) of chum salmon collected from a Japanese river where they migrated to spawn was integrated with δ15N isoscape in the northern North Pacific through a statistic model. The reconstructed migration route of chum salmon reproduced the known migration pattern133),134) with high spatial and temporal resolutions relative to previous studies with classical approaches (Fig. 12). After roaming the northwest Pacific Ocean for a year or two, the chum salmon gather at the later stage of their growth on the continental shelf in the eastern Bering Sea where they mature supported by high productivity. Coupling with the isoscape, amino acid isotope approach acquires versatility, potentially being applicable to most fish for tracking of their migration routes. Simultaneously, its drawback is that the amino acids are only minor components in the vertebral centra, limiting the sample amount, particularly for juvenile period and that after growth slows down due to sexual maturity.

Fig. 12.

The migration areas of chum salmon (Oncorhynchus keta) along growth stages (A–D) were estimated by matching δ15NPhe and δ15Nbulk recorded in centrifugal vertebral growth rings to the North Pacific δ15N isoscapes using a state-space model.130) The color gradients indicate presence probability of chum salmon in different growth stages. The δ15N isoscape was constructed with the analytical results of long-lived copepod.

Recently, eye lens has demonstrated as a potential archive for circumventing this issue.135)137) The eye lens composed mainly of proteins (i.e., amino acids), supply a plenty of sample and more suitable for this kind of study compared to the vertebral centra or otolith, even though its growth rate also slows down in the later life stage. The eye lens forms successive layers of cortical laminae, which can be manually separated layer by layer. We have successfully separated as many as 34 layers from a single eye lens of ∼5 mm i.d. from chub mackerel (Scomber japonicus).137) The amino acid isotopic results of some of these samples showed lifelong changes in the δ15N of baseline. In particular, early life histories of the fish is densely recorded, allowing to reconstruct the diet for juvenile, which is often important for conservation and aquaculture of the species.138)

The amino acids that form the eye lenses of all animals are little metabolized after accumulation. Our knowledge on the growth rate of the eyeball and metabolism of amino acids forming the eye lens still needs to be accumulated, but this method is potentially applicable not only to the marine organisms but also to the terrestrial animals and may be useful for reconstructing the detailed behavior analysis of a variety of organisms in the future.

Applications to paleodietary study

Another exciting application of the compound-specific isotopic approach is paleodiet reconstruction using mammalian bone collagen, a triple helix made of amino acid sequence.8),139)141) Since the early phase of the development of our investigation of this approach, various fossil bones of vertebrates from archaeological sites have been studied with the collaboration mainly with a research group at The University of Tokyo.

Once put aside its applications, the utility of the amino acid isotope approach to the archaeological studies can be partially validated by bone collagen samples from herbivorous animals which are reasonably expected to have TP of 2.0. Figure 13 summarizes the all available analytical results of bone collagens of various herbivores collected from archaeological sites in Japan, Syria, and France (3,000 to 8,000 years BP).142)149) Most plots (n = 51) fall along the trophocline of TP 2 (average: 2.0, standard deviation: 0.1), confirming that the method is effective for reconstructing at least herbivore TP using fossil bone collagen. Exceptions are results of deer and cattle from Kitakogane shell midden and Moyoro site (n = 6), both from coastal area of northern Japan with cold and snowy climate.150) These samples exhibit somewhat high TP ∼2.4 with relatively 15N-depleted phenylalanine. From the perspective of a top-down approach,151) these results are potentially explained by either microbial contribution in their diet or the consumption of marine-derived resources (δ15NPhe ranges 1 to 2‰ in these sites)142) such as seaweeds, as sometimes observed in cold regions when available food is limited for terrestrial herbivores in winter.152)

Fig. 13.

δ15NPhe15NGlu diagram for fossil bone collagen from herbivores collected from archaeological sites. Trophoclines in the diagram is drawn with 7.6‰ for (TDFGluTDFPhe) value and −8.4‰ for β value. The archaeological sites include Kitakogane shell midden,142) Moyoro site,143) Kitamura and Tochibara sites,154) Tell-el-Kerkh,146),147) and Noyen-sur-Sein.149) Plots enclosed by a broken circle are samples from Kitakogane shell midden and Moyoro site.

In the bulk isotope approach applied for decades, the ‘collagen’ samples with carbon-to-nitrogen (C/N) ratio outside of a certain range (e.g., 2.5–3.1 as weight ratio)153) have been interpreted to be ‘diagenetically altered’, and thus intentionally removed from the discussion. In the amino acid isotope approach, however, such a strict operational distinction may not be required. From both physicochemical and empirical perspectives, hydration rates (i.e., cleavage rate of peptide bond) in diagenetic processes of bone collagen are different among individual amino acids, modifying amino acid composition and thus δ15N of bulk collagen as diagenesis proceeds. Since δ15N values of individual amino acids widely vary especially for high-TP organisms, amino acid compositional change through diagenesis should substantially alter the δ15Nbulk record. The increase in the C/N ratio of the bulk collagen in the diagenetic samples may be primarily due to the removal of amino groups from the N-terminus of the peptide and amino acid residues (i.e., glutamine to glutamic acid, asparagine to aspartic acid, and arginine to glutamic acid).77) Also, the effect of contamination of exogenous amino acids cannot be entirely ruled out. In the future, some criteria for assessing the nature of amino acids from fossil bone would be required.

Because the amino acid approach focuses on nitrogen composing the protein backbone which is strongly connected with carbons for both sides, in theory, it is not exchangeable with ambient nitrogen such as ammonium during diagenesis. Unless exogenous amino acids are added secondarily, changes in amino acid composition during diagenesis do not alter the δ15N record of individual amino acids. This consideration is empirically supported by data shown in Fig. 14, which is a replot of the data previously reported.154) Both δ15NGlu and δ15NPhe values of collagen from human fossil bone, and hence estimated TP little vary (TP range: 2.6–2.8) even in the collagens with bulk isotope ratios of ∼4‰ variation and with C/N weight ratio up to ∼9. The observation does not necessarily guarantee the reliability of this approach, but it certainly warrants further study, including consistency with experimental observations that chemical hydrolysis of peptides at high temperatures leads to fractionation of nitrogen isotopes.155),156)

Fig. 14.

Evaluation of the diagenetic effects on nitrogen isotope ratios of amino acids. The samples shown here were collagens extracted from human bone samples collected at inland Tochibara and Kitamura sites, Japan (original data from Ref. 154). In the upper panel, estimated TP is plotted as a function of C/N (weight) ratio of the fossil collagen samples extracted from human bones. In the lower panel, nitrogen isotopic compositions of both glutamic acid and phenylalanine are plotted against C/N ratio of the collagen. The grey zone indicates the normal acceptable range of C/N weight ratio (2.5–3.1) in the bulk isotopic analysis.154)

Logically, the δ15N value of bulk collagen would strongly reflect that of glycine, which accounts for about one-third of amino acids in bone collagen. However, glycine metabolic pathways are multiple and complex, so the δ15N of glycine, and thus of bulk collagen is not explicitly related only to trophic processes,146) but could be regulated to some extent by the balance of the multiple metabolic pathways related to glycine.

The estimated TP of paleo-human from inland of Japan and Europe who are considered to little feed on aquatic resources range 2.2 to 2.9. Among them, Neanderthals from Spy Cave showed the higher TP compared to Homo sapiens, ranging from 2.7 to 2.9,145) supporting the previous thoughts of being primarily carnivorous.157)

Bridging bulk and amino acid approaches

In this section, I will consider how to connect the traditional bulk isotope data with the amino acid approach. Nitrogen isotopic ratio of bulk organism can be expressed as the weighted mean of various nitrogenous compounds composing the organism;   

\begin{equation} \delta_{\text{bulk}} = g_{1}\delta_{1} + g_{2}\delta_{2} + g_{3}\delta_{3} + \ldots + g_{\text{n}}\delta_{\text{n}}, \end{equation} [10]
where g represents the nitrogen fraction of the individual nitrogenous compounds in the total nitrogen, satisfying g1 + g2 + g3 + … + gn = 1. The subscripts 1 through n denote the respective nitrogenous compounds. The same can be applied to the TDF;   
\begin{align} \mathit{TDF}_{\text{bulk}} &= g_{1}\,\mathit{TDF}_{1} + g_{2}\,\mathit{TDF}_{2} \\ &\quad + g_{3}\,\mathit{TDF}_{3} + \ldots + g_{\text{n}}\,\mathit{TDF}_{\text{n}}. \end{align} [11]

The TDFs of individual amino acids appear to fall within certain ranges, especially in the lower trophic levels. Also, if the amino acid composition is relatively constant throughout the organisms (Fig. 3), it is reasonable to assume a relationship between TDFbulk and the mean TDF of amino acids. Figure 15 indicates the mean TDF at the lower TP organisms (i.e., not teleost) in the aquatic realm.29),68),120),158) Although the TDF values of cysteine, histidine, arginine, and tryptophan are not known, weighted mean of TDF values is ∼4.2‰ for these amino acids which account for ∼85% of nitrogen of proteinaceous amino acids. This value is significantly (∼0.8‰) greater than the known mean TDFbulk (∼3.4‰).19),22),59) The large mean TDF observed in amino acids must be offset somewhat by nitrogenous compounds with smaller TDF. The four amino acids and the amide nitrogen of glutamine and asparagine, which cannot be measured for analytical reasons, are not taken into account in the above calculations and may contribute to this difference. On the other hand, other nitrogen compounds such as nucleobases need to be taken into account. Quantitatively, proteinaceous amino acids are the major source of nitrogen in the organisms, accounting for 80–90% of biological nitrogen, whereas nitrogen from nucleobases (i.e., adenine, guanine, thymine, cytosine, and uracil) and other nitrogenous compounds roughly account for ∼10% and <5%, respectively.159),160)

Fig. 15.

The TDF (trophic discrimination factor) of individual amino acids for the lower TP (1 to 3) organisms in the aquatic environment (n.d.: no data).29),68),120),158) Weighted mean TDF is calculated based on the averaged distribution of amino acids shown in Fig. 3. TDF values of cysteine, histidine, arginine, and tryptophan are unknown and therefore not included in this calculation.

Like many proteinaceous amino acids, compounds with nitrogen as an amino group are potentially enriched for 15N by deamination. This means that they may have certain magnitudes of TDF. Nucleobases are principally included in the category of such compounds. Although there are only few reports on the nitrogen isotope ratios of nucleobases, those of the five major nucleobases from mammalian DNA rather suggest that the average TDF of the nucleobases is substantially (∼5‰) negative.161) The nucleobases may explain an important part of the difference between the amino acid TDF and TDFbulk. Furthermore, compounds that are primarily supplied through the diet and do not involve nitrogen in the metabolic process (in the same category as the source amino acids) would be expected to lower TDFbulk. Hemes and nitrogenous lipids are probably the compounds in this category.162) A recent analysis of heme B, the most abundant heme species in the natural environment, suggested that the TDF of this compound may be relatively small.163) The gap between the two approaches, amino acids and bulk, should narrow as more such data become available.

Conclusions

It is very significant that macroscale processes such as predation can be measured using 15N abundance of simple compounds that are regulated by statistical thermodynamics and physiology. However, it is obvious that much still remains to be done. As pointed out by the previous summary paper,34) accurate knowledge of TDF especially for higher-TP organisms is still lacking and, as is accurate knowledge of the isotope effects associated with the amino acid metabolism behind it. Even though the method still needs further evaluation, more amino acid isotope ratios would advance food web research, given that observation-based studies lag far behind theoretical studies in this field.164),165)

In this paper, I have focused on the nitrogen isotopic ratios of only three amino acids, glutamic acid, phenylalanine, and methionine, but in the future the isotopic information encoded in other amino acids may also be explored in depth for ecological applications. Although it has not been much of an issue among geochemists, knowledge of the mechanistic details of enzymatic reactions such as deamination and peptide cleavage will eventually be necessary to understand the magnitude of nitrogen isotope fractionation associated with amino acid metabolism.

Isotope ecology has been established and developed by isotopic compositions of bulk samples that accumulated over the last many decades. The bulk approach has the advantage over the amino acid approach providing a systematic perspective when the organism is considered as a box with diverse functions. Bulk information is something that can never be obtained by probing deeply into a single compound or compound group. However, at present the results of the bulk approach do not appear to be very well coupled with those of the amino acid isotopic approach which has been developed rapidly in recent years, as presented in this paper. Further growth of the amino acid approach is needed to link the two more systematically.

Acknowledgments

I am deeply indebted to E. Wada for encouraging me and my group, to pursue this topic over the decade. In writing this review, I would like to thank the three anonymous reviewers and E. Wada for constructive comments. In conducting the research presented here, I sincerely thank my marvelous colleagues working for many years, especially for Y. Chikaraishi, N. O. Ogawa, and Y. Takano who have enormously contributed what I described in this paper. I also wish to thank my excellent colleagues and collaborators around the world who have applied the method to various fields of science; They are T. Nagata, K. Maki, A. Wyatt, A. Kohzu, Y. Yamaguchi, Y. Takizawa, N.F. Ishikawa, H. Nomaki, M. Tsuchiya, K. Fujita, Y. Naito, N. Honch, Y. Itahashi, M. Yoneda, L. Rey, G. Goude, J. Matsubayashi, C. Yoshikawa, Y. Harada, I. Tayasu, C. Yoshimizu, K. Koba, T. Tsutaya, M. Nakabayashi, S.A. Steffan, B.R. Kruger, E.C. Minor, C. Campbell, D.M. Nelson, Q.H.S. Chan, J.T. Brenna, S. Miyachi, M.J. Miller, K. Tsukamoto, Y. Kashiyama, M. Ishikawa, A.S. Goto, S. Furota, T. Furota, K. Fujikura, Y. Fujiwara, A. Tame, T. Yoshida, K.H. Shin, H. Suga, H. Kitazato, N. Yoshida, N. Okada, Y. Miyake, B. Thibodeau, Y. Yokoyama, A. Urai, Y. Sun, Y. Isaji, and Y. Sasaki.

Notes

Edited by Eitaro WADA, M.J.A.

Correspondence should be addressed to: N. Ohkouchi, Biogeochemistry Research Center, Japan Agency for Marine-Earth Science and Technology, 2-15 Natsushima-cho, Yokosuka, Kanagawa 237-0061, Japan (e-mail: nohkouchi@jamstec.go.jp).

References
Non-standard abbreviation list

FCL

food chain length

GDH

glutamate dehydrogenase

TDF

trophic discrimination factor

TP

trophic position

Profile

In 1966, Naohiko Ohkouchi was born in Kyoto City. He graduated from The University of Tokyo and received his Ph.D. from the Department of Geology, The University of Tokyo in 1990 and 1995, respectively. After working at prestigious institutes such as the Kyoto University, Hokkaido University and Woods Hole Oceanographic Institution, U.S.A., he joined the Institute for Frontier Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology (JAMSTEC) in 2002. Subsequently, he studied the radiocarbon dating of organic compounds isolated from sediments and was the first to demonstrate that climate change affects the dynamics of deep sea particles. Furthermore, he and his colleagues at JAMSTEC also worked on the use of nitrogen isotopic ratios of amino acids as a tool for ecosystem studies and chlorophylls and porphyrins for understanding biogeochemical cycles in the photic zone of natural environments. Moreover, Dr. Ohkouchi is currently Director General of Research Institute for Marine Resources Utilization and Director of Biogeochemistry Research Center, JAMSTEC. He was visiting professors at the Department of Environmental Chemistry and Engineering, Tokyo Institute of Technology (2011–2019) and at the Department of Earth and Planetary Science, The University of Tokyo (2012–2014). Dr. Ohkouchi was also the recipient of the Oceanochemistry Award in 2013 and the Duke of Edinburgh Prize of the Japan Academy in 2022.

 
© 2023 The Author(s).

Published under the terms of the CC BY-NC license
https://creativecommons.org/licenses/by-nc/4.0/
feedback
Top