Pronounced Uptake and Metabolism of Organic Substrates by Diatoms Revealed by Pulse-Labeling Metabolomics

Diatoms contribute as a dominant group of microalgae to approximately 20% of the global carbon fixation. In the plankton, these photosynthetic algae are exposed to a plethora of metabolites, especially when competing algae are lysed. It is well established that diatoms can take up specific metabolites, such as vitamins, amino acids as nitrogen source, or dimethylsulfoniopropoionate to compensate for changes in water salinity. It is, however, unclear to which extent diatoms take up other organic resources and if these are incorporated into the cell´s metabolism. Here, we explore the general scope of uptake of metabolites from competitors. Using labeled metabolites released during lysis of algae grown under a 13CO2 atmosphere, we show that the cosmopolitan diatom Chaetoceros didymus takes up even dilute organic substrates from these lysates with little bias for molecular weight or polarity. This is reflected by a high degree of labeling in the metabolome of the exposed cells. The newly developed pulse label/mass spectrometry metabolomics approach reveals that polarity and molecular weight has no detectable influence on uptake. We further show that the taken-up metabolites are partly maintained without metabolic modification within the cells, but also a substantial part is subject to catabolic and anabolic transformation. One of the most dominant phytoplankton groups thus has the potential to compete with heterotrophs, suggesting that the observed osmotrophy may substantially impact organic material fluxes in the oceans. Our findings call for the refinement of our understanding of competition in the plankton.


INTRODUCTION
The traditional view on marine plankton distinguishes between phytoplankton as primary producers and zooplankton as consumers (Heesen, 1928). However, many planktonic eukaryotic organisms have been recognized as mixotrophs, which combine autotrophic photosynthesis with organic matter uptake (Selosse et al., 2017). Many microzooplankton grazers are mixotrophic and retain functional algal organelles or even algal endosymbionts. Also photosynthetically active organisms, such as phytoflagellates and dinoflagellates can engulf and consume prey organisms to acquire nutrients (Stoecker et al., 2017). As an additional strategy, the uptake of dissolved organic carbon termed absorbotrophic mixotrophy or osmotrophy can be observed in microalgae. This process seems to be ubiquitous but clearly less understood (Flynn et al., 2013a).
First investigations of osmotrophy in plankton focused on algae growth under extreme darkness in the presence of organic substrates (Neilson and Lewin, 1974). Administration and uptake studies of radiolabeled substrates deepened our mechanistic understanding. However, experiments were always limited to investigating one compound or a compound class, such as specific vitamins or amino acids. It is now clear that mixotrophy is widely distributed among planktonic eukaryotic organisms (Caron, 2016). However, the experiments with single compounds under limiting conditions conducted so far represent an oversimplification and do not reflect the situation in nature, where a cell is exposed to structurally more diverse metabolites. Consequently, the importance of osmotrophy for pelagic food webs and element cycling remains elusive and we are still far from quantitatively deciphering the trophic modes of primary producers forming the phytoplankton (Karlusich et al., 2020).
The availability of organic substrates for uptake will be highly variable. In the plankton mass occurrences of algae, so-called algal blooms can last over days to weeks before the population breaks down and is succeeded by other species that become dominant. Especially during the decay of such algal blooms, the surviving competitors will be exposed to the metabolites of the lysed algae. Also, lysis of specific phytoplankton members by pathogens, such as algicidal bacteria or viruses, results in situations where surviving resistant cells are exposed to the metabolomes released by the lysed species (Meyer et al., 2017). It is entirely unclear if and how all these compounds contribute to the metabolism of the phytoplankton and the potential ecological importance of phytoplankton as consumers of organic material is thus still poorly understood (Burkholder et al., 2008;Flynn et al., 2013b;Borowitzka, 2016;Ferroni et al., 2018). Modeling studies imply a leak-recovery system to be required to initiate a bloom, but the underlying mechanisms are not fully explored (Flynn et al., 2008) As critical primary producers, diatoms were initially classified as autotrophs, but the uptake capability of specific organic molecules has long been recognized (Raven et al., 2009;Flynn et al., 2013b). These include glucose, small polar organic acids such as acetate, succinate, fumarate, malate and lactate, amino acids, dipeptides and dimethylsulfoniopropionate (DMSP) (Lewin and Lewin, 1960;Hellebust and Guillard, 1967;Hellebust and Lewin, 1972;Vila-Costa et al., 2006;Cochlan et al., 2008;Mulholland and Lee, 2009;Spielmeyer et al., 2011;Petrou and Nielsen, 2018). Most experiments provided mechanistic insight but did not accurately reflect natural conditions, where the water in which the algae live harbors a diverse mixture of organic compounds. Under natural conditions, cells are exposed to these metabolites and also light is available to support photosynthesis. Thus, two competing mechanisms for carbon acquisition, phototrophy and heterotrophy will be active. Here we address osmotrophy in the diatom Chaetoceros didymus under non-limited conditions. The supply of nutrients and light in our study was non-limiting to allow efficient algal growth; organic metabolites were thus offered in addition to available inorganic sources.
We base our experimental setup on a well-investigated multipartner interaction involving an algicidal bacterium and two cosmopolitan diatom species. Kordia algicida is a marine Flavobacterium that possesses algicidal activity leading to cell lysis of several microalgal species, including the diatom Skeletonema costatum (Paul and Pohnert, 2011;Paul and Pohnert, 2013). C. didymus, in contrast, is a naturally cooccurring diatom that is resistant against K. algicida. The impact of K. algicida was recently shown in field experiments where it induces a population shift in a natural phytoplankton community towards resistant algae . We hypothesize that during bacterial lysis of S. costatum, resistant species can benefit by taking up metabolites of the lysate in an osmotrophic manner. Therefore, we exposed a culture of C. didymus to metabolites from a 50% diluted stationary culture of lysed S. costatum cells. We developed a novel analytical approach to test this hypothesis, including pulse labeling metabolomics and a novel data treatment routine to analyze isotopologue signatures. We show that highly universal uptake and metabolism of metabolites from the environment occur in the resistant alga.

