Mixed-layer lipidomes suggest offshore transport of energy-rich and essential lipids by cyclonic eddies
Introduction
Marine microbes are key players in global biogeochemical cycles including the carbon and nutrient cycles. About half of the global primary production occurs in the surface ocean via photosynthesis by phytoplankton1, which provides the basis for the transfer of organic matter and energy to marine ecosystems, supporting most marine food webs including fisheries production, and plays a critical role in ocean carbon sequestration2. Changes in the marine environment resulting from physical and chemical processes impact phytoplankton growth by providing access to light and nutrients. For example, upwelling of nutrient-rich subsurface water sustains high productivity in the ocean’s eastern boundary currents, while in the surface waters of the open-ocean subtropical gyres, productivity is limited by the supply of nutrients3. Besides these large-scale oceanographic features driving productivity patterns in the global ocean, smaller scale processes affecting areas of about 10–100 km, i.e. mesoscale processes, are important for shaping anomalies in productivity. One of the most prominent mesoscale features in the ocean are eddies, which are ubiquitous swirling currents that dominate local biogeochemical processes4. Eddies modulate light and nutrient conditions in the surface affecting abundance and activity of both phytoplankton and heterotrophic microbial communities5,6,7,8. Three types of eddies have been described: Cyclonic, anticyclonic and anticyclonic mode water eddies9. The three types show substantially different effects on community structure and food-web dynamics, with cyclones representing productive oases for marine life. By shoaling the pycnocline, cyclonic eddies can be sites of enhanced nutrient supply10,11 and light availability12 stimulating phytoplankton production. The stimulated phytoplankton activity has further been shown to attract zooplankton13 and enhance carbon export fluxes14,15. In all of the eastern boundary upwelling systems (EBUS), mesoscale eddies are frequently observed16,17. Besides affecting biogeochemical cycling of carbon and nutrients locally, eddies formed in EBUS might contribute to the redistribution of nutrients and carbon from the coast to the open oligotrophic ocean due to their westward propagation17,18,19,20. Based on satellite-derived measurements of particulate organic carbon (POC) in the California current system, trapping and transporting coastal water offshore by mesoscale eddies has been estimated to result in a significant offshore POC enrichment17. The connectivity between coastal upwelling and open-ocean regimes may thus contribute to sustaining the biological productivity in the oligotrophic ocean, which has been a long-standing conundrum21,22,23.
To date, biogeochemical studies of mesoscale eddies have focused primarily on bulk measurements of, for example, POC and dissolved organic carbon or primary production, while the effects of mesoscale eddies on changes in the biochemical composition of organic matter as well as the biogeochemical consequences of these for the ocean’s carbon cycle remain poorly constrained. One important component of the oceanic organic matter composition are microbial lipids24, which can make up to 20% of the surface ocean particulate organic matter pool25 and are produced and used by organisms from all three domains of life for energy storage, membrane structure, signaling and electron transport. Lipid biomarkers have mainly been used to study specific biological and chemical processes in the ocean, such as nutrient limitation26, but recent advances in lipidomics analysis using high resolution mass spectrometry27 enable the detection of hundreds to thousands of individual lipids in environmental samples. For example, a recent study by Holm et al.28 provided evidence for global-scale lipidome remodeling in surface ocean phytoplankton as a function of temperature. Lipids also carry coarse chemotaxonomic information that allows to distinguish the dominant microbial groups in the marine water column, such as different phytoplankton, heterotrophic bacteria, and archaea29,30,31. The observed effects of eddies on the distribution of microbes5,6 can thus be expected to be reflected in variations in the lipidome. Additionally, the lipidome of phytoplankton is one of the primary determinants of their nutritional quality32 and alterations in the types and amounts of lipids produced due to changes induced by mesoscale eddies will affect the activity and population dynamics of organisms at higher trophic levels and ultimately carbon export from the surface to the deep ocean.
Here, we investigated the surface microbial lipidome in and around a cyclonic eddy formed in the highly productive Mauritania EBUS in the Eastern Tropical North Atlantic (Fig. 1, see Supplementary Methods including Supplementary Figs. 1–3 for detailed eddy identification) where mesoscale eddies commonly occur16,33. Surface mixed-layer samples, where lipid abundances are typically highest25,30,34, were taken within the core, the edge of the core, and the periphery of the eddy. Samples outside the eddy were collected along a 18°N transect between the Cape Verde islands and the coast where the samples to the West of the eddy represent open-ocean sites. In previous studies, we showed that bulk microbial metabolic activities and abundances, and bulk pools of organic matter (POC and dissolved organic carbon) were significantly enhanced in the eddy studied with strong horizontal variability within the surface of the eddy7,35,36; this variability was covered by our sampling scheme for lipid analysis. The stimulated microbial growth in this cyclone might provide fresh organic matter to the open ocean during westward propagation. We sought to further elucidate microbial carbon cycling in this system through the lens of lipidomics and to examine the factors that determine ocean lipid composition at the mesoscale. Since mesoscale eddy activity is predicted to change as a result of climate change37, it is crucial to know how these features affect the lipid composition in order to better assess consequences for food web dynamics and carbon export.

