Figures

Abstract

Tracing the origin of nutrients is a fundamental goal of food web research but methodological issues associated with current research techniques such as using stable isotope ratios of bulk tissue can lead to confounding results. We investigated whether naturally occurring δ13C patterns among amino acids (δ13CAA) could distinguish between multiple aquatic and terrestrial primary production sources. We found that δ13CAA patterns in contrast to bulk δ13C values distinguished between carbon derived from algae, seagrass, terrestrial plants, bacteria and fungi. Furthermore, we showed for two aquatic producers that their δ13CAA patterns were largely unaffected by different environmental conditions despite substantial shifts in bulk δ13C values. The potential of assessing the major carbon sources at the base of the food web was demonstrated for freshwater, pelagic, and estuarine consumers; consumer δ13C patterns of essential amino acids largely matched those of the dominant primary producers in each system. Since amino acids make up about half of organismal carbon, source diagnostic isotope fingerprints can be used as a new complementary approach to overcome some of the limitations of variable source bulk isotope values commonly encountered in estuarine areas and other complex environments with mixed aquatic and terrestrial inputs.

Funding: The study was funded by the DFG-supported Cluster of Excellence “The Future Ocean” (EXC 80/1, CP0937) and Carlsbergfondet (2007_01_0301). TL was supported by Juan de la Cierva (JCI-2009-049330) and MV by the Spanish Government projects Fundalzoo (CGL2010-14841) and Invasivefish (427/2011). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

During the last 30 years, stable isotope analysis has emerged as one of most powerful tools for tracing organic carbon in food webs. Analyses of total organic matter (“bulk”) have become widespread due to the relative ease and low cost of sample preparation and analysis. However, the potential of variable environmental conditions to influence carbon isotopic ratios at the base of the food web (δ13C) is a serious drawback for disentangling aquatic and terrestrial sources [1], [2], [3], [4]. With the advent of new continuous flow technologies, compound specific isotope analysis is increasingly employed as a complementary tool for food web analysis. Isotope analysis of fatty acids in conjunction with their structural compositions are now widely used to investigate biosynthetic sources since their molecular structure are tied to their biosynthetic origins (e.g. [5]). However, fatty acids account for only a small fraction of total organic carbon fluxes, and they tend to undergo degradation and transformation during food web passage [6]. Amino acids (AAs), in contrast, account for the large majority of organic nitrogen, and about half of total carbon in most organisms [7] and therefore are among the major conduits of carbon through food chains. The main methodological drawback of the 20 protein AAs is that they are ubiquitous in all life forms. However, about half of the AAs can only be synthesized by bacteria, fungi and photoautotrophs, and are therefore essential or indispensable for animal diets. Since these essential AAs (EAAs) typically pass from food source to consumer without alteration to their carbon skeletons [8], [9], a method for tracking their origins and fluxes could greatly advance our understanding of nutrient cycling and trophic relationships.

Recent research has shown that naturally occurring δ13CAA patterns contain information of both biosynthetic origin and mode of carbon acquisition [10], and that the EAA group is particularly diagnostic of origin [11]. These δ13CAA patterns represent the sum of the isotopic fractionations associated with the individual biosynthetic pathways and associated branch points for each AA (e.g. [12]). Larsen et al. [11] found that the δ13CAA patterns of terrestrial plants, bacteria and fungi were distinct and consistent, and proposed these “stable isotope fingerprints” as a tool for tracing sources of organic matter in terrestrial ecosystems. Comparisons of δ13CEAA patterns between laboratory reared consumers and known diets have also indicated that stable isotope fingerprints are passed on to consumers [11], [13], [14]. In spite of these advances, it has until now been unresolved whether δ13CAA patterns can distinguish between aquatic and terrestrial primary producers [14] and to what extent variable growth conditions for aquatic producers may influence δ13CAA patterns. The factors that affect bulk and compound specific isotope patterns are fundamentally different. While bulk δ13C values for a given producer largely are determined by the ratio of carbon fixation to carbon flux into the cell [15], δ13CAA patterns are determined downstream of the Calvin cycle by AA biosynthetic pathways and associated branching points in the central metabolism. For bulk δ13C based studies, temporal and spatial variations of inorganic carbon sources and other environmental conditions therefore pose a challenge [16], [17]; particularly in the highly productive coastal benthic and estuarine habitats that are also the zones of most organic carbon deposition in modern biogeochemical cycles. A measure based on compound specific stable isotope analysis that could distinguish between terrestrial and aquatic sources of organic matter across varying ecological conditions would therefore have profound implications for biological research across a number of disciplines.