Experimental Design
Diatoms were grown under non-limited conditions and exposed to a labeled metabolome released by lysed cells of a competitor. The diatoms were extracted using optimized protocols for metabolomics sampling and extracts were analysed by liquid chromatography mass spectrometry to evaluate the uptake of labeled metabolites.

Algal Culturing
C. didymus was isolated by W. Kooistra from the Gulf of Naples (Stazione Zoologica Anton Dohrn, Naples, Italy) and S. costatum was obtained from the Roscoff Culture Collection (Roscoff, France). Both algae were cultivated in batch culture using nutrient-replete artificial seawater medium (Maier and Calenberg, 1994) in 50 mL Greiner Bio-One cell culture flasks at 11-13°C under a 14:10 h light: dark regime with an illumination of 20-25 μmol photons m -2 s -1 . Development of cultures was followed by in-vivo Chl a fluorescence using a Mithras LB 940 plate reader (excitation 430 nm, emission 665 nm).

Global 13 C-Labeling
For global 13 C-labeling of S. costatum we used autoclaved artificial seawater medium that was prepared without addition of NaHCO 3 . An aliquot of this medium was utilized to dissolve NaH 13 CO 3 (98 atom %, Sigma-Aldrich, Munich, Germany). This solution, containing sufficient NaH 13 CO 3 to reach a final concentration of 2.38 mM, was sterile filtered (0.2 μm pore size, Sarstedt Filtropur S) and transferred back to the medium bottle. Tissue culture flasks were filled to the neck in order to minimize the area for CO 2 exchange with the atmosphere and were inoculated with < 1% (v/v) of a stationary S. costatum culture. After growing to stationary phase, an aliquot was taken and transferred to fresh 13 C-enriched medium (< 1% (v/v)). After two of these cycles a plateau in the degree of labeling (verified by mass spectrometry as described below) was reached and the cultures were used for further experiments.

Bacterial Culturing
K. algicida (Sohn et al., 2004) was cultivated on marine broth agar at 30°C for 2 days. The bacterial lawn was removed with a sterile cotton swab and re-suspended in algal culturing medium to an OD 550 of 0.5 determined on a Genesys 10S UV-Vis spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).