Stations for lipidomics analysis collected during RV Meteor cruise M156 in the Eastern Tropical North Atlantic upwelling system off Mauritania. The gridded bathymetry map was obtained from http://www.gebco.net (General Bathymetric Chart of the Oceans, GEBCO 2024 Grid). The background in the zoom in shows the absolute dynamic topography (ADT) at the time of sampling obtained from https://www.aviso.altimetry.fr, last access: 23 July 2024, with the direction and speed of surface water geostrophic currents shown as arrows. Black circle indicates the eddy core and the dashed line indicates potential edge of the eddy-influenced area (eddy periphery). Colors of sampling stations indicate position with red colors representing eddy stations (red: eddy core, light red: edge of eddy core, dark red: eddy periphery) and blue colors outside stations (dark blue: coastal station, light blue: open-ocean stations). Station E3 appeared to be located just beyond the eddy periphery, but was shown to be within the eddy influenced area7 and thus assigned to be within the eddy periphery. While station E4 was located within the coastal-boundary current region, our lipid analysis revealed a high similarity to the open-ocean stations rather than the coastal station E5 and was assigned as open-ocean station.
Results and discussion
Lipidome variations in and around a cyclonic eddy formed in the EBUS off Mauritania
In total, almost 900 individual lipid species were identified in the mixed layer in and around the cyclonic eddy sampled during the RV Meteor cruise M156 (see Supplementary Table 1 for mixed-layer depths). Hierarchical clustering analysis using Euclidian distances of the lipidome (see Supplementary Figs. 4–6 for chemical structures and mass spectral identification) from all mixed-layer samples showed two main clusters of samples (Fig. 2). The first cluster contains the samples from all open-ocean stations and the coastal stations, while the second cluster only contained samples from inside the eddy as well as station E3, which was just beyond the eddy periphery and has been shown to be influenced by the cyclonic eddy in our earlier work7. The stations outside of the eddy further clustered into two subclusters: coastal and open-ocean stations. While station E4 was located within the coastal-boundary current region, it clustered together with the open-ocean stations rather than the coastal station E5. Within the eddy-influenced area, there are three subclusters: samples from the core of the eddy, the edge of the core, and the periphery. Thus, the eddy carries a distinct mixed-layer lipid signature with variability depending on the position within the eddy-influenced area.

Heatmap of log-transformed lipid concentrations (893 individual lipids) in mixed-layer samples collected during M156 in and around a cyclonic eddy formed in the upwelling system off Mauritania (see Supplementary Table 1 for mixed-layer depth of each station). Hierarchical clustering analysis of columns (samples) and rows (lipids) was performed using Euclidean distances. The samples clustered according to their relative location within the cyclonic eddy and the ambient waters. The locations of the stations are shown in Fig. 1. Each of the lipid (sub)clusters contained an array of molecularly diverse lipids from the four main categories: 1. intact polar lipids, 2. respiratory quinones, 3. pigments, and 4. triacylglycerols as described in main text. The cluster analysis with fully annotated lipids is provided in Supplementary Fig. 14.
Sources of lipid groups in the identified clusters
Our lipidome analysis covered lipids from four main categories: 1. intact polar lipids (IPLs), which are structural membrane components, 2. respiratory quinones, which are lipid-soluble electron and proton carriers, 3. pigments, which are involved in photosynthesis, and 4. triacylglycerols (TAGs), which are energy storage molecules. Most of these lipid groups showed significant differences in concentration between the eddy-influenced stations and stations outside of the eddy (one-way ANOVA, p < 0.05, Fig. 3). IPLs, quinones and pigments were most abundant in the eddy-influenced stations with the highest concentrations in the core and at the edge of the core. For example, concentrations of IPLs in the core of the eddy were found to be on average three times higher than in the open ocean, and pigment concentrations were even nine times higher, indicating enrichments in microbial and in particular photosynthetic biomass (Fig. 3a). TAGs were also more abundant in the eddy, but showed particularly high concentrations in samples from the eddy periphery (Fig. 3l). Here, TAGs were by far the most abundant lipid group in the eddy-influenced stations where total lipids contributed more than 6% to the POC pool in the eddy periphery, whereas the contribution was only ~2% for the coastal open-ocean stations (Supplementary Fig. 7). The large spread in TAG concentrations in the periphery is due to the extremely high concentrations in the sample from station EDZ1, which showed six times higher concentrations than the sample from the same depth from the other station in the eddy periphery (station E3, Supplementary Data File). TAGs are used by mainly by eukaryotic phytoplankton as energy reserves and TAG accumulation often occurs in growth limiting conditions such as dissolved inorganic nutrient (nitrogen) deprivations38. However, nutrient limitation was not observed in the periphery of the sampled cyclonic eddy as we could show based on dissolved inorganic nutrient profiles in our previous study7. Responses to the daily light cycle can also influence the TAG content in phytoplankton39 and some of the observed variability in TAG concentration may be related to the time of sampling; station EDZ-1 was sampled at the end of the day when TAG concentrations are highest, while station E3 was sampled at the end of the night when TAG concentrations are lowest (see Supplementary Table 1 for sampling information). However, all other stations within the eddy were also sampled at the end of the night or the early morning when daily TAG concentrations in phytoplankton are low39. We thus attribute the high TAG concentrations in the periphery of the eddy mainly to eukaryotic phytoplankton biomass that is particularly rich in TAGs. This is in line with our previous observation of the presence of large (>20 µm) eukaryotic phytoplankton in the periphery of this eddy based on combined chlorophyll a (Chl a) measurements, 14C-based primary production rates and flow cytometry analyses7.

