Factors Controlling the Stable Nitrogen Isotopic Composition (δ15N) of Lipids in Marine Animals.

Lipid extraction of biomass prior to stable isotope analysis is known to cause variable changes in the stable nitrogen isotopic composition (δ15N) of residual biomass. However, the underlying factors causing these changes are not yet clear. Here we address this issue by comparing the δ15N of bulk and residual biomass of several marine animal tissues (fish, crab, cockle, oyster, and polychaete), as well as the δ15N of the extracted lipids. As observed previously, lipid extraction led to a variable offset in δ15N of biomass (differences ranging from -2.3 to +1.8 ‰). Importantly, the total lipid extract (TLE) was highly depleted in 15N compared to bulk biomass, and also highly variable (differences ranging from -14 to +0.7 ‰). The TLE consisted mainly of phosphatidylcholines, a group of lipids with one nitrogen atom in the headgroup. To elucidate the cause for the 15N-depletion in the TLE, the δ15N of amino acids was determined, including serine because it is one of the main sources of nitrogen to N-containing lipids. Serine δ15N values differed by -7 to +2 ‰ from bulk biomass δ15N, and correlated well with the 15N depletion in TLEs. On average, serine was less depleted (-3‰) than the TLE (-7 ‰), possibly due to fractionation during biosynthesis of N-containing headgroups, or that other nitrogen-containing compounds, such as urea and choline, or recycled nitrogen contribute to the nitrogen isotopic composition of the TLE. The depletion in 15N of the TLE relative to biomass increased with the trophic level of the organisms.

Introduction

Stable isotopes of n class="Chemical">carbon and nitrogen (δ13C and δ15N) are routinely used in ecology to study a wide range of subjects such as trophic interactions, energy flow, diet composition, feeding habits, and migration (see e.g. [1,2]; and references therein). Stable carbon isotopes are generally used to distinguish between energy sources, such as terrestrial vs. aquatic, because differences in δ13C values are generated during primary production and largely conserved during heterotrophic processing (e.g. [3,4]). Nitrogen isotopes are mainly used to infer trophic transfers as each trophic step results in an increase in the δ15N signal of biomass (e.g. [5,6]). More recently, the δ15N of specific amino acids has been used to estimate trophic levels [7-9].

For application of stable isotopes in trophic ecology, n class="Chemical">lipids are sometimes removed from bulk biomass prior to stable isotope analysis. This is done because lipids are depleted in 13C compared to proteins and carbohydrates due to fractionation during lipid biosynthesis and because different tissues and organisms have variable lipid contents [10,11]. Interestingly, it has frequently been observed that lipid extraction also result in changes in the δ15N of the residual biomass [12-15]. The change in δ15N bulk biomass with lipid extraction (Δδ15Nresidue-bulk) can vary substantially (-2 to + 2.1 ‰ [1]), but unlike the change in δ13C from lipid extraction (e.g. [16-18]), no clear relationship between Δδ15Nresidue-bulk and parameters such as lipid content or C:N ratio of the organisms has been found. The relatively large range in Δδ15Nresidue-bulk is problematic as the change caused by lipid extraction is similar in magnitude to reported trophic fractionation (e.g. [6,19]) and diet-tissue discrimination factors (e.g. [6,15,20]). There is uncertainty about the cause for these changes. Some studies hypothesized that co-extraction of lipid-bound proteins leads to the removal of some amino acids (e.g. [13]), although this would imply that the extracted amino acids had a strongly different nitrogen isotopic composition compared to the remaining amino acids [14]. Another hypothesis is that lipid extraction leads to removal of cellular waste products (e.g. ammonia), which have quite different nitrogen isotopic compositions than that of organic nitrogen [21]. However, no experimental evidence has been provided to support these hypotheses.

In this study we investigate the cause of the changes seen in the δ15N following n class="Chemical">lipid extraction of tissues of several marine animals by determining the δ15N of bulk and residual, lipid-free, biomass as well as of the total lipid extract (TLE). In addition, we identify the intact polar lipids present in the lipid extract to elucidate the sources of nitrogen to the lipid extracts. Finally, we determine the δ15N of amino acids to show that δ15N of the total lipid extract relates to the δ15N of serine and the source amino acid phenylalanine.

Materials and Methods

2.1 Samples