Co-Culturing Experiment and Extraction of Released Metabolites
S. costatum cultures reared in medium prepared with NaHCO 3 (further referred to 12 C or 13 C medium were co-cultured in triplicates (185 mL each) with the K. algicida added to a final OD 550 of 0.01. After 6 days the lysed cultures were gently filtered (0.2 μm pore size) and the filtrate extracted for exometabolome analysis as follows. Solid phase extraction cartridges (Chromabond easy, Macherey-Nagel, Düren, Germany) were equilibrated with 4 mL of methanol (Chromasolv Plus, Sigma-Aldrich, Munich, Germany) and 4 mL of water (Chromasolv Plus, Sigma-Aldrich, Munich, Germany) before the filtrate (170 mL) was applied using vacuum with a flow rate < 1 L h -1 . The cartridge was washed with 4 mL water, air-dried and then extracted via gravity flow using 2 mL of methanol followed by 2 mL of methanol/tetrahydrofuran 1:1 (tetrahydrofuran HiPerSolv, VWR, Dresden, Germany). This extract was frozen until further chemical analysis.

Tests for Osmotrophy
Stationary cultures (45 mL) of S. costatum reared in 12 C or 13 C medium were centrifuged (500 x g, 15 min, 10°C) and washed three times by repeated addition of 45 mL of 12 C medium to the harvested pellets and centrifugation. After the third washing step the supernatant was collected as control medium (this was processed in parallel to the cells and later served as control for the natural 13 CO 2 ) and the cell pellet was re-suspended in 45 mL 12 C medium. To disrupt cells, the suspension was frozen at -20°C, thawed and treated in an ice-cold ultrasonic bath for 10 min. The lysate was filtered (1.2 μm pore size, GF/C, Whatman, GE Healthcare, Little Chalfont, United Kingdom), acidified to pH ≤ 1 for sterilization using 30% hydrochloric acid, incubated at 0°C for 10 min and subsequently neutralized under sterile conditions using a saturated sodium hydroxide solution. The solutions were stored at -20°C until use. After thawing, the solutions were diluted 1:1 with 12 C medium to yield the final medium for the investigation of osmotrophy ( 12 C organic & 12 C inorganic or 13 C organic & 12 C inorganic). 45 mL aliquots of both media were extracted as described above for exometabolomic analysis to determine the organic metabolites. To determine osmotrophy, 110 mL aliquots of these media were inoculated with 1% (v/v) stationary C. didymus culture in triplicates and cultivated for 3 days. Directly after inoculation and after 3 days of cultivation, samples (45 mL each) for intraand extracellular metabolomics were processed. Samples were filtered (1.2 μm, GF/C, Whatman, GE Healthcare, Little Chalfont, United Kingdom) and the flow-through was processed for exometabolomic analysis as described above (see co-culturing experiment). The cells were washed off the filter with an ice-cold freshly prepared mixture of methanol/ethanol/ chloroform (1:3:1) (ethanol LiChroSolv, Merck, Darmstadt, Germany; chloroform HiPerSolv, VWR, Dresden, Germany). Extracts were treated in an ultrasonic bath for 10 min, centrifuged (30,000 x g, 15 min, 4°C) and the supernatant was stored at -20°C. As a control to prove that no NaHCO 3 enrichment was present in mixotrophic media, the third wash supernatant (see above) was used as medium to cultivate C. didymus and metabolites were extracted as described. A graphical representation of the experimental setup can be found in Supplementary Figure 1.

Analysis of Exo-and Endometabolomes With LC-MS
Extracts from cells and media (see above) were dried in a nitrogen flow at room temperature and were resuspended in up to 200 μL methanol. Metabolites were separated on an UltiMate 3000 UHPLC (Thermo Fisher Scientific, Waltham, MA, USA) equipped with an Accucore C18 column (100x2.1 mm, 2.6 μm) at 25°C using water with 2% acetonitrile and 0.1% formic acid (A) and pure acetonitrile (B) as mobile phase. The gradient was as follows: 100% A for 0.2 min, linear gradient to 100% B in 7.8 min, 100% B for 3 min, linear gradient to 100% A in 0.1 min, 100% A for 0.9 min. The UHPLC was connected to a QEplus Orbitrap mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) equipped with heated electrospray ionization source (capillary temperature 360°C, sheath gas 60 nominal units, aux gas 20 nominal units, sweep gas 5 nominal units, aux gas temperature 400°C, spray voltage 3.3 kV, S-lens RF level 50) operated in positive or negative ion mode. Full scan measurements (m/z 100-1200, resolution 280k, AGC target 3 10 6 , maxIT 900 ms) were performed separately for positive and negative ion mode. All data are available at https://www.ebi.ac.uk/metabolights/MTBLS2333.

Isotopologue Detection
Full scan RAW files were converted to mzXML using ProteoWizard msConvert (Kessner et al., 2008) with the vendor's algorithm for peak picking. Isotopologue detection was achieved with R-based X 13 CMS (Huang et al., 2014). The R script can be found in the Supplementary Material. In brief, after peak-picking with centwave (3 ppm, peakwidth 5-20 s) and retention time alignment with orbiwarp, isotopologues with a mass difference of 1.00335 Da were assigned (RTwindow 10 s, 3 ppm). Either all isotopologue groups (a = 1) or only isotopologue groups significantly different from 12 C (a = 0.05) were reported. Afterwards, those significantly different isotopologue groups were manually curated with reference to the original spectrum in order to exclude groups that did not contain at least 3 consecutive isotopologues.

Calculation of the Degree of Labeling
To calculate the degree of labeling MS1 data from isotopologues are assumed to have a Bernoulli distribution (see Supplementary Material for formula) (Zachmann, 1994). Measured isotopologue intensities are compared to computed distributions and the squared coefficient of variation between both is minimized in an iterative process.