Average concentration of total intact polar lipids (a), chloroplast lipids (b), cell membrane lipids (c), archaeal lipids (d), total quinones (e), photosynthetic quinones (f), ubiquinones (g), menaquinones (h), total pigments (i), chloropigments (j), carotenoid pigments (k), and triacylglycerols (l) in the five sub-clusters of sampling stations identified by the hierarchical clustering analysis of surface mixed-layer samples collected during M156 off Mauritania. Boxes represent the interquartile ranges, the line within the box the median, and the whiskers extend to the minimum and maximum values within 1.5 times the interquartile range. Outliers are shown as circles. Asterisks indicate variables significantly different (one-way ANOVA, p < 0.05) from coastal (dark blue) and open-ocean (light blue) stations, respectively. The chloroplast lipids are the sum of glycolipids with diacylglycerol side chain, including the galactolipids mono- and digalactosyldiacylglycerol and the sulfolipid sulfoquinovosyldiacylglycerol. The cell membrane lipids include all phospho- and betaine lipids. The archaeal lipids contain all core and intact polar diether and tetraether lipids. The respiratory quinones are separated into photosynthetic quinones (plastoquinone and vitamin K1), ubiquinones and menaquinones. The pigments are divided into chloropigments, which is the sum of chlorophylls and related compounds (porphyrin ring-containing molecules) and carotenoid pigments, which include all xanthophylls (e.g., fucoxanthin and zeaxanthin) and β-carotene. For an overview of the chemical structures, the reader is referred to Supplementary Fig. 4.
The main lipid groups were further broken down into groups according to their sources and/or function. This revealed that the higher IPL concentrations in the eddy are mainly due to higher concentrations of chloroplast membrane lipids, i.e., the galactolipids monogalactosyl- (MGDG) and digalactosydiacylglycerol and the sulfolipid sulfoquinovosyldiacylglycerol (see Fig. 3b), which are of phytoplankton origin including both eukaryotic phytoplankton and cyanobacteria30. Similarly, within the pigments, we observed significantly higher concentrations of chloropigments and carotenoid pigments in the eddy-influenced stations compared to stations outside of the eddy (Fig. 3j, k), which are used by phytoplankton to perform photosynthesis. In contrast, cell membrane lipids, which include all phospholipids as well as betaine lipids (see Supplementary Fig. 4 for chemical structures), were only significantly higher in concentration in the eddy compared to the open ocean, but similar compared to the coastal samples (Fig. 3c). These lipids primarily originate from eukaryotic phytoplankton and to some extent of heterotrophic bacteria in surface ocean samples30. The higher concentration of chloroplast lipids together with the invariant concentration of cell membrane lipids in the eddy compared to the coastal station suggests higher abundances of cyanobacterial cells relative to eukaryotic phytoplankton in the core region of the cyclonic eddy, which is in line with our previous flow cytometry analysis7. This is further supported by the significantly higher concentration of photosynthetic quinones [plastoquinone and vitamin K1] in the core of the eddy compared to the other stations and the invariant ubiquinone (UQ) concentration (Fig. 3f, g) between those regions. While in eukaryotic algae UQ and photosynthetic quinones function as electron transporters in the electron transport chain of oxygenic photosynthesis and the aerobic respiratory chain, respectively, cyanobacteria only contain photosynthetic quinones and use plastoquinone in the electron transport pathways of both photosynthesis and respiration40. Higher abundances of cyanobacteria in the core of the eddy also explain the relatively low TAG, but high IPL and pigment concentrations in this part of the eddy, because naturally occurring cyanobacteria are not known to accumulate notable amounts of TAGs38. In several cultures of oceanic cyanobacteria that dominate in the study area, such as Synechococcus and Prochlorococcus41, even no TAGs were detected39, consistent with the lack of TAG biosynthesis in these groups42.
We additionally observed menaquinones (MKs) to have a similar concentration range as the photosynthetic quinones and UQs (Fig. 3h) and, similar to the photosynthetic quinones, MKs were significantly more abundant in the core region of the eddy compared to the coastal and open-ocean stations. The detected MKs were all fully unsaturated, i.e., they contained one double bond per isoprenoid unit and are thus of bacterial origin40. Typically, bacterial MKs are involved in anaerobic metabolism due to their low redox potential, while they cannot operate efficiently under aerobic conditions40. Accordingly, in the marine water column, bacterial MKs have so far only been found in the anoxic waters of the Black Sea43. In gram positive bacteria, however, MK biosynthesis, which is evolutionary the most ancient pathway of quinone biosynthesis44, is conserved45. Gram positive bacteria, in particular Actinobacteria, have been shown to be abundant in the Mauritania upwelling region41 making them a likely source of MKs that are particularly abundant in the core of the eddy.
Archaeal lipid concentrations were comparatively low, but showed highest values in the core and at the edge of the core of the eddy (Fig. 3d). The majority of the detected lipids were glycerol dibiphytanyl glycerol tetraethers (GDGTs) including intact polar and core crenarchaeol derivatives (see Supplementary Data File for all identified lipids), which are primarily derived from Nitrososphaerota (syn. Thaumarchaeota) in the marine water column46. For one station from the edge of the eddy core (station EDM5E), we obtained a depth profile down to 800 m, which showed that archaeal lipids had highest concentrations near the surface (20 m) within the surface mixed layer (Supplementary Fig. 8). This indicates that archaeal biomass is relatively enriched in the surface mixed layer at this station. Since the detected archaeal lipids mainly included core lipids and the intact polar lipids monohexose GDGTs (see Supplementary Data File), the differentiation between living and non-living archaeal biomass is difficult. It has been shown that archaeal IPLs, unlike bacterial and eukaryotic IPLs, may persist rapid degradation47,48,49. However, we have filtered only small volumes (1 L) of seawater and other archaeal IPLs might have simply been below the limit of detection. The dominance of monohexose GDGTs has further been observed in the marine water column from other regions and Nitrososphaerota were always suggested to be their dominant source29,31. Together with the distinct peak in concentration this points at least in part to the in-situ production of archaeal lipids in the surface mixed layer in and around the core of the cyclonic eddy. This was surprising to us, because Nitrososphaerota and their lipids are typically found to peak deeper in the water column50,51. The uplift of isopycnals in the cyclonic eddy brought nutrients including ammonia, which is the main substrate for Nitrososphaerota for energy generation52, into the mixed layer of the core of the eddy7 (Supplementary Fig. 9) and might explain their relatively high abundance at the shallow depth in our samples53.
Mixed-layer-integrated concentrations of specific lipid biomarkers
The distribution of individual biomarker lipids integrated over the surface mixed layer provided further insights into the distribution of microbes within the eddy-influenced area and manifests the observations from the lipid group concentrations. The distribution of Chl a (Fig. 4a), a generic biomarker for oxygenic photosynthetic organisms, clearly showed enhanced values in the eddy-influenced stations, but different groups of the phytoplankton community contributed differently to the Chl a signal in the different regions within the eddy as indicated by specific biomarkers. Long-chain, polyunsaturated chloroplast (MGDG 38:6, a monogalactosyldiacylglycerol with 38 acyl carbons and six acyl double bonds) and cell membrane (PC 38:8, a phosphocholine diacylglycerol with 38 acyl carbons and eight acyl double bonds) as well as energy storage (TAG 54:10, a triacylglycerol with 54 acyl carbon and ten acyl double bonds) lipid biomarkers, all indicative of eukaryotic phytoplankton such as diatoms and dinoflagellates54, were all particularly abundant in the stations from the periphery of the eddy (Fig. 4b–d). The high concentration of the carotenoid pigment fucoxanthin (Supplementary Fig. 7f) indicates that diatoms were the dominant eukaryotic phytoplankton group in this area of the eddy, which is in line with earlier studies55,56. In contrast, IPLs and pigments diagnostic for cyanobacteria, such as MGDG 30:130 and zeaxanthin57, were elevated in the core and at the edge of the core of the eddy (Fig. 4e, f). As expected, divinyl Chl a, a unique biomarker for the cyanobacterium Prochlorococcus58, occurred mainly in the open-ocean stations (Fig. 4g), consistent with their predominance in (sub)tropical oligotrophic waters59. This suggests that other cyanobacteria, e.g., Synechococcus, likely dominated the photosynthetic biomass in the core of the eddy, which is consistent with our previous flow cytometry analysis7 and which has also been shown for other cyclonic eddies60. Similarly to the cyanobacterial lipids in the core, a short-chain, monounsaturated phosphatidylethanolamine with 30 acyl carbons (PE 30:1), which is derived primarily from heterotrophic bacteria in marine samples30, was dominant in the surface mixed layer of the core of the eddy (Fig. 4h). The increased concentrations of MK7:7 (Fig. 4i), which is likely derived mainly from gram positive bacteria45 that may be involved in the degradation of complex polymers61, is another indicator for abundant heterotrophic bacterial biomass and matches heterotrophic bacterial cell abundances quantified by flow cytometry from our previous study of this eddy7. The mixed-layer-integrated concentration of monohexose crenarchaeol (Fig. 4j), a diagnostic biomarker for ammonia oxidizing Nitrososphaerota62, confirms their predominance in the core of the eddy core as outlined above. An example of a biomarker for the open-ocean stations is 19’Hexanoyloxy-4-ketofucoxanthin (Fig. 4k), which is indicative of haptophytes, in particular cocclithophores63. This is consistent with the general observation of increasing coccolithophore abundance from coastal to open-ocean regimes64.

