The fatty acid trophic marker approach, i.e., the assimilation of fatty acids as entire molecule and transfer into consumer tissue with no or minor modification, is hampered by knowledge gaps in fatty acid metabolism of small soil invertebrates. Isotopologue profiling is provided as a valuable tool to disentangle trophic interactions.
Fatty acids (FAs) are useful biomarkers in food web ecology because they are typically assimilated as a complete molecule and transferred into consumer tissue with minor or no modification, allowing the dietary routing between different trophic levels. However, the FA trophic marker approach is still hampered by the limited knowledge in lipid metabolism of the soil fauna. This study used entirely labelled palmitic acid (13C16:0, 99 atom%) as a tracer in fatty acid metabolism pathways of two widespread soil Collembola, Protaphorura fimata and Heteromurus nitidus. In order to investigate the fate and metabolic modifications of this precursor, a method of isotopologue profiling is presented, performed by mass spectrometry using single ion monitoring. Moreover, the upstream laboratory feeding experiment is described, as well as the extraction and methylation of dominant lipid fractions (neutral lipids, phospholipids) and the related formula and calculations. Isotopologue profiling does not only yield the overall 13C enrichment in fatty acids derived from the 13C labeled precursor but also produces the pattern of isotopologues exceeding the mass of the parent ion (i.e., the FA molecular ion M+) of each labeled FA by one or more mass units (M+1, M+2, M+3, etc.). This knowledge allows conclusions on the ratio of dietary routing of an entirely consumed FA in comparison to de novo biosynthesis. The isotopologue profiling is suggested as a useful tool for evaluation of fatty acid metabolism in soil animals to disentangle trophic interactions.
In a cryptic habitat such as soil, trophic relationships are difficult to address and are further restricted by the small size of the fauna. The last decade has seen advances in biochemical ecology, particularly in the use of fatty acids as biomarkers for defining feeding strategies of the soil fauna under field conditions1,2,3. This is based on the fact that fatty acids from resources can be incorporated in consumer tissue as entire molecules, a process termed dietary routing4. Transfer of fatty acids has been reported over three trophic levels, i.e., from fungi to nematodes to Collembola5. Recently, the predatory fauna was considered6,7 and the first reviews on fatty acids as trophic markers in soil food webs have been published8,9.
More detailed information on trophic interactions is attained by fatty acid stable isotope probing (FA-SIP). The determination of 13C/12C ratios in fatty acids in diets and consumers can ascribe binary links and estimate the associated carbon flow, and has been employed in terrestrial, fresh water, and marine food webs10,11,12,13. The basic assumption is that dietary routed fatty acids are not subject to enzymatic processes; therefore, their 13C signal, i.e., the 13C/12C ratio of the fatty acid, in the consumer is similar to that in the diet1. However, a gradual depletion of the 13C signature up the food chain has been reported in aquatic systems, thereby hindering broad application of FA-SIP in trophic studies14,15,16. Moreover, knowledge in the lipid metabolism in most invertebrates in terrestrial food webs is still limited.
An understanding of the lipid metabolism pathways in consumers is essential for the usage of trophic marker fatty acids as means for the determination of the quantitative carbon flow in food web ecology. With this in mind, 13C-isotopologue profiling, which in principle can be applied for investigation of the carbon metabolism of any biological system17, is a promising method. Following the introduction of a 13C-labelled carbon substrate, the distribution of the 13C in the metabolic network is traceable since the generated metabolic products in the consumer show a specific isotopologue distribution. This can be assessed by quantitative nuclear metabolic resonance spectroscopy18,19 or mass spectrometry20,21, with the latter favored in biological samples with low biomass due to its higher sensitivity.
Although isotopologue profiling has been successfully applied to amino acids and provided insight into the in vivo carbon metabolism of bacterial pathogens17,22,23, its implementation in fatty acids has lagged behind. The first detailed analysis on the fate of a stable isotope labelled precursor fatty acid, its dietary routing or degradation via β-oxidation, in soil invertebrate consumers, was recently performed by Menzel et al.24. Here, the methodological basics for incorporation experiments with 13C labelled fatty acids followed by isotopologue analysis of key descendants in frequent soil invertebrates, the Collembola, are provided. These microarthropods are a good model group as they form important components of the soil food web and are well investigated for their trophic marker fatty acids8,25.
An understanding of the lipid metabolism pathways in consumers is essential for the usage of trophic marker fatty acids as means for the determination of the quantitative carbon flow in food web ecology. The present protocol gives the design and set up for a laboratory feeding experiment, and the biochemical procedures for extraction and methylation of dominant lipids fractions (neutral lipids, phospholipids) from Collembola. It demonstrates how the isotopologue composition of fatty acids is analyzed by mass spectrometry and describes the related formula and calculations. This procedure results in: (i) the ratios of isotopologues exceeding the mass of the parent ion (i.e., the fatty acid molecular ion M+) by one or more mass units (M+1, M+2, M+3, etc.) and (ii) the overall 13C enrichment in fatty acids derived from the 13C labelled precursor. Although used for Collembola, this approach can generally be applied to any other predator-prey interaction on the premise that these are culturable in sufficient quantity under controlled conditions to ensure a successful label uptake and subsequent verification.
The described protocol does not fall under the competence of Animal Ethics. However, when people adapt the described protocols to higher animals, take care that the institutional Animal Ethics committee approved the protocol for animal handling.
1. Cultivation of Animals
NOTE: All explained experimental steps are based on well-established protocols26,27,28. Biotests in the laboratory need a continuous supply of easily culturable organisms. Here, the Collembola species Protaphorura fimata (Gisin, 1952) and Hetermurus nitidus (Templeton, 1835) have been used. Both species are simple to maintain as productive laboratory cultures fed with baker's yeast.
Figure 1: Cultivation of Collembola. (A) Microcosm filled with breeding substrate, a dried mixture of plaster of Paris, activated charcoal and distilled water. (B) and (C) Representative specimen of a Protaphorura fimata culture; note the small nuggets of dry Baker's yeast used as the food source and also as holes in the breeding substrate (black arrow) (B) as well as two eggs (white arrow) (C). Please click here to view a larger version of this figure.
2. Labelling Diet, Harvest, and Sample Handling
3. Lipid Extraction from Animal Tissue and Methanolysis
4. Quantification of Fatty Acids by GC-FID
5. 13C Analysis by Isotopologue Profiling
6. Calculations of 13C Enrichment
Fresh weight and lipid content of Collembola
In the course of the described experiment, the content in NLFAs and PLFAs did not change significantly over time, whereas the fresh weight of specimens increased slightly but not significantly24. Both parameters indicate a good level of physical fitness of the Collembola specimens. Be aware to investigate Collembola's fresh weight and lipid content throughout the experiment corresponding to the sampling days for fatty acid and isotope analysis. Note that a loss of weight and/or a decrease of lipid content during the experimental period indicate a reduced fitness of the test organism and the derived data drastically lose significance.
Isotopologue detection
The special significance of isotopologue profiling projects lies in the individual detection and quantification of both the molecular ion and all isotopologues in a specific fatty acid. For example, in case of 16:0 the ions range from 270, i.e., molecular ion M+ (C skeleton mainly composed of 12C atoms) to 286 (isotopologue M+16 – carbon chain completely substituted with the heavy 13C). Figure 2A presents a SIM-MS spectrum of pure 16:0, palmitic acid. The natural abundance of 13C in this compound becomes detectable by the presence of M+1 and M+2 isotopologues in addition to the molecular ion (M+). By comparison, Figure 2B shows the spectrogram of the entirely labeled palmitic acid (99 atom% 13C) used in this study. Here, the sole occurrence of the ion M+16 reflects the high purity of this synthetically labeled compound.
Figure 2: Representative SIM scans of pure (A) unlabeled and (B) entirely labeled palmitic acid. Note the natural abundance of M+1 and M+2 isotopologues in addition to the molecular ion (M+) in the unlabeled C 16:0 (A) but the exclusive presence of the M+16 isotopologue in the entirely labeled fatty acid (99 atom% 13C). Please click here to view a larger version of this figure.
Elongation/desaturation versus de novo synthesis
After harvesting the test organisms, lipid extraction and derivatization, the generated FAMEs are ready to be analyzed by GC and SIM-MS. By examination of isotopologue profiling data, one can now distinguish between de novo synthesis and desaturation/elongation events of the entire labeled marker fatty acid. The former can take place, at least partially, based on 13C Acetyl-CoA as product of the degradation of the labelled precursor via beta-oxidation. Figure 3 shows representative examples of SIM scans derived from four dominant fatty acids in the PLFA fraction of H. nitidus at the first sampling day after labelling. The assimilation of the labelled marker molecule into palmitic acid becomes visible by the abundance of the M+16 – isotopologue, as it represents the entirely labeled fatty acid (Figure 3A – 16:0). Further, desaturation of 16:0 (Figure 3B – 16:1ω7) or elongation plus desaturation (Figure 3C – 18:1ω9, Figure 3D – 20:4ω6) can be assigned. Thereby, the detection of isotopologues bigger than M+16, e.g. M+17 and M+18 demonstrate chain elongation of the labelled precursor 16:0 by the use of 13C2-fragments. De novo biosynthesis based on 13C Acetyl-CoA of fatty acids is indicated if isotopologues close to the molecular ion, i.e., M+1 to M+2 get significantly more frequent than in the control represented by day 0 sampling.
Figure 3: Representative SIM (selected ion monitoring) scans of fragment ions in fatty acids of Heteromurus nitidus (PLFA fraction, 1 day after labelling). (A) Palmitic acid (16:0) (B) Palmitoleic acid (16:1ω7) (C) Oleic acid (18:1ω9), and (D) Arachidonic acid (20:4ω6). Please click here to view a larger version of this figure.
Further illustrating this, Figure 4 compares the initial isotopologue pattern of the five most heavily labeled fatty acids (analyzed both in PLFAs and NLFAs); compared are the data from the Collembola species P. fimata on day 0 and day 1. On day 1, the 13C signal in C16 and C18 FAs was almost entirely based on the occurrence of the M+16molecule ion, resulting from entirely labeled 13C 16:0, supplemented to the food. As mentioned above, the observed traces of M+18 of 18:0 and 18:1ω9 suggest the introduction of a minor portion of 13C Acetyl-CoA via chain elongation from 16:0 to 18:0. For 18:1ω9 this was more pronounced in the PLFA faction. Such 13C acetyl-CoA derives from β-oxidation of labeled 13C 16:0 molecules. Incorporation of 13C acetyl-CoA also occurred in the elongation step from chain length C18 to C20, as assigned by the 20:4ω6 of the PLFA samples. Thereby, the 13C enrichment in M+16 and M+18 significantly increased with time until day 14 (Figure 4). Moreover, M+2 of this PLFA 20:4ω6 at day 14 increased when compared to day 0 or 1. This 13C allocation in the different isotopologues indicates that the formation of 20:4ω6 is based both on de novo synthesis with inclusion of acetyl-CoA labeled with 13C or at natural 13C abundance and the elongation/desaturation of the entire 13C16:0 precursor molecule. In contrast, distinctly labeled 20:4ω6 did not appear in the NLFA fraction (Figure 4).
Figure 4: Isotopologue pattern of the 13C incorporation by SIM (selected ion monitoring). Presented are the relative portions of the ions M+1 to M+3 and M+16 to M+20 of the five most heavily labeled NLFAs and PLFAs from Protaphorura fimata, estimated on day 0 and day 1. Only for 20:4ω6 the day 14 data are presented in addition. P<0.05 (Dunnett's test with day 0 data used as reference). This figure is a reprint of Figure 5 published in Menzel et al.24 Please click here to view a larger version of this figure.
Isotopologue profiling
A detailed analysis of the quantitative aspects in 13C distribution in FAs needs cutting-edge technology to assign carbon partitioning in food webs. The present work employed isotopologue profiling to assess the 13C/12C ratios in common FA biomarkers for tropic interactions. This method is well established for amino acid analysis by liquid chromatography (LC-MS) and was applied for investigations of carbon metabolism in pathogenic bacteria17,23. Only recently, was isotopologue profiling further developed as a tool to study dietary routing of lipids in small soil invertebrates by Menzel et al.24 Since this method is not hampered by the "carry over" of the 13C signal, it is therefore advantageous for the evaluation of organic tissue highly enriched with a specific isotope compared to GC-C-IRMS32.
With this approach, the mass-to-charge (m/z) values in the molecular ion and its isotopologues is determined by GC-MS using EI and SIM mode. In comparison to the full scan, the sensitivity in SIM can be increased by about factor 100, which even outstrips the GC-FID quantification range33. Further, by acquisition of a specified group of ions in a given time window, the analyte is measured reliably despite a background of co-eluting peaks30. Thurnhofer and Vetter34 proved that SIM is suitable for both FA identification and quantification in food samples even with low amounts of lipids available.
When using SIM for 13C isotopologue profiling at least two methodological constraints have to be considered. Firstly, molecular ions are of relatively low abundance and their isotopologues are even less. With a dirty ion source, these become less abundant and as a consequence the SIM method will only be semi-quantitative33. Secondly, the increased occurrence of 13C at one or more carbon positions of the FA decreases the 13C natural abundance of associated carbon atoms. These changes with labelling are non-linear, a phenomenon called "skew" in natural abundance, which can lead to underestimation of enrichment at these masses20. This skew effect must be accounted for either by using a correction matrix based on FAME standards (see Fernandez et al.) or by including natural controls, i.e., FAs from experimental organisms without label. The latter was performed in the present work comparing the same FA in labelled (day 1 and later) and non-labelled Collembola (day 0).
Dietary routing and 13C flux
While the methodological advantages of isotopologue profiling for 13C determination in FAs are obvious, in vivo metabolism of organisms can be hindered by biological means. Isotopic and metabolic steady state are difficult to achieve when, e.g., diet usage or the physiological status of consumers are not known, hampering obtaining absolute carbon fluxes. However, by comparing 13C isotopologues profiles of FAs in the diet and the consumer, the relative proportion of FA assimilation and metabolic pathways can be obtained. This approach allows following the specific fate of a marker FA based on their 13C label by tracing its dietary routing. MS-SIM demonstrated that the marker FA was mainly routed into the NFLA fractions of both Collembola species analyzed24. Moreover, the marker FA was stored in the neutral lipids with only minor or no modifications. Only elongation of the dietary 13C16:0 to 18:0 and desaturation to the monoenoic equivalents were detected (see Figure 4). The presence of the ion M+16 in the C16 and C18 FAs of day 1 samples confirms that these FAs are descendants of the supplemented entirely labeled palmitic acid. This finding underpins earlier studies that were based on the FA pattern of consumer and diet solely5,35 and suggests that marker FAs are incorporated as entire molecules into consumer tissue (for a review, see Ruess and Chamberlain8).
Meanwhile, the isotopologue 13C pattern of PLFAs revealed that the observed increase in 13C20:46 over time is due to both elongation/desaturation of entirely labeled precursor FAs and de novo synthesis with the latter also including 13C labeled acetyl-CoA as indicated by the presence of M+18 ion traces in C20 and also C18 FAs (see Figure 4). Figure 5 summarizes the dietary routing of supplemented 13C16:0 within PLFAs and NLFAs of the observed Collembola species. Interestingly, only 13C palmitic acid of PLFAs is subject to stronger changes, whereas the metabolic transformation in the NLFAs remained limited to only one elongation and two desaturation steps.
Figure 5: Fate of dietary 13C palmitic acid in consumer lipids. The flow-charts show the proposed fate of incorporated 13C16:0 in the investigated Collembola species. Arrows symbolize the different levels of metabolic turnover with different thickness. Please click here to view a larger version of this figure.
To summarize, isotopologue profiling by MS-SIM clearly shows that FA use as trophic marker is a robust and reliable method in the ecology of food webs. Although the experimental design did not allow the precise quantification of the 13C flux (e.g., as atom% 13C of the supplemented FA) for the different steps in the pathways of FA metabolism, isotopologue profiling is a useful tool for all ecological studies using fatty acids as trophic markers. A goal for the future is to obtain a quantitative consumer FA carbon budget by using closed microcosm systems and feeding consumers defined amounts of labeled precursors of fatty acids.
The financial support of R. Menzel and L. Ruess by the Deutsche Forschungsgemeinschaft (RU RU780/11-1) is gratefully acknowledged. R. Nehring was funded by RU 780/10-1. Finally, we are extremely thankful to Dr. Hazel Ruvimbo Maboreke for proofreading our manuscript.
neoLab-Round jars | neoLab | 2-1506 | 69 x 40 mm, 10 pacs/pack |
Charcoal activated | Carl Roth | X865.1 | p.a., powder, CAS No. 7440-44-0 |
Alabaster Dental | RÖHRICH-GIPSE | — | http://www.roehrich-gipse.de/dentalgipse.php |
Chloroform | Carl Roth | 7331.1 | HPLC ≥ 99,9 % |
Methanol | Carl Roth | P717.1 | HPLC ≥ 99,9 % |
Hexan | Carl Roth | 7339.1 | HPLC ≥ 98 % |
tert-Butyl methyl ether (MTBE) | Carl Roth | T175.1 | HPLC ≥ 99,5 % |
Aceton | Carl Roth | 7328.2 | HPLC ≥ 99,9 % |
NaOH | Carl Roth | 6771.1 | p.a. ≥99 %, in pellets |
di-Natriumhydrogenphosphat | Carl Roth | P030.1 | p.a. ≥99 % , water free |
Na-dihydrogenphosphat Dihydrat | Carl Roth | T879.1 | p.a. ≥99 % |
Hypochloric acid (6 N) | VWR International | 26,115,000 | AVS TITRINORM vol. solution |
Bond Elut (Columns) | Agilent Tech. | 14102037 | HF Bond Elut-SI, 500 mg, 3 mL, 50/PK |
Präparatengläser Duran | Glasgerätebau Ochs | 135215 | Ø 16 x 100 mm, plus screw cap with handy knurl and integrated PTFE/silicone gasket |
Supelco 37 Component FAME Mix | Sigma-Aldrich | 47885-U Supelco | 10 mg/mL in methylene chloride, analytical standard |
FlowMesh | Carl Roth | 2796.1 | Polypropylene mesh, approximately 0.3 mm thick, with 1 mm strand spacing |
Bacterial Acid Methyl Ester (BAME) Mix | Sigma-Aldrich | 47080-U Supelco | 10 mg/mL in methyl caproate, analytical standard |
Methyl nonadecanoate | Sigma-Aldrich | 74208 | analytical standard ≥ 98.0 % |
Hexadecanoic acid-1-13C (Palmitic) | Larodan Fine Chemicals | 78-1600 | GC ≥ 98.0 % (13C: 99.0 %) |
RVC 2-25 CDplus | Martin Christ Gefrier-trocknungsanlagen | Compact benchtop midi concentrator | |
Alpha 2-4 LDplus | Martin Christ Gefrier-trocknungsanlagen | Drying manifold | |
MZ 2C NT | Vacuubrand GMBH | Vacuum pump | |
Roto-Shake Genie | Scientific Industries | Combined rocking and rotating device | |
XP64 Micro Comparator | Mettler Toledo | Super high precision balance | |
GC-System 7890A | Agilent Tech. | Gas chromatograph | |
7000 GC/MS Triple Quad | Agilent Tech. | Triple Quad mass spectrometer | |
7683B Series Injector | Agilent Tech. | Sample injector | |
Heraeus Multifuge 3SR+ | Thermo Scientific | Centrifuge with 10 ml tube rotor |