Generation and Evaluation of a Labeled Metabolome for Uptake Experiments
A complex medium containing the entire labeled metabolome of lysed diatom cells could be generated from 13 C-labeled S. costatum. We grew S. costatum in a medium containing NaHCO 3 as a sole carbon source to obtain a labeled metabolome. With a repeated exchange of the medium, we reached up to 65% labeling in the algal metabolites, sufficient for our analytical protocol. Mechanical lysis of the labeled culture and removing the cell debris gave an axenic medium, rich in organic metabolites (Supplementary Figure 1).
We hypothesized that the resistant C. didymus encounters these metabolites also in the field when algicidal bacteria lyse its competitors. To test this hypothesis, labeled S. costatum cultures were infected with K. algicida, which resulted in lysis of more than half of the diatom cells within six days. Ultrahigh-pressure liquid chromatography-high resolution mass spectrometry (UHPLC-HRMS) of the exometabolome revealed similar metabolic profiles in mechanically lysed S. costatum cells and those lysed by the bacteria which is also reflected in Table 1 where two thirds of the identified metabolites were found in both sample types.
Before incubation with C. didymus, the metabolites from lysed S. costatum were sterilized and diluted with a medium containing inorganic carbon with normal isotope distribution (1.1% 13 C, Supplementary Figure 1). Thereby, we could ensure that 13 C labeled organic metabolites taken up from the medium can be distinguished from de novo synthesized compounds. We exposed C. didymus to a three day pulse of labeled metabolites released from lysed cells of S. costatum and then extracted the culture. Pulse labeling for three days (i.e. ca. two cell divisions) instead of prolonged incubation prevents scrambling of isotopes detected in the mass spectra due to excessive cellular metabolism during incubation.
Following the same procedure, a control medium was generated using S. costatum cells grown in a medium with natural isotope distribution that could be used to generate mass spectra for structure elucidation. Using fragmentation trees and database comparison of the unlabeled metabolites and subsequent comparison to authentic standards, we identified several of the labeled metabolites in C. didymus (Table 1).
Our study is limited to the fraction of metabolites that can be extracted and separated in UHPLC (and extracted with SPE in the case of exometabolites) using our standardized protocols. Thus, a broad range of metabolites is investigated, but we cannot claim to cover the entire metabolome (which would be technically not feasible in any metabolomics approach to date). Since we did not use a biasing extraction protocol for the generation of labeled medium, but rather used the cell lysate and since extraction success is identical for isotopologes and the unlabeled compounds we can conclude that our results will not be distorted due to extraction selectivity. We did not aim for a quantification of the respective metabolites in this metabolomics approach but rather surveyed the physiochemical properties of metabolites taken up during exposure to organic metabolites released from a lysed culture.

Evaluation of Uptake
For uptake experiments, C. didymus was cultivated in a medium containing the sterilized labeled S. costatum endometabolome or control medium for three days (n=3). We selected the concentration of added metabolites to be equivalent to 50% of those released by a lysed stationary culture. Cells were then recovered by filtration and washed extensively. After being exposed to this medium under otherwise optimum growth conditions, UHPLC-HRMS analysis of the C. didymus endometabolome revealed that the alga took up substantial amounts of labeled compounds from various metabolic classes. Quantitative analysis of labeling proved to be challenging in terms of chemoinformatic data treatment. Therefore, ions of the same metabolite, only differing in their number of incorporated 13 C were summarized in an isotopologue group using the software X 13 CMS (Huang et al., 2014). Of 5587 isotopologue groups (positive and negative ionization mode) detectable in the endometabolome 2381 were significantly labeled with 13 C ( Figure 1A). After manual curation, 548 isotopologue groups were categorized according to their labeling pattern and analyzed regarding their retention time and mass to charge ratio. The degree of labeling was estimated using the probability mass function for Bernoulli trials (random experiments with precisely two possible outcomes, "success i.e. incorporation of 13 C" and "failure i.e. incorporation of 12 C ", in which the probability of success is the same every time the experiment is conducted): The degree of labeling p for a metabolite with n carbon atoms is estimated from the intensities of a set of isotopologues I(m) with m 13 C atoms (Pohnert and Jung, 2003).

