News

May 01, 2023

Fungal parasitism on diatoms alters formation and bio

Communications Biology volume

Communications Biology volume 6, Article number: 206 (2023) Cite this article

1381 Accesses

1 Citations

3 Altmetric

Metrics details

Phytoplankton forms the base of aquatic food webs and element cycling in diverse aquatic systems. The fate of phytoplankton-derived organic matter, however, often remains unresolved as it is controlled by complex, interlinked remineralization and sedimentation processes. We here investigate a rarely considered control mechanism on sinking organic matter fluxes: fungal parasites infecting phytoplankton. We demonstrate that bacterial colonization is promoted 3.5-fold on fungal-infected phytoplankton cells in comparison to non-infected cells in a cultured model pathosystem (diatom Synedra, fungal microparasite Zygophlyctis, and co-growing bacteria), and even ≥17-fold in field-sampled populations (Planktothrix, Synedra, and Fragilaria). Additional data obtained using the Synedra–Zygophlyctis model system reveals that fungal infections reduce the formation of aggregates. Moreover, carbon respiration is 2-fold higher and settling velocities are 11–48% lower for similar-sized fungal-infected vs. non-infected aggregates. Our data imply that parasites can effectively control the fate of phytoplankton-derived organic matter on a single-cell to single-aggregate scale, potentially enhancing remineralization and reducing sedimentation in freshwater and coastal systems.

Phytoplankton grows fast, renewing its biomass once per week and leaving behind several gigatons of decaying organic matter in the water column of the aquatic biosphere1,2. Most of this decaying organic matter is remineralized within the euphotic zone, but a critical fraction gets exported to deeper layers as sinking particulate material in diverse systems, including freshwater3,4 and coastal environments5,6. Particulate organic matter is thereby transported through the water column towards sediments—a mechanism that drives the bentho–pelagic coupling, i.e., the exchange of energy, mass, and nutrients between pelagic and benthic habitats. Sinking organic matter thus sustains life below the euphotic zone7, while, on the other hand, it can also lead to expanding oxygen-deficient waters worldwide in areas where oxygen consumption through organic matter remineralization exceeds the available oxygen supply8,9,10. The fate of phytoplankton-derived organic matter is, therefore, key to aquatic food webs and biogeochemical cycles.

After blooming, phytoplankton cells together with fecal pellets and other detritus can coagulate as aggregates—referred to as lake or marine snow—which sink rapidly with velocities of up to several hundred meters per day11. In this regard, sinking aggregates are primary vehicles of (in-)organic matter to depth12. The formation, remineralization, and sedimentation of aggregates are commonly explained by physical and biological processes, the latter being intimately linked to zooplankton grazing, bacterial degradation, and viral lysis12,13,14,15,16. Recent data, however, suggests that we are still missing some crucial links in the network of biological control mechanisms on vertical organic matter fluxes17. Aquatic fungi, for instance, are hardly considered in this context although they are abundant and active across various aquatic systems18,19, reaching from high-mountain lakes20 over coastal areas21 to the deep sea22. For instance, members of the fungal division Chytridiomycota, referred to as chytrids, can thrive as microparasites on phytoplankton cells. These chytrids are traditionally well documented in lakes23, especially during bloom events when host abundances are high24,25. More recently, they have also gained increasing attention in coastal systems after being observed via microscopy, for instance, in the highly productive upwelling region off Chile26, in the Arctic Ocean27,28, and during harmful algae blooms in the Mediterranean Sea29. Using high-throughput sequencing methods, Chytridiomycota have also been detected in other marine regions30,31,32,33,34, and DNA-based abundances of Chytridiomycota have been correlated to phytoplankton blooms or chlorophyll in coastal regions33,35,36,37,38. Chytrids thus occur frequently in freshwater and coastal regions where a strong benthic–pelagic coupling is commonly assumed7,39, whereas in open ocean regions, their distribution is less evident with rare observations so far40.

Parasitic chytrids infect up to 1–44%24 or even 80–100% of their specific phytoplankton host population41,42,43, including diatoms, dinoflagellates, and cyanobacteria. They thereby alter the dynamics of phytoplankton populations44 and are presumed to impact organic matter remineralization and sedimentation45. In a natural lake, chytrid-infected phytoplankton has been shown to mainly decompose in surface waters, whereas non-infected cells settled to depth41. It has further been estimated that 20–25% of the photosynthetically derived carbon from the infected host population is channeled to parasitic chytrids (through the fungal shunt)46 and further to zooplankton (through the mycoloop)47. Chytrids have therefore been proposed to be an integral part of aquatic food webs and element cycling23,48. Yet, their impact on the fate of phytoplankton biomass in terms of aggregate formation and vertical organic matter fluxes remains unknown.

Diatoms, which form fast-sinking aggregates49,50 and often mediate high proportions of the particle export51,52, are highly susceptible to fungal infections in freshwater and coastal systems26,28,44. We thus used one of the few available diatom–parasite model pathosystems (freshwater diatom Synedra and chytrid Zygophlyctis)53 to investigate the impact of fungal infections on the formation and characteristics of sinking diatom aggregates and associated organic matter fluxes. We hereby considered cell aggregation, mass density, settling velocity, sticky polymers, bacterial colonization, and respiration—parameters that are all intimately linked to organic matter remineralization and sedimentation processes. Moreover, we complemented these culture-based investigations with analyses of bacterial abundances and infection prevalence in aggregates that were formed from field-sampled phytoplankton communities. Our results demonstrate that epidemics of fungal microparasites need to be appreciated as an episodically impactful control mechanism on organic matter fluxes in the aquatic biosphere.

Data obtained from the cultured diatom Synedra sp. and chytrid Zygophlyctis planktonica pathosystem are indicated as model system and data from the field-sampled populations as natural systems.

We incubated non-infected and chytrid-infected diatom co-cultures—hereafter non-infected and fungal-infected treatment, respectively. Both culture treatments included the freshwater pennate diatom Synedra and co-growing bacteria. The designated infected treatment additionally included the parasitic chytrid Zygophlyctis, which developed Synedra-associated sporangia and free-living zoospores, both being part of the life cycle of this fungal microparasite53. One set of both treatments (each in triplicates) was grown and sub-sampled over a growth period of nine days, to monitor the culture development (referred to as culturing). An additional set was grown for six days and thereafter transferred into rotating cylinders, to analyze the formation and characteristics of aggregates (referred to as rotating cylinders). Symbols, terms, and abbreviations used in the following are summarized in Supplementary Box S1.

For cell enumeration, we distinguished non-infected, early-infected, maturely-infected, post-infected, and decaying Synedra cells (Fig. 1). Non-infected Synedra were recognized as healthy, intact cells with chlorophyll autofluorescence. Early-infected Synedra showed an encysted zoospore, which grew into a mature sporangium (incl. visible zoospores) on maturely-infected cells. Post-infected Synedra cells displayed remains of chitinous cell walls as a sign of previous infections after zoospore discharge. Decaying cells displayed no chlorophyll autofluorescence and also no signs of previous infections, and thus, they had most likely not undergone any fungal infection46.

a Non-infected (non-inf) healthy diatom cells. b Infected diatom with encysted zoospore (encyst, early-infected), which developed into a c mature sporangium (mat sp, maturely-infected) with new zoospores inside. The zoospores are eventually discharged into the ambient water, leaving behind an empty transparent sporangium (empty sp, post-infected) on a non-viable host cell (d). Shown are brightfield and epifluorescence (Fluorescence) images. Sporangia are displayed in blue (stained with Calcofluor White) and diatom chloroplasts in red (chlorophyll autofluorescence). Scale bar is 20 µm (applies to all panels).

During the 9-day growth period, total Synedra abundances increased from approximately 0.8 × 104 to 3 × 104 cells mL−1 (Fig. 2a, b), with similar growth rates for non-infected and infected Synedra populations during exponential growth (day 0–4, 0.26 ± 0.02 and 0.26 ± 0.02 d−1, respectively, P = 0.64, t-test, N = 3 incubation flasks, see also Supplementary Fig. S1). In the infected treatment, the infection prevalence (including early-infected, maturely-infected, and post-infected Synedra cells) reached 47% on day 6 (mostly maturely-infected cells) and remained at the same percentage until day 9 (mostly post-infected cells, Fig. 2b). Time series data of chlorophyll a were significantly different between infected and non-infected cultures (P < 0.0001, generalized linear mixed model GLMM, F4,16 = 217.63), with similar chlorophyll a contents in both treatments during day 0–2 (P > 0.05, t-test, N = 3 flasks) but significantly lower chlorophyll a contents in the infected vs. non-infected cultures during day 4–9 (e.g., day 9: 9.5 ± 1.2 and 30.1 ± 1.7 ng chl a mL−1, respectively, P < 0.001, t-test, N = 3 flasks, Fig. 2c). Nutrient concentrations (analyzed at the end of the incubations on day 9) were similar in both culture treatments (NO3− + NO2−: 225 ± 2 and 243 ± 9 µmol L−1, P = 0.07; soluble reactive phosphorous: 31 ± 2 and 34 ± 3 µmol L−1, P = 0.32; and dissolved organic carbon: 332 ± 24 and 299 ± 11 µmol L−1 in the non-infected and infected treatment, respectively, P = 0.12, t-test, N = 3 flasks).

Diatom abundances and infection prevalence in the non-infected (a) and fungal-infected culture (b), c chlorophyll a, as well as d TEP and e CSP concentrations, monitored over nine days. Stacked bars a, b display mean values (bar height) with triplicate data points (circles) and standard deviations (error bars). Data in c–e are shown as single data points, while lines connect mean values for consecutive time points (symbols for mean values are not displayed). P-values, shown in the upper right corner of each graph, indicate statistically significant differences between the entire time series of each parameter in both treatments (GLMM). Statistically significant differences between treatments at single time points are indicated with asterisks (*P < 0.05, ***P < 0.001). a–e N = 3 incubation flasks. c–e Blue and red items represent data from the non-infected and infected cultures, respectively. TEP Transparent Exopolymer Particles, CSP Coomassie Blue Stainable Particles, XG equ. Xanthan Gum equivalents, BSA equ. Bovine Serum Albumin equivalents.