Four species of benthic invertebrates (Common cockle, n class="Species">Cerastoderma edule; Pacific oyster, Crassostrea gigas; Green shore crab, Carcinus maenas; and lugworm, Arenicola marina) and three fish species (Atlantic herring, Clupea harengus; Brown trout, Salmo trutta; and Twait shad, Alosa fallax) were collected in the Dutch Wadden Sea in 2011. Fishes were collected using a passive fishing device (kom-fyke net; [22]) in the Marsdiep area near the southern part of Texel (Netherlands), and at the entrance to the western Dutch Wadden Sea, in spring (April to June) and autumn (September to October) of 2011. Individuals were either processed or frozen immediately after capture. All fish were measured to the nearest cm total length and weighed, and gill tissue and white muscle were sampled (see [18]). Muscle tissue was taken from below the dorsal fin and cleaned from skin and scales. Invertebrates were collected through the SIBES (Synoptic Intertidal Benthic Survey) program (https://www.nioz.nl/sibes). Animals from sediment cores were collected, cleaned, and identified. Depending on the type and size of the animal, different types of tissue samples were collected from each individual (e.g. muscle, head, or whole), which were then frozen, followed by freeze-drying for 72 h in glass vials, and stored at -20°C until further processing. Before lipid extraction and stable isotope analysis, samples were homogenized using a mortar and pestle or a ball mill grinder (Retsch, Düsseldorf, Germany).

2.2 Lipid extraction (TLE)

Total lipid extracts (n class="Disease">TLE) were prepared as described in Svensson et al. [18]. In short, samples were extracted four times using dichloromethane (DCM) and methanol (MeOH) (2:1 v/v) and ultrasonication (1x10 min plus 3x5 min) and centrifuged at 1000xg, 2.5 min. Organic solvents were pipetted off after each extraction and combined as the total lipid extract (TLE). Residual biomass (lipid free) and TLEs were evaporated to dryness under a gentle stream of N2 at room temperature.

2.3 Intact polar lipid analysis

Intact polar lipids (IPLs) were analyzed using HPLC/ESI/MS according to Sturt et al. [23] with some modifications as described in Schouten et al. [24] and Bale et al. [25]. In short, n class="Chemical">lipid extracts were re-dissolved in hexane:2-propanol:water (72:27:1, v/v/v) at a concentration of 2 mg mL-1 and filtered through a 0.45 μm regenerated cellulose (RC) filter (Alltech Associates Inc., Deerfield, IL) prior to injection. An Agilent 1200 series LC (Agilent, San Jose, CA), equipped with thermostatted auto-injector and column oven, and coupled to a Thermo LTQ XL linear ion trap with Ion Max source with electrospray ionization (ESI) probe (Thermo Scientific, Waltham, MA), was used. Separation was achieved on a LiChrospher diol column (250 x 2.1 μm, 5 μm particles; Alltech Associates Inc., Deerfield, IL) maintained at 30°C. The following elution program was used with a flow rate of 0.2 mL min-1: 100% A for 1 min, followed by a linear gradient to 66% A: 34% B in 17 min, maintained for 12 min, followed by a linear gradient to 35% A: 65% B in 15 min, where A = hexane:2-propanol:formic acid:NH3aq (14.8M) (79:20:0.12:0.04, v/v/v/v) and B = 2-propanol:water:formic acid:NH3aq (14.8M) (88:10:0.12:0.04, v/v/v/v). Total run time was 60 min with a re-equilibration period of 20 min in between runs. The lipid extracts were analyzed by an MS routine where a positive ion scan (m/z 400–2000) was followed by a data dependent MS2 experiment where the base peak of the mass spectrum was fragmented (normalized collision energy 25, isolation width 5.0, activation Q 0.175). This was followed by a data dependent MS3 experiment where the base peak of the MS2 spectrum was fragmented under identical fragmentation conditions. This process was repeated on the 2nd to 4th most abundant ions of the initial mass spectrum. Major IPL classes were identified as described in Brandsma et al. [26].

2.4 Stable isotope analysis

The stable nitrogen isotopic ratio was determined on bulk biomass (δn class="Chemical">15Nbulk) and on residual biomass after extraction (δ15Nresidue), as well as on the total lipid extracts (δ15NTLE). For δ15Nbulk and δ15Nresidue, ca. 0.4–0.8 mg of freeze dried, homogenized sample was weighed into tin cups. These samples were analyzed for δ15N and percent total organic carbon (%TOC) and percent total nitrogen (%TN) content in duplicate using isotope ratio monitoring mass spectrometry (IRMS) with a Delta V Advantage-IRMS coupled to a Flash 2000 elemental analyzer (Thermo Scientific). Total lipid extracts were dissolved in ethyl acetate and pipetted into tin cups for a final weight of ca. 0.4 mg for δ15N determination. Ethyl acetate was allowed to evaporate completely at room temperature (minimum 6 h) before folding the cups for analysis. Due to the low amount of nitrogen compared to carbon in the lipid extracts, the TLE fraction was analyzed for δ15N on a Delta XL isotope ratio MS (Thermo Finnigan) coupled to a Flash 1112 Series elemental analyzer (CE Instruments) equipped with a liquid nitrogen trap to remove CO2 from the sample stream.

Stable isotope ratios are expressed using the δ notation in units n class="Chemical">per mil according to: where R = 15N/14N, and expressed in per mil versus air. An acetanilide standard with a δ15N value of 1.18‰ (standard obtained from Arndt Schimmelmann, Indiana University; [27]), and known %TOC and %TN content, was used for calibration. The average repeatability of δ15N determination was 0.2‰ based on repeated analysis of the acetanilide standard. The pooled standard deviation of replicate measurements (n = 2–5) were ≤0.6 ‰ for bulk and lipid-free biomass, and ≤0.9‰ for total lipid extract.

2.5 Amino acid nitrogen isotope analysis

The method used for amino acid isotope anan class="Chemical">lysis was a modified version of the method from Metges et al.[28] and Chikaraishi et al. [29]. In short, tissue samples were hydrolyzed in 6M HCl at 110°C for ca. 18 h. Hydrolyzates were then filtered (GHP Nanosep centrifugal filtration devices, Pall Co), and defatted using n-hexane:DCM (3:2, v/v). Samples were dried under a flow of N2 with repeated additions of methanol. Carboxylic acid groups of amino acids were esterified by addition of 200 μl 2-propanol:acetyl chloride (4:1, v/v) and heating at 110°C (2 h). Amine groups were subsequently acetylated using DCM:pivaloyl chloride (200 μl; 4:1,v/v) at 110°C (2 h). After each derivatization step, any leftover reagents were removed by addition of DCM and drying under a gentle stream of N2 (2x). Derivatized amino acids were dissolved in ca. 200 μl bidistilled water and extracted using n-hexane:DCM (3:2, v/v), dried over MgSO4, diluted to a suitable concentration with dried (MgSO4) and de-gassed DCM, and stored at -20°C in amber vials until analysis.

The N-pivaloyl/n class="Chemical">2-propyl derivatives of amino acids were analyzed by gas chromatography/combustion/isotope ratio mass spectrometry (GC/C/IRMS) using a Thermo Delta V Advantage connected to an Agilent 6890 GC. The GC and IRMS were interfaced via a combustion furnace (980°C), reduction furnace (650°C), and a liquid nitrogen cold trap to remove CO2. Separation of the derivatized amino acids was achieved on a DB-5ms column (Agilent J&W, 60 m x 0.32 mm i.d., 0.50 μm film thickness; Agilent Technologies), using the following temperature program: Initial temperature 70°C for 1 min; ramp up to 140°C at 10°C min-1, dwell for 5 min; ramp up to 190°C at 2°C min-1; ramp up to 300°C at 10°C min-1, hold for 10 min. Carrier gas was helium at a continuous flow of 2 ml min-1 (29 cm s-1). Injection volumes ranged from 0.4–2 μl. An in-house standard mixture consisting of five amino acids (glycine, norleucine, glutamic acid, phenylalanine, and tyrosine) with known δ15N values (determined offline) was used to evaluate daily system performance. Long-term reproducibility, based on the standard deviation of multiple injections (n = 52) of the in-house standard mixture was 0.7 (glycine), 0.9 (norleucine), and 1.2‰ (phenylalanine, glutamic acid, and tyrosine).

2.6 Data analysis

Correlations were evaluated using n class="Chemical">Pearson correlation analysis. The non-parametric Kruskal-Wallis test was used to evaluate differences in isotopic compositions, because of the relative small dataset. For the few cases with larger sample size, a Student t-test was used and the results were similar. Data were non-transformed and evaluated at a 5% significance level. Statistical tests were done using XL-Stat version 2015.4.01.20780.

Data used in this study are also available at http://doi.pangaea.de/10.1594/PANGAEA.855456

2.7 Ethical statement

This study was done with permission from the Dutch Fisheries Insn class="Chemical">pection of the Ministry of Agriculture, Nature and Food Quality (Visserijinspectie) and reviewed and approved by the Animal Experiments Committee (Dierenexperimentencommisie, DEC) under DEC protocol NIOZ 2010.03.

Results and Discussion

3.1 δ15N contents of extracted biomass and lipid extracts

Lipids were extracted from biomass of several sn class="Chemical">pecies of aquatic animals and the δ15N values were determined both before (bulk biomass) and after lipid extraction (residual, lipid-free biomass), as well as that of the total lipid extract (TLE). Bulk biomass δ15N values for different species and tissue types ranged, on average, from 6.7 to 17‰ (Table 1). The δ15N of residual biomass after lipid extraction differed by -0.9 to + 1.8 ‰ compared to bulk biomass (Δδ15Nresidue-bulk; Fig 1 and S1 Table) which is consistent with previous studies [30,31].

Table 1

Ranges of δ15N values and C:N ratios for bulk biomass, residual (lipid-free) biomass and lipids (total lipid extract) per species and tissue type.

n = number of individuals analyzed. Standard deviation of δ15N values from replicate measurements were for bulk and lipid-free biomass ≤0.6 ‰, for total lipid extract ≤0.9.

		Bulk biomass				Residual biomass			Total lipid extract
Species	Tissue	n	%lipids	δ15N	C:N	n	δ15N	C:N	n	δ15N	C:N
Atlantic herring	Gill	6	19–32	10.6–14.4a	4.8–8.0a	6	10.7–13.7	3.0–3.3	4	0.4–9.2	62.9
	Muscle	6	7–40	11.4–16.2a	3.3–6.1a	6	11.4–16.3	2.4–3.1	5	-0.5–3.1	14.6–23.2
Brown trout	Gill	6	4–11	13.9–16.3a	3.7–4.0a	6	14.2–16.7	3.2–3.5	3	7.2–11.0	19.5–42.2
	Muscle	7	5–25	13.7–16.5a	3.2–5.1a	7	14.4–17.0	3.0–3.2	5	4.3–9.7	21.3
Twait shad	Gill	2	6–9	14.7–16.7a	4.2–4.3a	2	14.8–17.0	3.3–3.6		n.d	n.d
	Muscle	2	6–7	15.8–17.1a	3.2–3.2a	2	16.7–17.9	3.1–3.1	1	10.3	14.1
Green shore crab	Muscle	3	2–6	13.2–15.8	3.7–5.3	3	13.4–16.5	3.3–4.4	1	8.2	n.d
Common cockle	Muscle	5	4–5	11.3–12.6	4.1–4.9	5	11.7–12.7	3.7–4.3	2	6.9–7.2	n.d
Pacific oyster	Muscle	2	2–6	12.3–12.3	3.0–3.3	2	13.3–13.5	3.2–3.3	2	6.0–7.0	n.d
Lugworm	Head	2	3–5	11.7–12.4	4.2–4.7	2	10.8–11.3	3.7–4.0	2	8.4–8.5	n.d
	Whole	1	1	6.7	3.1	1	6.6	3.5	1	7.4	n.d

a Data from Svensson et al. [18].

n.d. = not determined.

Fig 1

δ15N values of residual biomass, the total lipid extract and serine normalized to bulk biomass.

Median values (symbols) and ranges of differences in δ15N of residual biomass, the total lipid extract (TLE) and serine compared to bulk biomass (Δδ15Nfraction-bulk) for the different animals.

n = number of individuals analyzed. Standard deviation of δ15N n class="Chemical">values from replicate measurements were for bulk and lipid-free biomass ≤0.6 ‰, for total lipid extract ≤0.9.

Median values (symbols) and ranges of differences in δn class="Chemical">15N of residual biomass, the total lipid extract (TLE) and serine compared to bulk biomass (Δδ15Nfraction-bulk) for the different animals.

The TLEs were depleted in n class="Chemical">15N compared to bulk biomass in all samples except one (whole lugworm, Δδ15NTLE-bulk = +0.7‰), with Δδ15NTLE-bulk values ranging from -14 to +0.7‰ (Fig 1, Table 1, and S1 Table). The lipid content of the investigated tissues ranged from 1–40% with the majority having lipid contents ≤ 10% and lipids thus form a relative small portion of the total biomass. Furthermore, the relative amount of N in the lipid extract was small in most cases (average C:N ratio ranging from 14 to 63; Table 1) compared to bulk biomass (C:N ratio ranging from 3–8; Table 1) and highly variable. No correlation was observed between Δδ15Nresidue-bulk and either the C:N ratio or %lipid of bulk biomass (Pearson correlation, r = -0.130 [P = 0.417], n = 40, and 0.239 [P = 0.132], n = 40, respectively).

Although there is a large variability in 15N of the n class="Chemical">lipid extracts, on average they were significantly depleted in 15N compared to biomass by approx. -7‰ (Kruskal-Wallis test, p<0.001, n = 26). The Δδ15Nresidue-bulk of lipid-free biomass also varies, but the residue was on average significantly enriched in 15N by 0.4‰ (Kruskal-Wallis test, p<0.001, n = 41) compared to the original biomass. This suggests that, in general, lipid extracts contain 15N depleted nitrogen, which induces small but significant changes in 15N in the residual biomass. Remarkably, however, there is a large variability in both δ15N and C:N ratio of lipid extracts as well as Δδ15Nresidue-bulk (Table 1). This may suggest that there is a mixture of different nitrogen sources in the lipid extract, such as lipids, lipid-bound protein and urea [12-14], all present in variable ratios and with different δ15N values. Below we investigate the δ15N of lipids as a potential source for the variable depletion in 15N in the lipid extract.

3.2 Sources of 15N depleted nitrogen in lipid extracts

To investigate the origin of the 15N depleted n class="Chemical">nitrogen in the TLE we investigated its lipid composition by HPLC/ESI/MS. The majority of identified lipids in the TLEs comprised the nitrogen-containing phosphatidylcholines (PC; Fig 2 and S2 Table), which is a common lipid class in animal tissue (e.g. [32]). Taurine conjugated lipids as well as phosphatidylethanolamines (PE), betaines, and phosphatidylinositols (PI) were also detected in some species. With the exception of PI, all of these lipids contain one nitrogen atom in the lipid headgroup (Fig 3). The dominance of nitrogen-containing lipids suggests that they form an important source of nitrogen to the lipid extracts, although a contribution of nitrogen from other sources, such as extractable proteins and/or nitrogen containing waste products, cannot be excluded. Some evidence for the latter comes from the C:N ratio of the lipid extracts: the C:N ratio of some lipid extracts (Table 1) were lower than the theoretical C:N ratios of the identified lipids (25–50) suggesting the presence of additional nitrogen.

Fig 2

LC–MS chromatogram showing identified intact polar lipids in brown trout (Salmo trutta) gill tissue.

Base peak LC–MS chromatogram (Gaussian smoothed) of MS1 of intact polar lipids (IPLs) in brown trout (Salmo trutta) gill tissue showing the prevalence of the nitrogen-containing IPL phosphatidylcholine. PC = phosphatidylcholine.

Fig 3

Structures of identified intact polar lipids in lipid extracts of animal tissues.

PC = phosphocholine.

Base peak LC–MS chromatogram (Gaussian smoothed) of MS1 of intact polar n class="Chemical">lipids (IPLs) in brown trout (Salmo trutta) gill tissue showing the prevalence of the nitrogen-containing IPL phosphatidylcholine. PC = phosphatidylcholine.

PC = n class="Chemical">phosphocholine.

To assess whether the N-containing headgroups, in particular n class="Chemical">PC, were 15N-depleted compared to bulk biomass, as observed for the lipid extracts, the biosynthetic source of nitrogen for these headgroups was considered. Headgroups like PC and PE derive their nitrogen either from amino acids, choline or from recycled nitrogen within the cell [33,34]. The amino acid serine in particular is interesting as it is one of two main sources of cellular nitrogen to many nitrogen-containing lipids (the other being choline), while for PE it is the only nitrogen source. If serine contributes a large fraction of nitrogen to the lipids during biosynthesis, the δ15N values of the precursor serine should be reflected in the lipid product. We, therefore, determined the δ15N of serine in a selection of tissues (Table 2). Serine was indeed depleted in 15N compared to bulk biomass (Fig 1) with Δδ15Nser-bulk ranging from -1 to -8‰. This agrees with previous observations in animals that serine generally is 15N depleted compared to total biomass (e.g. [35-37]). The range of depletion in 15N of serine was highly variable and resembled the variation in Δδ15NTLE-bulk of the total lipid extract. Indeed, the Δδ15Nser-bulk of serine is strongly and significantly correlated with Δδ15NTLE-bulk of TLE (Pearson correlation, r = 0.94, [P<0.001], n = 8). The strong correlation between Δδ15Nser-bulk and Δδ15NTLE-bulk is in line with the biochemical evidence that serine is indeed an important source of nitrogen to the lipid pool and shows as well that the δ15N of lipid extracts were mainly determined by the δ15N of N-containing lipids. It should be noted, however, that the extent of depletion of serine was less than for the lipids by 3–4‰ in some tissues (Fig 1). These differences in δ15N between serine and lipid extracts may be due to one or a combination of several factors: (i) a (large) contribution of nitrogen recycled within the cell or originating from dietary choline, (ii) fractionation during biosynthesis of the headgroups from serine or choline. Further research using e.g. compound specific analysis of headgroups of lipids may shed further light on this.

Table 2

δ15N values of amino acids.

Avg and s.d. = average and standard deviation of n injections. ALA = alanine; ASP = aspartic acid; GLU = glutamic acid; GLY = glycine; ILE = isoleucine; LEU = leucine; LYS = lysine; OH-PRO = hydroxy-proline; PHE = phenylalanine; PRO = proline; SER = serine; THR = threonine; TYR = tyrosine; VAL = valine.

	Brown trout muscle			Brown trout gill			Atlantic herring muscle			Atlantic herring gill			Green shore crab			Pacific oyster			Lugworm head			Lugworm whole
	n	Avg	s.d.	n	Avg	s.d.	n	Avg	s.d.	n	Avg	s.d.	n	Avg	s.d.	N	Avg	s.d.	n	Avg	s.d.	n	Avg	s.d.
ALA	6	29.3	0.9	8	27.5	0.6	4	26.3	0.9	3	22.5	0.5	5	23.3	2.2	4	23.0	0.7	3	20.9	0.9	4	18.1	1.4
ASX	5	22.6	1.1	8	21.4	0.7	3	21.7	0.7	3	18.0	0.9	5	20.2	1.3	4	19.9	0.3	3	20.0	1.7	4	18.2	1.2
GLU	6	26.5	0.7	8	27.2	0.5	4	26.6	1.0	3	23.9	0.5	5	24.1	0.9	4	22.7	0.2	3	19.4	1.1	4	17.4	0.9
GLY	6	7.6	0.8	8	8.8	0.3	4	5.0	0.8	3	5.2	0.9	5	10.0	1.4	4	9.8	0.6	3	9.9	0.9	4	9.1	1.0
ILE	5	26.0	0.8	5	25.7	0.6	3	25.7	0.1	n.d.			5	18.3	1.8	4	19.1	0.7	3	16.7	0.7	2	16.2	1.4
LEU	6	26.3	0.7	8	26.9	0.8	4	25.8	0.7	3	21.9	1.1	5	19.8	1.9	4	19.1	0.8	3	17.7	0.5	5	16.7	0.8
LYS	6	4.0	0.8	8	2.8	0.3	4	5.0	0.4	n.d.			n.d.			4	8.0	0.8	3	3.6	0.5	2	2.1	0.2
OH-PRO	n.d.			8	22.9	0.9		n.d.		3	20.5	1.8	n.d.			n.d.			n.d.			n.d.
PHE	6	8.6	0.9	6	9.9	0.8	3	7.6	1.1	3	6.0	2.4	5	8.9	1.0	4	11.7	0.7	3	9.3	0.8	2	11.0	0.5
PRO	n.d.			6	27.6	0.7	n.d.			n.d.			n.d.			n.d.			n.d.			n.d.
SER	3	7.7	0.5	6	8.4	0.7	3	5.9	1.0	2	8.3	0.5	2	10.7	1.6	4	10.3	0.9	3	8.4	1.5	2	8.8	1.0
THR	5	-15.3	3.2	2	-22.2	0.0	3	-17.5	2.4		n.d.		1	-4.7		1	4.6		2	4.7	1.6	2	3.0	0.2
TYR	6	15.5	1.5	8	15.5	0.7	4	14.0	0.9	3	8.7	0.8	5	11.4	0.5	4	15.0	1.1	3	11.3	0.3	2	11.6	0.9
VAL	7	27.6	0.9	8	28.5	1.4	4	27.2	1.2	3	23.2	1.1	5	21.2	2.3	4	22.4	0.5	3	19.6	1.2	4	18.5	2.3

n.d. = not detected.

Avg and s.d. = average and standard deviation of n injections. ALA = n class="Chemical">alanine; ASP = aspartic acid; GLU = glutamic acid; GLY = glycine; ILE = isoleucine; LEU = leucine; LYS = lysine; OH-PRO = hydroxy-proline; PHE = phenylalanine; PRO = proline; SER = serine; THR = threonine; TYR = tyrosine; VAL = valine.

The cause for the observed large variability in the n class="Chemical">15N depletion of N-containing lipids is not clear. However, it is noticeable that the organisms with higher trophic levels, such as fish, show a larger depletion in 15N than those associated with lower trophic levels, such as lugworms. Indeed, when we calculate trophic levels for the different organisms based on the δ15N of glutamic acid and phenylalanine, using the equation of Chikaraishi et al. [29] for marine food webs, we observe an increasing depletion in 15N of the lipid extract relative to biomass with increasing trophic level (Fig 4). Further studies are needed to explore the cause for the large variability in the 15N depletion of N-containing lipids and whether the observed correlation with trophic level is causal.

Fig 4

Difference in δ15N of lipid extracts (TLE) and bulk biomass (Δδ15NTLE-bulk) plotted against trophic levels.

Δδ15NTLE-bulk data points represent averages of several individuals with the error bar reflecting the standard deviation of multiple individuals (S1 Table). Trophic level data points are plotted as averages (with error) of replicate analysis of amino acids of a single individual. Trophic levels were calculated according to Chikaraishi et al [29] using δ15N values of the amino acids phenylalanine and glutamic acid.

Δδ15Nn class="Disease">TLE-bulk data points represent averages of several individuals with the error bar reflecting the standard deviation of multiple individuals (S1 Table). Trophic level data points are plotted as averages (with error) of replicate analysis of amino acids of a single individual. Trophic levels were calculated according to Chikaraishi et al [29] using δ15N values of the amino acids phenylalanine and glutamic acid.

Conclusions

Our study shows that lipid extracts are generally highly depleted in n class="Chemical">15N compared to bulk biomass. The majority of lipids in the lipid extracts from biomass in this study consisted of glycerophospholipids with a phosphatidylcholine headgroup, a nitrogen-containing lipid which is commonly found in animal tissues. The 15N depletion in the lipids is likely a reflection of the δ15N of the biosynthetic nitrogen sources such as the amino acid serine, which is also depleted in 15N compared to bulk biomass. The depletion of δ15N of lipid extracts compared to bulk biomass seems to be higher with organisms of a higher trophic level, the reasons for which need to be further explored.

Stable nitrogen isotope values (δ15N), C:N ratios, and lipid content of bulk and residual biomass, and total lipid extracts of samples used in this study.

Samples highlighted in yellow indicate those used for amino acid analysis.

(XLSX)

Click here for additional data file.

Fatty acid composition, indicated by number of carbons and degree of unsaturation, of most abundant intact polar lipids (IPLs) in different aquatic species in this study.

n.d. = not detected. PC: n class="Chemical">Phosphatidylcholine; PE: Phosphatidylethanolamine; MMPE: monomethyl–PE; PI: Phosphatidylinositol; Plasmalogen: Fatty acid with vinyl linkage to glycerol backbone; Lyso: IPL where one fatty acid has been removed (lysed); Sphingo: Sphingobase. See Fig 3 for structures.

(DOCX)

Click here for additional data file.