Classification of Taken-Up and Processed Metabolites
We categorized the metabolites according to their degree of labeling in the receiving cells ( Figures 1B, C). 73% of the isotopologue groups were lowly labeled with a degree of labeling < 5% indicative of compounds mainly synthesized de novo in C. didymus. The low degree of labeling that still exceeds the natural 13 C content of 1.1% can be explained by general utilization of taken-up metabolites in the metabolism: metabolites that are assimilated are catabolized to metabolic building blocks that are used together with the autotrophic metabolic pool and used for anabolism again. Seven percent of the isotopologue groups were highly labeled. They contained a degree of labeling similar to that of metabolites in the medium (ca. 65%). The cellular content of these highly labeled metabolites taken up from lysed S. costatum shows that C. didymus assimilates metabolites that it does not (or only to a minimal extent) produce itself. Certain compounds with a high degree of labeling can be found in C. didymus but not in the DOM labelenriched medium. These could arise from the metabolic transformation of more complex metabolites released by S. costatum. 15% of the isotopologue groups were labeled in a mixed manner. In the mass spectrum both, signals from a lowly labeled and a highly labeled proportion of the respective compound can be detected. This pattern can be explained by metabolites that are biosynthesized by C. didymus and also acquired from the medium. This compound class includes a wide range of natural products, from small charged molecules like choline and carnitine to non-polar lipids and fatty acids. A few signals (5%) had a complex labeling pattern that Bernoulli statistics could not describe. These can be interpreted as compounds that result from metabolites taken up and further metabolized using the pool of de novo synthesized metabolites ( Figures 1B, C). For example, intermediate metabolic products like fatty acids can be utilized in the anabolism of more complex ones like lipids.

Properties of Taken Up Metabolites
Labeled metabolites span a wide range of polarity from charged small molecules to non-polar fatty acids. Most labeled products are found in the non-polar region of the chromatogram ( Figure 1D, Table 1 and Supplementary Figure 2). The identified polar metabolites that are efficiently taken up include the amino acids leucine and/or isoleucine and several small  Figure 2; Uptake in % (the total amount of the respective metabolite is defined as 100%); presence of the metabolite in the medium enriched with the S. costatum metabolome; presence in medium after lysis of S. costatum by K. algicida; level of identification according to (Sumner et al. 2007), level 1 identified compound, level 3 putatively characterized compound class. Uptake values in () are estimates, since the complex labeling pattern does not allow a proper Bernoulli fit. *released from S. costatum during lysis by K. algicida.