Here we investigated whether δ13CAA patterns based approaches could transcend some of the limitations associated with bulk δ13C, by characterizing δ13CAA values for a large set of different algae and vascular plants. We tested whether these two groups have different δ13CAA patterns, and compared these with δ13CAA patterns from heterotrophic bacteria and fungi. Further, to test the potential for δ13CAA source information to be independent of variation in bulk δ13C values, we analyzed two marine primary producers sampled from a range of different environments within their natural habitats. Finally, to assess the practical relevance of using δ13CAA patterns as diagnostic and quantitative biomarkers in actual ecosystems we analyzed consumers from freshwater, pelagic, and estuarine systems.

Materials and Methods

Sampling Design

To test whether freshwater, marine and terrestrial primary producers have different δ13CAA patterns, we collected and cultured samples from wide variety of primary producers (micro- and macroalgae, and terrestrial plants). In the field, we only collected fresh and newly emerged thallus and leaves. For microbial reference samples we obtained axenically cultured bacteria and fungi. For testing the potential for δ13CAA source information to be independent of variation in bulk δ13C values among aquatic producers, we acquired multiple samples of the seagrass Posidonia oceanica and the giant kelp Macrocystis pyrifera collected from a range of different environments within their natural habitats. For P. oceanica, the maximum photosynthetically active radiation above their canopies ranged from 50 to 310 mmol m−2 s−1 between the sampling locations and the coverage of leaf necrosis ranged from 0 to 37.5%. M. pyrifera samples were collected either in spring or late fall, periods of contrasting ocean conditions that result in widely divergent bulk isotopic values [18]. Note that our sampling was not designed for investigating the influence of geographical region on terrestrial plants since a previous study [14] found no systematic differences in δ13CAA patterns from greenhouse plants and those collected in boreal and mangrove ecosystems. Finally, to assess the practical relevance of using δ13CAA patterns as diagnostic biomarkers in actual ecosystems we analyzed consumers from three well-studied ecosystems: from a marine pelagic ecosystem in the central North Pacific Ocean the carnivorous fish species opah (Lampris guttatus), common dolphinfish (Coryphaena hippurus) and broadbill swordfish (Xiphias gladius), from a littoral marine system (estuarine) the California mussel (Mytilus californianus), and from oligotrophic arctic lakes the water flea Daphinia sp. and seston.

Sample Acquisition and Preparation

A detailed list of all our field samples and their locations is provided in Table S1. Here follows a general description of sampling locations, protocols and permits. All macroalgae (22 species), two seagrasses and two mussels were collected by the Californian shore. The macroalgae except M. pyrifera (see below) were collected on state tidelands for which collection without a permit is allowed for less than ten pounds fresh weight. The seagrass samples were provided under a permit to Joseph M. Long Marine Laboratory, Santa Cruz. Subsamples of mussels were obtained from recreationally harvested mussels under a permit to Natasha Vokhshoori. The three pelagic fish samples were collected as in Choy et al. [19] with NOAA longline observers, and locations for the fish samples are approximate and are reported as the centers of 5×5 degree cells in accordance with NOAA confidentiality policies. Soils and the majority of terrestrial plants (10 out of 12 species) were collected adjacent to five Alaskan tundra lakes in which we also collected seston (5–80 µm size particles) and Daphnia. The Alaskan samples were collected on land owned by the Alaskan state or Bureau of Land Management where no permission was required because we sampled on day trips by foot or rafts leaving no permanent marks; soil samples were collected between 2 and 8 cm below the soil surface and amounted to <50 g dry weight for each location; each plant sample amounted <1 g dry weight. No permit was required for sacrificing Daphnia, which are not protected or endangered; they were sampled with a zooplankton net (150–200 µm) and kept alive until sorting at Toolik Field Station. The terrestrial plant collection was supplemented with two species from public owned land in Denmark where no permission was required for sampling. We obtained five P. oceanica samples from the Catalonian shore in the Mediterranean Sea under a permit to Teresa Alcoverro by the Catalan Water Agency, and the five M. pyrifera samples from the Californian shore were sampled under a permit to Melissa M. Foley [18] by NOAA’s Monterey Bay National Marine Sanctuary. All field samples were kept on crushed ice in coolers and stored between 1 and 3 days before thoroughly washing them in milliQ water (except for animal and seston samples). Seagrass leaves were further cleaned by scraping off epiphytes with a razor blade. The seston samples that during sampling had been collected into 500 ml bottles were in the laboratory concentrated onto GF/C glass microfiber filters (Whatman). For analysis of animal samples, we prepared whole tissue samples of Daphnia and muscle tissue for mussels and fish. All samples were dried at 50°C except for animal, seston and soil samples, which were freeze-dried. After drying the samples were homogenized with a mortar.

