Modeling the elemental stoichiometry and silicon accumulation in diatoms

Diatoms are important microorganisms involved in global primary production, nutrient cycling, and carbon sequestration. A unique feature of diatoms is their silica frustules, which impact sinking speed, defense against predators and viruses, and growth cycling. Thus, frustules are inherently linked to their role in ecosystems and biogeochemical cycles. However, constraints on cellular silicon levels remain unclear and few existing models resolve diatom elemental stoichiometry to specifically include variable silicon levels. Here, we use a coarse-grained model of the diatom, Thalassiosira pseudonana, compared with laboratory results to illustrate the relationship of silicon uptake with elemental stoichiometry of other nutrients. The model-data comparison suggests the balance between growth rate and silicon uptake constrains the amount of cellular silicon. Additionally, it expresses relationships between silicon, nitrogen, phosphorus, and carbon to changing growth rates in nitrogen-limited and phosphorus-limited regimes. First, our model-data comparison suggests Si uptake hits a maximum cellular quota at low growth rates and below this maximum there is independent Si uptake. In each nutrient regime, Si:N, Si:P, and Si:C decrease exponentially with growth rate when Si is below the maximum limit. This is explained by independent Si uptake and increased loss of Si to new cells. These results provide predictions of diatom stoichiometry and allocation, which can be used in ecosystem models to differentiate phytoplankton types to better represent diatoms' contribution to global biogeochemical cycles and ecosystems.