Analysis of Selected Metabolites
Detailed analyses of the isotopic pattern enabled us to determine the ratio of heterotrophic uptake to de novo biosynthesis and look for evidence of mixed strategies. For example, the isotopologues of 6Z,9Z,12Z,15-hexadecatetraenoic acid ( Figure 2A) originate from two distinct pools ( Figure 2D), a lowly labeled pool from de novo biosynthesis and a highly labeled pool from uptake. Modeling with Bernoulli statistics showed that the de novo biosynthesis pool had a degree of labeling of 2.8%, slightly higher than what would be expected from photosynthesis using exclusively natural inorganic carbon ( Figure 2B). The highly labeled pool has a very similar isotopologue distribution compared to the fatty acid from the lysed labeled alga ( Figure 2C). Modeling a mixture of these two pools showed that ca. 20% of the metabolite in the algae results from uptake, while 80% are synthesized de novo. This complementation of de novo synthesized products with externally available metabolites is observed for several metabolites with variable proportions of the two sources from only a few percent to nearly 80% of the metabolite acquired by uptake ( Table 1).
The isotopologue patterns of metabolites that can be explained by an uptake of resources from the medium followed by transformations within the cell using the pool of de novo synthesized precursors are more complex. They do not follow the Bernoulli statistics since different resources can be utilized in d i ff e r e n t r e l a t i v e a m o u n t s . A n e x a m p l e i s t h e lysophosphatidylcholine shown in Figure 3A. The isotope pattern of the lipid ( Figure 3B) cannot be interpreted with the model described above, but tandem MS experiments allow for dissecting the lipid. This reveals a unique labeling pattern for those parts of the molecule that are derived from different biosynthetic pathways. The isotopologue pattern of oleic acid in the lysophosphatidylcholine ( Figure 3C) can also not be fitted with a Bernoulli statistic. The fatty acid is thus assembled using resources that were taken up as well as de novo produced. This labeling pattern in the fatty acid moiety is also observed in free oleic acid and 11Z-eicosenoic acid. Lipid assembly thus does not discriminate between acquired and de novo synthesized resources. In contrast, the choline fragment detected in the same substance shows that most choline is highly labeled and therefore taken up ( Figure 3D). The glycerol moiety of the molecule is not giving charged fragments, but its labeling could be established indirectly. Therefore, we conducted fragmentation of the M+8 isotopologue of the lysophosphatidylcholine. The MS/MS of this ion gave rise to an oleic acid fragment with isotopologues containing down to zero 13 C ( Figure 3E). The remaining eight carbon atoms in the uncharged C8-fragment have thus to be labeled in different degrees. The fragments can thus only be derived from a precursor with labeled, thus acquired glycerol.
These exemplarily discussed mass spectra stand for several hundreds of labeled peaks in the chromatogram of the endometabolome of exposed C. didymus (see Figure 1 and Supporting information). Using the combined results, we can now draw a picture of the osmotrophic metabolism in C. didymus (Figure 4). The uptake introduced here for C. didymus is not limited to this one species but broader distributed in diatoms. When the diatom Thalassiosira weissflogii was raised in the labeled medium as described for C. didymus, labeling patterns similar to those described above were detected (data not shown). shows that the lowly labeled pool contains 2.8% 13 C and thus likely derives from anabolism. (C) The highly labeled pool is taken up from the medium, black bars represent measured data of cellular hexadecatetraenoic acid, white bars measured data of hexadecatetraenoic acid in medium. (D) The measured mass spectrum in D can be explained by 20% hexadecatetraenoic acid derived from uptake and 80% from de novo synthesis (grey bars). Relative intensity refers to the sum of isotopolog ions of the molecular peak. All data are mean ± SD from biological triplicates.  In order to study potential osmotrophy in the diatom C. didymus a medium rich in 13 C labeled organic metabolites was prepared by mechanical lysis of a globally labeled S. costatum culture and removal of the cell debris by filtration. We diluted this filtrate to reach concentrations for all individual compounds in the incubation that correspond to 50% of those found in a stationary, lysed culture. Given the fact that bacterial lysis of phytoplankton in the lab and field happens within hours and that entire blooms of algae can be lysed, we conclude that the concentrations employed are in the range of those found during lysis of a dense bloom . The employed concentrations will nevertheless be higher than under natural non-bloom conditions. However, we observed uptake of metabolites irrespective of their relative concentration in the extract. Not only abundant metabolites were found to be labeled and therefore acquired compounds. Also, minor and trace signals in the chromatograms indicated the incorporation of labels. Given the sensitivity of the mass spectrometer with detection limits in the femptomolar range, this indicates that the uptake mechanism is not driven by the concentration of the metabolites in the water and does thus not result from a "force-feeding" of C. didymus.
The resulting metabolome enriched medium showed a similar (but labeled) metabolic profile as the medium of an infection experiment where the lytic bacterium K. algicida lysed S. costatum. We thus conclude that K. algicida is a sloppy feeder not utilizing the entire algal metabolome but instead leaving substantial organic resources that it does not require or recover in the water. Surviving algae, such as the resistant competitor C. didymus will be exposed to these metabolites. Earlier, we observed that such exposure to metabolites released from lysed algae supports the growth of C. didymus if administered at low concentrations. However, no information about the underlying mechanism was available . The ability to take up released organic molecules may counterbalance the metabolic costs to maintain resistance mechanisms and would be highly advantageous, for example, during the collapse of a competing algal bloom as simulated in this study (Sohn et al., 2004;Meyer et al., 2017).
In addition to the release of organic material upon cell lysis, intact algae can also be a source of dissolved metabolites. Field studies indicate that natural phytoplankton communities can release 20 to 40% of newly assimilated dissolved inorganic nitrogen in form of dissolved organic nitrogen into the environment (Flynn and Berry, 1999). Also, secretion of dissolved organic carbon e.g., in the form of mucus-forming saccharides is well-documented (Alderkamp et al., 2007). If the metabolic uptake also functions in such situations or phases where dilute phytoplankton abundance results in low concentrations of dissolved metabolites would have to be verified in follow-up studies. Future work focusing on the quantification of selected metabolites could support uptake mechanisms and kinetics. We have, however, no indication from the analysis of labeling that there would be a bias for the uptake of higher concentrated metabolites.
We reached up to 65% labeling in the S. costatum metabolome using repeated exchange of medium containing NaHCO 3 as an exclusive carbon source. Higher degree of labeling might be reached in gas-tight culturing devices after more generation of culturing, but was not required for the present study. Analysis of the C. didymus endometabolome after being exposed to the labeled metabolome of S. costatum for 3 days under otherwise optimum growth conditions revealed that the alga took up substantial amounts of labeled compounds from various metabolic classes. Given that analysis of every single isotopologue of a labeled metabolite would generate highly complex and partially redundant information, we reverted to a statistical treatment assuming that labeling results in an isotopologue distribution in MS1 data that can be described with a Bernoulli statistic. In the Bernoulli trial a random experiment with exactly two possible outcomes is evaluated. In this case "success" would represent the incorporation of 13 C and "failure" that of 12 C. We estimate that excess substrate is present, thus the probability of success is the same every time a carbon is incorporated. This data treatment allowed to match labeling patterns to predicted spectra that would result from specific degrees of labeling. Thereby the average 65% labeling of the S. costatum metabolome could be determined.
Labeling Patterns of Metabolites in C. didymus After Exposure to 13 C-Labeled S. costatum Extract Analysis of the endometabolome of the resistant alga C. didymus after exposure to metabolites from 13 C-labeled S. costatum revealed unlabeled metabolites, metabolites with unaltered full labeling and those with more complex mass spectra that could be assigned to different metabolic processing (Figure 1). The methodology described here thus not only shows the uptake of one metabolite, but allows simultaneous quantification of the uptake and analysis of the metabolic fate of all taken up metabolites in the receiving alga. This is a valuable expansion of the classical fluxomics approach, where feeding of one single labeled metabolite to heterotrophs can only reveal its uptake and metabolism. Our experimental approach reflects the situation in the plankton with highly complex microbial communities and complex exometabolomes (Barofsky et al., 2010). It allows evaluation of the uptake capacity and incorporation of the taken-up metabolites under ecologically relevant conditions. General Patterns in the Uptake of Organic Metabolites by the Receiving Diatom C. didymus Detailed analysis of the mass spectra allows to draw a picture of osmotrophy in C. didymus (Figure 4). More than a quarter of all detected metabolites in C. didymus were labeled to different degrees ( Figure 1). Nearly 10% showed the identical labeling pattern as those in the 13 C-labeled S. costatum metabolome. These metabolites are not (or only to a minor extent) synthesized by the receiver but taken up and maintained in the cells or the uptake could down-regulate synthesis in the receiver. The major part of taken-up metabolites showed mixed labeling, indicative for the uptake of a metabolite that is also synthesized by the receiver. This includes a wide range of natural products from small charged molecules like choline and carnitine to non-polar lipids and fatty acids. Metabolites acquired from the outside water can thus be metabolized in the same way as the de novo produced compounds, indicating that no compartmentation of assimilated material occurs. Instead, the osmo(hetero)trophic and the photo(auto)trophic pool are used for anabolism and catabolism. For example, the presence of complex labeling patterns can be explained by the use of intermediate fatty acids in the anabolism of more complex lipids. The complex isotopologue pattern of oleic acid and 11Z-eiosenoic acid, for example, might reflect a dynamic system with rounds of beta oxidation releasing labeled acetate and subsequent fatty acid reassembly from this labeled and the own unlabeled acetate pool. Since the isotopologues of these fatty acids are not Bernoullidistributed the relative amount of labeled acetate in the total acetate pool used for fatty acid biosynthesis might vary over time.