Microalgae comprised of five diatoms, two cyanobacteria, two chlorophytes, three haptophytes, and two crysophytes obtained from existing axenic cultures at GEOMAR, Kiel, Germany or the culture collection of algae of Goettingen University, Germany. See Table S2 for detailed list of laboratory samples. The microalgae were cultured between April 2010 and January 2011 at GEOMAR in sterile 225 cm2 tissue culture flasks with vented cap in brackish water (13.9 psu) collected from Kiel Fjord or seawater (31.2 psu) collected at Multimar-Wattforum by the German North Sea. Water mixed with added nutrients was sterile filtered (Whatman celluloseacetate 0.2 µm filter) (see Table S2 for light, temperature and nutrient regimes). In addition we collected three non-sterile filtered water samples from Kiel Fjord February 2012 for the use of culturing of natural assemblages of microalgae in culture tanks. Visual inspection of these samples revealed between 60–80% dominance of diatoms. All microalgae were harvested after 5 to 21 days during exponential growth or right after the onset of the lag phase on either a Durapore PVDF 0.22 µm pore size filter (Sigma-Aldrich, Germany), GF/C glass microfiber filters (Whatman) or by centrifugation in 50 ml vials at 3300 g. The algae were subsequently freeze dried, and surplus substrates with live cultures were autoclaved.

Laboratory grown fungi and bacteria samples were either obtained from a previous study [11] or isolated and cultured at the Institute of Arctic Biology, University of Alaska Fairbanks from soil or water samples collected during the Alaskan field sampling mentioned above (see Table S2 for details). For isolating fungal and bacterial strains, the growth media were treated with either antibiotics (Streptomycin and Tetracycline hydrochloride) or fungicides (Amphotericin and Nystatin, all chemicals were from Sigma-Aldrich, St. Louis, Missouri, USA). After isolation under sterile conditions, all bacteria samples were grown on solid media (Bacto Agar BD, Sparks, Maryland, USA) and fungi were grown in liquid media flasks in a shaking incubator (Innova 4230, New Brunswick Scientific, Edison, New Jersey, USA) and harvested after one to three weeks. We used amino acid free nutrient media containing 15 g L−1 of one of the following nutrient mixes ‘Czapek’, ‘MAG’ and ‘MMN’ (See Table S2 for temperature and nutrient regimes). Microbial samples in liquid media were harvested by centrifugation in 50 ml tubes (2200 g) and freeze dried after harvesting. All substrates with viable cultures were autoclaved after harvesting.

Elemental and Isotope Analysis

Elemental content and bulk isotope ratios of plants, bacteria, fungi, macroalgae and animals were measured at the UCSC Stable Light Isotope Facility. Approximately 1 mg of sample was pelletized into tin capsules and analyzed on a Carlo Erba 1108 linked to a Thermo Finningan DeltaPlus XP mass spectrometer with an analytical standard deviation of typically<±0.15‰ (n = 3). Elemental content and bulk isotope ratios of microalgae were determined on 2 mg samples pelletized into tin capsules at the UC Davis Stable Isotope Facility using a PDZ Europa ANCA-GSL elemental analyzer interfaced to a PDZ Europa 20–20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK). Isotope data are expressed in delta (δ) notation as ((Rsample/Rstandard) –1) × 1000‰, where R is the ratio of heavy to light isotope; and the standard is Vienna Pee Dee Belemnite (VPDB) for carbon and air for nitrogen.