Concentrations of individual lipids representative of different microbial groups integrated over the surface mixed layer in samples collected in and around a cyclonic eddy off Mauritania. a chlorophyll a (Chl a) as biomarker for oxygenic photosynthetic organisms, b–d polyunsaturated fatty acid-containing intact polar and energy storage lipids indicative of eukaryotic algae, e, f a monounsaturated glycolipid (MGDG) and the carotenoid pigment zeaxanthin indicative of cyanobacteria, g divinyl (DV) Chl a diagnostic for Prochlorococcus, h a monounsaturated phospholipid (PG) indicative of heterotrophic bacteria, i a menaquinone (MK) indicative of gram-positive bacteria, j an intact Crenarchaeol diagnostic for Nitrososphaerota, k a carotenoid pigment (Hex-kfuco) indicative of haptophytes, and l a ubiquinone (UQ) as generic biomarker for aerobic respiration. The black circle indicates the eddy core and dashed line indicates potential edge of the eddy-influenced area (eddy periphery). For interpolation between data points, weighted-average gridding was used, which is integrated in Ocean Data View103. MGDG monogalactosyldiacylglycerol, PC phosphatidylcholine, PE phosphatidylethanolamine, TAG triacylglycerol, 1G monohexose, Hex-kfuco 19’Hexanoyloxy-4-ketofucoxanthin. Nomenclature of lipids: For intact polar diacylglycerol lipids and triacylglycerols the sum composition of acyl double bonds and acyl carbon atoms in each compound is provided as C a:b, where C indicates the lipid group, a the number of acyl carbons and b the number of acyl double bonds. For respiratory quinones, the nomenclature after Elling et al.104 is used (Qm:n) where Q indicates head group type, m the number of isoprenoid units in the side chain and n the number of double bonds.
Variations in acyl side chain structure of the microbial lipids
We were interested in the detailed structure of the side chains of the acyl containing lipids, i.e., chloroplast lipids, cell membrane lipids and triacylglycerols, across the identified clusters. On the one hand this provides additional structural information and on the other hand unsaturated fatty acids, which make up the acyl chains, have nutritional properties important for higher trophic levels54. For example, phytoplankton are the main producer of essential fatty acids, such as eicosapentaenoic acid, in the marine environment. Essential fatty acids are important nutrients that higher trophic levels, such as zooplankton, cannot synthesize and must obtain from their diet65. Since our mass spectrometry analysis only allows glycerolipid identification to the lipid species level for some lipid groups (i.e., sum of carbon atoms and sum of double bonds on all fatty acyl moieties due to the lack of diagnostic MS/MS fragment ions, see Supplementary Figs. 5 and 6), the average number of total acyl carbon and total acyl double bonds (unsaturations) were calculated for each lipid group. We observed significant differences for both the number of acyl carbons and the degree of unsaturation for the chloroplast lipids in the periphery of the eddy compared to all other stations with the higher values occurring in the periphery (Fig. 5; Supplementary Fig. 10). This indicates a higher content of long-chain polyunsaturated fatty acids (PUFA) in this region of the eddy. Indeed, the mass spectral analysis of the dominant long-chain, polyunsaturated compounds (MGDG 36:9 and MGDG 36:10) showed that the essential fatty acid eicosapentaenoic acid (20:5n3) is the dominant fatty acid in C36 and C38 diacylglycerols (Supplementary Fig. 11). The second fatty acid in MGDG 36:9 was a 16:4 fatty acid, which is used as a marker for diatoms66, further supporting that larger eukaryotic phytoplankton were the dominant source of these lipids in the eddy periphery. The difference in fatty acid composition in the periphery is likely a result of the differences in phytoplankton community composition. While the chloroplast lipids were of predominantly cyanobacterial origin in the core of the eddy as well as the open ocean, larger eukaryotic phytoplankton were the dominant source of these lipids in the eddy periphery, which contain high abundances of PUFAs32,54. For cell membrane lipids and TAGs, differences only occurred between the eddy and the open ocean, but none was observed between the eddy-influenced stations and the coastal one (Fig. 5; Supplementary Fig. 12). This suggests invariant PUFA content in these lipid groups between the eukaryotic phytoplankton in the coast and the cyclonic eddy, but much lower PUFA content in open-ocean eukaryotic phytoplankton. Changes in temperature are well known to be one of the main drivers for changes in membrane unsaturation in a wide range of marine phytoplankton cultures67,68,69. Additionally, on the global scale, fatty acid unsaturation has recently been shown to be driven by temperature irrespective of the plankton community28. However, temperature variations were small in our local data set (Supplementary Fig. 2). Thus, our results indicate that shifts in phytoplankton group dominance due to mesoscale eddy activity can be the dominant driver for fatty acid unsaturation on local scales.

