Contrasting size evolution in marine and freshwater diatoms

Diatoms are key players in the global carbon cycle and most aquatic ecosystems. Their cell sizes impact carbon sequestration and energy transfer to higher trophic levels. We report fundamental differences in size distributions of marine and freshwater diatoms, with marine diatoms significantly larger than freshwater species. An evolutionary game theoretical model with empirical allometries of growth and nutrient uptake shows that these differences can be explained by nitrogen versus phosphorus limitation, nutrient fluctuations and mixed layer depth differences. Constant and pulsed phosphorus supply select for small sizes, as does constant nitrogen supply. In contrast, intermediate frequency nitrogen pulses common in the ocean select for large sizes or the evolutionarily stable coexistence of large and small sizes. Size-dependent sinking interacts with mixed layer depth (MLD) to further modulate optimal sizes, with smaller sizes selected for by strong sinking and shallow MLD. In freshwaters, widespread phosphorus limitation, together with strong sinking and shallow MLD produce size distributions with smaller range, means and upper values, compared with the ocean. Shifting patterns of nutrient limitation and mixing may alter diatom size distributions, affecting global carbon cycle and the structure and functioning of aquatic ecosystems.