For δ13C analysis of individual AAs (δ13CAA) we transferred between 1.5 and 7 mg of sample to Pyrex culture tubes (13×100 mm). Samples were flushed with N2 gas, sealed, and hydrolyzed in 1–2 ml 6 N HCl (37% HCl diluted with Milli-Q water, Merck, Darmstadt, Germany) at 110°C in a heating block for 20 h. After hydrolysis, samples collected on GF/C filters and coralline algae were purified with Dowex 50WX8 cation exchange resin according to Amelung & Zhang [20] and He et al [21]. In the remaining samples we removed lipophilic compounds by adding 2 ml n-hexane/DCM (6∶5, v/v) to the Pyrex tubes that were flushed shortly with N2 gas and sealed before vortexing for 30 s. The aqueous phase was then filtered through a Pasteur pipette lined with glass wool that had been pretreated at 450°C. All samples were transferred into 4 ml dram vials before evaporating the samples to dryness under a steam of N2 gas for 30 minutes at 110°C in a heating block. The samples were stored at −18°C. To volatize the AAs, we followed the derivatization procedure by Corr et al [22] methylating the dried samples with acidified methanol and subsequently acetylating them with a mixture of acetic anhydride, triethylamine and acetone (NACME: N-acetyl methyl ester derivatives). As a precautionary measure to reduce oxidation of amino acids during derivatization, we flushed and sealed reaction vials with N2 gas prior to the methylation and acetylation reactions. To account for carbon added during derivatization [23] and variability of isotope fractionation during analysis, we also derivatized and analyzed pure amino acids with known δ13C values. Nor-leucine was used as an internal standard. Amino acid derivatives were injected with an autosampler into an Agilent Single Taper Ultra Inert Liner (#5190-2293) that was held at 280°C for 2 min. The compounds were separated on a Thermo TraceGOLD TG-200MS GC column (60 m×0.32 mm×0.25 um) installed on an Agilent 6890N gas chromatograph (GC). The oven temperature of the GC started at 50°C and heated at 15°C min−1 to 140°C, followed by 3°C min−1 to 152°C and held for 4 min, then 10°C min−1 to 245°C and held for 10 min, and finally 5°C min−1 to 290°C and held for 5 min. The GC was interfaced with a MAT 253 isotope ratio mass spectrometer (IRMS) via a GC-III combustion (C) interface (Thermo-Finnigan Corporation). All samples were analyzed in triplicate. The average reproducibility for the internal standard nor-leucine (Nle) was ±0.4‰ (n = 3) and the amino acid standards ranged from ±0.1‰ for Phe to ±0.6‰ for Thr (n = 12–15 for each batch). Of the amino acids we were able to analyze (see Fig. S1 for GC-C-IRMS chromatogram) the following were defined as non-essential for animals: alanine (Ala), asparagine/aspartic acid (Asx), glutamine/glutamic acid (Glx), glycine (Gly), and tyrosine (Tyr). The following were defined as essential: histidine (His), isoleucine (Ile), leucine (Leu), lysine (Lys), methionine (Met), phenylalanine (Phe), threonine (Thr), and valine (Val).

Calculations and Statistical Analyses

All statistical analyses except the mixing modeling were performed in R version 2.12.1 [24] with RStudio interface version 0.96.330. All values in the text are given as mean ± standard deviation. To explore patterns and group memberships in our dataset we performed Ward’s hierarchical clustering (R-package cluster) and principal component analysis (PCA, R-package vegan) on δ13C values of amino acids that had been normalized to their respective sample means denoted as δ13CAAn. Prior to applying statistical analysis the data were tested for univariate normality by visually checking whether there were departures from normality on Q-Q plots. His and Met were excluded from the analyses due to missing measurements caused by concentrations below detection limits. Differences in each amino acid between different producer groups were tested with ANOVA with Tukey HSD post-hoc tests. To examine combinations of independent variables (i.e. δ13CAA values) that best explained differences between the categorical variables (i.e. the groups defined by the PCA and one-way ANOVA tests) and to construct models for predicting membership of unknown samples, we performed linear discriminant function analysis (LDA, R package MASS [25]) on δ13CAA values. For calculating the probability of group membership of the classifier samples we used a leave-one-out cross-validation approach. To test the null hypothesis that there was no difference in classification among the groups we applied Pillai’s trace (MANOVA). Relative contributions of EAAs from diets to consumers was estimated in the software FRUITS (version 0.1, http://sourceforge.net/projects/fruits) [26] with normalized isotope values. FRUITS also considers the biochemical composition of sources and which sources are most likely to contribute the most. FRUITS is executed with BUGS, which is a software package for performing “Bayesian inference Using Gibbs Sampling” that includes an expert system for determining an appropriate Markov chain Monte Carlo scheme based on the Gibbs sampling.

Results

To assess differences in δ13CAA patterns we first carried out a PCA with δ13CAAn values that in contrast to a LDA does not require categorical variables. With this analysis we found that the samples clustered according to their major phylogenetic associations: algae, bacteria, fungi and vascular plants (Fig. 1). The first principal component accounted for 33% of the variation, and separated the photoautotrophs from the microorganisms. The second principal component accounted for 27% of the variation and separated algae from vascular plants, and fungi from bacteria. The AAs in the PCA grouped largely according to their biosynthetic families, with the exception that Ile grouped with the pyruvate AAs (Ala, Leu and Val), but biosynthetically belongs to the oxaloacetate family (Asx, Lys and Thr). The third group consisted of the aromatic AAs, Tyr and Phe. See Table S3 for the δ13C values for each AA measured, and Table S4 for detailed PCA results. We then tested which of the δ13CAAn values were significantly different between algae, terrestrial plants, bacteria and fungi with ANOVA (Table S5). Seven amino acids (Ala, Phe, Tyr, Thr, Val, Leu and Lys) were significantly different between algae and terrestrial plants, all AAs except Glx were significantly different between algae and bacteria, and six (Glx, Gly, Ile, Leu, Lys and Val) were significantly different between algae and fungi (see Table S5 for direction of the differences). The microalgal groups (chlorophytes, chrysophytes, cyanobacteria, diatoms, and haptophytes) only displayed subtle differences between δ13CAAn values (Table S5). For this reason we grouped them to a single group and tested them against the two macroalgal groups. We found the most notable difference between brown algae and microalgae with five significantly different AAs (Ala, Asx, Ile, Lys, Tyr).

Figure 1. The principal component analysis of δ13CAAn values of different producers show a range of different isotope patterns between bacteria, fungi, vascular plant and algae.

None of the microalgal or macroalgal group clustered separately from one another. Values in parentheses are the percentage variation accounted by each axes. The first axis separates the photoautotrophs from the microbes, and the second axis separates vascular plants from algae, and fungi from bacteria. The fairly similar vector lengths show that almost all amino acids were important for the variations of the two first ordination components. See Table S4 for analytical details.

We then applied LDA to our δ13CEAA data to identify which of the six EAAs (Ile, Leu, Lys, Phe, Thr, Val) were most important for distinguishing between algae, bacteria, fungi and terrestrial plants. The seagrass samples were omitted from this analysis due to a small number of samples (the number of samples did not exceed the number of EAAs). Bacterial, fungal, and terrestrial plant samples classified with 99.1±2.6% posterior probability within their own groups (Fig. 2, Table S6). The microalgal samples classified with 97.3±8.1% probability as algae with Melosira varians and Isochrysis galbana having the lowest probabilities (84% and 61%, respectively). All brown algae samples (Phaeophyceae) classified with 100% probability as algae, and eight out of nine red algae (Rhodophyta) classified with 98.6±3.8% probability as algae. The ninth red algae, Osmundea spectabilis, classified as a bacterium (54% probability) rather than an alga (46% probability). The most important linear discriminants for separating the four categorical variables were Lys, Phe, Leu and Val. We created a second LDA model based on the five most informative EAAs (Ile, Leu, Lys, Phe, Val) to assess to what extent it would be possible to separate the three algal groups (microalgae, brown algae and red algae) and seagrass samples from each other. Of all 27 microalgal samples, 24 samples classified as microalgae, 10 out of 12 brown algae classified as brown algae, 7 out of 9 red algae classified as red algae, and all 7 seagrass samples classified as seagrass (Fig. 3, Table S7). None of the algal samples classified as seagrass.

In the left figure (a) displaying the scores of the first two discriminant axes, fungi and terrestrial plants each cluster separately from algae and bacteria. In the right figure (b) displaying the second and third discriminant axes, bacteria are separated apart from the algae, fungi and terrestrial plants. The dotted lines represent confidence ranges at P = 0.5. See Table S6 for details.

It separates all seagrass samples from the three algal groups. The majority of the algal samples classified correctly within their own groups (Table S7). The dotted lines represent confidence ranges at P = 0.5; confidence ranges are only displayed in the left figure (a) because the third linear discriminant in right figure (b) only explained 14%.

The second major question asked by our study was to what extent δ13CAA patterns are affected by environmental conditions. To answer this question we analyzed seagrass (P. oceanica) and giant kelp samples (M. pyrifera) [16] across a variety of growth conditions (see Table S1 for details). For both species the range in δ13C values was five- to ten-fold greater for bulk than δ13CAAn values (Fig. 4). Individual δ13CAAn values typically spanned between 0.4 to 0.6‰ compared to 2.6% and 5.2‰ for bulk δ13C values of P. oceanica and M. pyrifera, respectively.

Finally, we investigated how δ13CEAA patterns of animals from three different aquatic ecosystems resembled the main primary production sources in their respective environments. For Arctic shallow lakes in Northern Alaska, we used bacteria, microalgae and terrestrial plants as the most likely end members for Daphnia. The LDA model classified all the categorical variables correctly, and the δ13CEAA patterns of the five Daphnia samples resembled microalgae with 84.0±16.9% probability and bacteria with 11.9±12.3% probability (Fig. 5a, Table S8). The seston sample resembled microalgae with 96.8% probability. While three out of four soil samples resembled plants with >97% probability, the remaining sample resembled plants with 69% probability. In an open pelagic system, we used microalgae, bacteria and fungi as end members for three predatory fish species (C. hippurus, L. guttatus and X. gladius) from the central Pacific. The three categorical variables (algae, bacteria and fungi) were distinctly different and the δ13CEAA patterns of all three fish samples matched those of microalgae with 100% probability (Fig. 5b, Table S8). In the estuarine system we selected microalgae, giant kelp and bacteria as the most likely particulate organic matter sources for the California mussel (Mytilus californianus). The LDA model classified all the categorical variables correctly, and the δ13CEAA patterns of the two mussels resembled microalgae with ≥99.7% probability and brown algae with 0.3% probability (Fig. 5c, Table S8). We used the California mussel samples to exemplify how EAA isotope values can be used to obtain relative proportions of food sources for a consumer. The mixing model was based on δ13Cn values of the three most informative EAAs for separating microalgae, brown algae and bacteria: Leu, Val and Lys. For brown algae, we only considered kelp since it is the most dominant brown alga in the mussels’ habitat. We also included information about the relative proportion of the three EAAs (Leu, Lys and Val) in the food sources [27]–[29], and the likelihood that the mussels would consume one food sources over another (See Appendix S1). We found that the mussels obtained about two-third of their EAAs from microalgae (66.7±13.3%), about a quarter from kelp (26.0±11.4%) and the remaining fraction from bacteria (7.3±6.0%) (Fig. 6, Appendix S1).

Figure 5. Application of source diagnostic δ13CEAA patterns in food web studies across three different ecosystems.

(a) In oligotrophic arctic lakes in Alaska, Daphinia sp. and seston cluster closely to each other, and their EAAs appear to derive predominantly from microalgae although a part of their EAAs may have come from foods reworked by bacteria or from allochtonous sources (i.e soils). (b) In the central North Pacific Ocean the EAAs of the carnivorous fish species (opah; Lampris guttatus, common dolphinfish; Coryphaena hippurus, broadbill swordfish; Xiphias gladius) resembled microalgae rather than EAAs from bacteria and fungi. (c) In a complex littoral marine system by the Californian shore, the δ13CEAA fingerprints of California mussel (Mytilus californianus) resemble microalgae and not bacteria or brown algae, i.e. kelp. In the figure legend, ‘Pr’ signifies predicted samples. See Table S8 for analytical details.

The boxes provide a 68% confidence interval (corresponding to the 16th and 84th percentiles) and the whiskers provide a 95% confidence interval. The horizontal continuous line indicates the average while the horizontal discontinuous line indicates the median (50th percentile). See Appendix S1 for detailed information.

Discussion

Our results show that δ13CAA patterns can be used as powerful and ubiquitous tracers for discerning carbon origins in both terrestrial and marine settings (Fig. 1). We found that δ13C values for seven out of eleven AAs were significantly different between algae and terrestrial plants, which enabled us to create a classification model that determined whether the AAs originated from terrestrial or marine sources (Figs. 2 and 3). The fact that variations in δ13CAA patterns among algae were sufficiently constrained to make a clear distinction between aquatic and terrestrial primary producers seems remarkable, considering that our algal samples encompassed a large number of species from both freshwater and marine environments. Within the algal groups that comprised two macroalgal and five microalgal domains, brown algae stood out as having the most distinct δ13CEAA patterns (Fig. 5c, Table S5). We also found that δ13CEAA patterns of terrestrial and aquatic primary producers were different from bacteria and fungi. While the mechanistic reasons why algae vs. vascular plants have such unique δ13CAA patterns are not currently clear as discussed below, the strong diagnostic potential we observe for these major classes of primary producers is consistent with broad ability of δ13CAA patterns to distinguish between other major organism groups.

Our results also indicate that δ13CEAA patterns applied as source diagnostic isotope fingerprints may offer a partial solution to one of the major issues for the application of stable isotopes in ecological and biogeochemical research: confounding variable source bulk isotope values. It has remained a persistent challenge for isotope ecologists to disentangle confounding δ13C values caused by variations in inorganic carbon sources and other environmental parameters. Here we show that despite substantial shifts in bulk δ13C values for the seagrass P. oceanica and the giant kelp M. pyrifera linked to season or growth conditions, δ13CAAn values were constant within a 0.5‰ standard deviation (Fig. 4). These results are consistent with the notion that δ13CAA patterns are mostly determined by major evolutionary AA metabolic pathways of an organism [11], [12], [30] rather than the factors affecting bulk δ13C values such as carbon availability, growth rates, and cell surface area [3]–[4]. For application in food web studies it is particularly encouraging that δ13CAAn values only shifted by 0.5‰ compared to the >5‰ shift in bulk δ13C values of giant kelp. Another important observation is that light attenuation and leaf necrosis for the seagrass samples did not affect δ13CAAn values notably. Light intensity and necrosis are important factors for seagrass growth and often associated with changes in leaf composition of phenolic compounds, carbohydrates and chlorophyll [31]–[33].

We directly tested the ability of δ13CAA patterns to transfer information about major primary producer sources through food webs by examining δ13CAA values in consumers from diverse habitats (freshwater arctic lakes, subtropical pelagic ocean, estuarine marine). In every habitat, the results were consistent with the hypothesis that δ13CEAA source patterns based are conserved in passage through food webs. In the freshwater ecosystem, the resemblance of Daphnia and seston AAs with microalgae rather than bacteria and terrestrial plants agrees with limnological food web studies based on bulk isotope data [34] (Fig. 5a). However, it is interesting that neither Daphnia nor seston were projected directly on top of the algal samples pointing to a possible influence of microbial reworking or allochtonous input, as has been previously suggested in Arctic lakes [35]. In the open ocean almost all primary production derives from single-celled algae [36], so δ13CAA patterns would be expected to align with algal sources in oceanic consumers. The δ13CEAA patterns exactly followed this prediction (Fig. 5b), which also implies that δ13CEAA diagnostic information was not substantially altered by microbial reworking [37]–[40], or transformed during digestive processes [41]. Our estuarine system is much more complex, with microalgae, giant kelp, and reworked organic matter as possible contributors to suspended particular organic matter (POM) [18]. For the estuarine classification of California mussels, we included bacteria in addition to microalgae and brown algae as end members in our model, and found that microalgae had a much higher probability being an AA source for the mussels than kelp (Fig. 5c). Since classification models are not suited to estimate relative proportions of food sources, we applied a mixing model to the study with estuarine mussels. Based on the three most informative EAAs for separating the mussel’s most likely food sources, we found that the mussels obtained about two-third of their EAAs from mussels, about a quarter from kelp and the remaining fraction from bacteria (Fig. 6, Appendix S1). Thus, our findings are consistent with the expectation that mussels mainly feed on POM derived from fresh phytoplankton, and that microbially reworked POM is a minor source [42].

If algae and vascular plants share similar amino acid biosynthetic pathways [43], [44], it raises the question why they have consistently different δ13CAA patterns. It has been proposed by Hayes [30] that growth rates potentially could affect intramolecular isotope distribution. During the stationary growth phase the removal of carbon from the tricarboxylic acid cycle would be slower than in the exponential phase. This would lead to the accumulation of 13C enriched compounds at the ends of the biosynthetic pathways in turn delivering 13C depleted precursors to the first steps of the pathways [30]. While we did not test this hypothesis, the rather uniform 13CAA patterns among microalgae with contrasting growth rates suggest that this influence is rather small. In support of this view, we found that the Mediterranean seagrass samples had similar δ13CAA patterns in spite of the different lighting regimes and hence growth rates in their natural habitats. Isotope fractionation at metabolic branch points may also explain some of the observed differences in δ13CAA patterns [11], [12,]. We found that algae were 13C enriched relative to plants for the AAs belonging to the pyruvate group (Ala, Val, Leu). A similar case of greater 13C enrichment in algae than plants was found by Chikaraishi and Naraoka [44] for n-alkanes that like most other lipids have pyruvate as a precursor [45]. In contrast to the pyruvate AAs, the aromatic AAs (Tyr, Phe) were 13C depleted in algae relative to plants. From a carbon mass balance point of view it is possible that 13C enrichment of pyruvate AAs in algae lead to depletion of the aromatic AAs because they are less coupled to lipid synthesis by having phosphoenolpyruvate and erythrose-4-phosphate as precursors. The aromatic AAs also serve as precursors for the synthesis of numerous primary and secondary metabolites such as alkaloids and lignins [46], which also could have influenced 13CAA fractionation of the aromatic AAs. For deepening our understanding of the biochemical processes leading to the 13CAA fractionation patterns further inquiry is needed into the overall carbon mass balance between the most abundant hydrocarbon groups and AAs.

Taken together, our results show that δ13CAA patterns can transcend bulk isotope analyses across diverse ecological environments, and be used to understand carbon sources and transfer in ecological research as exemplified with our classification and mixing modeling approaches. The strong diagnostic potential for algae and plants may be particularly powerful for applications in complex estuarine and coastal systems, where mixed aquatic and terrestrial inputs occur, or in freshwater environments with strong allochtonous influence. While our findings indicate that source-specificity of δ13CEAA patterns is conserved across environmental gradients, further controlled physiological studies are also warranted to better understand under what circumstances these patterns may be altered. For a broad understanding of food web cycling of nutrients, we also stress that it is also important to consider other major biochemical classes, such as lipids and carbohydrates. Finally, in addition to investigating ecosystem level transfer of carbon and nitrogen, we suggest that this fingerprinting method can help assessing symbiotic contributions of AAs from bacteria to animal hosts. There is mounting evidence that protein-nitrogen assimilation is possible in the lower gut of some animals during digestion [47]–[49]. Our results indicate that δ13CEAA patterns may offer a direct way to assess the importance of such microbial AA contributions, not only in the specific animals where this may occur, but more broadly up food chains.

9.
O’Brien DM, Fogel ML, Boggs CL (2002) Renewable and nonrenewable resources: Amino acid turnover and allocation to reproduction in lepidoptera. Proceedings of the National Academy of Sciences of the United States of America 99: 4413–4418.

19.
Choy CA, Popp BN, Kaneko JJ, Drazen JC (2009) The influence of depth on mercury levels in pelagic fishes and their prey. Proceedings of the National Academy of Sciences of the United States of America 106: 13865–13869.