Transparent Exopolymer Particles (TEP) and Coomassie Blue Stainable Particles (CSP) were determined microscopically and spectrophotometrically54,55,56. Photometrically-analyzed TEP concentrations increased until day 6 and decreased thereafter in both treatments (range: 0.07–0.75 µg XG equ. mL−1, Fig. 2d). The time-dependent development of TEP concentrations, however, was significantly different between treatments (P = 0.005, GLMM, F4,16 = 5.618). TEP concentrations were significantly higher (P < 0.001) on day 2, but significantly lower on day 9 (P = 0.04, pair-wise comparison) in the infected treatment in comparison to the non-infected treatment. On the other sampling days, no significant differences were detected between both treatments (P > 0.22). Similar to TEP concentrations, CSP concentrations increased until day 6 and decreased thereafter (range: 0.07–0.39 µg BSA equ. mL−1, Fig. 2e), but no statistically significant differences were detected between treatments (P = 0.34 for GLMM, F4,16 = 1.2145 and P > 0.09 for pair-wise comparison between single sampling days). Microscopy-derived numbers of TEP and CSP were highly variable, ranging from 119 to 3672 TEP mL−1 (dTEP = 2–4 µm, A = 3–12 µm2) and from 966 to 21,196 CSP mL−1 (dCSP = 3–6 µm, A = 9–32 µm2) in both culture treatments (ranges represent mean values during day 0–9, data are listed in Supplementary Data 6). No statistically significant differences in TEP and CSP sizes and areas between treatments were detected (P > 0.13).

Using fluorescence microscopy, we enumerated (i) free-living bacteria as well as diatom-associated bacteria based on their association with (ii) non-infected, (iii) early-infected, (iv) maturely-infected, (v) post-infected, and (vi) decaying Synedra cells. On day 9, the abundance of diatom-associated bacteria was lowest in conjunction with non-infected Synedra cells (6.5 ± 4.2 bacteria diatom−1) and increased with increasing infection stage, from early-infected to maturely-infected, post-infected, and decaying Synedra cells (9.9 ± 10.5, 9.6 ± 7.2, 16.1 ± 12.0, and 24.2 ± 13.7 bacteria diatom−1, respectively, Fig. 3a). The same cell stage of Synedra was colonized by fewer bacteria in the non-infected culture than in the infected culture (3.1 ± 2.2 and 6.5 ± 4.2 bacteria per non-infected Synedra cell, respectively, and 17.8 ± 15.3 and 24.2 ± 13.7 bacteria per decaying Synedra cell, respectively). Consequently, diatom-associated bacteria were in total 3-times more abundant in the infected culture (Fig. 3b). Abundances of free-living bacteria increased with time, reaching significantly higher abundances on day 4 (P = 0.0004) and day 6 (P = 0.0001) in the infected culture. On day 9, however, the abundances of free-living bacteria were similar in both treatments (P = 0.59, t-test, N = 3 incubation flasks, Fig. 3c).

Abundances of diatom-associated bacteria are shown per individual Synedra cell (a) and per culture volume (b) on day 9. Bacteria were grouped based on their association with individual Synedra cells of different infection stages. The letters a–d in (a) denote significantly different groups (Kruskal–Wallis, P < 0.05, N = numbers of analyzed Synedra cells). Data are shown as single data points (white-filled circles), the 25th, 50th, and 75th percentile (grey boxes), mean values (blue and red-filled circles), and distribution curves (kernel smoothing). Data from the previous sampling days 0–6 are listed in Supplementary Data 2. Stacked bars in b represent individual sample replicates. c Abundances of free-living bacteria. Single data points are plotted as circles and lines connect mean values over time (symbols of mean values are not displayed). Time series data of free-living bacterial abundances were significantly different between both treatments (P = 0.012, GLMM, F4,16 = 4.591). Blue and red items represent data from the non-infected and fungal-infected cultures, respectively. Statistically significant differences at single time points are indicated with asterisks (***P < 0.001). b, c N = 3 incubation flasks.