Utilization of Assimilated Metabolites
The fact that uptake has no apparent preference for nitrogen containing metabolites contradicts the assumption that diatoms use specific heterotrophic mechanisms to acquire reduced nitrogen (Mulholland and Lee, 2009). Since we worked in nonlimiting conditions this lack of specificity might also reflect that diatoms do not need additional nitrogen. Flynn and Berry (1999) suggested that these microbes leak metabolites from internal mM levels and that they thus must have some recovery (uptake) potential, else they would 'bleed to death'. Our findings will likely also hold true for such a recovery of metabolites released by the cells themselves. Thus, such a recovery mechanism would be in place and could be quantified using a modified labeling approach.
Notably, many metabolites with essential physiological functions are thus not exclusively produced de novo but taken up in high proportions. Control of physiological concentrations will thus have to include biosynthesis, catabolism as well as uptake. A universal uptake could be highly advantageous in the event of a lysis of a competing algal bloom (Sohn et al., 2004;Meyer et al., 2017). But even under regular conditions in plankton, algae might encounter metabolites released by other members of the phytoplankton that they can take up and benefit from (Sommer et al., 2012).
The observation of universal uptake of metabolites requires re-thinking of the interpretation of incubation studies with single labeled substrates. It is well documented that specific primary metabolites can be taken up by the cells under limiting conditions. Thus, under darkness, uptake of organic acids like lactate and malate was observed (Hellebust and Lewin, 1972). Amino acid usage is prompted by both N and/or C stress (Flynn and Wright, 1986) and uptake of the essential metabolite dimethylsulfoniopropionate (DMSP) was discussed as a way to compensate for the lack of own DMSP biosynthetic capabilities (Spielmeyer et al., 2011). We now expand this view and introduce that the diatoms C. didymus and T. weissflogii can take up a structural variety of metabolites released from competitors even under non-limiting growth conditions. Thus, a universal complementation of the own metabolome occurs and not a selected uptake of potentially valuable metabolites. This unspecific uptake of metabolites over a wide range of polarities and masses newly defines the investigated diatoms as a general sink of organic carbon in the sea and might be generalized to the entire group of these algae.