Total average number of acyl carbon atoms and acyl unsaturations in chloroplast lipids (a, d), cell membrane lipids (b, e) and triacylglycerols (c, f) across the five clusters shown in Fig. 2. Boxes represent the interquartile ranges, the line within the box the median, and the whiskers extend to the minimum and maximum values within 1.5 times the interquartile range. Outliers are shown as circles. The chloroplast lipids include the glycolipids with diacylglycerol side chain (the galactolipids mono- and digalactosyldiacylglycerol and the sulfolipid sulfoquinovosyldiacylglycerol). The cell membrane lipids include all detected phospho- and betaine lipids. Asterisks indicate variables significantly different (one-way ANOVA, p < 0.05) from coastal (dark blue) and open-ocean (light blue) stations, respectively.
Drivers of surface ocean lipidome changes in and around a cyclonic eddy
Our results revealed that the cyclonic eddy formed in the EBUS off Mauritania carried a distinct lipidome signature indicative of enhanced abundance of eukaryotic phytoplankton, bacterial and archaeal biomass compared to the surrounding coastal and open-ocean stations. Since the eddy formed near the coast when upwelling still occurred, but it was sampled during the following relaxation season (see Material and Methods; Supplementary Fig. 3), the lipid composition could reflect the composition found near the coast at the time of eddy formation. Eddy trapping of coastal waters together with the supply of nutrients into the euphotic zone by shoaling the pycnocline likely sustains productivity of the coastal community during westward propagation. There was strong spatial variability in lipid distribution within the eddy in line with our previous observations of a spatial mosaic of microbial processes within this eddy7. This variability can be associated with eddy dynamics. Satellite measurements of upper-ocean chlorophyll revealed that horizontal advection of chlorophyll by the rotational velocity within the interior of eddies results in dipoles with extrema outside of the eddy cores, rather than monopoles of positive or negative chlorophyll anomalies trapped at the eddy cores70. In westward-propagating cyclones that rotate anticlockwise in the Northern Hemisphere, anomalously high and low chlorophyll concentrations typically occur in the southwest and northeast quadrants of the eddy, respectively, but upwelling of nutrient-rich water can still result in high chlorophyll concentrations in the core through enhanced phytoplankton production10. Nutrients are also transported laterally within eddies via advection and diffusion71,72,73 affecting phytoplankton distributions on the submesoscale74. Thus, both vertical and lateral processes on the submesoscale affect the distribution of plankton in eddies leading to patchy phytoplankton distribution and activity7,75,76 and are likely the driver for the observed variations in lipid composition within the eddy. Interestingly, the core of the eddy appears to consolidate important members of the microbial community involved in different processes of the carbon and nutrient cycles through the increased presence of not only photoautotrophic cyanobacteria and heterotrophic bacteria, but also chemoautotrophic ammonia oxidizing Nitrososphaerota in the surface mixed layer as indicated by specific lipid biomarkers for these microbial groups. Lipids provide information mainly about biomass, but not necessarily activity. However, a recent study of a cyclonic eddy in the subtropical Pacific showed enhanced ammonia and nitrite oxidation rates in the euphotic zone indicating that mesoscale mechanisms can affect the distribution of marine nitrifiers53. Future studies combining lipidomics with 16S and 18S diversity data will shed further light on the relationship between changes in microbial community composition and lipid diversity in mesoscale eddies. Additionally, the study of other omics, such as metatranscriptomic surveys of lipid biosynthesis genes, will be useful for linking microbial community dynamics within eddies to lipid distributions.
Effects of lipid distributional changes for the nutritional quality of phytoplankton
The profound impact of the eddy on the microbial lipidome and with this its biochemical composition has further implications for the nutritional quality of the microbes and in particular the phytoplankton. The strong enrichment in TAG storage lipids together with the higher concentration of PUFA in the periphery of eddy suggest that eukaryotic phytoplankton in this region, which are mainly diatoms, are a food source of superior quality for zooplankton compared to non-cyclonic eddy influenced areas. Lipids provide the densest form of energy in the marine environment and contain more than twice the energy per gram compared to carbohydrates or proteins32,77. Additionally, PUFAs are essential fatty acids that are provided by the phytoplankton to higher trophic levels that cannot synthesize these molecules and must obtain them through their diet32,54. The enrichment of these nutritious lipids is thus a plausible driver of the observed changes in zooplankton feeding behavior and population dynamics in the periphery of cyclonic eddies78,79,80,81,82 that ultimately contribute to fisheries resources. It has recently also been shown that human fishing activities are linked to mesoscale eddy activity with highest fishing activities occurring at the margins of cyclones83, which we associate to the high nutritional quality of the phytoplankton. Variability in phytoplankton community composition has further been shown to result in strongly enhanced export of carbon at the periphery of eddies due to high diatom abundances82. Mesoscale eddies can thus efficiently transfer CO2 from the atmosphere to deep waters via the biological carbon pump, but with horizontal variability depending on the position within the eddy. Although the sampling of eddies is difficult due to their transient nature, our lipidomics results support the view that considering spatial heterogeneity within these features will be important for more accurate predictions of the effects of eddies on oceanic carbon cycling75.
Lateral export of lipid carbon from the Mauritanian upwelling area to the open ocean
During westward propagation, mesoscale eddies formed in EBUS can travel hundreds of kilometers into the subtropical gyres, where they may be involved in supporting the relatively high productivities – elevated biological productivity has been shown to be sustained across the extensive ocean subtropical gyres, but the underlying mechanisms have been under debate for decades20,21,22,23. Quantitatively important amounts of POC have been shown to be exported from the California coastal upwelling region to the gyres by mesoscale eddies with an eddy-induced enrichment of more than 20 Gg year−1 17 potentially providing labile carbon for the open-ocean microbes and other organisms. Similarly, model simulations suggest lateral POC export from the Mauritanian upwelling system to the open ocean to be in the Gg to Tg C year−1 range84. Based on our quantitative lipid analysis from the cyclonic eddy, we estimated the amount of lipids carried in the surface mixed layer through the sustained productivity from the Mauritanian upwelling region to the subtropical gyre within cyclonic eddies using the modeling results from Schütte et al.16 for average cyclonic eddy dimensions (radius of 51 km) and the yearly number of cyclonic eddies that live long enough to propagate into the open ocean (3.9 eddies year−1). With our observed mixed-layer-integrated lipid concentration averaged for all eddy stations of 203.6 ± 39.3 mg m−2 (Supplementary Fig. 13), this adds up to 12.9 ± 2.5 Gg lipid year−1, which translates to 9.7 ± 2.0 Gg lipid C year−1. Considering that all major EBUS are hotspots for eddy generation85,86, cyclonic eddies may play a quantitatively important role in the lateral export of labile carbon in the form of energy and nutrient-rich lipids to the open ocean. Indeed, using the number of long-lived cyclonic eddies generated in the four major EBUS based on 10-year records of satellite-altimetry data and Argo float profiles of temperature and salinity from Pegliasco et al.85 (between 40 and 50 cyclonic eddies in each EBUS), 1.8 Tg lipid C year−1 globally may be transported within the mixed layer from coastal upwelling regions to the open ocean. The fate of the energy-rich and essential lipid-producing phytoplankton upon decay of the eddy remains unknown, but will be an important target for future studies as the different agents of mortality differ in their impact on the aquatic food web and carbon cycling. While phytoplankton lysis (autolysis and viral lysis) results in the production of DOM available for heterotrophic microbes, potentially contributing to the often observed net heterotrophy in the open ocean21, zooplankton grazing mortality transfers the carbon to higher trophic levels and enhances export of sinking organic matter from the upper ocean to depth87. Additionally, biological activity and particulate organic matter characteristics have been shown to be affected by the intensity and the age of eddies88,89, which will have consequences for the quality of organic matter delivered to the open ocean. For example, a recent study in the California Current System showed that both the physicochemical and diatom community structure within cyclonic eddies evolve during propagation offshore with a depletion of nutrients leading to lower biological activities over time8. Modeling results further suggest that the carbon content of coastal cyclonic eddies decreases while propagating westward84. Our estimate of surface ocean lipid transport by cyclonic eddies thus likely represents the upper limit as the eddy was still in its intensification stage35. Understanding changes in lipid abundance and composition in the same eddy over time including its vertical structure as well as variability between eddies formed during different seasons will help to constrain the effect of the supplied carbon for the indigenous open-ocean community and the carbon cycle. To this end, future studies that follow eddies during their westward propagation from formation at the coast until decay in the open ocean and concomitant comprehensive biogeochemical analysis will be essential, especially since eddy activity is projected to change significantly in the future90.
Material and methods
Sample collection
Seawater samples for lipidomics were collected during RV Meteor cruise M15691 in the Eastern Tropical North Atlantic off Mauritania (see Fig. 1) using standard Niskin-type bottles attached to a CTD rosette. After collection, the samples (1 L) were filtered using vacuum filtration (ca. −200 mm Hg) onto 47 mm diameter 0.2 µm hydrophilic Durapore filters (Millipore) and immediately flash-frozen and stored at −80 °C until processing. Samples were collected from eleven stations, of which six were located within the zonal eddy passage at 18°N between the Mauritanian Shelf and the Cape Verde Islands. One of these stations (E3) was just beyond the southwestern periphery of a cyclonic eddy, which was identified during the cruise based on satellite imagery (sea surface height) and ship-based Acoustic Doppler Current Profiler analysis. Four additional stations were located inside of the eddy-influenced waters with two of them (EDM4 and EDM5E) at the edge of the eddy core, one inside of the core (EDM4E) and one in the southwestern periphery (EDZ1). As a reference site, the Cape Blanc (CB) time series station (21°10.037’N, 20°54.983’W) was sampled. The exact location of the eddy center was validated and refined by Acoustic Doppler Current Profiler (ADCP) reconstructions and assuming an axis-symmetric vortex analysis (see Supplementary Methods). For details about the eddy identification the reader is also referred Devresse et al.7. We have further used an angular momentum eddy detection and tracking algorithm (AMEDA)92 to estimate the age and origin of the eddy. This analysis revealed that the eddy formed near the coast, was 1.5 months old and still in the intensification stage during sampling7 (see Supplementary Fig. 3a). While the upwelling season is typically from December to April93, satellite-derived chlorophyll concentrations together with the reconstructed ADT using AMEDA from the time of eddy formation in May indicates that coastal upwelling has still occurred at that time and that the sampling in July took place during the relaxation phase (see Supplementary Fig. 3b).
Lipid extraction and analysis
For lipidomics analysis, samples from the surface mixed layer only were selected (Supplementary Table 1) in this study. Lipids were ultrasonically extracted following a modified Bligh & Dyer protocol94 using a monophasic mixture of methanol, dichloromethane, and 50 mM phosphate buffer at pH 7.4 (2:1:0.8, v:v:v). 1,2 dihenarachidoyl-sn-glycero-3-phosphocholine (Avanti Polar Lipids, Alabaster, AL, USA) was used as internal standard. The total lipid extracts (TLE) were dried under a stream of N2 and stored at −20 °C until measurement.
Lipids were analyzed by injecting TLE aliquots dissolved in methanol:dichloromethane (9:1, v:v) on a Dionex Ultimate 3000 high performance liquid chromatography (HPLC) system connected to a Bruker maXis Ultra-High Resolution quadrupole time-of-flight tandem mass spectrometer (MS) equipped with an ESI ion source operating in positive mode (Bruker Daltonik, Bremen, Germany). The MS was set to a resolving power of 27,000 at m/z 1222 and every analysis was mass-calibrated by loop injections of a calibration standard and correction by lock mass, leading to a mass accuracy of better than 1–3 ppm95.
Analyte separation was achieved using reversed phase HPLC on an Acquity UPLC BEH C18 column (1.7 µm, 2.1 × 150 mm, Waters, Eschborn, Germany) maintained at 65 °C as described by Wörmer et al.96. Analytes were eluted at a flow rate of 0.4 mL min−1 using linear gradients of methanol:water (85:15, v:v, eluent A) to methanol:isopropanol (50:50, v:v, eluent B) both with 0.04% formic acid and 0.1% NH3. The initial condition was 100% A held for 2 min, followed by a gradient to 15% B in 0.1 min and a gradient to 85% B in 18 min. The column was then washed with 100% B for 8 min.
The samples were additionally analyzed on the same system under different chromatographic conditions using normal phase HPLC on an Acquity UPLC BEH Amide column (1.7 µm, 2.1 × 150 mm; Waters Corporation, Eschborn, Germany) maintained at 40 °C as described by Wörmer et al.96. Aliquots of the TLE were dissolved in dichloromethane:methanol (9:1, v:v) and analytes were eluted at a flow rate of 0.5 ml min−1 with 99 % eluent A (acetonitrile:dichloromethane, 75:25, with 0.01 % formic acid and NH3, v:v:v:v) and 1% eluent B (methanol:water, 50:50, with 0.4 % formic acid and NH3) for 2.5 min, increasing B to 5% at 4 min, to 20% B at 22.5 min and 40% B at 26.5 min. The column was then flushed with 40% B for 1 min.
Lipids were identified by retention time as well as accurate molecular mass and isotope pattern match of proposed sum formulas in full scan mode and MS2 fragment spectra39,94,97,98 (see Supplementary Figs. 5 and 6 for representative MS2 spectra of each compound group). The normal phase HPLC-MS analysis was used for bacterial and eukaryotic intact polar lipid detection, while reversed phase HPLC-MS analysis was used for archaeal lipids, pigments, respiratory quinones and triacylglycerols (TAGs). Integration of peaks was performed on extracted ion chromatograms of ±10 mDa width and included the [M + H]+, [M + NH4]+ and [M+Na]+ ions using Data Analysis™ and Quant Analysis™ software (Bruker Daltonik, Bremen, Germany). Where applicable, doubly charged ions were included in the integration. Lipid quantification was achieved by comparison of parent ion responses relative to known amounts of the internal standard. Lipid concentrations were corrected for response factors of commercially available standards in external standard curves. For phosphoethanolamine (PE), phosphocholine (PC), and phosphatidylglycerol (PG), 1,2-dipalmitoyl-sn-glycero-3-PE, -PC, and -PG, standards were used, respectively (Avanti Polar Lipids, Alabaster, AL, USA). Monogalactosyldiacylglycerol (MGDG), digalactosyldiacylglycerol (DGDG), sulfoquinovosyldiacylglycerol (SQDG) standards from Avanti Polar Lipids were used to correct MGDG, DGDG and SQDG, respectively. For the betaine lipids Diacylglyceryl-carboxyhydroxymethylcholine, Diacylglycerylhydroxy-methyltrimethylalanine and Diacylglyceryltrimethylhomoserine, a Diacylglyceryltrimethylhomoserine standard also from Avanti Polar Lipids was used. Quinone concentrations were corrected for the relative response of menaquinone (MK; MK4:4 for MKs and vitamin K1) and ubiquinone (UQ; UQ10:10 for UQs, and plastoquinones) standards (Sigma-Aldrich, St. Louis, MO, USA). For chlorophylls and their degradation products, a chlorophyll a standard and for carotenoids an astaxanthin standard also from Sigma-Aldrich was used. TAGs were quantified using a series of individual standards from Nu-Check-Prep, (Elysian, MN, USA). TAG response factors were based on the equivalent carbon number of each TAG after Holčapek et al.99. Purified archaeal lipid standards were obtained from extracts of Archaeoglobus fulgidus as described by Elling et al.100. The concentrations of monohexose glycerol dibiphytanyl glycerol tetraethers (GDGTs) were corrected for the response of a purified acyclic monohexose GDGT standard. The concentration of GDGT core lipids and glycerol dibiphytanol diethers were corrected for the response factors of purified acyclic GDGT as described previously100, while the concentration of archaeol was corrected for the response factor of the respective purified standard. The lower limit of detection for lipids was <1 pg on column. Isotopically labeled standards are not needed for samples from the surface ocean39. For normal phase HPLC-MS analysis, peak areas were additionally corrected for isotopic overlying effects from co-eluting molecular species containing an additional double bond by subtracting the intensity contributed from the natural abundance of isotopologues (e.g., the 2x13C containing isotopologues). For example, PC 30:1 containing 2x13C atoms has a molecular mass of 706.5288 and PC 30:0 with only 12C atoms as a mass of 706.5381, which cannot be distinguished with the used methods. The natural abundance of 2x13C PC 30:1 is 10.5% compared to an only 12C-containing PC 30:1. Thus, 10.5% of the peak area of PC 30:1 was subtracted from the peak area of PC 30:0. Higher isotopologues were neglected due to their low natural abundance (<1%).
Weighted mean calculations
A weighted arithmetic mean of lipid species concentrations was used to describe the average number of acyl carbon atoms and the average number of unsaturations in the acyl-containing lipid groups. They were calculated as follows:
where C is the total number of acyl chain carbons and %lipid the percent abundance of all lipids with C amount of carbon.
where U is the total number of acyl chain unsaturations and %lipid the percent abundance of all lipids with U amount of unsaturations.
Particulate organic carbon analysis
For analyses of particulate organic carbon, water samples (0.5–1 L) were filtered onto pre–combusted (450 °C for 5 h) Whatman GF/F filters (25 mm, 0.7 µm) under low pressure (<200 mbar). Filters were stored frozen at −20 °C until analysis. Prior to analysis, filters were acid–fumed (37% HCl for 24 h) in order to remove inorganic carbon and dried at 40 °C for 24 h. Subsequently, filters were wrapped in tin cups (8 × 8 × 15 mm) and analyzed according to Sharp101 using a EuroEA Elemental Analyzer (HEKAtech, Wegberg, Germany). Blank filters were used to correct for contaminations.
Hierarchical clustering analysis
The hierarchical clustering analysis was carried out using R version 4.3.0 on lipid concentrations (ng L−1) converted into log-space to account for the large value range. First, the Euclidean distance matrix was calculated using the dist function. Then, the hierarchical clustering analysis was performed using the hclust function with default parameters and the average linking method. The functions dist and hclust are built-in functions in R. Clustered heatmaps were generated using the R package pheatmap102.
Responses