Aggregate formation was reduced substantially in the presence of the fungal parasite, that is, aggregates were less abundant and smaller in fungal-infected vs. non-infected Synedra treatments (Fig. 4a–d). In detail, during the initial 12 h of aggregate formation, small aggregates (d = 0.4–1.3 mm) were 38 ± 14-times less abundant (range: 20–58-times, # L−1) and comprised 46 ± 24-times less cumulative aggregate volume (range: 13–78-times, mm3 L−1, N = 5 time points) in the fungal-infected Synedra treatment. Thereafter, also large aggregates (>1.3–5.6 mm) became abundant, and both size classes (d = 0.4–1.3 mm and >1.3–5.6 mm) were 6 ± 4-times less abundant (range: 2–14-times) and comprised 9 ± 7-times less cumulative aggregate volume (range: 3–25-times, N = 10 time points) during fungal infections. Aggregates within the smallest size bin (d = 0.35–0.46 mm) were most abundant (up to 80 and 26 aggregates L−1 in the non-infected and fungal-infected treatment, respectively, Fig. 4e, f), while large aggregates (d = 2.2–3.8 mm) accounted for most of the cumulative aggregate volume (up to 20 and 2 mm3 L−1 in the non-infected and fungal-infected treatments, respectively, Fig. 4g, h).

a, b Size spectra of aggregates during the initial 30 h of cylinder rotation. Size spectra show the number of aggregates per water volume normalized to the width of each size bin. Log-transformation was used to visualize the wide size spectrum range. Values of −1.0 indicate that no aggregates in the respective size bin were present (i.e., log10(0) was set to −1.0). Small black dots mark the sampled time points and size bins. c, d Exemplary particle images recorded after 30 h. e–h Aggregate concentration and volume spectra for different size bins are shown for aggregates imaged after 6 h and 30 h. The diameter represents the equivalent circular diameter (ECD). Data in e–h are shown as single data points (triplicates) and lines connect mean values for consecutive time points (symbols for mean values are not displayed, but error bars represent standard deviations). Empty and filled circles show data after 6 and 30 h. The imaging setup is shown in Supplementary Fig. S2. Supplementary Fig. S3 shows exemplary aggregates imaged after sampling and the aggregation status after 45 h in the infected treatment.

We sampled individual aggregates from the rotating cylinders to analyze their bio–physical properties. Those aggregates were similar in size, with mean diameters of 2.4 ± 1.3 mm for non-infected aggregates and 2.4 ± 1.0 mm for infected aggregates (P = 0.93, t-test, Table 1). Bacteria were 5-times more abundant on fungal-infected aggregates than on non-infected aggregates (12.3 ± 0.8 × 104 and 2.7 ± 0.5 × 104 bacteria agg−1, P < 0.0001, t-test, Fig. 5a) due to two aspects. Firstly, abundances of diatom-associated bacteria were lowest on non-infected Synedra, and increased consecutively on early-infected, maturely-infected, post-infected, and decaying Synedra, (0.8–18.8 bacteria diatom−1, Fig. 5b). And secondly, fungal-infected cells were 1.7-times enriched within aggregates as compared to the ambient water (71 ± 3% and 42 ± 3% infection prevalence, respectively, P < 0.0001, t-test, Fig. 5c). This enrichment was mostly due to post-infected cells, which were 2.6-times more abundant within the aggregates as compared to the ambient water (51 ± 6% and 19 ± 3% relative abundance, P < 0.0001, t-test, Fig. 5d).

Bacterial abundances on whole aggregates (a) and single Synedra cells (b), infection prevalence (c), and relative abundances of different Synedra cell types within aggregates and in the ambient water (d). Bacteria in a were grouped based on their association with individual Synedra cell types within single aggregates. The letters a–e in b denote significantly different groups (Kruskal–Wallis, P < 0.05). N indicates the number of analyzed cells. Data in b are shown as single data points (white circles), the 25th, 50th, and 75th percentile (grey boxes), mean values (blue and red-filled circles), and distribution curves (kernel smoothing). Stacked bars in a, d represent individual replicates. Mass density (e), settling velocities (f), and respiration rates (g, h) of single aggregates. f Settling velocities were positively correlated with the aggregate diameter. Fitted curves (power function) were significantly different between both treatments (P = 0.04). Blue and red items represent data from the non-infected and fungal-infected aggregates, respectively. g Respiration rates given as nmol O2 per µmol agg-POC × h define the oxygen consumption per mol aggregate-POC and hour. Rates of 0.1 d−1 indicate that 10% of the aggregate-POC were respired per day. h Respiration rates per meter settled, estimated for differently-sized aggregates. Rates of 0.1% m−1 indicate that 0.1% of the aggregate-POC were respired per meter settled. Data in c, e, g are shown as single data points (white circles), the 25th, 50th, and 75th percentile (white boxes), and mean values (black circles).

The mass density ρs-agg of aggregates, i.e., the density of the particulate material in hydrated aggregates (also termed solid hydrated density)57, was significantly higher for non-infected aggregates than for infected aggregates (1.284 ± 0.003 and 1.258 ± 0.004 g cm−3, respectively, P < 0.0001, t-test, Fig. 5e, Table 1). Similarly, mass densities of single cells ρs-cells were significantly higher for non-infected Synedra cells than for infected Synedra cells (1.269 ± 0.001 and 1.251 ± 0.004 g cm−3, respectively, P = 0.002, t-test, Table 1). Excess densities, defined as the mass densities of aggregates in excess relative to the ambient water, were also significantly higher for non-infected aggregates in comparison to fungal-infected aggregates (0.8 ± 0.3 and 0.6 ± 0.1 mg cm−3, Welch-test, P = 0.009, see Supplementary Fig. S4a for the size~excess density curve fit). Settling velocities ranged from 52–534 m d−1, and increased with increasing aggregate size. The size~settling velocity curve fit (power function as used by58) was significantly different between non-infected and infected aggregates (P = 0.04, F-test comparing the curve fit, F1,22 = 3.651, Fig. 5f). The fractal dimension D3, which describes the fractal geometry of aggregates, was significantly lower for non-infected aggregates than for fungal-infected aggregates (2.37 ± 0.30 and 2.64 ± 0.24, P = 0.005, Table 1), indicating a more compact, less porous structure of infected aggregates (see Supplementary Fig. S4b for the curve fit used to derive D3). This indication, however, was not confirmed by similar porosities of non-infected and infected aggregates (Table 1) and faster settling velocities for non-infected vs. infected aggregates (Fig. 5f). Carbon contents were similar for both aggregate types of similar sizes (6.7 ± 2.8 and 8.0 ± 5.3 µg POC agg−1, P = 0.52, t-test), but carbon-specific respiration rates were on average 2-fold higher for fungal-infected aggregates as compared to non-infected aggregates (0.16 ± 0.08 and 0.34 ± 0.08 d−1, P = 0.0002, t-test, Fig. 5g). These respiration rates are high compared to previously measured rates of laboratory-formed and field-sampled aggregates (approx. range: 0.06–0.13 d−1)6,57,59,60,61,62,63. Due to the slower settling velocities and higher respiration rates, infected aggregates were estimated to respire 2.4‒3.7-times more carbon per meter settled than non-infected aggregates (diameter range: 1‒5 mm, Fig. 5h).

We sampled a natural phytoplankton community in a temperate lake (Lake Stechlin), which comprised mostly filamentous cyanobacteria (Planktothrix, 398 ± 74 filaments mL−1 and Aphanizomenon/Pseudanabaena, combined 443 ± 11 filaments mL−1) and diatoms (Synedra 280 ± 63 cell mL−1 and Stephanodiscus/Fragilaria, combined 1382 ± 345 cells mL−1). Chytrid infections were highest on Planktothrix, Synedra, and Fragilaria cells, and thus, we inspected these taxa in more detail (~2% infection prevalence during our sampling, while prevalences of up to 44% on Synedra and Fragilaria have been observed previously at the same location)24. Bacterial abundances were ≥17-times higher on infected cells than on non-infected cells (Planktothrix 0.1 ± 1.2 and 13.4 ± 14.8 bacteria 100 µm-filament−1, Synedra: 0.9 ± 2.2 and 17.5 ± 20.9 bacteria diatom−1, and Fragilaria: 0.5 ± 1.5 and 8.8 ± 6.3 bacteria diatom−1, respectively, Fig. 6a). For Planktothrix and Fragilaria, fungal-infected cells were relatively more abundant in aggregates than in the ambient water (Planktothrix: 3.4 ± 2.7% and 0.2 ± 0.4% infection prevalence, P < 0.001, t-test, N = 30 and 29, respectively; and Fragilaria: 2.1 ± 1.6% and 0.5 ± 1.3%, P < 0.001, t-test, N = 30 and 29, respectively, Fig. 6b). In contrast, relative abundances of infected Synedra cells were similar in both aggregates and the ambient water (2.6 ± 3.1% and 2.7 ± 3.5% prevalence, P = 0.47, t-test, N = 30 and 29, respectively).

a Bacterial abundances on non-infected and fungal-infected phytoplankton cells. The numbers of associated bacteria are given per 100 µm filament length (Planktothrix) or per diatom cell (Synedra and Fragilaria). b Infection prevalence in the ambient water and aggregates, after cell aggregation in rotating cylinders. N indicates the number of analyzed cells a or the number of aggregates and ambient water samples b. Data are shown as single data points (white circles), the 25th, 50th, and 75th percentile (grey boxes), mean values (blue and red circles), and distribution curves (kernel smoothing).

Understanding the sinking flux of particulate organic matter marks a long-standing frontier in aquatic science, while over the last two decades, the consideration of the multiple mechanisms that control those fluxes has expanded substantially14,64,65,66,67. We suggest that epidemics of eukaryotic microparasites, such as the herein-studied parasitic fungi, should be regarded as an additional biological mechanism that can attenuate vertical organic matter fluxes, next to the commonly considered bacterial degradation, viral lysis, and zooplankton grazing64,68,69. We build upon this suggestion in the following discussion. Our data, however, derive mostly from one model system and from the freshwater environment. We, therefore, encourage future studies on eukaryotic parasites in diverse aquatic systems to uncover their full potential to modulate sinking organic matter fluxes. Moreover, parasitic infections induce cascading effects on phytoplankton community composition70, zooplankton production71,72, and carbon recycling46, all of which are known to alter vertical mass fluxes6,50,73,74,75,76,77. The weighted impact of parasite epidemics on sinking organic matter fluxes thus likely varies in multifarious plankton networks—variations that we are only starting to disentangle.

The photosynthetic biomass of the diatom host population was substantially diminished during fungal infections, as implied by up to 3-fold lower chlorophyll a at 47% infection prevalence in comparison to the non-infected culture (Fig. 2c). Whether this decrease in chlorophyll results in an actual decrease in photosynthetic activity has not been quantified to date, but in seaweeds, chytrid infections were shown to reduce the photosynthetic efficiency78. The reduction in photosynthetic pigments may thus also reduce the photosynthetic biomass that could get exported to depth. Bacterial abundances, especially those of phytoplankton-associated bacteria, were promoted up to 3-fold on the cultured Synedra cells during fungal infections (Fig. 3a, b) in line with previous laboratory-based studies46,71,79,80. We here confirmed this pattern also for field-sampled populations, which demonstrated an even higher—at least one order of magnitude—increase in bacterial colonization on fungal-infected vs. non-infected phytoplankton cells (Fig. 6a). This observation may be explained by phytoplankton–bacteria interactions that are more essential for cell growth under nutrient-depleted conditions in nature, as opposed to nutrient-replete conditions in culture. Bacterial colonization on infected phytoplankton was possibly promoted by bacterial chemotaxis as bacteria are commonly attracted by decomposing, nutrient-leaking phytoplankton cells81,82,83. Additionally, healthy diatom cells can actively repel opportunistic bacteria through chemical signaling84, an ability that fungal-infected host cells might have lost and consequently have become overgrown by bacteria in our incubations.

Parasitic infections reduced the formation of sinking diatom aggregates substantially. For instance, aggregates larger than 0.5 mm, which can dominate the downward flux in natural systems85,86, reached abundances of 5 L−1 after 6 h when parasitic fungi were absent, whereas it took 12 h to reach comparable abundances when the parasites were present (Fig. 4). Aggregate formation in rotating cylinders generally depends on cell concentrations and cell-to-cell stickiness12,13. In our culture incubations, cell concentrations were similar in both treatments (~5000 Synedra L−1, resembling Synedra blooms87), but the stickiness, i.e., the probability that two cells coagulate upon encounter, was presumably reduced during fungal infections. The stickiness of cells and particles can be increased by the presence of exopolymers like TEP and CSP, which are broadly classified as acidic polysaccharides and alkaline amino acids, respectively88,89. TEP concentrations showed significantly different dynamics over time in our non-infected vs. infected Synedra culture (Fig. 2d). However, differences in absolute TEP concentrations on each sampling day were minor, and we presume that fungal infections changed the molecular composition rather than concentration of TEP. Concentrations of CSP were similar in both culture treatments, supporting previous suggestions that CSP are less involved in aggregate formation than TEP90.

Interestingly, infected cells preferentially aggregated over non-infected cells in the fungal-infected population, and yet, overall aggregation was reduced in comparison to non-infected populations. To explain this observation, we consider mechanical and chemical stickiness91. Concerning mechanical stickiness, the development of fungal sporangia on the outside of infected Synedra cells increased the effective size of cells, potentially increasing the interparticle encounter rate and entangling of cells (see Supplementary Fig. S5 for illustration). Concerning chemical stickiness, we may need to distinguish between cell-associated and freely-suspended polymers, as shown recently for diatoms92. We suggest that fungal infections shifted polysaccharides towards (1) fewer freely-suspended polymers, reducing the overall aggregation, but (2) more adhesive polymers at the cell surface of infected Synedra cells, enhancing their cell-to-cell stickiness. Indeed, the observed faster growth of bacterial populations in the infected treatments (Fig. 3) indicates a modification in the organic carbon pool during fungal infections, as shown previously79. Moreover, bacterial activities, which are impacted by fungal infections46, are known to either promote or inhibit aggregate formation93,94. We, therefore, hypothesize that freely-suspended polymers were degraded more effectively by free-living bacteria during fungal infections, whereas cell-associated polymers accumulated at the surface of post-infected diatoms. To test this hypothesis, we recommend investigating the molecular composition and stickiness of both cell-associated and freely-suspended polymers during parasitic infections.

Synedra cells in culture, as well as field-sampled Planktothrix and Fragilaria, aggregated preferentially when being infected as compared to their non-infected counterparts. In contrast, infected Synedra cells sampled in the field did not show such preferential aggregation since non-infected and infected cells were proportionally equally distributed in aggregates and the ambient water. Our Synedra cell isolate originated from a different lake (Lake Haussee) than the field-sampled cells (Lake Stechlin), and their morphologies were different (e.g., field-sampled cells were ca. 4-times longer than the cell isolate). The different aggregation patterns between cultured and field-sampled Synedra may thus have resulted from genera-specific differences (involving diatom-specific bacteria interactions) or even strain-specific differences, as known from viral infections on phytoplankton95. Moreover, in culture, post-infected Synedra cells showed a higher aggregation potential than early-infected and maturely-infected cells (Fig. 5). The infection stage is therefore crucial when examining aggregation dynamics. In the field-sampled populations, we did not differentiate between different infection stages (cells were solely classified as non-infected or infected), but considering the low infection prevalence, we suggest that the field-sampled Synedra population was at an early epidemic stage with mostly early-infected cells.

The development of sporangia on the outside of the host cell expanded the biovolume of diatom cells by 16 ± 5% (N = 20 cells, Table 1) with mostly organic material with, according to literature values, rather low density (cytoplasm: 1.10 g cm−3, chitin: 1.425 g cm−3) as compared to siliceous frustules (1.82 g cm−3)96,97. The lower mass density of infected Synedra cells and aggregates in comparison to non-infected cells and aggregates may further be explained by more cell-associated polymers (such as TEP), fewer cells per aggregate, or thinner silica frustules as bacteria are known to dissolve silica98. Following Stokes’ law (Eq. 5), lower mass densities result in slower settling, which was confirmed by two lines of evidence in our data set. Firstly, based on the size~settling velocity regression (Fig. 5f), we estimated that infected aggregates (71% prevalence, d = 1–5 mm) sank 11–48% slower than non-infected aggregates. And secondly, the herein-measured lower mass densities (measured in a density gradient) and lower excess densities (Eq. 4) of infected vs. non-infected aggregates would both result in ca. 30% slower settling velocities of infected aggregates (following Stokes’ law), similar to the size~settling velocity regression. We thus conclude that host-associated sporangia increased the cell size of diatom cells but reduced their mass per volume (g cm−3), slowing down their settling. Reduced aggregate formation and hence, smaller aggregates during fungal epidemics may have an even stronger impact on sinking velocities. For example, maximum aggregate diameters were 1 mm in the non-infected treatment but only 0.5 mm in the infected treatment after 6 h of aggregate formation. Such a 2-fold decrease in aggregate diameter led to 3-fold slower sinking velocities (according to the size~settling velocity regression in Fig. 5f).

Aggregates that comprised 71% fungal-infected cells respired twice as much of their carbon content per day as compared to non-infected aggregates with similar size and POC content (Fig. 5g). We presume that those higher respiration rates were due to the above-mentioned increase in bacterial abundances during fungal infections and due to the respiration of the parasitizing fungi. Chytrid fungi grow as host-associated sporangia, which after maturation, discharge free-swimming zoospores to seek out a new host. Especially those zoospores are expected to have a high oxygen demand since they are rich in sterols and fatty acids99,100 to support their intensive swimming and seeking behavior. Considering the herein determined carbon respiration rates and settling velocities (Fig. 5h), we roughly estimated that infected aggregates (d = 1 mm) would have lost 37% and non-infected aggregates only 12% (i.e., 3-times less) of their initial carbon content after settling 50 m, indicating a substantial reduction in export efficiencies due to fungal infections.

Taken together, our incubations revealed that fungal infections decreased the chlorophyll content, cell mass density, cell aggregation, aggregate size, and settling velocity, but increased bacterial colonization and carbon respiration on a single-cell to single-aggregate scale (Fig. 7). As a result, fungal infections may be expected to shift organic matter fluxes towards decreased sedimentation and increased remineralization through several mechanisms, with implications for ecosystem-wide bentho–pelagic coupling. Given that fungal parasites are widespread in aquatic systems across different climate zones globally27,28,29,44, including productive upwelling regions26, commercial mass cultures101, and areas impacted by harmful algal blooms29, their epidemics thus have the potential to diminish the strength and efficiency of vertical organic matter fluxes in natural and engineered aquatic environments.

The blue and red text and arrows indicate the fungal-induced increase and decrease of each property, respectively. Diatom cells a and aggregates b are illustrated in light grey, diatom chlorophyll in green, fungal sporangium and zoospores in orange, and bacteria in dark grey.

The model pathosystem comprised the pennate diatom host Synedra sp. Ehrenberg, 1830 (also referred to as Ulnaria102, strain HS-SYN2, l = 87.6 ± 4.1, w = 3.3 ± 0.5, and h = 5.0 ± 0.8 µm, N = 50 cells) and the parasitic chytrid Zygophlyctis planktonicum (strain SVdW-SYN-CHY1), including host-associated sporangia (7.5 ± 0.8 µm, N = 20) and zoospores (d = 3 µm, measured in53) isolated from lakes in Northern Germany53. Bacteria were co-isolated with the diatom and chytrid and maintained in co-culture in a nutrient-replete medium for several months, which likely selected for copiotrophic taxa46. Chytrid infections on the diatom host were initiated by free-swimming zoospores, which attached to a host cell, encysted, and developed into epibiotic sporangia while penetrating and digesting the host's interior through a rhizoidal system53. Within the (zoo-)sporangium, the next generation of zoospores was formed and finally released through the rupture of the sporangium. One infection cycle took 1–2 days. Each infection was lethal and prohibited further reproduction of the host.

Batch cultures were grown in CHU-10 medium (Supplementary Table S1) at constant temperature (17 °C). The light regime was 16:8 h, providing 40 µE s−1m−2 during the 16-h light phase. Synedra was grown as 12 × 650 mL in 1 L Erlenmeyer flasks until reaching ca. 7500 cells mL−1. Half of those flasks (6 × 650 mL) were thereafter inoculated with a Synedra–Zygophlyctis co-culture (40 mL with ca. 98% prevalence), including mature fungal sporangia, which were projected to discharge new zoospores within the next hours (infected treatment). The resulting infection prevalence was 9% in the infected treatment on day 0. The other 6 × 650 mL flasks remained without Zygophlyctis (non-infected treatment). To monitor the culture development over time, 6 × 650 mL (three per treatment) were grown and sub-sampled for nine days (see culture characteristics). The remaining 6 × 650 mL (three per treatment) were grown for six days and subsequently transferred into rotating cylinders to follow the formation of aggregates (see formation and characterization of sinking aggregates). Erlenmeyer flasks were gently shaken by hand once per day for cell resuspension and water mixing.

Sub-samples were taken from the non-infected and fungal-infected diatom cultures (in triplicates) after 0, 2, 4, 6, and 9 days. CHU-10 medium served as a blank for all analyses.

For chlorophyll a and nutrients analyses, 3 mL were transferred into 15 mL Falcon tubes and centrifuged at 3000 rpm for 20 min at 4 °C. The supernatant was gently removed and the remaining pellet was extracted in 90% acetone, frozen at −20 °C, and analyzed after one month using a fluorometer (436/680 nm, Hitachi Fluorescence Spectrometer F-7000). Concentrations of chlorophyll a were calculated and corrected for phaeopigments, measured after acidification103. Calibration was done with a chlorophyll a standard (Anacystis nidulans, Sigma C6144, 0.3–300 ng mL−1, precision ±1%). Concentrations of nitrate/nitrite and soluble reactive phosphorous were determined from 0.45  µm-filtered water by flow injection analysis (FIA) and spectrometric detection (ISO-13395-D28 and ISO/DIS-15681-2, respectively). Concentrations of dissolved organic carbon (DOC) were analyzed from GF/75-filtered water on a TOC-V CPH with a nondispersive infrared sensor (Shimadzu, Kyoto, Japan) using combustion catalytic oxidation.

Diatom abundances and infection prevalence were analyzed from 2 mL sub-samples, which were transferred into microcentrifuge tubes, preserved with Lugol (50 µL per 100 mL sample, see Supplementary Table S1 for Lugol recipe), and stored at 4 °C in the darkness. Synedra cells were counted in Utermöhl plankton chambers (Hydrobios, Germany) under an inverted epifluorescence microscope (Nikon Ti2-E, x400 magnification). The chitinous cell walls of sporangia were stained with Calcofluor White (CFW, 5 µg mL−1, excitation wavelength Ex 387/11, emission wavelength Em 442/46 nm, Merck F3543)104. Synedra cells were differentiated as (i) non-infected, (ii) early-infected, (iii) maturely-infected, (iv) post-infected, and (v) decaying cells (see Fig. 1 for cell type i–iv). The chlorophyll autofluorescence was inspected at Ex 635/18 nm and Em 680/42 nm. We counted at least 50 cells per group (i–iv), or half of the chamber if fewer cells were present.

For bacterial abundance analyses, 0.5–1 mL were preserved with paraformaldehyde (PFA, 1.5% final concentration), filtered onto polycarbonate filters (PC, 0.2 µm, 25 mm, Whatman), and stored at −20 °C. Prior analyses, filters were stained with 4′,6-diamidino-2-phenylindole (DAPI, 1 µg mL−1) and Wheat Germ Agglutinin (WGA conjugated to Alexa Fluor™ 488, 5 µg mL−1, binding to chitinous cell walls, Thermo Fisher Scientific W11261)104. Bacteria were counted under a fluorescence microscope (Leitz Leica DMRB) at x1000 magnification. For each group and replicate, 20 Synedra cells (for attached bacteria) or 20 counting grids (for free-living bacteria) were examined, to reach representative mean values (standard error ≤5%).

Transparent Exopolymer Particles (TEP) and Coomassie Blue Stainable Particles (CSP) were determined microscopically and spectrophotometrically54,55,56. Volumes of 22 mL for spectrophotometric and 2 mL for microscopic analyses were gently filtered (<150 mbar) onto PC filters (0.4 µm, 25 mm, Nucleopore, Whatman). TEP filters were stained with Alcian Blue (AB, 0.02%, pH = 2.5) for 5 s, and CSP filters with Coomassie Brilliant Blue G (CBB-G, 0.04%, pH = 7.4) for 30 s. After staining, filters were rinsed with MilliQ to remove excess dye and stored at −20 °C. Sterile-filtered water (0.2 µm) from the cultures was used as blank. For spectrophotometric analyses, AB was extracted in 80% H2SO4 and CBB-G in 3% SDS in 50% isopropyl alcohol, as specified in Supplementary Methods S1. TEP concentrations are reported relative to xanthan gum (µg XG equivalents L−1) and CSP concentrations relative to bovine serum albumin (µg BSA equivalents L−1).

For microscopic analyses, the stained filters were placed onto CytoClear slides and inspected under a microscope (brightfield, x200 magnification). For each filter, 40 images were taken randomly along two transects across the entire filter area. Care was taken to use the same microscope and camera settings during the entire image acquisition. Image analyses were done in ImageJ 1.51p105. Images were corrected for uneven light intensities in the background and split into separate RGB channels. The blue channel was subtracted from the red channel, to remove chlorophyll-related areas and outlines of diatoms cells. The obtained 8-bit greyscale images were reversed in color and a global threshold was applied. Automated particle analyses in ImageJ provided the number of particles per filter area (used to calculate the number of particles per filtered volume, # mL−1) and the cross-sectional area A of each particle (used to calculate the equivalent circular diameter ECD, assuming spherical geometry, denoted as diameter). Particles with ECD < 3 µm were not included. Since the automatic TEP and CSP recognition was not perfect, the processed images were compared visually with the original images, and falsely outlined TEP or CSP were manually removed from the data set. Cell outlines and the infection stages were poorly visible on our images, and thus, we removed all cell-associated TEP and CSP outlines to not introduce any bias. Only freely-suspended TEP and CSP were therefore included in the microscopy analyses, whereas the spectrophotometric analyses included freely-suspended and cell-associated TEP and CSP.

To analyze the mass density of diatom cells (on day 6), 10 mL were transferred into 15 mL Falcon tubes and centrifuged (3000 rpm, 5 min). The supernatant was gently removed, and cells were re-suspended in 1 mL CHU-10 medium and transferred onto a density gradient, consisting of six viscous layers with increasing density from top to bottom (1.18–1.38 g cm−3). The viscous layers were prepared with mixtures of Ludox TM colloidal silica (50 wt. % suspension in H2O, Merck 420778), sucrose, and distilled water49. After 12 h and an additional centrifugation run (3000 rpm, 5 min), cells had settled into the density layer that was equivalent to the mass density ρs-cells of their particulate material when hydrated in liquid (also termed solid hydrated density)57. The intact density layers (visible to the naked eye) were sub-sampled into 2 mL tubes and stored at 4 °C. Half of each sample was analyzed for its mass density ρs-cells (g cm−3) using a density meter (Anton Paar DMA 38). The other half was used to enumerate non-infected and infected diatom cells (after CFW staining) under an inverted fluorescence microscope with combined bright-field and UV excitation (x400 magnification, Olympus CKX41, Japan).

Triplicates of non-infected and fungal-infected cultures were grown for six days until the infection prevalence had reached 59 ± 3% (N = 3 flasks) in the infected treatment. Thereafter, the cultures were diluted with CHU-10 medium to approx. 5000 cells mL−1 (resembling cell abundances during bloom scenarios87) and transferred into rotating cylinders (r = 8.5 cm, h = 10 cm, V = 2.3 L, bubble-free) to mimic the aggregation of cells due to their stickiness and differential settling behavior (setup is shown in Supplementary Fig. S2). This technique allows forming aggregates that are commonly larger and sink faster than natural aggregates (e.g., compare with106) but their physical properties such as density, porosity, and composition are comparable to natural aggregates107,108. Initial diatom abundances were similar in both treatments (5363 ± 495 cells mL−1 and 4765 ± 367 cells mL−1, respectively, P = 0.17, t-test, N = 3 cylinders). Cylinders were rotated at 1–2 rpm at 17 °C in darkness on a roller table. The aggregate formation was recorded over time using a Mini Deep-focus Particle Imager (MDPI, Bellamare, La Jolla, CA, USA). The MDPI is a ‘shadowgraph imager’ using a near-infrared LED light source (A007 Indus star, LED dynamics, Randolph, VT, USA) behind a pinhole and a set of identical plano-convex collimator lenses to create parallel light beams between the LED light pod and the camera pod (Supplementary Fig. S2). Using this parallel light, i.e., telecentric optics, ensured that all aggregates were in focus, independent of their position between the lenses109. We recorded two image sequences per cylinder every 2–9 hours. Each sequence lasted for at least one rotation (50 images). The sampling volume of two image sequences covered 2 × 0.35 L, equal to 30% of the entire cylinder volume.

Image sequences were processed similarly to TEP and CSP images (see above). Aggregates were counted, measured, and binned into size classes, in which the upper diameter was 1.3-times the lower size class. The resulting eleven size classes ranged from 0.4 to 5.6 mm (aggregates with an ECD < 0.4 mm were excluded from the data set, and aggregates with ECD > 5.6 mm were not present). We counted the number of aggregates per size class, reported as aggregate number concentration N(d) (# L−1). To describe the number concentration as a function of size, the aggregate size spectrum n(d) (# L−1 mm−1) was calculated as

where ∆d (mm) is the width of each size class (d equals ECD)110. The volume spectrum nVd was calculated by normalizing the volume distribution to the aggregate size

where Vagg (mm3) is the average volume of single aggregates in each size class111. Vagg and d (equal to ECD) were calculated from A derived from the image analyses and assuming spherical geometry. After stopping the rotation, individual aggregates were picked with a wide-bore glass pipette for various analyses. To ensure that aggregate sizes were similar for non-infected and infected aggregates, aggregates were sampled after 30 h from the non-infected treatment and after 45 h from the non-infected treatment (see Supplementary Fig. S3 for details). For subsequent analyses, aggregates were grouped into five sets (set I–V, each with 4–12 replicate aggregates) since not all parameters could be measured on the same aggregate.

Set I was used to determining the size, porosity, excess density, settling velocity, and fractal dimension. Aggregates were imaged directly after sampling using a compact camera (Panasonic DMC-TZ22, 14 MP). The cross-sectional area A was determined in ImageJ and used to calculate the ECD and volume V, assuming spherical geometry. The porosity φ was calculated as

where ∆ρ is the excess density of aggregates and ρs-agg is the mass density of the particulate material in the aggregates when hydrated in liquid (solid hydrated density)57. We measured the settling velocity U of aggregates in a settling column (d = 11 cm, h = 40 cm). The column was filled with the same water as the rotating cylinders and kept at the same temperature (17 °C). It was double-walled, with a 2 cm interspace filled with water and a sealed top, except for a 2 cm inlet to introduce the aggregates. This setup minimized convective currents in the column during settling experiments112. Each aggregate was allowed to freely settle from the wide-bore pipette into the water. The settling was timed along 15 cm using a stopwatch. Excess densities ∆ρ of aggregates were calculated using Stokes’ law, which applies to ballasted phytoplankton aggregates113

where ν is the kinematic viscosity (1.085 × 10−2 cm2 s−1 at 17 °C) and g the gravitational acceleration (981 cm s−2). The drag coefficient CD and Reynolds number Re derived from

We determined the fractal structure of aggregates, using their three-dimensional fractal number D3. D3 was derived from the slope of d and U after log–log transformation (Supplementary Fig. S4b)112.

Set II was used for oxygen respiration and particulate organic carbon (POC) analyses. Single aggregates were transferred into 5.9 mL gas-tight Exetainer® vials with 0.2 µm-filtered water from the rotating cylinders. Oxygen concentrations were measured with a Clark-type oxygen microsensor (OX-500, Unisense A/S, Denmark) immediately after aggregates were added and after 24 h-incubations at 17 °C in darkness. During incubations, the Exetainer rotated to keep the aggregates freely suspended and to minimize diffusion-limited gas exchange at the aggregate surface114. Oxygen consumption rates were converted to carbon-specific (C-specific) respiration rates by assuming a respiratory quotient of 1.2 mol O2 to 1 mol CO2115 and by normalization to the aggregate-specific POC contents. Shifts in oxygen concentrations due to, e.g., temperature changes, were tested in control vials with 0.2 µm-filtered water but without aggregates. The observed shift in oxygen concentration was significantly lower in the control vials than in the vials with aggregates (2 ± 1% and 11 ± 6%, N = 6 and 18 vials, respectively, P < 0.001). After incubations, aggregates were filtered individually onto pre-combusted GF/F filters (25 mm, Whatman) and stored at −20 °C for later POC analyses. After storage, filters were fumed over HCl and dried at 50 °C overnight. POC was analyzed after combustion using a carbon analyzer (Surface C-800, Eltra Analysers, 0.01 mg C standard, 2% precision). Carbon respiration rates are given in the units d−1 and % m−1 (Fig. 5g, h), which define the fractional carbon respiration within single aggregates per day and per meter settled, respectively. Respiration rates per meter settled, also referred to as remineralization length scale L6, were calculated as

for aggregates with d = 1–5 mm, based on the size~settling velocity curve fit (Fig. 5f) and mean respiration rates per day (Fig. 5g).

Set III and IV served to determine the infection prevalence and bacterial abundance. Single aggregates and 1 mL of the ambient water were sampled and inspected under the microscope, as described above (see section Culture characteristics). In short, the infection prevalence in Synedra populations was analyzed in Lugol-preserved samples (x400 magnification), whereas abundances of Synedra-associated bacteria were analyzed in PFA-preserved samples (x1000 magnification). Aggregates were disaggregated via gentle shaking in 2 mL microcentrifuge tubes prior to filtration and/or counting, to be able to inspect individual Synedra cells under the microscope. The number of aggregate-associated bacteria (per aggregate) was calculated by multiplying the number of associated bacteria per Synedra-cell type with the abundance of each Synedra-cell type per aggregate. We assumed 20,000 Synedra cells per aggregate since we counted on average 22,843 ± 34,443 Synedra cells per aggregate with d = 2.1 ± 0.8 mm (N = 18).

Set IV was used to derive the mass density of single aggregates ρs-agg. Single aggregates were gently transferred on top of a density gradient. The further procedure resembled the protocol used for ρs-cell, i.e., after centrifugation and settling, the layer containing the aggregate was sampled and analyzed for its mass density (see above, last paragraph in the section Culture characteristics).

A plankton community was sampled from Lake Stechlin (Northern Germany) during the spring bloom in April 2018. Cells were concentrated from 0–10 m using a hand-held plankton net (25 µm mesh size, Hydro-Bios), resuspended in 10 L of in situ surface water, and transferred into three rotating cylinders. Cylinders rotated at 7.5 °C (in situ temperature, 1.5 rpm) for 12 h in darkness until macroscopic aggregates (d = 2–4 mm) were formed. Individual aggregates were sampled with a wide-bore glass pipette, and sub-volumes of the ambient water were sampled with a plastic syringe. Ten replicates (10x aggregates and 10 × 2 mL ambient water) were sampled from each cylinder and split in half. One half was preserved for determining the infection prevalence on single phytoplankton taxa, and the other half for counting cell-associated bacteria. Sample preservation and microscopy analyses were done as described above (see section Culture characteristics).

Pairwise comparisons between mean values were calculated with the Mann–Whitney test for non-normally distributed data, the Welsh-test for normally distributed data with non-equal variance, and the t-test for normally distributed data with equal variance (agricolae v. 1.3.1 in R). Normal distribution was tested using the Shapiro-test and data variance with the F-test. Statistical differences between multiple groups were determined with the Kruskal–Wallis test (if normal distribution was rejected, with Bonferroni correction for P-value adjustment). Statistically significant differences between time series data were calculated based on generalized linear mixed models (GLMM, lmerTest v. 3.1.3 in R). Data were log-transformed if the assumption of homogeneity of variance was rejected (inspected visually on residuals versus fits plots). Pair-wise comparisons in time series data were run with emmeans (v. 1.7.2 in R). Statistical tests and plotting were done in R 3.3.0116 and Origin2021. Statistical tests were run two-sided. The significance level was 0.05. P values at very highly significant levels are given as <0.001 or <0.0001. Uncertainties (±sd) that derived from combined uncertainties of single variables were calculated following the laws of error propagation. The number of replicates is indicated for each mean value in the results section.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data plotted in the figures and mentioned in the text are listed in the supplementary data files (Supplementary Data 1–6). Images used for particle analyses and the image processing code can be downloaded on Dryad https://doi.org/10.5061/dryad.2rbnzs7s7. Further data requests should be directed to the corresponding author.

Falkowski, P. The power of plankton. Nature 483, 17–20 (2012).

Article Google Scholar

Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 281, 237 (1998).

Article CAS PubMed Google Scholar

Grossart, H. P. & Simon, M. Significance of limnetic organic aggregates (lake snow) for the sinking flux of particulate organic matter in a large lake. Aquat. Microb. Ecol. 15, 115–125 (1998).

Article Google Scholar

Weyhenmeyer, G. A. & Bloesch, J. The pattern of particle flux variability in Swedish and Swiss lakes. Sci. Total Environ. 266, 69–78 (2001).

Article CAS PubMed Google Scholar

Fender, C. K. et al. Investigating particle size-flux relationships and the biological pump across a range of plankton ecosystem states from coastal to oligotrophic. Front. Marine Sci. 6, https://doi.org/10.3389/fmars.2019.00603 (2019).

Iversen, M. H., Nowald, N., Ploug, H., Jackson, G. A. & Fischer, G. High resolution profiles of vertical particulate organic matter export off Cape Blanc, Mauritania: Degradation processes and ballasting effects. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 57, 771–784 (2010).

Article CAS Google Scholar

Griffiths, J. R. et al. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Change Biol. 23, 2179–2196 (2017).

Article Google Scholar

Diaz, R. J. & Rosenberg, R. Spreading dead zones and consequences for marine ecosystems. Science 321, 926–929 (2008).

Article CAS PubMed Google Scholar

Jenny, J. P. et al. Global spread of hypoxia in freshwater ecosystems during the last three centuries is caused by rising local human pressure. Glob. Chang Biol. 22, 1481–1489 (2016).

Article PubMed Google Scholar

Carstensen, J., Andersen, J. H., Gustafsson, B. G. & Conley, D. J. Deoxygenation of the Baltic Sea during the last century. Proc. Natl Acad. Sci. 111, 5628–5633 (2014).

Article CAS PubMed PubMed Central Google Scholar

Simon, M., Grossart, H. P., Schweitzer, B. & Ploug, H. Microbial ecology of organic aggregates in aquatic ecosystems. Aquat. Microb. Ecol. 28, 175–211 (2002).

Article Google Scholar

Burd, A. B. & Jackson, G. A. Particle aggregation. Ann. Rev. Mar. Sci. 1, 65–90 (2009).

Article PubMed Google Scholar

Kiørboe, T., Lundsgaard, C., Olesen, M. & Hansen, J. L. S. Aggregation and sedimentation processes during a spring phytoplankton bloom: A field experiment to test coagulation theory. J. Mar. Res. 52, 297–323 (1994).

Article Google Scholar

Boyd, P. W. & Trull, T. W. Understanding the export of biogenic particles in oceanic waters: Is there consensus? Prog. Oceanogr. 72, 276–312 (2007).

Article Google Scholar

Legendre, L. & Rivkin, R. B. Fluxes of carbon in the upper ocean: regulation by food-web control nodes. Mar. Ecol. Prog. Ser. 242, 95–109 (2002).

Article Google Scholar

Kaneko, H. et al. Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean. iScience 24, 102002 (2021).

Article CAS PubMed Google Scholar

Guidi, L. et al. Plankton networks driving carbon export in the oligotrophic ocean. Nature 532, 465 (2016).

Article CAS PubMed PubMed Central Google Scholar

Grossart, H.-P. et al. Fungi in aquatic ecosystems. Nat. Rev. Microbiol. 17, 339–354 (2019).

Article CAS PubMed Google Scholar

Amend, A. et al. Fungi in the marine environment: Open questions and unsolved problems. mBio 10, e01189–01118 (2019).

Article CAS PubMed PubMed Central Google Scholar

Ortiz-Álvarez, R., Triadó-Margarit, X., Camarero, L., Casamayor, E. O. & Catalan, J. High planktonic diversity in mountain lakes contains similar contributions of autotrophic, heterotrophic and parasitic eukaryotic life forms. Sci. Rep. 8, 4457 (2018).

Article PubMed PubMed Central Google Scholar

Gutiérrez, M. H., Pantoja, S., Tejos, E. & Quiñones, R. A. The role of fungi in processing marine organic matter in the upwelling ecosystem off Chile. Mar. Biol. 158, 205–219 (2011).

Article Google Scholar

Edgcomb, V. P., Beaudoin, D., Gast, R., Biddle, J. F. & Teske, A. Marine subsurface eukaryotes: The fungal majority. Environ. Microbiol. 13, 172–183 (2011).

Article CAS PubMed Google Scholar

Frenken, T. et al. Integrating chytrid fungal parasites into plankton ecology: research gaps and needs. Environ. Microbiol. 19, 3802–3822 (2017).

Article PubMed Google Scholar

Van den Wyngaert, S. et al. Seasonality of parasitic and saprotrophic zoosporic fungi: linking sequence data to ecological traits. ISME J. 16, 2242–2254 (2022).

Article PubMed PubMed Central Google Scholar

Gsell, A. S. et al. Long-term trends and seasonal variation in host density, temperature, and nutrients differentially affect chytrid fungi parasitising lake phytoplankton. Freshwat. Biol. https://doi.org/10.1111/fwb.13958 (2022).

Gutiérrez, M. H., Jara, A. M. & Pantoja, S. Fungal parasites infect marine diatoms in the upwelling ecosystem of the Humboldt current system off central Chile. Environ. Microbiol. 18, 1646–1653 (2016).

Article PubMed Google Scholar

Kilias, E. S. et al. Chytrid fungi distribution and co-occurrence with diatoms correlate with sea ice melt in the Arctic Ocean. Commun. Biol. 3, 183 (2020).

Article CAS PubMed PubMed Central Google Scholar

Hassett, B. T., Ducluzeau, A. L. L., Collins, R. E. & Gradinger, R. Spatial distribution of aquatic marine fungi across the western Arctic and sub-arctic. Environ. Microbiol. 19, 475–484 (2017).

Article CAS PubMed Google Scholar

Lepelletier, F. et al. Dinomyces arenysensis gen. et sp. nov. (Rhizophydiales, Dinomycetaceae fam. nov.), a chytrid infecting marine dinoflagellates. Protist 165, 230–244 (2014).

Article PubMed Google Scholar

Hassett, B. T. & Gradinger, R. Chytrids dominate arctic marine fungal communities. Environ. Microbiol. 18, 2001–2009 (2016).

Article CAS PubMed Google Scholar

Le Calvez, T., Burgaud, G., Mahé, S., Barbier, G. & Vandenkoornhuyse, P. Fungal diversity in deep-sea hydrothermal ecosystems. Appl. Environ. Microbiol. 75, 6415–6421 (2009).

Article PubMed PubMed Central Google Scholar

Richards, T. A. et al. Molecular diversity and distribution of marine fungi across 130 european environmental samples. Proc. R. Soc. B Biol. Sci. 282, 20152243 (2015).

Article Google Scholar

Taylor, J. D. & Cunliffe, M. Multi-year assessment of coastal planktonic fungi reveals environmental drivers of diversity and abundance. ISME J. 10, 2118–2128 (2016).

Article CAS PubMed PubMed Central Google Scholar

Comeau, A. M., Vincent, W. F., Bernier, L. & Lovejoy, C. Novel chytrid lineages dominate fungal sequences in diverse marine and freshwater habitats. Sci. Rep. 6, 30120 (2016).

Article CAS PubMed PubMed Central Google Scholar

Wang, Y., Sen, B., He, Y., Xie, N. & Wang, G. Spatiotemporal distribution and assemblages of planktonic fungi in the coastal waters of the Bohai Sea. Front. Microbiol. 9, 584 (2018).

Article PubMed PubMed Central Google Scholar

Gao, Z., Johnson, Z. I. & Wang, G. Molecular characterization of the spatial diversity and novel lineages of mycoplankton in Hawaiian coastal waters. ISME J. 4, 111–120 (2009).

Article PubMed Google Scholar

Duan, Y. et al. A high-resolution time series reveals distinct seasonal patterns of planktonic fungi at a temperate coastal ocean site (Beaufort, North Carolina, USA). Appl. Environ. Microbiol. 84, e00967–00918 (2018).

Article CAS PubMed PubMed Central Google Scholar

Cleary, A. C., Søreide, J. E., Freese, D., Niehoff, B. & Gabrielsen, T. M. Feeding by Calanus glacialis in a high arctic fjord: Potential seasonal importance of alternative prey. ICES J. Mar. Sci. 74, 1937–1946 (2017).

Article Google Scholar

Renaud, P. E., Morata, N., Carroll, M. L., Denisenko, S. G. & Reigstad, M. Pelagic–benthic coupling in the western Barents Sea: Processes and time scales. Deep Sea Res. Part II: Topical Stud. Oceanogr. 55, 2372–2380 (2008).

Article CAS Google Scholar

Lepère, C., Ostrowski, M., Hartmann, M., Zubkov, M. V. & Scanlan, D. J. In situ associations between marine photosynthetic picoeukaryotes and potential parasites – a role for fungi? Environ. Microbiol. Rep. 8, 445–451 (2016).

Article PubMed Google Scholar

Kagami, M., Gurung, T. B., Yoshida, T. & Urabe, J. To sink or to be lysed? Contrasting fate of two large phytoplankton species in Lake Biwa. Limnol. Oceanogr. 51, 2775–2786 (2006).

Article Google Scholar

Gerphagnon, M., Colombet, J., Latour, D. & Sime-Ngando, T. Spatial and temporal changes of parasitic chytrids of cyanobacteria. Sci. Rep. 7, 6056 (2017).

Article PubMed PubMed Central Google Scholar

Ibelings, B. W. et al. Chytrid infections and diatom spring blooms: Paradoxical effects of climate warming on fungal epidemics in lakes. Freshwat. Biol. 56, 754–766 (2011).

Article Google Scholar

Gsell, A. S. et al. Spatiotemporal variation in the distribution of chytrid parasites in diatom host populations. Freshwat. Biol. 58, 523–537 (2013).

Article Google Scholar

Grami, B. et al. Functional effects of parasites on food web properties during the spring diatom bloom in Lake Pavin: A linear inverse modeling analysis. PLOS ONE. 6, e23273 (2011).

Article CAS PubMed PubMed Central Google Scholar

Klawonn, I. et al. Characterizing the "fungal shunt": Parasitic fungi on diatoms affect carbon flow and bacterial communities in aquatic microbial food webs. Proc. Natl Acad. Sci. 118, e2102225118 (2021).

Article CAS PubMed PubMed Central Google Scholar

Kagami, M., Miki, T. & Takimoto, G. Mycoloop: Chytrids in aquatic food webs. Front. Microbiol. 5, 166 (2014).

Article PubMed PubMed Central Google Scholar

Laundon, D. & Cunliffe, M. A call for a better understanding of aquatic chytrid biology. Front. Fungal Biol. 2, https://doi.org/10.3389/ffunb.2021.708813 (2021).

Ploug, H., Iversen, M. H. & Fischer, G. Ballast, sinking velocity, and apparent diffusivity within marine snow and zooplankton fecal pellets: Implications for substrate turnover by attached bacteria. Limnol. Oceanogr. 53, 1878–1886 (2008).

Article Google Scholar

Laurenceau-Cornec, E. C., Trull, T. W., Davies, D. M., De La Rocha, C. L. & Blain, S. Phytoplankton morphology controls on marine snow sinking velocity. Mar. Ecol. Prog. Ser. 520, 35–56 (2015).

Article Google Scholar

Tréguer, P. et al. Influence of diatom diversity on the ocean biological carbon pump. Nat. Geosci. 11, 27–37 (2018).

Article Google Scholar

Alldredge, A. L., Gotschalk, C., Passow, U. & Riebesell, U. Mass aggregation of diatom blooms: Insights from a mesocosm study. Deep Sea Res. Part II: Topical Stud. Oceanogr. 42, 9–27 (1995).

Article CAS Google Scholar

Seto, K., Van den Wyngaert, S., Degawa, Y. & Kagami, M. Taxonomic revision of the genus Zygorhizidium: Zygorhizidiales and Zygophlyctidales ord. nov. (Chytridiomycetes, Chytridiomycota). Fungal Syst. Evol. 5, 17–38 (2020).

Article CAS PubMed Google Scholar

Engel, A. in Practical Guidelines for the Analysis of Seawater (eds Wurl O & Raton B) (CRC Press, 2009).

Cisternas-Novoa, C., Lee, C. & Engel, A. A semi-quantitative spectrophotometric, dye-binding assay for determination of Coomassie Blue stainable particles. Limnol. Oceanogr. Methods. 12, 604–616 (2014).

Article Google Scholar

Passow, U. & Alldredge, A. L. A dye-binding assay for the spectrophotometric measurement of transparent exopolymer particles (TEP). Limnol. Oceanogr. 40, 1326–1335 (1995).

Article CAS Google Scholar

Iversen, M. H. & Ploug, H. Temperature effects on carbon-specific respiration rate and sinking velocity of diatom aggregates - potential implications for deep ocean export processes. Biogeosciences 10, 4073–4085 (2013).

Article Google Scholar

van der Jagt, H., Friese, C., Stuut, J.-B. W., Fischer, G. & Iversen, M. H. The ballasting effect of Saharan dust deposition on aggregate dynamics and carbon export: Aggregation, settling, and scavenging potential of marine snow. Limnol. Oceanogr. 63, 1386–1394 (2018).

Article Google Scholar

Grossart, H. P. & Ploug, H. Microbial degradation of organic carbon and nitrogen on diatom aggregates. Limnol. Oceanogr. 46, 267–277 (2001).

Article CAS Google Scholar

Iversen, M. H. & Ploug, H. Ballast minerals and the sinking carbon flux in the ocean: Carbon-specific respiration rates and sinking velocity of marine snow aggregates. Biogeosciences 7, 2613–2624 (2010).

Article CAS Google Scholar

Ploug, H. & Grossart, H. P. Bacterial growth and grazing on diatom aggregates: Respiratory carbon turnover as a function of aggregate size and sinking velocity. Limnol. Oceanogr. 45, 1467–1475 (2000).

Article CAS Google Scholar

Belcher, A. et al. Depth-resolved particle-associated microbial respiration in the northeast Atlantic. Biogeosciences 13, 4927–4943 (2016).

Article Google Scholar

Ploug, H., Grossart, H. P., Azam, F. & Jørgensen, B. B. Photosynthesis, respiration, and carbon turnover in sinking marine snow from surface waters of Southern California Bight: Implications for the carbon cycle in the ocean. Mar. Ecol. Prog. Ser. 179, 1–11 (1999).

Article CAS Google Scholar

Turner, J. T. Zooplankton fecal pellets, marine snow, phytodetritus and the ocean's biological pump. Prog. Oceanogr. 130, 205–248 (2015).

Article Google Scholar

Nguyen, T. T. H. et al. Microbes contribute to setting the ocean carbon flux by altering the fate of sinking particulates. Nat. Commun. 13, 1657 (2022).

Article CAS PubMed PubMed Central Google Scholar

Zakem, E. J., Cael, B. B. & Levine, N. M. A unified theory for organic matter accumulation. Proc. Natl Acad. Sci. 118, e2016896118 (2021).

Article CAS PubMed PubMed Central Google Scholar

Alcolombri, U. et al. Sinking enhances the degradation of organic particles by marine bacteria. Nat. Geosci. 14, 775–780 (2021).

Article CAS Google Scholar

Buesseler, K. O. & Boyd, P. W. Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean. Limnol. Oceanogr. 54, 1210–1232 (2009).

Article CAS Google Scholar

Henson, S., Le Moigne, F. & Giering, S. Drivers of carbon export efficiency in the global ocean. Glob. Biogeochem. Cycles. 33, 891–903 (2019).

Article CAS Google Scholar

Gsell, A. S., De Senerpont Domis, L. N., Verhoeven, K. J. F., Van Donk, E. & Ibelings, B. W. Chytrid epidemics may increase genetic diversity of a diatom spring-bloom. ISME J. 7, 2057–2059 (2013).

Article PubMed PubMed Central Google Scholar

Agha, R., Saebelfeld, M., Manthey, C., Rohrlack, T. & Wolinska, J. Chytrid parasitism facilitates trophic transfer between bloom-forming cyanobacteria and zooplankton (Daphnia). Sci. Rep. 6, 35039 (2016).

Article CAS PubMed PubMed Central Google Scholar

Rasconi, S. et al. Parasitic chytrids upgrade and convey primary produced carbon during inedible algae proliferation. Protist 171, 125768 (2020).

Article CAS PubMed Google Scholar

Guidi, L. et al. Effects of phytoplankton community on production, size, and export of large aggregates: A world-ocean analysis. Limnol. Oceanogr. 54, 1951–1963 (2009).

Article Google Scholar

Boyd, P. W. & Newton, P. P. Does planktonic community structure determine downward particulate organic carbon flux in different oceanic provinces?. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 46, 63–91 (1999).

Article CAS Google Scholar

van der Jagt, H., Wiedmann, I., Hildebrandt, N., Niehoff, B. & Iversen, M. H. Aggregate feeding by the copepods Calanus and Pseudocalanus controls carbon flux attenuation in the arctic shelf sea during the productive period. Front. Mar. Sci. 7, 543124 (2020).

Article Google Scholar

Steinberg, D. K. et al. Bacterial vs. zooplankton control of sinking particle flux in the ocean's twilight zone. Limnol. Oceanogr. 53, 1327–1338 (2008).

Article Google Scholar

Cavan, E. L., Henson, S. A., Belcher, A. & Sanders, R. Role of zooplankton in determining the efficiency of the biological carbon pump. Biogeosciences 14, 177–186 (2017).

Article CAS Google Scholar

Gachon, C. M. M., Küpper, H., Küpper, F. C. & Šetlík, I. Single-cell chlorophyll fluorescence kinetic microscopy of Pylaiella littoralis (Phaeophyceae) infected by Chytridium polysiphoniae (Chytridiomycota). Eur. J. Phycol. 41, 395–403 (2006).

Article Google Scholar

Senga, Y., Yabe, S., Nakamura, T. & Kagami, M. Influence of parasitic chytrids on the quantity and quality of algal dissolved organic matter (AOM). Water Res. 145, 346––353 (2018).

Article PubMed Google Scholar

Roberts, C., Allen, R., Bird, K. E. & Cunliffe, M. Chytrid fungi shape bacterial communities on model particulate organic matter. Biol. Lett. 16, 20200368 (2020).

Article CAS PubMed PubMed Central Google Scholar

Blackburn, N., Fenchel, T. & Mitchell, J. Microscale nutrient patches in planktonic habitats shown by chemotactic bacteria. Science 282, 2254–2256 (1998).

Article CAS PubMed Google Scholar

Smriga, S., Fernandez, V. I., Mitchell, J. G. & Stocker, R. Chemotaxis toward phytoplankton drives organic matter partitioning among marine bacteria. Proc. Natl Acad. Sci. 113, 1576–1581 (2016).

Article CAS PubMed PubMed Central Google Scholar

Buchan, A., LeCleir, G. R., Gulvik, C. A. & González, J. M. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698 (2014).

Article CAS PubMed Google Scholar

Shibl, A. A. et al. Diatom modulation of select bacteria through use of two unique secondary metabolites. Proc. Natl Acad. Sci. 117, 27445–27455 (2020).

Article CAS PubMed PubMed Central Google Scholar

Guidi, L. et al. Relationship between particle size distribution and flux in the mesopelagic zone. Deep-Sea Res. Part I Oceanogr. Res. Papers. 55, 1364–1374 (2008).

Article CAS Google Scholar

Jackson, G. A. et al. Particle size spectra between 1 μm and 1 cm at Monterey Bay determined using multiple instruments. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 44, 1739–1767 (1997).

Article Google Scholar

Frenken, T. et al. Warming accelerates termination of a phytoplankton spring bloom by fungal parasites. Glob. Change Biol. 22, 299–309 (2016).

Article Google Scholar

Mari, X., Passow, U., Migon, C., Burd, A. B. & Legendre, L. Transparent exopolymer particles: Effects on carbon cycling in the ocean. Prog. Oceanogr. 151, 13–37 (2017).

Article Google Scholar

Passow, U. Transparent exopolymer particles (TEP) in aquatic environments. Prog. Oceanogr. 55, 287–333 (2002).

Article Google Scholar

Prieto, L. et al. Scales and processes in the aggregation of diatom blooms: high time resolution and wide size range records in a mesocosm study. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 49, 1233–1253 (2002).

Article Google Scholar

Kiørboe, T., Andersen, K. P. & Dam, H. G. Coagulation efficiency and aggregate formation in marine phytoplankton. Mar. Biol. 107, 235–245 (1990).

Article Google Scholar

Vidal-Melgosa, S. et al. Diatom fucan polysaccharide precipitates carbon during algal blooms. Nat. Commun. 12, 1150 (2021).

Article CAS PubMed PubMed Central Google Scholar

Gärdes, A., Iversen, M. H., Grossart, H. P., Passow, U. & Ullrich, M. S. Diatom-associated bacteria are required for aggregation of Thalassiosira weissflogii. ISME J. 5, 436–445 (2011).

Article PubMed Google Scholar

Grossart, H. P. & Simon, M. Interactions of planktonic algae and bacteria: Effects on algal growth and organic matter dynamics. Aquat. Microb. Ecol. 47, 163–176 (2007).

Article Google Scholar

Short, S. M. The ecology of viruses that infect eukaryotic algae. Environ. Microbiol. 14, 2253–2271 (2012).

Article PubMed Google Scholar

Carlström, D. The crystal structure of α-chitin (Poly-N-acetyl-d-glucosamine). J. Biophysical Biochemical Cytol. 3, 669–683 (1957).

Article Google Scholar

Miklasz, K. A. & Denny, M. W. Diatom sinkings speeds: Improved predictions and insight from a modified Stokes’ law. Limnol. Oceanogr. 55, 2513–2525 (2010).

Article Google Scholar

Bidle, K. D. & Azam, F. Accelerated dissolution of diatom silica by marine bacterial assemblages. Nature 397, 508 (1999).

Article CAS Google Scholar

Gerphagnon, M. et al. Comparison of sterol and fatty acid profiles of chytrids and their hosts reveals trophic upgrading of nutritionally inadequate phytoplankton by fungal parasites. Environ. Microbiol. 21, 949–958 (2019).

Article CAS PubMed Google Scholar

Kagami, M., Von Elert, E., Ibelings, B. W., De Bruin, A. & Van Donk, E. The parasitic chytrid, Zygorhizidium, facilitates the growth of the cladoceran zooplankter, Daphnia, in cultures of the inedible alga, Asterionella. Proc. R. Soc. B Biol. Sci. 274, 1561–1566 (2007).

Article Google Scholar

Carney, L. T. & Lane, T. W. Parasites in algae mass culture. Front. Microbiol. 5, 1–8 (2014).

Article Google Scholar

Williams, D. M. Synedra, Ulnaria: definitions and descriptions – a partial resolution. Diatom Res. 26, 149–153 (2011).

Article Google Scholar

Arar, E. J. & Collins, G. B. Method 445.0: In vitro determination of chlorophyll and phaeophytin a in marine and freshwater algae by fluorescence. U.S. Environemental Protection Agency, Cinncinnati, Ohio Revision 1.2, 1–22 (1997).

Klawonn, I., Dunker, S., Kagami, M., Grossart, H.-P., Van den Wyngaert, S. Intercomparison of two fluorescent dyes to visualize parasitic fungi (Chytridiomycota) on phytoplankton. Microb. Ecol. 85, 9–23 (2023).

Article CAS PubMed Google Scholar

Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).

Article CAS PubMed PubMed Central Google Scholar

Alldredge, A. L. & Gotschalk, C. In situ settling behavior of marine snow. Limnol. Oceanogr. 33, 351 (1988).

Article Google Scholar

Jackson, G. A. Coagulation in a rotating cylinder. Limnol. Oceanogr. Methods. 13, e10018 (2015).

Article Google Scholar

Shanks, A. L. & Edmondson, E. W. Laboratory-made artificial marine snow: a biological model of the real thing. Mar. Biol. 101, 463–470 (1989).

Article Google Scholar

Cowen, R. K. & Guigand, C. M. In situ ichthyoplankton imaging system (ISIIS): system design and preliminary results. Limnol. Oceanogr. Methods. 6, 126–132 (2008).

Article Google Scholar

Jackson, G. A. & Burd, A. B. Simulating aggregate dynamics in ocean biogeochemical models. Prog. Oceanogr. 133, 55–65 (2015).

Article Google Scholar

Petrik, C. M., Jackson, G. A. & Checkley, D. M. Aggregates and their distributions determined from LOPC observations made using an autonomous profiling float. Deep-Sea Res. Pt I Oceanogr. Res. Papers. 74, 64–81 (2013).

Article Google Scholar

Johnson, C. P., Li, X. & Logan, B. E. Settling velocities of fractal aggregates. Environ. Sci. Technol. 30, 1911–1918 (1996).

Article CAS Google Scholar

Laurenceau-Cornec, E. C. et al. New guidelines for the application of Stokes’ models to the sinking velocity of marine aggregates. Limnol. Oceanogr. 65, 1264–1285 (2020).

Article CAS Google Scholar

Ploug, H. & Grossart, H. P. Bacterial production and respiration in suspended aggregates - A matter of the incubation method. Aquat. Microb. Ecol. 20, 21–29 (1999).

Article Google Scholar

Berggren, M., Lapierre, J.-F. & del Giorgio, P. A. Magnitude and regulation of bacterioplankton respiratory quotient across freshwater environmental gradients. ISME J. 6, 984–993 (2012).

Article CAS PubMed Google Scholar

R.CoreTeam. R: A language and environment for statistical computing. Vienna, Austria. Retrieved from https://www.R-project.org/ (2016).

Download references

We are grateful to Hannah Bachmann, Matthias Bodenlos, Monika Degebrodt, Hannah Gebhardt, Maren Lentz, Uta Mallok, Monika Papke, Solvig Pinnow, Ignacio Rodanes Ajamil, Michael Sachtleben, and Kensuke Seto for their fantastic support during experiments and sample analyses. We especially thank Christiane Hassenrück for statistical advice. Funding was provided by the German Science Foundation DFG grant IV 124/3-1 (to H.P.G. and M.H.I. to fund the project and support I.K.) and the Emmy Noether Project KL 3332/1-1 (to I.K. to fund the project and support I.K.). Further funding was provided by the DFG (GR1540-33-1 to H.P.G. to support the culture maintenance), the Helmholtz Association Young Investigator Group SeaPump VH-NG-1000 (to M.H.I. to support C.M.F. and M.H.I.), the AWI Strategy Fund project "EcoPump" (to C.M.F. to support C.M.F.), the Japan Society for the Promotion of Science (grants JSPS KAKENHI 15KK0026, 16H02943, and 19H05667 to M.K. to support M.K.), and the DFG-Research Center of Excellence "The Ocean Floor – Earth's Uncharted Interface": EXC-2077-390741603 and framework of the HGF Infrastructure Program FRAM of the Alfred-Wegener-Institute Helmholtz Center for Polar and Marine Research (to M.H.I. to support M.H.I.). The MDPI was financed by the German Federal Ministry of Education and Research, BMBF (033W034A to J.C.N.).

Open Access funding enabled and organized by Projekt DEAL.

Silke Van den Wyngaert

Present address: Department of Biology, University of Turku, 20014, Turku, Finland

Clara M. Flintrop

Present address: The Inter-University Institute for Marine Sciences in Eilat, Eilat, 8810302, Israel

Carolina Cisternas-Novoa

Present address: Ocean Sciences Centre, Memorial University of Newfoundland, St. John's, NL, A1C 5S7, Canada

Department of Plankton and Microbial Ecology, Leibniz-Institute of Freshwater Ecology and Inland Fisheries (IGB), 16775, Stechlin, Germany

Isabell Klawonn, Silke Van den Wyngaert, Tim J. W. Walles, Jens C. Nejstgaard & Hans-Peter Grossart

Department of Biological Oceanography, Leibniz Institute for Baltic Sea Research Warnemünde (IOW), 18119, Rostock, Germany

Isabell Klawonn

Alfred Wegener Institute (AWI), Helmholtz Centre for Polar and Marine Research, 27570, Bremerhaven, Germany

Morten H. Iversen & Clara M. Flintrop

Centre for Marine Environmental Sciences (MARUM) and University of Bremen, 28359, Bremen, Germany

Morten H. Iversen & Clara M. Flintrop

Helmholtz Centre for Ocean Research (GEOMAR), 24148, Kiel, Germany

Carolina Cisternas-Novoa

Faculty of Science, Toho University, Funabashi, Chiba, 274‑8510, Japan

Maiko Kagami

Faculty of Environment and Information Sciences, Yokohama National University, Yokohama, Kanagawa, 240‑8502, Japan

Maiko Kagami

Institute of Biochemistry and Biology, Potsdam University, 14469, Potsdam, Germany

Hans-Peter Grossart

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Conceptualization: I.K. and S.Vd.W. Methodology: I.K., S.Vd.W, M.H.I., T.J.W.W., C.C.N., J.C.N., M.K., and H.P.G. Investigation: I.K., S.Vd.W., C.M.F., and M.K. Visualization: I.K. Supervision: I.K., M.H.I., and H.P.G. Writing original draft: IK. Proofreading and editing: I.K., S.Vd.W., M.H.I., T.J.W.W., C.M.F., C.C.N., J.C.N., M.K., and H.P.G.

Correspondence to Isabell Klawonn.

The authors declare no competing interests.

Communications Biology thanks Trang Nguyen, Kyle Mayers and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Linn Hoffmann and David Favero. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Klawonn, I., Van den Wyngaert, S., Iversen, M.H. et al. Fungal parasitism on diatoms alters formation and bio–physical properties of sinking aggregates. Commun Biol 6, 206 (2023). https://doi.org/10.1038/s42003-023-04453-6

Download citation

Received: 19 October 2022

Accepted: 10 January 2023

Published: 21 February 2023

DOI: https://doi.org/10.1038/s42003-023-04453-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.