Physiological Function of Assimilated Metabolites
Of the multiple labeled metabolites, two groups will be exemplarily highlighted here, polar nitrogen containing metabolites and fatty acids. Polar metabolites that are taken up include glycine betaine, b-alanine betaine, carnitine and choline, with diverse physiological functions. Glycine betaine is a known osmolyte and uptake in diatoms has already been described under N-limited conditions (Keller et al., 1999;Gebser and Pohnert, 2013). Choline is a biosynthetic precursor of glycine betaine (Kageyama et al., 2018) and an important building block of phospholipids (Borowitzka, 2016). In this study, it has been found as a free form and in lysophosphatidylcholine ( Table 1). However, the labeling pattern of the choline fragment of lysophosphatidylcholine differs from the free form ( Supplementary Figure 3), indicating different origins. balanine betaine is a known osmoprotectant in plants and is biosynthesized via the methylation of b-alanine, a building block of coenzyme A (Rathinasabapathi et al., 2001). It has also been found in marine algae (Blunden et al., 1986). Carnitine is a central metabolite in energy metabolism of all eukaryotic cells. It plays an essential role for the transport of fatty acids across the mitochondrial membrane. The utilization of carnitine by diatoms has been reported previously by measuring oxidation rates in a biofilm-forming freshwater diatom (Tuchman et al., 2 0 0 6 ). The presence of labeled N -m e t hy l g ro u p s (Supplementary Figure 4) excludes lysine degradation (Paik et al., 1977) as a metabolic origin. The high proportion of heterotrophically acquired carnitine and the presence of many fatty acids among the taken-up metabolites is striking. It suggests a highly active transport mechanism for high-energy fatty acids and the N-containing carnitine. It would be interesting to investigate whether an acylcarnitine-type transport system might facilitate the uptake of fatty acids across the cellular membrane. The presence of a protein homologous to a class I carnitine/acylcarnitine translocase in a diatom cell wall proteome supports this hypothesis (Frigeri et al., 2006).
Several labeled fatty acids and derivatives have also been identified ( Table 1). Eicosapentaenoic acid is one of the dominant fatty acids in diatoms and precursor of many bioactive oxylipins (Pohnert and Boland, 2002;Borowitzka, 2016;Rettner et al., 2018;Zulu et al., 2018). 6Z,9Z,12Zhexadecatrienoic acid is the precursor of octadienal, an allelopathic polyunsaturated aldehyde (D'ippolito et al., 2003). Also, oleic acid, 11Z-eicosenoic acid, and 6Z,9Z,12Z,15hexadecatetraenoic acid are common in diatoms (Ackman et al., 1968;Dunstan et al., 1994;Pohnert et al., 2004;Chen, 2012;Zulu et al., 2018). Two other C18-and C20 polyunsaturated fatty acids were identified according to retention time and mass spectra. Comparison with synthetic standards showed that these fatty acids were isomers of linoleic and arachidonic acid. Fatty acids and their derivatives have multiple physiological functions. The fact that they are present in labeled and unlabeled forms after the pulse labeling experiments indicates an unbiased utilization of de novo and assimilated compounds. This might serve as a mechanism to avoid costly biosynthesis of fatty acids, which is supported by the observation that the fatty acids are not only internalized or adsorbed to lipidic structures due to their physicochemical properties but rather incorporated into the primary metabolism of the diatom.
The broad range of polarity and molecular weight of assimilated compounds raises the question of the uptake mechanism. Candidate systems for fatty acid shuttling include specific transporters and genome data shows homolog candidate sequences (Frigeri et al., 2006;Indiveri et al., 2011). But indeed, also unspecific osmotrophy and incorporation mechanisms that do not require transporters will most likely be involved. Nevertheless, the question of how diatoms acquire exogenous metabolites is still open.

Concluding Remarks
The diatom C. didymus takes up and incorporates metabolites from lysed competitors with surprising little bias. Preliminary results suggest that the same is true for the diatom T. weissflogii (data not shown), which indicates a rather universal mechanism. The multitude of cellular functions in which these metabolites are involved suggests that uptake complements the internal metabolic pool from many different biosynthetic pathways.
The universal osmotrophic lifestyle of members of one of the most abundant algal classes in the oceans substantially changes our view of the metabolic shuttling in phytoplankton communities. These algae take up more metabolites than just a few highly polar compounds for e.g. nitrogen supply. They rather complement their metabolism quite universally with resources from the surrounding seawater. This process occurs even under illumination and is not related to the complementation of lacking photosynthate in the dark. It will have to be verified if this is a universal property of diatoms and if also other members of the phytoplankton utilize similar strategies, as it is already indicated by work on Emiliania huxleyi (Godrijan et al., 2021). In addition to bacteria, diatoms and other phytoplankton compete for the dissolved organic carbon in the plankton (Berge et al., 2012;Georges et al., 2014). Our study thus has consequences for element cycling in the oceans and ecosystem dynamics that will have to be addressed in the future.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: MetaboLights [accession: MTBLS2333].

AUTHOR CONTRIBUTIONS
NM developed the pipeline and analyzed the data. AR acquired and analyzed data. NM and GP conceived the study, directed all experiments, and wrote the manuscript with contributions from the co-author. All authors approved the manuscript.