# A Mechanistic Model of Macromolecular Allocation, Elemental Stoichiometry, and Growth Rate in Phytoplankton

^{1}School of Oceanography, University of Washington, Seattle, WA, United States^{2}Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA, United States^{3}Department of Microbiology, University of Tennessee, Knoxville, Knoxville, TN, United States^{4}National Herbarium of New South Wales, The Royal Botanic Gardens and Domain Trust, Sydney, NSW, Australia

We present a model of the growth rate and elemental stoichiometry of phytoplankton as a function of resource allocation between and within broad macromolecular pools under a variety of resource supply conditions. The model is based on four, empirically-supported, cornerstone assumptions: that there is a saturating relationship between light and photosynthesis, a linear relationship between RNA/protein and growth rate, a linear relationship between biosynthetic proteins and growth rate, and a constant macromolecular composition of the light-harvesting machinery. We combine these assumptions with statements of conservation of carbon, nitrogen, phosphorus, and energy. The model can be solved algebraically for steady state conditions and constrained with data on elemental stoichiometry from published laboratory chemostat studies. It interprets the relationships between macromolecular and elemental stoichiometry and also provides quantitative predictions of the maximum growth rate at given light intensity and nutrient supply rates. The model is compatible with data sets from several laboratory studies characterizing both prokaryotic and eukaryotic phytoplankton from marine and freshwater environments. It is conceptually simple, yet mechanistic and quantitative. Here, the model is constrained only by elemental stoichiometry, but makes predictions about allocation to measurable macromolecular pools, which could be tested in the laboratory.

## Introduction

Phytoplankton are responsible for the majority of photosynthesis in the ocean (Field et al., 1998) and more than half in lakes (Vadeboncoeur et al., 2002). The elemental stoichiometry of phytoplankton varies significantly through acclimation and adaptation (Quigg et al., 2003, 2011; Finkel et al., 2016), modulates fitness in different environments (Deutsch and Weber, 2012), global ocean carbon storage (Galbraith and Martiny, 2015), and the nutrition of higher trophic levels (Mitra et al., 2007). Population growth rates of phytoplankton depend on resource availability (Caperon and Meyer, 1972a,b; Paasche, 1973; Laws and Bannister, 1980; Pedersen and Borum, 1996; Xu et al., 2010) and also vary through acclimation and adaptation (Falkowski and Owens, 1980; Levasseur et al., 1993; Litchman et al., 2002, 2003; Collos et al., 2005; Litchman and Klausmeier, 2008; Van Mooy et al., 2009; Lewis et al., 2019). The environmentally dependent growth rate of a population is an important component of its fitness and significant for ecological and biogeochemical modeling.

The elemental stoichiometry and growth rate of phytoplankton are not independent. Robust qualitative relationships between growth rate, elemental stoichiometry, and resource availability are evident in controlled laboratory cultures spanning wide taxonomic and allometric ranges. We illustrate this in Figure 1, Supplementary Figure 1 with data compiled from published, continuous culture laboratory studies of 12 species, including marine, freshwater, prokaryotic, and eukaryotic phytoplankton. In all cases, at a fixed irradiance, Chl:C (chlorophyll per carbon) increases linearly with growth rate, μ (Laws and Bannister, 1980; Healey, 1985; Sakshaug and Andersen, 1989; Chalup and Laws, 1990), and both the slope and the intercept increase as the irradiance declines (Figures 1A–C) (quantitative fits with *R*^{2} values in Supplementary Table 1). Similarly, cellular N:C (nitrogen:carbon) increases linearly with growth rate (Caperon and Meyer, 1972a; Laws and Caperon, 1976; Laws and Bannister, 1980; Healey, 1985; Sakshaug and Andersen, 1989; Chalup and Laws, 1990; Figures 1D–G) and its slope and intercept both increase with decreasing photon flux (Healey, 1985; Sakshaug and Andersen, 1989; Chalup and Laws, 1990; Figures 1D,F, *R*^{2} values in Supplementary Table 2). In contrast, the cellular P:C (phosphorus:carbon) increases non-linearly with growth rate (Figure 1H).

**Figure 1**. Compiled laboratory data of growth rate and light dependence of chlorophyll and elemental stoichiometry, and light dependence of ${\mu}_{max}^{I}$: nutrient replete growth rate. Data are all from chemostat cultures where the growth rate, μ, is controlled by the dilution rate, *D*. Here, illustrated N:C data were all under N limitation and P:C were all under P limitation. Points represent the original data and curves represent regression lines. **(A**–**C)** Chlorophyll:C for various light intensities for *Pavlova lutheri* (Chalup and Laws, 1990), *Skeletonema costatum* (Sakshaug and Andersen, 1989), and *Synechococcus linearis* (Healey, 1985). Legends indicate light intensities in μmol m^{−2} s^{−1}. **(D**–**F)** N:C for various light intensities for the same organisms as **(A**–**C)**, respectively. Legend values are light intensities as in **(A**–**C). (G)** N:C for other organisms. (a) *Synechococcus* (WH8102) (Garcia et al., 2016). (b) *Synechococcus* (WH7803) (Liu et al., 1999). (c) *Coccochloris stagnina* (Caperon and Meyer, 1972a). (d) *Thalassiosira pseudonana* (Claquin et al., 2002). (e) *Dunaliella tertiolecta* (Caperon and Meyer, 1972a). (f) *Monochrysis lutheri* (Caperon and Meyer, 1972a). **(H)** P:C for three organisms. (g) *Synechococcus* (WH8102) (Garcia et al., 2016). (h) *Selenastrum minutum* (Elrifi and Turpin, 1985). (i) *Synechococcus linearis* (Healey, 1985). In **(G,H)** a light intensity varied between experiments. **(I)** ${\mu}_{max}^{I}$- light relationships for three organisms as in **(A**–**C)**. (j) *Pavlova lutheri* (Chalup and Laws, 1990). (k) *Skeletonema costatum* (Sakshaug and Andersen, 1989). (l) *Synechococcus linearis* (Healey, 1985). For *Synechococcus linearis*, only data points where the limiting nutrient was fully consumed are plotted (this applies to other figures).

Here we define ${\mu}_{max}^{I}$ as maximum growth rate for a given light intensity. ${\mu}_{max}^{I}$ can be also considered as a nutrient replete growth rate. In the chemostat culture, as the dilution rate increases, the rate of nutrient input increases. Despite the increased rate, when the dilution rate is above ${\mu}_{max}^{I}$, the cells are flushed away, since cellular growth cannot increase further. Thus, ${\mu}_{max}^{I}$ is indicated by the termination of the linear increase in Chl:C and N:C, and the termination of the non-linear increase in P:C with μ (Figures 1A–H, and indicated schematically by the dotted line in Figures 2A–C). In Figure 1I, we plot ${\mu}_{max}^{I}$ as a function of light intensity revealing the typical saturation of growth rate at high light intensities (Healey, 1985) (we note that none of the illustrated experiments were in a regime of photo-inhibition). Figures 1A–I, Supplementary Figure 1 thus reveal a set of robust qualitative relationships between light intensity, growth rate, and the elemental stoichiometry of diverse phytoplankton under steady-state growth conditions (summarized schematically in Figure 2).

**Figure 2**. Illustration of general trends in laboratory data in chemostat culture studies. Growth rate (μ) and light dependence of **(A)** Chl:C, **(B)** N:C under N limitation, and **(C)** P:C under P limitation, respectively, and **(D)** ${\mu}_{max}^{I}$–light relationships. Whereas, Chl:C and N:C are described by linear curves, P:C has non-linear relationships with μ. Black dotted lines represent ${\mu}_{max}^{I}$.

The common patterns in Figure 1 reflect the fact that, despite the diversity of species represented, there are shared physiological underpinnings. The elemental stoichiometry of a cell depends on the relative abundances of the macromolecules from which it is composed and they, in turn, are linked to environment and physiological state (Sterner and Elser, 2002). The C:N:P stoichiometry of phytoplankton can be largely accounted for by the sum of contributions from a handful of major macromolecular components: protein, pigment, carbohydrate, lipid, DNA, RNA, and storage molecules (Liefer et al., 2019) each of which has a distinct elemental stoichiometry (see Table 1). For example, proteins are relatively rich in nitrogen so increasing the cellular allocation to protein typically raises cellular N:C (Sterner and Elser, 2002; Klausmeier et al., 2004). The broad-brush response of the macromolecular allocation of phytoplankton to changes in environmental factors is common across broad taxonomic groupings; for example in laboratory studies of nitrogen starvation in four marine species (Liefer et al., 2019) and with changing temperature and growth rate amongst a wide variety of freshwater phytoplankton (Fanesi et al., 2017, 2019).

Models of phytoplankton physiology have sought to relate growth rate (related to fitness) and elemental stoichiometry (related to biogeochemical impacts) to external resource availability (Riley, 1946; Monod, 1949), internal stores of resources (Caperon, 1968; Droop, 1968), and the internal allocation between functional pools and storage molecules (Shuter, 1979; Geider et al., 1998; Kooijman, 2010). Recent models also explicitly represent trade-offs associated with allocation of the resource and proteome (Bonachela et al., 2013; Burnap, 2015; Smith et al., 2016; Reimers et al., 2017; Chen and Smith, 2018; Faizi et al., 2018; Jahn et al., 2018; Faizi and Steuer, 2019). We provide a more comprehensive review of published physiological models in Supplementary Text. The model presented here also aims to explicitly connect growth rate, elemental stoichiometry, and environmental conditions. It is based on the allocation of resources between and within the major macromolecular pools. We seek to frame the model in terms of measurable (rather than abstracted) pools, to provide interpretations of observed laboratory relationships, and to keep the model efficient and simple for practical applications.

## Model Description

In Figure 3A we sketch the broad-brush allocation of C, N, and P to key macromolecular pools in the phytoplankton model (Cell Flux Model of Phytoplankton: CFM-Phyto). Cells also allocate resources within the macromolecular pools. For example, lipids incorporate lipid membranes and lipids storage molecules (Shifrin and Chisholm, 1981; Lengeler et al., 1999). The protein pool includes enzymes devoted to hundreds of reactions which may be coarse-grained into several major categories including those related to light-harvesting and electron transport, and those related to biosynthesis, growth, and reproduction (Figure 3A). Recent proteomic analyses are quantifying broad-brush protein allocation (McKew et al., 2013, 2015; Christie-Oleza et al., 2017; Jahn et al., 2018; Zavřel et al., 2019) in ways which connect to such coarse-grained models (Scott et al., 2010; Burnap, 2015; Reimers et al., 2017; Faizi et al., 2018; Faizi and Steuer, 2019). Recent studies have revealed that a large and highly variable fraction of phytoplankton proteome is devoted to light-harvesting and electron transport (Jahn et al., 2018; Zavřel et al., 2019) and this investment increases as light intensity decreases. The macromolecular pools identified in Figure 3A are potentially measurable (Scott et al., 2010; McKew et al., 2013, 2015; Jahn et al., 2018; Fanesi et al., 2019; Zavřel et al., 2019), which is useful for testing and calibrating. At the same time, it is useful to recast these macromolecular pools into “functional units.” For example, in low light, cyanobacteria allocate to increase light harvesting proteins, but also pigments and lipids in the thylakoid membrane. We illustrate this re-organization in Figure 3B. Allocation-based models of phytoplankton populations have abstracted the system at this level (Shuter, 1979; see Supplementary Text), providing mechanistic representations of trade-offs but which are more difficult to directly constrain from observations. The relationship between Figures 3A,B indicates how the two perspectives are compatible. The model we outline below is developed in terms of the measurable pools indicated in Figure 3A, but interpreted in terms of their aggregated, functional allocation as depicted in Figure 3B.

**Figure 3**. Schematic of the two different views of the model: CFM-Phyto. **(A)** Allocation of C, N, and P to key macromolecular pools. Orange outline, cell membrane layers; green background is cytoplasm. Arrows represent elemental fluxes; yellow, C; red, N; blue, P. Colors in boxes represent elements that are influenced by each molecule most. Yellow: C. Red: N. Blue: P. Black dotted lines indicate possible intramolecular associations of C and P. Orange outline, cell membrane layers; green background, cytoplasm; *Chl*, chlorophyll; *thyla*, thylakoid; *CH*, carbohydrate. **(B)** Simple view of macromolecular allocation grouped into four different functions; *Bio*, biosynthetic molecules; *Photo*, photosynthetic molecules; *Other*, other constant molecules; *Store*, storage molecules. *Bio* affects growth rates and *Photo* affects photosynthesis rates. Essential in *Other* indicates essential lipids and carbohydrates. Black dashed arrows indicate positive influences. Red dashed arrow represents light.

In the following sections, we outline an idealized, allocation-based model of phytoplankton physiology and growth rate under a range of resource conditions (N, P, light). We show that the observed relationships between Chl:C, N:C, P:C growth rate and light (Figure 1, Supplementary Figure 1; summarized schematically in Figure 2) can be quantitatively modeled by understanding carbon allocation between and within the major macromolecular pools. By considering the allocation of nitrogen, we find an interpretation for the linear relationship of N:C with growth rate and its dependence on light intensity. By relating allocation of phosphorus to the rate of biosynthesis, we model and interpret the non-linear relationship between P:C and growth rate. The model is developed with particular reference to Healey's study of *Synechococcus linearis*, a freshwater cyanobacterium (Healey, 1985; Figures 1C,F,i,l), which provides a comprehensive set of constraints on elemental stoichiometry at multiple growth rates and light intensities, under both N and P limitation. In the experimental data, there are no direct constraints on macromolecular allocation, so we infer the latter through combination of observed elemental stoichiometry and model structure and discuss the inferred macromolecular allocation with reference to other published studies. While the model has been developed by exploiting the comprehensive data set of Healey (1985), the physiology of allocation at this level is common across taxa and it can be fitted to data from other phytoplankton, as we also illustrate.

## Modeled Macromolecular Composition of the Cell

### Carbon Allocation

Laboratory studies have shown that almost all the cellular carbon in phytoplankton is accounted for by the major macromolecular pools: proteins (*C*_{Pro}), chlorophyll and other pigments (*C*_{Chl}), nucleic acids (*C*_{Nuc}), carbohydrates (*C*_{Carb}), and lipids (*C*_{Lip}) (Anderson, 1995; Liefer et al., 2019). The carbon quota of a phytoplankton cell, *C*_{Cell} (mol C cell^{−1}) is thus defined as the sum of these components along with carbon associated with nitrogen storage molecules, *C*_{Nsto}:

While chlorophyll is a relatively minor contribution in this regard, it provides a constraint on the light harvesting capacity of the cells and is routinely measured. Here we neglect the contribution from intra-cellular, dissolved metabolites which are typically minor [e.g., ~4% of cellular dry weight in *E. Coli* (Lengeler et al., 1999) and predicted to be only ~1% as inorganic ions in *Synechococcus elongatus* (PCC 7942) (Reimers et al., 2017)].

Elemental variations referenced to carbon (C:C, N:C, P:C) present clearer relationships with μ and light than cellular quotas (Caperon and Meyer, 1972a; Laws and Bannister, 1980; Elrifi and Turpin, 1985), so we define carbon normalized macromolecular composition by dividing both sides of Equation (1) by *C*_{Cell}:

where ${Q}_{C}^{i}$ are ratios with units of (mol C mol C^{−1}). Proteins account for a large fraction of the carbon and nitrogen in a phytoplankton cell (Anderson, 1995; Geider and La Roche, 2002).

Recent quantitative proteomics studies revealed a coarse-grained reorganization of the proteome of *Synechocystis* in response to a changing light environment: an increase in photon flux drove an increase in growth rate with an associated downregulation of light harvesting machinery and upregulation of translational machinery (McKew et al., 2013; Jahn et al., 2018; Zavřel et al., 2019). This motivates the resolution of protein pools related to biosynthesis, ${Q}_{C}^{Pro-Bio},$ and photosynthesis, ${Q}_{C}^{Pro-Pho}$, the latter including contributions from light absorbing antennas, as well as proteins for photosystems and electron transport. Proteomic studies have shown that light-harvesting proteins contribute as much as 38% of the proteome of *Synechococcus* sp. WH7803 in culture (Christie-Oleza et al., 2017) and can vary considerably (Jahn et al., 2018; Zavřel et al., 2019). We seek to exploit these observed proteomic trade-offs and so resolve the allocation of protein into three functional pools:

We also resolve a fixed-size pool of “essential” proteins, ${Q}_{C}^{Pro-Other}$ which is necessary to close the cellular budget and notionally includes enzymes associated with essential metabolism (Jahn et al., 2018; Zavřel et al., 2019) and structure.

Nucleic acids include contributions from DNA and RNA:

where the contribution from RNA is significantly more variable and related to growth rate; discussed in more detail below. Intracellular dissolved pools are not resolved since they generally represent <5% of the total cellular mass (Lengeler et al., 1999).

The lipid pool can be separated into three components. A large fraction of thylakoid membrane is lipid (~30%) (Kirchhoff, 2014), and we resolve a phospholipid fraction of it, ${Q}_{C}^{Plip-Thy}$, which also contributes significantly to the cellular phosphorus budget. ${Q}_{C}^{Lip-Sto}$ is a flexible pool of storage molecules (Werner, 1977; Shifrin and Chisholm, 1981), and ${Q}_{C}^{Lip-Other}\text{}$represents essential structural components of the cell membrane (Neidhardt et al., 1990; Lengeler et al., 1999), which we consider as a non-flexible pool. Hence, we resolve three lipid pools:

We represent the total cellular carbohydrate pool as the sum of two contributions: a flexible component, ${Q}_{C}^{Carb-Sto}$, representing storage (Shifrin and Chisholm, 1981; Deschamps et al., 2008; Dron et al., 2012) and the pool of essential carbohydrate metabolites (Lengeler et al., 1999; Michal, 1999), along with an “essential,” fixed carbohydrate contribution, ${Q}_{C}^{Carb-Other}$ (Harrison et al., 1990; Anderson, 1995; Biersmith and Benner, 1998):

### Nitrogen Allocation

Cellular nitrogen (*N*_{Cell}) is mostly associated with protein (*N*_{Pro}) (Anderson, 1995; Liefer et al., 2019), along with contributions from RNA (*N*_{RNA}), DNA (*N*_{DNA}), chlorophyll (*N*_{Chl}), and nitrogen storage (*N*_{Sto}):

Since N:C presents clearer relationships with μ and light than *N*_{Cell} (Caperon and Meyer, 1972a; Laws and Bannister, 1980; Elrifi and Turpin, 1985), we divide both sides of Equation (7) by *C*_{Cell}:

where ${Q}_{N}^{i}$ are ratios with units of (mol N mol C^{−1}). Each of the macromolecular pools has a distinct elemental stoichiometry (see Table 1), and carbon-normalized nitrogen quotas (${Q}_{N}^{i}$, mol N mol C^{−1}) are constructed accordingly. For example, the total nitrogen content of cellular protein, ${Q}_{N}^{Pro}={Q}_{C}^{Pro}{Y}_{Pro}^{N:C}$, where ${Y}_{Pro}^{N:C}$ is the empirically-informed, average N:C of protein (see Table 1).

### Phosphorus Allocation

Nucleic acids (*P*_{RNA} and *P*_{DNA}), phospholipids in the thylakoid membrane (*P*_{Thy}), and storage compounds including polyphosphate (*P*_{Sto}) are observed to account for most of the cellular phosphorus in phytoplankton and bacteria (Anderson, 1995; Lengeler et al., 1999; Table 1). Phosphorus may also be distributed in non-photosynthetic phospholipids and associated with other molecules (e.g., phosphorylation; Lengeler et al., 1999) which are here assumed to be in a fixed pool, *P*_{Other0}. We also account for the flexible part of non-Thylakoid P-lipid in *P*_{Sto}. Here we resolve phosphorus allocation to these distinct pools:

We note that a full accounting for the phosphorus in phytoplankton has not been experimentally characterized to date (Moreno and Martiny, 2018; Liefer et al., 2019). As we have done for *C*_{Cell} and *N*_{Cell}, we divide both sides of Equation (9) and obtain P:C:

where ${Q}_{P}^{i}$ are ratios with units of (mol P mol C^{−1}). As for N:C, ${Q}_{P}^{i}$ is constructed according to the distinct elemental stoichiometry of each macromolecule (Table 1).

### Representing Relationships Between Macromolecular Pools and Rates

In addition to these statements of mass conservation and allocation, we must connect macromolecular allocation to rates. We do this assuming four mathematical representations of which three are well-supported by laboratory observations:

(i) *The per chlorophyll gross rate of photosynthesis, v*_{I} *(mol C (mol C in Chl)*^{−1} *d*^{−1}*) is a saturating function of irradiance I (*μ*mol m*^{−2} *s*^{−1}*)*. Following established models rooted in empirical observations and target theory (Cullen, 1990; Geider et al., 1998), we model photosynthesis as a function of light intensity:

Here ${v}_{I}^{max}$ is the maximum photosynthesis rate per chlorophyll and *A*_{I} is a coefficient characterizing the absorption cross-section and turnover time of the photosynthetic unit (Cullen, 1990).

(ii) *The components of the photosynthetic machinery, namely chlorophyll, light-related proteins, and the thylakoid phospholipids co-vary linearly*. In other words, the composition of the thylakoid apparatus remains constant but its amount per cell is varied with acclimation. Hence, the allocations to photosynthetic protein and thylakoid phospholipids are assumed linearly proportional to cellular chlorophyll content:

and

where *A*_{Pho} and ${A}_{Pho}^{P:Chl}$ are constants of proportionality. It is observed that the size of thylakoid membranes increases under low light in phytoplankton (Geider et al., 1996) and the chloroplasts of plants (Lichtenthaler et al., 1982). The thylakoid membranes are generally highly crowded with proteins (Folea et al., 2008; Kirchhoff et al., 2008; Kirchhoff, 2014) but how the fraction of proteins might change with growth conditions is less clear. Thus, we choose the simplest model and assume a fixed composition of the light harvesting apparatus.

(iii) *Investment in biosynthetic protein is proportional to growth rate:*

This is consistent with the observed linear increase in the investment in ribosomal proteins with growth rate in multiple cultures of *Synechocystis* (Jahn et al., 2018; Zavřel et al., 2019). Cultures of *Scenedesmus* sp. (Rhee, 1978), also show a near-linear relationship between protein-based nitrogen and growth rate under constant light. We note that ${Q}_{C}^{Pro-Bio}$ represents not only ribosomal proteins but includes those involved in synthesis of lipid and nucleic acids, C metabolism, and cell division.

(iv) *The investment in RNA*, ${Q}_{P}^{RNA}$, *varies linearly with total protein content and with growth rate*. This relationship is derived from the observation that the RNA:protein ratio is linear with growth rate in phytoplankton (Nicklisch and Steinberg, 2009; Liefer et al., 2019), as is the case for heterotrophic bacteria (Bremer and Dennis, 1996; Scott et al., 2010). Thus, we model investment in RNA as

where ${Q}_{P,\text{}min}^{RNA}$ is the minimum RNA (mol P mol C^{−1}), which occurs at zero growth rate. This relationship says that cells need more RNA to divide faster and/or to reproduce a higher cellular protein quota.

Relationships in (i), (iii), and (iv) above are directly supported by empirical data in the associated citations. The relationship between components of the light harvesting and photosynthesis machinery in (ii) is logical and simple, but unconfirmed by direct empirical data to our knowledge.

Using the above statements of mass conservation (Equations 12, 13) and representations of key relationships between fluxes and pools (Equations 14, 15), we model the observed dependencies of cellular stoichiometry (i.e., Chl:C, N:C, and P:C) on growth rate, light intensity and limiting factor, as well as the variation of maximum growth rate, ${\mu}_{max}^{I}$, under different light intensities. In the following sections we outline the model, making some approximations to the full equations presented above, in order to provide an illustrative and instructive discussion. A complete approach to solution of the model is presented in the Methods section and the solutions presented in all figures were generated using the un-approximated forms. Final equations of the un-approximated model are also summarized in tabular form in Supplementary Table 3. We frame our discussion of the model and its application by seeking to explain the trends identified in Figures 1, 2.

### Model Representation and Analysis

#### Why Does Chl:C Vary Linearly With Growth Rate?

Consider the rate of change of the cellular carbon quota, which is increased by photosynthesis, and reduced by division with population growth rate μ (d^{−1}) and maintenance respiration rate *m* (d^{−1}) (e.g., Geider et al., 1998; Pahlow and Oschlies, 2009):

where ${Q}_{C}^{Cell}$ is cellular chlorophyll to cellular carbon ratio [(mol C in Chl) mol C^{−1}], and *v*_{I} is the per chlorophyll rate of photosynthesis as defined in Equation (11). *E* is the respiratory cost of synthesis (moles C respired per mole C synthesized). *E* is estimated based on the production of biomass with stoichiometry of C_{5}H_{7}O_{2}N_{1}P_{1/30} using nitrate as the nitrogen source with energy transfer efficiency of 0.6 (Rittmann and McCarty, 2001). The assumed elemental stoichiometry is based on the suggested values of C:H:O:N (Rittmann and McCarty, 2001) and within the range observed in the laboratory experiments (Healey, 1985). Idealized models have typically assumed that *E* is proportional to cellular nitrogen content, assuming associated costs with nitrate reduction and protein synthesis (Laws and Wong, 1978; Geider et al., 1998; Pahlow and Oschlies, 2009). However, many aspects of metabolism consume ATP, including the synthesis of lipid and carbohydrate that contain little nitrogen (Lengeler et al., 1999; Michal, 1999) so here we assume the respiratory cost of biosynthesis is proportional to cellular carbon (i.e., *E**C*_{Cell} in Equation 16).

In steady-state, the solution of Equation (16) anticipates the observed linear relationship between the cellular chlorophyll to carbon ratio (${Q}_{C}^{Chl}=Chl\text{}:C$) and growth rate for any given photon flux, *I* (a similar relationship as in Laws and Bannister, 1980):

where *A*_{Chl}(*I*) = (1 + *E*)/*v*_{I}(*I*) and *B*_{Chl}(*I*) = *m*/*v*_{I}(*I*). The predicted linear trend in Equation (17) is consistent with the data in Figure 1. We fitted Equation (17) to the data of Healey (1985) using the Metropolis-Hastings procedure (Metropolis et al., 1953; Hastings, 1970; Omta et al., 2017), a Markov Chain Monte Carlo method which is described in detail in Methods. The values of the parameters *m*, ${v}_{I}^{max},$ and *A*_{I} were optimized in the model-data fit, with single values for ${v}_{I}^{max}$ and *A*_{I} for all light levels.

Equation (17), indicates that the cell must maintain pigments to sustain maintenance respiration even at a net zero growth rate (y-intercept). It predicts a linear relationship between Chl:C and μ, as well as an increase in both the slope and intercept of the Chl:C ratio with decreasing irradiance (Figure 4). The linear trend reflects the increased investment in light-harvesting machinery to maintain the same growth rate or maintenance costs at different light levels. The form of the model is also consistent with data from cultures of other phytoplankton (*Pavlova lutheri* and *Skeletonema costatum*; Sakshaug and Andersen, 1989; Chalup and Laws, 1990), as illustrated by model simulations in Supplementary Figure 2. Hence the model is not specific to a single organism, though each species requires different parameter values, representing inter-species variations in traits. We discuss the plausibility of inferred parameter values at the end of the Results section (see the subsection “Plausibility of the predicted macromolecular allocation”). We note that Equation (17) does not predict ${\mu}_{max}^{I}$, which requires additional constraints brought to bear through N and/or P limitation when the chlorophyll, C, N, and P conservation equations are coupled (see the subsection “What is the maximum growth rate for a given light intensity, ${\mu}_{max}^{I}$?”).

**Figure 4**. Model-data comparison of chlorophyll per C of *Synechococcus linearis* for various growth rates and light intensities. Curves: model solution (Equation 17). Points: data (Healey, 1985; circles, N limited; diamonds, P limited). Dotted line represents ${\mu}_{max}^{I}$ at various light intensities; high ${\mu}_{max}^{I}$ for higher light intensity. μ in the *x* axis represents growth rates (d^{−1}). Legend values are light intensities (μmol m^{−2} s^{−1}).

#### Why Does N:C Vary Linearly With Growth Rate, and Why Does It Change With Photon Flux?

Consider the case of nitrogen limitation, when allocation to nitrogen storage is small. In this case, to a first approximation, the cellular quota of nitrogen is dominated by that of protein (Liefer et al., 2019).

Nucleic acids account for <7% of cellular dry weight in phytoplankton (Parsons et al., 1984; Anderson, 1995), small relative to that in heterotrophic bacteria. Thus, for explanatory purposes we will use an approximate form of Equation (8) that considers only the contribution from protein (a full, un-approximated solution is provided in the Methods section and is used in the figures of model solutions). In the case where most nitrogen is associated with protein, Equation (8) becomes:

Here nitrogen and carbon are linked by the constant elemental ratio for protein, ${Y}_{Pro}^{N:C}$ (Table 1). Other macromolecules have different elemental stoichiometries. ${Q}_{C}^{Pro-Other}$ represents the fixed, minimum complement of protein essential for the cell. The cellular investment in the photosynthetic protein is assumed to vary linearly with chlorophyll, and the investment in biosynthetic protein pools is assumed to vary linearly with growth rate using Equations (12, 14), respectively, as discussed above (in “Representing relationships between macromolecular pools and rates”). Combining Equations (12, 14, 17, 18) leads to an expression that describes the relationship between N:C of the population and growth rate under N-limiting conditions:

Equation (19) predicts a linear relationship between N:C and growth rate which has decreasing slope and intercept as with photon flux, qualitatively consistent with the observed data in Figures 1D–F, 2B, 5A, Supplementary Figures 1, 2. Extending the model to include N storage allows predictions of N:C under P-limitation (Supplementary Figure 4) which is discussed later.

**Figure 5**. Model-data comparison of N:C and model prediction of macromolecular allocation of *Synechococcus linearis*. **(A)** N:C under N limitation under different light intensities; light intensities (μmol m^{−2} s^{−1}) are in the legend. Curves are model results and points are data (Healey, 1985). Dotted lines represent ${\mu}_{max}^{I}$ at various light intensities; high ${\mu}_{max}^{I}$ for higher light intensity. **(B)** Macromolecular allocation in N, normalized by cellular C under the light intensity of 144 (μmol m^{−2} s^{−1}). Black points are data for total values under the same light intensity (Healey, 1985). See the legend for color definitions: *Bio*, biosynthetic protein + RNA; *Photo*, chlorophyll + photosynthetic protein; *Other*, other molecules. *Bio* and *Photo* is dominated by biosynthetic protein and photosynthetic protein, respectively. Detailed macromolecule allocation in Supplementary Figure 3 and N:C under P limitation in Supplementary Figure 4.

The model suggests that N:C increases with growth rate at a fixed light intensity because there is a linear increase in the investment in both biosynthetic protein (Equation 14) and photosynthetic protein, latter being the associated with the linear increase in Chl:C with growth rate at fixed light intensity (Equations 12, 17). Likewise, a reduction of light intensity at a fixed growth rate also demands a higher investment in both chlorophyll and photosynthetic proteins, hence the slope of N:C increases with decreasing light intensity.

Equation (19) can be fit to the data on *S. linearis* also using the Metropolis-Hastings procedure (see Methods). The values of parameters *m*, ${v}_{I}^{max},$ and *A*_{I} (and thus *A*_{Chl} and *B*_{Chl}) were solved by fitting Equation (17) above so the N:C vs. μ data provide constraints on *A*_{Pho} and *A*_{Bio} (the parameters which scale photosynthetic protein to chlorophyll and biosynthetic protein to growth rate) as well as the fixed pool of “essential” protein, ${Q}_{C}^{Pro-Other}$. In Figure 5A we illustrate fitted solutions of the of the un-approximated version of Equation (19) (where the minor contribution to N:C from RNA and chlorophyll are also resolved). Similar simulations of *Pavlova lutheri* and *Skeletonema costatum* are shown in Supplementary Figure 2. The ability of the model to fit the data suggests that the model captures major processes.

The model result shows a similar relative increase in the investment in photosynthetic and biosynthetic protein at moderate to high light levels (Figure 5B, Supplementary Figures 2–4). At lower photon fluxes N:C increases more rapidly (Figure 5A) because there is a much higher demand for investment in photosynthetic machinery to achieve the same growth rate (Supplementary Figures 2–4). The model suggests that N:C increases linearly with growth rate at fixed light because investment in both biosynthetic and photosynthetic protein must increase linearly and protein investment dominates the N:C ratio. The un-approximated model (Equation 33 in Methods) suggests a non-linearity due to investment in RNA, but its overall contribution to the cellular nitrogen budget means that the non-linearity is very small and Equation (19) is a good approximation.

#### Why Does P:C Increase Non-linearly With Growth Rate?

A significant fraction of cellular phosphorus is present in nucleic acids, lipid membranes, and storage compounds such as polyphosphate. Consider the case for P-limited culture in which luxury storage is small and the cellular quota of phosphorus is approximated by the sum of the three pools:

where ${Q}_{P}^{Other}$ groups relatively stable pools of P: ${Q}_{P}^{Other}={Q}_{P}^{DNA}+{Q}_{P}^{Other0}$, required even in the absence of net growth.

Here we invoke two of the fundamental relationships discussed earlier: the investment in thylakoid phospholipid, ${Q}_{P}^{Thy}$, is assumed linearly proportional to chlorophyll (Equation 13) and the investment in RNA is modeled as linearly proportional to total protein content and growth rate (Equation 15). (These models, and the evidence for them, are discussed in the section entitled “*Representing relationships between macromolecular pools and rates*”). Substituting Equations (13, 15, 3, 12, 14, 17) (in this order) into Equation (20), yielding a quadratic relationship between cellular P:C with growth rate:

where

Equation (21) predicts a quadratic relationship of P:C with growth rate, μ, is qualitatively consistent with the non-linear relationship in the P-limited cultures of *S. linearis* (Figure 1H), marine *Synechococcus* (WH8102) and *Selenastrum minutum* (Figure 1H). The qualitative fit enabled the optimization of parameters to match the *S. linearis* data (using the Metropolis Hastings algorithm, see Methods) and the resulting fit is shown in Figure 6A. Using the optimized parameters, the model provides a prediction of the allocation of phosphorus, shown as a function of growth rate at a single light intensity in Figure 6B. ${Q}_{P}^{Thy}$ increases linearly with growth rate for a fixed photon flux (Figure 6B) in concert with ${Q}_{C}^{Chl}$ and the photosynthetic proteins (Figures 4, 5A). The inference is that the non-linear relationship of P:C with growth rate (Figure 6B) is due to the investment in phosphorus-rich RNA (Figure 6B) which increases in proportion to both the growth rate and the cellular quota of protein, which also increases with growth rate due to investment in biosynthesis and light harvesting (Figure 5B). A similar non-linear relationship of P:C vs. μ was shown to be consistent in the culture of *Selenastrum minutum* (Elrifi and Turpin, 1985; Ågren, 2004).

**Figure 6**. Model-data comparison of P:C and model prediction of macromolecular allocation in P for *Synechococcus linearis*. **(A)** P:C under P limitation under different light intensities; light intensities (μmol m^{−2} s^{−1}) are in the legend. Curves are model results and points are data (Healey, 1985). Dotted line represents P:C at ${\mu}_{max}^{I}$ for various light intensities; high ${\mu}_{max}^{I}$ for higher light intensity. **(B)** Macromolecular allocation in P, normalized by cellular C under the light intensity of 62 (μmol m^{−2} s^{−1}). Black points are data for total values under the same light intensity (Healey, 1985). See the legend for color definitions: *P thyla*, P in thylakoid membranes; *Other*: other molecules. Detailed macromolecule allocation in Supplementary Figure 3 and P:C under N limitation in Supplementary Figure 4.

#### What Is the Maximum Growth Rate for a Given Light Intensity, ${\mu}_{max}^{I}$?

Expanding the cellular carbon quota in terms of the macromolecular components as described by Equations (2)–(6), and accounting for only the most quantitatively influential molecules in order to provide an analytic solution (the full, un-approximated model is described in Methods), we describe carbon allocation in the cell as

Using Equation (3), we can further resolve the proteomic contributions into photosynthetic, biosynthetic, structural/other:

where ${Q}_{C}^{Other}\approx {Q}_{C}^{Pro-Other}+{Q}_{C}^{Carb-Other}+{Q}_{C}^{Lip-Other}$ and ${Q}_{C}^{Csto}=\text{}{Q}_{C}^{Carb-Sto}+\text{}{Q}_{C}^{Lip-Sto}$. Using Equations (12, 14, 17), we can solve Equation (23) for the population growth rate, μ;

This equation indicates that as the investment in carbon storage ${Q}_{C}^{Sto}$ decreases, growth rate μ increases.

This inference is logical: the maximum growth rate for a given light intensity, ${\mu}_{max}^{I}$, should occur when as much biomass can be allocated to growth related macromolecules as possible; in other words when carbon storage is minimal and ${Q}_{C}^{Csto}$ approaches 0. In Figure 7 we illustrate this in terms of carbon allocation as a function of growth rate in model solutions where the parameters were fitted for *S. linearis*. Solutions are shown for two light levels: at low light, the rapidly increasing allocation to photosynthetic machinery as a function of growth rate means that cellular allocation to storage becomes small (and allocation to functional machinery becomes large) at quite a low growth rate. In contrast, at high light, the lower demand for photosynthetic apparatus allows a greater investment in biosynthesis and higher maximum growth rate (Figure 7, Supplementary Figures 5A, 6).

**Figure 7**. Simulated carbon allocation and ${\mu}_{max}^{I}$ of *Synechococcus linearis* under N limitation. The light intensities are **(A)** 12 and **(B)** 62 (μmol m^{−2} s^{−1}). See **(A)** for color definitions: *C store*, C storage; *Bio*, Biosynthetic protein + RNA; *Photo*, Chlorophyll + Photosynthetic protein + P-lipid in thylakoid membranes; *Other*, other molecules with constant cellular investment. *Bio* and *Photo* are dominated by biosynthetic and photosynthetic proteins, respectively. C allocation and ${\mu}_{max}^{I}$ of *Synechococcus linearis* under P limitation in Supplementary Figure 5. C allocations of other species (*Pavlova lutheri* and *Skeletonema costatum*) in Supplementary Figure 6.

The limit of the model, which occurs when ${Q}_{C}^{Csto}$ approaches zero in Equation (24), reproduces the observed ${\mu}_{max}^{I}$–light relation of *S. linearis* and two marine phytoplankton (*Pavlova lutheri* and *Skeletonema costatum*; Figure 8, Supplementary Figure 5B). The ${\mu}_{max}^{I}$ curve ultimately saturates because photosynthesis per chlorophyll *v*_{I} saturates. When ${\mu}_{max}^{I}$ increases, due to an increase in light, there is a decreased investment in light-harvesting proteins which is traded off against increased biosynthetic proteins and a higher maximum growth rate. This is qualitatively consistent with recent proteomic studies (McKew et al., 2013; Jahn et al., 2018; Zavřel et al., 2019).

**Figure 8**. Simulated light dependence of maximum growth rate (${\mu}_{max}^{I}$) and photosynthesis rate (*v*_{I}) under N limitation for three different phytoplankton. **(A)** *Synechococcus linearis*, **(B)** *Pavlova lutheri*, and **(C)** *Skeletonema costatum*. Modeled ${\mu}_{max}^{I}$ is compared to data (Healey, 1985; Sakshaug and Andersen, 1989; Chalup and Laws, 1990). There is one data point at the light intensity of 1,203 (μmol m^{−2} s^{−1}) for *Skeletonema costatum* not included or considered since the cells were likely photodamaged under such high light intensity, negatively altering ${\mu}_{max}^{I}$. Light dependence of ${\mu}_{max}^{I}$ of *Synechococcus linearis* under P limitation in Supplementary Figure 5, where ${\mu}_{max}^{I}$ is almost identical as **(A)**, similarly reproducing the data (Healey, 1985).

#### Differences Between ${\mu}_{max}^{I}$ and *v*_{I}

Growth rates and photosynthesis are often used interchangeably and the relationships for photosynthesis and light have often been applied to growth rates in ecosystem models (Moore et al., 2004; Buitenhuis et al., 2013; Dutkiewicz et al., 2015; Coles et al., 2017). However, photosynthesis and growth (biosynthesis) are metabolically distinct and need not be equivalent. The two rates are qualitatively similar because both share a saturating dependence on light (Figure 8)_{.} Growth rate has a non-zero intercept on the light axis, which represents the minimum light intensity required for cellular maintenance. Notably, ${\mu}_{max}^{I}$ approaches the saturated value at a lower light intensity than *v*_{I}. This can be seen in Figure 8 where the model is fit to data sets for *S. linearis* (Healey, 1985) and two photosynthetic algae (Sakshaug and Andersen, 1989; Chalup and Laws, 1990). As photon flux decreases, investment in photosynthetic apparatus increases rapidly (Figure 7, Supplementary Figures 5A, 6) at the expense of biosynthetic machinery and such that ${\mu}_{max}^{I}$ saturates at a lower light intensity than *v*_{I}. This highlights the high cost of the photosynthetic apparatus.

#### What Does N:P Depend Upon?

The N:P ratio of plankton has been a topic of interest going back to Redfield (1934, 1958). Could the framework presented above be used to predict and interpret the N:P of primary producers? The answer is not immediately clear: Equations (19, 21) represent the N:C ratio under N limitation and P:C under P limitation. Simply dividing the two does not provide an accurate prediction of N:P because, typically, one resource is limiting and the other in excess in the environment, and non-limiting resources accumulate in intracellular storage pools.

Using the comprehensive data set of Healey (1985), it is possible to quantify the storage of N and P when each is the non-limiting resource. For stoichiometric purposes, we assume that P-storage is in the form of polyphosphate which has no carbon content and that N-storage is in the form of cyanophycin (Table 1). The storage capacity of the cells cannot be predicted a priori so the model allows cells to take up and store the non-limiting resource (N or P) until a maximum storage is reached, constrained by the observed elemental stoichiometry. Thus, we introduce two new, empirically constrained parameters; the maximum storage capacities for N and P (see the subsection “Evaluating cellular C concentration and N and P storage” in Methods for details). Using this approach, we can model N:P under both N and P limitation (illustrated in Figure 9), as well as N:C under P limitation and P:C under N limitation (Supplementary Figure 4). In Figure 9, model parameter values are constrained (as discussed above) with data for *S. linearis*, including the maximum storage capacities.

**Figure 9**. Model-data comparison of N:P. **(A)** Under N limitation. **(B)** Under P limitation. Curves: model. Points: data of *Synechococcus linearis* (Healey, 1985). Dotted black line indicates values at ${\mu}_{max}^{I}$ for various light intensities. Numbers in the legend show light intensities (μmol m^{−2} s^{−1}).

The model qualitatively captures the variations in N:P with growth rate under both N and P limitation. Under N-limitation, N:P (Figure 9A) has a linearly increasing trend with growth rate, following N:C (Figure 5) because of the increasing investment in N-rich proteins with growth rate. In this case P:C is relatively constant regulated by the contribution of P storage (Supplementary Figure 4B). Under P-limitation, N:P declines rapidly with growth rate (Figure 9B) because the increase in RNA with growth rate is quadratic (Figure 6), while that of protein is linear (Supplementary Figure 4A), so the investment in P increases more rapidly. N storage is a relatively moderate contribution to the cellular N quota (Supplementary Figure 4A). The model with the N and P storages captures and interprets the trends of N:P observed in the laboratory study by simply allowing storage up to an empirically informed limit. Further basic study of the dynamics of, and limits to, nutrient storage pools would be necessary to inform a model rooted in first principles.

#### Plausibility of the Predicted Macromolecular Allocation

We have used laboratory data on elemental stoichiometry to constrain a model which resolves macromolecular allocation. As such, the model makes testable predictions. Our estimated photosynthetic parameters (${v}_{I}^{max}$ and *A*_{I}) sit within the range of observation (Platt et al., 1980; Cullen, 1990; Moore and Chisholm, 1999) for the given range of chlorophyll vs. μ (Figure 4, Supplementary Figure 2). While there are not direct macromolecular or proteomic data available for the particular laboratory studies which we simulated, some recent culture studies have resolved macromolecular and proteomic allocation. The predictions of changes in macromolecular allocation (e.g., Proteins, Lipids+ Carbohydrate and RNA) with growth rate and light are qualitatively similar to those observed in a recent laboratory studies (Liefer et al., 2019). Also, in Figures 10A,B, we compare the inferred allocation to protein, carbohydrates and lipids as a function of growth rate for *S. linearis* with measurements of these bulk macromolecular pools in chemostat cultures of *Prochlorococcus marinus* (PCC 9511) (Felcmanová et al., 2017). The model qualitatively captures the observed trends in allocation and the general magnitude of the observed pools, though the specific values differ due to either inter-species differences or model limitations. In Figures 10C,D, we also compare the predictions for broad-scale protein allocation with light-dependent growth rate (i.e., maximum growth rate for a given light intensity) from the constrained *Synechococcus* model, and direct observations of similar proteomic categories from turbidostat cultures of *Synechocystis* (Jahn et al., 2018). Again, there is a qualitative agreement between the predicted trends: as the light-limited growth rate increases, the investment in light-harvesting proteins declines while the investment in biosynthesis increases. The inferred allocation to photosynthetic proteins tends to be rather high relative to the direct proteomic study, likely reflecting inter-species differences. Direct laboratory studies where elemental stoichiometry, macromolecular composition, and proteomics are all concurrently measured are possible and would allow a more strenuous test and calibration of such models.

**Figure 10**. Model-data comparison of macromolecular allocation. **(A)** Modeled C allocation of *Synechococcus linearis* at the light intensity of 50 μmol m^{−2} s^{−1} under N limitation. *Carb*.: carbohydrate. **(B)** Measured C allocation of *Prochlorococcus marinus* (PCC 9511) at the light intensity of ~50 μmol m^{−2} s^{−1} under N limitation (Felcmanová et al., 2017). C allocation (%) for protein and carbohydrate (*Carb*.) are on the left axis and that for lipid is on the right axis. Dotted lines are linear interpolation and error bars are the standard deviation. **(C)** Modeled protein allocation of *Synechococcus linearis* at ${\mu}_{max}^{I}$, nutrient replete growth rate. Here changes in ${\mu}_{max}^{I}$ are caused by changes in light intensities. *Photo*, photosynthetic proteins; *Bio*, biosynthetic proteins; *Other*, other proteins. **(D)** Observed protein allocation in *Synechocystis* sp. (PCC 6803) (Jahn et al., 2018). *LHC*, light-harvesting complex; *RIB*, ribosome and protein production; *CBM*, proteins for C uptake, fixation and metabolism. Compare the general trends in *Photo* with *LHC* and *Bio* + *Other* with *RIB* and *CBM*.

## Discussion

### A Model of the Elemental Stoichiometry of Phytoplankton in Relation to Light Intensity, Growth Rate, and Limiting Resource

The model presented above provides a conceptually simple, yet quantitative description of the relationship between the elemental stoichiometry of phytoplankton, their growth rate, resource availability, and macromolecular allocation (Figure 11). It is based on a straightforward accounting of the allocation between and within major pools of macromolecules, along with four representations of the relationships between pools and fluxes. These are (i) a saturating relationship between light intensity and photosynthetic efficiency, (ii) a constant ratio of chlorophyll to other light harvesting and photosynthesis apparatus, (iii) a linear relationship between allocation of biosynthetic protein and growth rate, and (iv) a linear relationship between RNA:protein and growth rate. Representations (i), (iii), and (iv) are empirically driven, while (ii) is hypothetical, though simple and logical. The ability of the model to fit laboratory data for diverse phytoplankton taxa indicates that the model framework is generally applicable.

**Figure 11**. How macromolecular allocation of phytoplankton responds to light and nutrient, constraining the growth rate based on CFM-Phyto. *Bio*, biosynthetic molecules; *Photo*, photosynthetic molecules; *Other*, other constant molecules; *Store*, storage molecules. Red dashed arrows represent light intensities. Here the growth rate μ is proportional to the ratio of *Bio*. When light and nutrient are sufficient, there is a high ratio of *Bio* with high ${\mu}_{max}^{I}$ (the maximum viable growth rate at a certain light intensity) **(A)**. Light limitation increases the fraction of *Photo* limiting that of *Bio* leading to lower ${\mu}_{max}^{I}$ **(B)**. Nutrient limitation leads to a smaller allocation to photosynthetic and growth-related molecules (*Photo, Bio*) leading to a growth rate μ lower than ${\mu}_{max}^{I}$ **(C)**. Under both light and nutrient limitation, these effects are combined, further limiting the ratio of *Bio*, thus limiting the growth rate **(D)**. C storage lowers N:C of phytoplankton.

The framework of the model is conceptually simple and steady-state solutions can be solved algebraically and parameters optimized to match empirical data. We have used it to model and interpret laboratory data relating the elemental stoichiometry of *S. linearis* to growth rate and light intensity under both N- and P-limiting conditions. The stoichiometric data provide indirect constraint on macromolecular allocation, which imply a common set of allocation strategies amongst phytoplankton. Below, we discuss some of the limitations and simplifications of the approach.

By explicitly resolving the macromolecular allocation, the model captures and provides a simple interpretation for the contrast between cellular nitrogen quota (or N:C ratio) which varies linearly with growth rate under N-limitation, and the phosphorus quota (or P:C) which varies non-linearly with growth rate under P-limitation. At fixed light intensity and under nitrogen limitation, cellular protein increases linearly with growth rate (Felcmanová et al., 2017), driving the linear trend in N:C vs. μ captured in Equation (19). In contrasts, the observed linear relationship between RNA:protein (Nicklisch and Steinberg, 2009) combines with the increase in protein with growth rate (Rhee et al., 1981; Liefer et al., 2019) to drive the non-linear relationship of P:C with growth rate, captured in Equation (21). This contrast was also captured in the more abstract model of Ågren (2004). The explicit macromolecular resolution also allows the model to capture the observed decrease of N:P with growth rate under P limitation (Perry, 1976; Elrifi and Turpin, 1985; Healey, 1985; Garcia et al., 2016; Figure 9B) which is driven by the increasing RNA/protein ratio with μ (Bremer and Dennis, 1996; Nicklisch and Steinberg, 2009; Scott et al., 2010).

We suggest that the explicit resolution of macromolecular reservoirs provides an important advantage over more idealized frameworks. It allows the exploitation of key observed relationships (e.g., RNA:protein vs. μ) and it explicitly couples the dynamics of N, P and C providing a comprehensive framework (c.f. N:C only in the case of Geider et al., 1998, for example). The more idealized internal-stores model (Droop, 1968) treats all resources alike and implicitly does not capture this contrast in N:C and P:C variations with μ, unless significantly different parameters are applied between N:C and P:C. The chain model (Pahlow and Oschlies, 2009) differentiates the dynamics of N and P, but still remains abstracted from the macromolecular foundations. It predicts non-linear relationships for both N:C and P:C with μ, which also leads to compensation and predictions of a rather constant N:P under P limitation (Pahlow and Oschlies, 2009), whereas as the data generally shows decreasing N:P with μ (Perry, 1976; Elrifi and Turpin, 1985; Healey, 1985; Garcia et al., 2016) as our model predicts (Figure 9B).

We suggest that the explicit macromolecular representation also has some interpretive and predictive advantages. It has the potential to be directly compared with direct proteomic and macromolecular observations (McKew et al., 2013; Felcmanová et al., 2017; Jahn et al., 2018; Liefer et al., 2019; Zavřel et al., 2019), which will leverage new data sets and technologies more directly. For example, the calibrated model for *Synechoccocus linearis* presented above indicates a strong relationship between investment in photosystem proteins and light intensity. Under fixed light intensity, both biosynthetic and photosynthetic proteins increase with growth rate. However, under light co-limitation, when varying the growth rate (i.e., ${\mu}_{max}^{I}$), by changing light intensity, the model predicts a reduced requirement for photo-proteins with increasing light accompanied by an increase in biosynthetic protein, consistent with data from proteomic studies (McKew et al., 2013; Jahn et al., 2018; Zavřel et al., 2019). This mechanism stabilizes the amount of total protein and explains relatively stable N:C with simultaneously varied light and growth rate (Geider et al., 1985).

The model illustrates the relationship between the maximum growth rate at a given light intensity and storage. In order to increase the growth rate, cells invest in protein at the expense of storage compounds. The maximum growth rate for a given light intensity occurs when storage is minimized and functional allocation is maximized. In some circumstances, maximizing growth rate will be the best measure of fitness, but in others storage is likely to be advantageous. For example, if we consider phytoplankton-bloom conditions, maximizing growth rate may be more important since in such situations, phytoplankton with faster growth can outcompete others and dominate the region (Dutkiewicz et al., 2009). However, in environments where the nutrient level rapidly fluctuates (e.g., with time scale of days), phytoplankton with high storage capacity might be advantageous by being able to grow under nutrient depletion with stored nutrients before another pulse of high nutrient occurs (Tozzi et al., 2004; Grover, 2009, 2011).

### Generality of the Model

While we have focused our development and discussion around the data set of Healey (1985) for *S. linearis*, the framework is sufficiently coarse-grained and rooted in basic, common physiology that it is qualitatively compatible with data from numerous phytoplankton, spanning a wide range of cell size and taxonomic groups (Figure 1). We have fit the same framework to several of these data sets (e.g., Figures 4–6, 8, Supplementary Figure 2; Healey, 1985; Sakshaug and Andersen, 1989; Chalup and Laws, 1990). While different parameter values are required, reflecting different allocation strategies or traits, the basic framework is general, predicting the common trends from the laboratory studies (Figures 1, 2, Supplementary Figure 1). However, the allocation strategies (and parameter values) differ between species. The use of different nitrogen substrates (e.g., nitrate vs. ammonium) could be represented by changes to the respiratory cost of synthesis, *E* (Rittmann and McCarty, 2001).

### Model Simplifications and Limitations

The model presented here represents an attempt to provide a minimal, transparent and biologically meaningful framework which relates allocation between and within the major macromolecular pools to elemental ratios and growth rates under diverse environmental conditions. It is framed so that the internal allocation is, in principle, in terms of measurable quantities (though not all were available in the data sets studied here). These measurable pools can be mapped into categories which are not directly measurable, but which are grouped by function (Figures 3B, 11), in the spirit of allocation models (Shuter, 1979; Scott et al., 2010). This functional mapping enables simple interpretations of the relationships of interest.

As with any quantitative model, there is a trade-off between realism, data constraints and insight, and we have not resolved a number of potentially important factors. We have assumed a fixed composition of thylakoid membranes which may vary in reality. In particular, the fraction of light harvesting machinery might change relative to other components, which would alter the chlorophyll to protein or lipid ratio. The model could be improved with an additional layer of detail, separating the light harvesting and other components. To constrain the model, combined measurements of chlorophyll and proteomics (e.g., McKew et al., 2013, 2015; Zavřel et al., 2019) as well as thylakoid lipid would be useful. Photo-inhibition has not been addressed here and presumably would demand the resolution of photoprotective proteins (Geider et al., 2009). Healey's (1985) data set focused on lower photon fluxes so this case was not addressed here. We have also constrained a single set of photosynthesis-vs.-light parameters to simulate experiments at all light levels, though we are aware that acclimation would likely modify them, however we found no significant improvement in model fits when allowing this extra degree of freedom. We have not addressed the potential for variations in allocation and elemental ratios as a function of temperature, though it is known that there are indeed sensitivities (Thrane et al., 2017) and this could be an interesting and important extension to the study. We have also not explicitly resolved allocation to nutrient uptake transporters, which varies with nutrient concentration (McKew et al., 2015; Lin et al., 2016). However, in overall elemental stoichiometry, the influence appears limited since proteomic studies suggest that allocation to the nutrient transporter proteins is modest (e.g., on the order of 10^{−1}% of total spectral counts in a recent proteomic study; McKew et al., 2015) relative to investments in light harvesting and biosynthesis (McKew et al., 2013; Jahn et al., 2018; Zavřel et al., 2019). Also, investment in transporters increases under low nutrients (McKew et al., 2015), so if the transporter is a dominant part of protein, it would lead to high N:C at low growth rate, but the data show otherwise (Figure 1, Supplementary Figure 1). Hence, we suggest that investment in transporters is a next-order effect which would have a small impact for elemental stoichiometry in the phytoplankton addressed here. Resolution of transporter allocation would also introduce further unconstrained parameters. Other studies have placed more emphasis on the allocation to transporters showing it to be an important factor in expressing phenotypic diversity and acclimation to different nutrient and light regimes (Smith et al., 2009, 2016; Bonachela et al., 2013; Garcia et al., 2016; Chen et al., 2019).

We have also neglected the substitution of non-P-lipids for P-lipids under low P concentrations (Van Mooy et al., 2009), which has a significant impact on the P budget of the cell, and focused on the major macro-nutrient elements (C, N, P) though trace metal allocation is also of significance (Ho et al., 2003; Saito et al., 2011). The allocation of trace metals to specific protein groups would provide a way to link them in such a model. As discussed above, storage of non-limiting elements is important for consideration of the N:P ratio and the applicability of the Growth Rate Hypothesis (specifically under N-limitation). The limits to storage and maximum quotas are not clearly defined at present and worthy of further work. Such simplifications and omissions could ultimately be addressed with coordinated laboratory and modeling studies.

### Perspective and Outlook

Despite the limitations of the study listed above, we have shown that a conceptually simple model rooted in mass balance and a few basic, empirically sound representations can capture the relationships between growth rate and elemental stoichiometry under a variety of environmental conditions accurately. We suggest that the explicit representation of measurable macromolecular pools allows an advantage over more abstracted forms rooted in elemental quotas. It allows the exploitation of key physiological observations such as the changes in RNA:protein with μ, as well as testable predictions regarding macromolecular allocation. Parameters controlling rates and allocation can be calibrated with laboratory data, either inverted from stoichiometric data as we have done here, or directly measured (McKew et al., 2013, 2015; Felcmanová et al., 2017; Jahn et al., 2018; Zavřel et al., 2019) though this is not yet routinely the case.

Physiological models of “intermediate complexity” such as this have a role to play in ecological and biogeochemical studies. While Monod (1949) and Droop (1968) kinetics provide much simpler frameworks which have fewer parameters and are mathematically convenient, they lack some important biological detail, especially if one wishes to relate elemental stoichiometry to growth rate and environment. The approach presented here, while still idealized, is economical and, with some modifications, could be efficiently employed in biogeochemical and ecological simulations. While modern Flux Balance Analysis approaches now allow genome-scale representations of microbial physiology (Orth et al., 2010) they are typically subject to an imposed macromolecular composition (the “biomass function”) which is generally empirically determined and invariant and so do not address the elemental stoichiometry of the cell prognostically. Laboratory studies reveal visible changes in biomass function (i.e., the relative allocation to different macromolecules) over reasonable ranges of environmental conditions (Rhee, 1978; McKew et al., 2013, 2015; Felcmanová et al., 2017; Jahn et al., 2018; Liefer et al., 2019; Zavřel et al., 2019). Thus, models of the type presented here complement, and could potentially couple to, more detailed genome-scale simulations. We suggest that integrated laboratory and modeling studies in which a comprehensive set of physiological measurements (i.e., elemental stoichiometry, proteome, transcriptome) and a hierarchy of models (coarse-grained and genome-scale) would be valuable.

## Methods

### Full Model Description and Parameter Estimation

Here we provide a complete version of the model: CFM-Phyto (Figure 3). We first detail the organization of macromolecular components into four functional classes (*Photo, Bio, Store*, and *Other*). Then we discuss how maximum growth rate, ${\mu}_{max}^{I}$, can be predicted from macromolecular allocation. Then we provide details on how storage, population density, elemental stoichiometry and carbon biomass density are evaluated. Finally, we describe how model parameters are estimated.

### Re-framing the Model According to Functional Allocation

Equations (2)–(6) lead to the following accounting for total cellular carbon in various macromolecular pools:

where ${Q}_{C}^{Other0}={Q}_{C}^{Carb-Other}+{Q}_{C}^{Lip-Other}$. The cellular pools, defined in the main text, are described in carbon units, relative to total cellular carbon (mol C mol^{−1} C). The cellular components can be re-arranged and gathered into four functional classes, as depicted in Figure 3B:

where

${Q}_{C}^{Pho}$ includes all components in the thylakoid membranes, which are allocated according to the light intensity and growth rate. The model assumes that all components of ${Q}_{C}^{Pho}$ are adjusted in concert; i.e., the relative proportions of the components of ${Q}_{C}^{Pho}$ are fixed and independent of variations in the magnitude of ${Q}_{C}^{Pho}$. In other words, the makeup of the photosynthetic machinery is invariant. ${Q}_{C}^{Bio}$ contains all components whose allocation depends mostly on growth rate. ${Q}_{C}^{Other}$ includes components which are assumed to represent fixed fractions of the cell: both ${Q}_{C}^{Other}$ and its components are independent of light intensity and growth rate (Figure 7, Supplementary Figures 5A, 6).

### Evaluating Allocation to Photosynthetic Apparatus

We first obtain ${Q}_{C}^{Chl}$ with Equation (17), and then ${Q}_{C}^{Pro-Pho}$ and ${Q}_{P}^{Thy}$ from Equations (12) and (13), respectively. ${Q}_{P}^{Thy}$is stoichiometrically related to ${Q}_{C}^{Plip-Thy}$:

where ${Y}_{Plip}^{C:P}$ is C:P of phospholipids.

### Evaluating Allocation to Biosynthetic Apparatus

To compute biosynthetic apparatus, we first compute ${Q}_{C}^{Pro-Bio}$ and ${Q}_{P}^{RNA}$ from Equations (14) and (15), respectively. ${Q}_{P}^{RNA}$ is stoichiometrically related to ${Q}_{C}^{RNA}$:

where ${Y}_{RNA}^{C:P}$ is the C:P of RNA.

### Evaluating Maximum Growth Rate, ${\mu}_{max}^{I}$

${\mu}_{max}^{I}$ is the value of growth rate, μ, when all of the flexible component of cellular carbon has been allocated to the growth-related apparatus (${Q}_{C}^{Pho}$ and ${Q}_{C}^{Bio}$) and allocation to C or N storage is negligible. Phosphorus storage is assumed to be polyphosphate, thus not contributing to the carbon budget. Given ${Q}_{C}^{Sto}$ = 0, and substituting Equations (27, 13, 28, 15, 3, 14, 12, 17) (in this order) into Equation (25) leads to the following quadratic equation in μ:

where

Here, the positive solution for μ*-* equals ${\mu}_{max}^{I}$:

### Obtaining N:C

N:C is represented by the sum of N from N-containing molecules normalized by cellular C quota:

Here, we define ${Y}_{Chl}^{N:C}$, ${Y}_{Pro}^{N:C}$, ${Y}_{DNA,}^{N:C}$ and ${Y}_{RNA}^{N:P}$ as N:C of chlorophyll, protein and DNA and N:P of RNA, respectively (value in Table 1). Using these conversion terms:

Then, by substituting Equations (15, 3, 14, 12, 17) (in this order) into Equation (32), we obtain

where

When the nitrogen content of RNA is accounted for, we predict a quadratic relationship between N:C and growth rate. However, since the contribution from RNA is small, that from protein dominates and the linear approximation of Equation (19) works well (as seen in the data and un-approximated solution shown in Figure 5A). We define *Q*_{N} as N:C and ${Q}_{N}^{NonSto}$ as N:C without N storage:

### Obtaining P:C

P:C is represented by the sum of N from N-containing molecules normalized by cellular C quota:

We define ${Y}_{DNA}^{P:C}$ as P:C of DNA (value in Table 1), which leads to

By substituting Equations (15, 3, 14, 13, 12, 17) (in this order) into Equation (36), we obtain

where

We define *Q*_{P} as P:C and ${Q}_{P}^{NonSto}$ as P:C without P storage:

Since RNA is the dominant contribution to cellular phosphorus, the relationship between P:C and growth rate is non-linear (Figures 1H, 6A).

### Obtaining N:P

Once we obtain N:C and P:C, N:P can be obtained as follows:

### Evaluating Cellular C Concentration and N and P Storage

There are three types of storage: C in carbohydrates and lipids, N storage assumed to be cyanophycin, and P storage assumed to be polyphosphate. Only C and N storage affect the carbon budget. To compute N and P storage, which we assume accumulate only when each element *is not* limiting, we must first determine which nutrient *is* limiting. To do that, we first compute the carbon-based biomass in the culture under N or P limitation, [*C*_{Cell}]_{i} (mol C m^{−3}), where *i* is *N* or *P*, respectively. Carbon biomass, [*C*_{Cell}]_{i}, is by definition the product of the cellular carbon quota, *C*_{Cell} (mol C cell^{−1}), and the cell density, *X*_{i} (cell m^{−3}):

Under N limitation, the time variation of dissolved inorganic N (or NO${}_{{3}^{-}}$) [*N*] (mol N m^{−3}) in the culture is based on the balance between dilution and uptake:

where *D* is the dilution rate (d^{−1}), [*N*]_{in} (mol N m^{−3}) is the concentration of dissolved inorganic N (or NO${}_{3}^{-}$) in the incoming medium, *V*_{N} is the N uptake rate per cellular C (mol N mol C^{−1} d^{−1}), *X*_{N} (cell m^{−3}) is the cell density in the culture under N limitation. We also consider the time variation of *X*_{N}:

At steady state (i.e., d[N]/dt = 0), Equation (41) suggests that

where [*C*_{Cell}]_{N} (mol C m^{−3}) (= *C*_{Cell}*X*_{N}) is the carbon biomass in the culture under N limitation. The steady state of Equation (42) leads to the following well-known relation for a chemostat at steady state:

To relate *V*_{N} to known parameters, we further consider the balance of *Q*_{N} :

The steady state of this equation leads to a simple relation between the uptake and consumption of N:

as *Q*_{N} = ${Q}_{N}^{NonSto}$ under N limitation (given N storage would be small). By assuming that the amount of limiting nutrient is small relative to that in the incoming medium (here [*N*]_{in} >> [*N*]) as in previous chemostat simulations (Inomura et al., 2017, 2018) and as justified by laboratory observations (Laws and Bannister, 1980; Healey, 1985; Bühler et al., 1987), and by substituting Equations (44, 46) into Equation (43), we obtain the simple expression for [*C*_{Cell}]_{N}:

We follow the same procedures above (Equations 41–47) by replacing N with P to obtain an expression for the carbon biomass in the culture under P limitation [_{CCell]P} (mol C m^{−3}):

Here [*P*]_{in} (mol m^{−3}) is the concentration of dissolved inorganic P (PO${}_{4}^{3-}$) in the incoming medium.

We assume that the limiting resource is that which gives the smallest cellular C concentration in the culture. For example, when [*C*_{Cell}]_{N} < [*C*_{Cell}]_{P}, the culture is limited by N since this relationship with Equations (47) and (48) leads to the following equation:

showing that the input N:P (left hand side) is lower than required N:P (right hand side). On the other hand, when [*C*_{Cell}]_{N} > [*C*_{Cell}]_{P}, input N:P is higher than required N:P indicating excess N and thus P limitation. Once the nutrient limitation is determined, we define the actual cellular C concentration in the culture [*C*_{Cell}] (cell^{−1} m^{−3}):

With this equation, we have simulated the relationship between biomass, [*C*_{Cell}], and growth rate, μ, which reveals a decreasing trend with dilution rates, capturing the observation (Healey, 1985; Supplementary Figure 7).

If the culture is N limited, there is an excess of P, which could be stored in the cell. To determine the potential level of cellular P based on the P availability (${Q}_{P}^{Pot}$) (mol P mol C^{−1}), we follow the same procedure as when determining [*C*_{Cell}]_{N} from Equations (41)–(47), by using P instead of N, except for using *X*_{N} and [*C*_{Cell}]_{N} obtained previously:

Then we compare this potential quota with the maximum capacity of cellular P (${Q}_{P}^{max}$) (mol P mol C^{−1}) (mol P mol C^{−1}) and define *Q*_{P}:

and by rearranging Equation (38) we determine ${Q}_{P}^{Sto}$:

Under P limitation, N storage is accumulated since excess N is available. By following the same steps as above with reversed N and P and by rearranging Equation (34), we obtain the following relations:

where ${Q}_{N}^{Pot}$ (mol N mol C^{−1}) is the potential cellular quota of N normalized by C based on the N availability and the cellular C determined by P limitation and ${Q}_{N}^{max}$ (mol N mol C^{−1}) is the maximum cellular capacity of N normalized by cellular C. We note that in most cases, ${Q}_{N}^{Pot}$ > ${Q}_{N}^{max}$ (or ${Q}_{P}^{Pot}$ > ${Q}_{P}^{max}$). However, when the N:P ratio without storage molecules (i.e., ${Q}_{N}^{NonSto}$/${Q}_{P}^{NonSto}$) is close to the N:P ratio of the resource (i.e., [*N*]_{in}/[*P*]_{in}), the excess nutrient is relatively small and ${Q}_{N}^{Pot}$ < ${Q}_{N}^{max}$ (or ${Q}_{P}^{Pot}$ < ${Q}_{P}^{max}$) can occur (e.g., at the light intensity of 12 μmol m^{−2} s^{−1} in Supplementary Figure 4A).

The data and model together reveal that N and P “storage” work differently (Supplementary Figure 4) in *S. linearis* (Healey, 1985). Under nitrogen limitation, the total phosphorus quota per carbon appears relatively constant while under P-limitation, N storage appears to be relatively constant. Hence, to model the storage contributions, we have imposed maximum total P quota per carbon (${Q}_{P}^{max}$) (mol P mol C^{−1}) and maximum N storage per carbon (${Q}_{N}^{Sto,max}$) (mol N mol C^{−1}); thus ${Q}_{N}^{max}={Q}_{N}^{Sto,max}+{Q}_{N}^{NonSto}$. In other words, maximum P storage depends on the level of other P molecules while maximum N storage is independent from the level other N molecules. These simple assumptions allow us reproduce N:C under P limitation and P:C under N limitation (Supplementary Figure 4) as well as N:P (Figure 9).

This model of the storage pools is simple and logical, yet still somewhat *ad hoc* and empirically driven. Because of the importance of storage of the non-limiting element, prediction of the N:P ratio depends on this. Clearer understanding of the dynamics of the storage pools will be necessary to provide more mechanistic models. Laboratory data which resolves the macromolecular pools in sufficient detail would aid this effort. Interpretations of the Redfieldian N:P ratio as a homeostatic protein:RNA ratio (Loladze and Elser, 2011), while revealing the central controls, do not necessarily reflect these important “storage” dynamics (Supplementary Figure 4).

N storage has an associated C (e.g., cyanophycin), which can be obtained from a given elemental ratio of N storage:

### Evaluating C Storage

The difference between the total cellular C and computed sum of macromolecular C is assumed to be C storage:

### Parameter Estimation

Elemental stoichiometry of each molecule and some parameters are assumed based on available information (Supplementary Table 4). As a result, there are 11 parameters that need to be estimated from the data: *m*, ${v}_{I}^{max}$, *A*_{I}, *A*_{Pho}, *A*_{Bio}, ${Q}_{C}^{Pro-Other}$, ${A}_{RNA}^{P}$, ${Y}_{Pho}^{P:Chl}$, ${Q}_{P}^{Other0}$, ${Q}_{N}^{Sto,max}$, ${Q}_{C}^{Other0}$ (values in Supplementary Table 5). To estimate parameter values that best fit experimental data, we adapt the Metropolis-Hastings algorithm (Metropolis et al., 1953; Hastings, 1970; Omta et al., 2017), a Markov Chain Monte Carlo method. The implementation of the algorithm is described in detail below. We begin by using the Metropolis-Hastings algorithm to fit Equation (17) to the chlorophyll:C vs. μ data for *S. linearis* (Healey, 1985; Figure 4), estimating the values for 3 chlorophyll related parameters (*m*, ${v}_{I}^{max}$, *A*_{I}). We estimate these independently of the other parameters, which do not influence chlorophyll. Next, we estimate nitrogen related parameters (*A*_{Pho}, *A*_{Bio}, ${Q}_{C}^{Pro-Other}$, ${Q}_{N}^{Sto,max}$) and the parameter ${Q}_{C}^{Other0}$, by fitting Equations (33, 50, 30) to the observed data for N:C vs. μ under both N and P limitations (Figure 5A, Supplementary Figure 4Ai) for C concentration vs. μ under N limitation (Supplementary Figure 7A) and for ${\mu}_{max}^{I}$ vs. *I* under both N and P limitations (Figure 8A, Supplementary Figure 5B), respectively. Finally, we estimate the remaining three parameters (${A}_{RNA}^{P}$, ${Y}_{Pho}^{P:Chl}$, ${Q}_{P}^{Other0}$) by fitting Equations (50, 37) to data for the concentration of cellular C (Supplementary Figure 7B) and for P:C (Figure 6A) under P limitation.

For *Pavlova lutheri* (Chalup and Laws, 1990) and *Skeletonema costatum* (Sakshaug and Andersen, 1989), we use the data of chlorophyll (Supplementary Figure 2i) to estimate *m*, ${v}_{I}^{max}$, *A*_{I}, and N:C (Supplementary Figure 2ii) and ${\mu}_{max}^{I}$ (Figures 8B,C) to estimate *A*_{Pho}, *A*_{Bio}, ${Q}_{C}^{Pro-Other}$, ${Q}_{C}^{Other0}$. For these cases, certain parameters are adopted from *S. linearis* since the experimental information is less comprehensive and appropriate data is not available to constrain them; these other parameters have limited influence on the specific model results illustrated here (i.e., Chl:C, N:C, ${\mu}_{max}^{I}$).

#### Algorithm

The model is first solved with an initial set of parameters, which are determined by manually tuning values until model solutions are reasonably consistent with the data. The algorithm then proceeds in a series of steps (a “chain” of steps) that introduce random perturbations to the parameter values. It finds a set of values that provide a good fit to the data by keeping new parameter values that fit the data well (the new parameters become the “current” state of the parameters), and usually discarding others. That is, if the new parameters fit the data much more poorly than the “current” state, there is a high probability they will be rejected. In aggregate, the algorithm evaluates many combinations of parameters in the search for globally optimal solutions.

#### Evaluating the Fit Between Model and Data

For a given set of parameters (*Pset*), beginning with the initial set, we compute a measure of the fit between the model and data. We use the sum of squared errors for each data set between the model (with parameters *Pset*) and data points, given by

where *i* and *k* indicate the *i*th measurement of *k*th data set, *j* is the *j*th iteration (initial step, *j* = 1), $Data({\mu}_{i}^{k},{I}_{i}^{k})$ are data for different growth rates μ and light intensities *I*, $Model({\mu}_{i}^{k},{I}_{i}^{k},Pse{t}_{j})$is the model estimated for the same μ and *I* with a parameter set *Pset*_{j}, and σ_{k} is an estimate of the measurement error of the *k*th dataset. We estimate these values based on the magnitude of scatter among measurements that are made under similar experimental conditions. Once we obtain an error value for each dataset, we normalize the data with the number of data *n*_{k} and add them up to obtain the error covering all the datasets:

This normalization by *n*_{k} is intended to give similar weight to data sets with different numbers of observations and resulted in slightly improved model-data fit for ${\mu}_{max}^{I}$ at high light intensities while keeping other model outputs visually unchanged.

#### Iteration

At each step in the chain, we generate a new parameter set with small random perturbations of the previous set (in this study, within the range of ±20%). In *Pset*_{j} we only accept positive values and chlorophyll related parameters less than certain values (mostly ~5 times of the estimated values), since values outside of these ranges are less likely. We then compare the fit of the *j*th proposed parameter set to the current set, based on the likelihood ratio, given by

Once we obtain the likelihood ratio, we generate a uniform random number (*Random*) between 0 and 1 and compare it with *Ratio*_{j}. If *Ratio*_{j} > *Random*, we update the current parameter set to be the *j*th set. Therefore, if *Error*_{j} is smaller than *Error*_{Current} (i.e., the model with *Pset*_{j} fits better than the current model), the current parameters are updated to *Pset*_{j}. However, if *Error*_{Current} is smaller than *Error*_{j}, the *j*th state will be accepted with a probability that declines as a function of the difference between the error terms. This means good parameter sets tend to be kept, but the acceptance of poorer sets provides a mechanism to get out of the local maxima. After many steps (10^{6} steps), we identified the parameter set that gives the smallest errors between model and data.

## Data Availability Statement

The model code for this study can be found in GitHub https://github.com/ag105020/Phyto1 (10.5281/zenodo.1203718).

## Author Contributions

KI, JB, and MF designed the study. KI, AO, DT, and MF gathered data. KI developed a model and led the project. KI parameterized the model with help of AO, DT, and JB. KI, CD, and MF acquired funding. CD supervised KI. MF advised KI, AO, and DT. All the authors wrote the manuscript.

## Funding

This research was supported by the Japan Student Service Organization (L11171020001, KI), the Gordon and Betty Moore Foundation (GBMF 3775, CD; GBMF 3778, MF), US National Science Foundation (OCE-1558702, MF), and the Simons Foundation (Simons Postdoctoral Fellowship in Marine Microbial Ecology, Award 544338, KI; Simons Collaboration on Ocean Processes and Ecology, Award 329108, MF; Simons Collaboration on Computational Biogeochemical Modeling of Marine Ecosystems, CBIOMES, Award 549931, MF).

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

We thank Rogier Braakman, Stephanie Dutkiewicz, Takako Masuda, Jodi N. Young, and Sallie W. Chisholm for useful discussions. We are grateful to Hedy Kling and Christine Sherratt for sharing the information about *Synechococcus linearis* and members of Writing and Communication Center at MIT for their advice on writing. We appreciate thoughtful and helpful comments and suggestions from the reviewers, which have greatly improved the manuscript.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00086/full#supplementary-material

## References

Ågren, G. I. (2004). The C:N:P stoichiometry of autotrophs - theory and observations. *Ecol. Lett.* 7, 185–191. doi: 10.1111/j.1461-0248.2004.00567.x

Anderson, L. A. (1995). On the hydrogen and oxygen content of marine phytoplankton. *Deep Sea Res. I* 42, 1675–1680. doi: 10.1016/0967-0637(95)00072-E

Biersmith, A., and Benner, R. (1998). Carbohydrates in phytoplankton and freshly produced dissolved organic matter. *Mar. Chem.* 63, 131–144. doi: 10.1016/S0304-4203(98)00057-7

Bonachela, J. A., Allison, S. D., Martiny, A. C., and Levin, S. A. (2013). A model for variable phytoplankton stoichiometry based on cell protein regulation. *Biogeosciences* 10, 4341–4356. doi: 10.5194/bg-10-4341-2013

Bremer, H., and Dennis, P. (1996). “Modulation of chemical composition and other parameters of the cell by growth rate,” in *Escherichia coli and Salmonella typhimurium*, ed F. Neidhardt (Washington, DC: American Society for Microbiology, 1553–1569.

Brown, M. R. R. (1991). The amino-acid and sugar composition of 16 species of microalgae used in mariculture. *J. Exp. Mar. Biol. Ecol.* 145, 79–99. doi: 10.1016/0022-0981(91)90007-J

Bühler, T., Monter, U., Sann, R., Kuhla, J., Dingier, C., and Oelze, J. (1987). Control of respiration and growth yield in ammonium-assimilating cultures of *Azotobacter vinelandii*. *Arch. Microbiol.* 148, 242–246. doi: 10.1007/BF00414819

Buitenhuis, E. T., Hashioka, T., and Le Quéré, C. (2013). Combined constraints on global ocean primary production using observations and models. *Global Biogeochem. Cycles* 27, 847–858. doi: 10.1002/gbc.20074

Burnap, R. L. (2015). Systems and photosystems: cellular limits of autotrophic productivity in cyanobacteria. *Front. Bioeng. Biotechnol.* 3:1. doi: 10.3389/fbioe.2015.00001

Caperon, J. (1968). Population growth response of *Isochrysis galbana* to nitrate variation at limiting concentrations. *Ecology* 49, 866–872. doi: 10.2307/1936538

Caperon, J., and Meyer, J. (1972a). Nitrogen-limited growth of marine phytoplankton-I. Changes in population characteristics with steady-state growth rate*. Deep Sea Res.* 19, 601–618. doi: 10.1016/0011-7471(72)90089-7

Caperon, J., and Meyer, J. (1972b). Nitrogen-limited growth of marine phytoplankton-II. Uptake kinetics and their role in nutrient limited growth of phytoplankton*. Deep Sea Res. Oceanogr. Abstr.* 19, 619–632. doi: 10.1016/0011-7471(72)90090-3

Chalup, M. S., and Laws, E. A. (1990). A test of the assumptions and predictions of recent microalgal growth models with the marine phytoplankter *Pavlova lutheri*. *Limnol. Oceanogr.* 35, 583–596. doi: 10.4319/lo.1990.35.3.0583

Chen, B., and Smith, S. L. (2018). Optimality-based approach for computationally efficient modeling of phytoplankton growth, chlorophyll-to-carbon, and nitrogen-to-carbon ratios. *Ecol. Modell.* 385, 197–212. doi: 10.1016/j.ecolmodel.2018.08.001

Chen, B., Smith, S. L., and Wirtz, K. W. (2019). Effect of phytoplankton size diversity on primary productivity in the North Pacific: trait distributions under environmental variability. *Ecol. Lett.* 22, 56–66. doi: 10.1111/ele.13167

Christie-Oleza, J. A., Sousoni, D., Lloyd, M., Armengaud, J., and Scanlan, D. J. (2017). Nutrient recycling facilitates long-term stability of marine microbial phototroph-heterotroph interactions. *Nat. Microbiol.* 2:17100. doi: 10.1038/nmicrobiol.2017.100

Claquin, P., Martin-Jézéquel, V., Kromkamp, J. C., Veldhuis, M. J. W., and Kraay, G. W. (2002). Uncoupling of silicon compared with carbon and nitrogen metabolisms and the role of the cell cycle in continuous cultures of *Thalassiosira pseudonana* (Bacillariophyceae) under light, nitrogen, and phosphorus control. *J. Phycol.* 38, 922–930. doi: 10.1046/j.1529-8817.2002.t01-1-01220.x

Coles, V. J., Stukel, M. R., Brooks, M. T., Burd, A., Crump, B. C., Moran, M. A., et al. (2017). Ocean biogeochemistry modeled with emergent trait-based genomics. *Science* 1154, 1149–1154. doi: 10.1126/science.aan5712

Collos, Y., Vaquer, A., and Souchu, P. (2005). Acclimation of nitrate uptake by phytoplankton to high substrate levels. *J. Phycol.* 41, 466–478. doi: 10.1111/j.1529-8817.2005.00067.x

Cullen, J. J. (1990). On models of growth and photosynthesis in phytoplankton. *Deep Sea Res.* 37, 667–683. doi: 10.1016/0198-0149(90)90097-F

Deschamps, P., Colleoni, C., Nakamura, Y., Suzuki, E., Putaux, J. L., Buléon, A., et al. (2008). Metabolic symbiosis and the birth of the plant kingdom. *Mol. Biol. Evol.* 25, 536–548. doi: 10.1093/molbev/msm280

Deutsch, C., and Weber, T. (2012). Nutrient ratios as a tracer and driver of ocean biogeochemistry. *Annu. Rev. Mar. Sci.* 4, 113–141. doi: 10.1146/annurev-marine-120709-142821

Dron, A., Rabouille, S., Claquin, P., Chang, P., Raimbault, V., Talec, A., et al. (2012). Light:dark (12:12 h) quantification of carbohydrate fluxes in *Crocosphaera watsonii*. *Aquat. Microb. Ecol.* 68, 43–55. doi: 10.3354/ame01600

Droop, M. R. (1968). Vitamin B_{12} and marine ecology. IV. The kinetics of uptake, growth and inhibition in *Monochrysis lutheri. J. Mar. Biol. Assoc. U.K.* 48, 689–733. doi: 10.1017/S0025315400019238

Dutkiewicz, S., Follows, M. J., and Bragg, J. G. (2009). Modeling the coupling of ocean ecology and biogeochemistry. *Global Biogeochem. Cycles* 23:GB4017. doi: 10.1029/2008GB003405

Dutkiewicz, S., Hickman, A. E., Jahn, O., Gregg, W. W., Mouw, C. B., and Follows, M. J. (2015). Capturing optically important constituents and properties in a marine biogeochemical and ecosystem model. *Biogeosciences* 12, 4447–4481. doi: 10.5194/bg-12-4447-2015

Elrifi, I. R., and Turpin, D. H. (1985). Steady-state luxury consumption and the concept of optimum nutrient ratios: a study with phosphate nitrate limited *Selenastrum minutum* (Chlorophyta). *J. Phycol.* 21, 592–602. doi: 10.1111/j.0022-3646.1985.00592.x

Faizi, M., and Steuer, R. (2019). Optimal proteome allocation strategies for phototrophic growth in a light-limited chemostat. *Microb. Cell Fact.* 18:165. doi: 10.1186/s12934-019-1209-7

Faizi, M., Zavřel, T., Loureiro, C., Červený, J., and Steuer, R. (2018). A model of optimal protein allocation during phototrophic growth. *Biosystems* 166, 26–36. doi: 10.1016/j.biosystems.2018.02.004

Falkowski, P. G., and Owens, T. G. (1980). Light—shade adaptation. *Plant Physiol.* 66, 592–595. doi: 10.1104/pp.66.4.592

Fanesi, A., Wagner, H., Birarda, G., Vaccari, L., and Wilhelm, C. (2019). Quantitative macromolecular patterns in phytoplankton communities resolved at the taxonomical level by single-cell Synchrotron FTIR-spectroscopy. *BMC Plant Biol.* 19:142. doi: 10.1186/s12870-019-1736-8

Fanesi, A., Wagner, H., and Wilhelm, C. (2017). Phytoplankton growth rate modelling: can spectroscopic cell chemotyping be superior to physiological predictors? *Proc. R. Soc.* 284:20161956. doi: 10.1098/rspb.2016.1956

Felcmanová, K., Lukeš, M., Kotabová, E., Lawrenz, E., Halsey, K. H., and Prášil, O. (2017). Carbon use efficiencies and allocation strategies in *Prochlorococcus marinus* strain PCC 9511 during nitrogen-limited growth. *Photosynth. Res.* 134, 71–82. doi: 10.1007/s11120-017-0418-3

Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P. (1998). Primary production of the biosphere: integrating terrestrial and oceanic components. *Science* 281, 237–240. doi: 10.1126/science.281.5374.237

Finkel, Z. V., Follows, M. J., Liefer, J. D., Brown, C. M., Benner, I., and Irwin, A. J. (2016). Phylogenetic diversity in the macromolecular composition of microalgae. *PLoS ONE* 11:e0155977. doi: 10.1371/journal.pone.0155977

Folea, I. M., Zhang, P., Nowaczyk, M. M., Ogawa, T., Aro, E. M., and Boekema, E. J. (2008). Single particle analysis of thylakoid proteins from *Thermosynechococcus elongatus* and *Synechocystis* 6803: localization of the CupA subunit of NDH-1. *FEBS Lett.* 582, 249–254. doi: 10.1016/j.febslet.2007.12.012

Galbraith, E. D., and Martiny, A. C. (2015). A simple nutrient-dependence mechanism for predicting the stoichiometry of marine ecosystems. *Proc. Natl. Acad. Sci. U.S.A*. 112:201423917. doi: 10.1073/pnas.1423917112

Garcia, N. S., Bonachela, J. A., and Martiny, A. C. (2016). Growth-dependent cell size controls interactions between nutrient supply and cellular elemental stoichiometry of marine *Synechococcus*. *ISME J.* 10, 2715–2724. doi: 10.1038/ismej.2016.50

Geider, R. J., and La Roche, J. (2002). Redfield revisited: variability of C : N : P in marine microalgae and its biochemical basis. *Eur. J. Phycol.* 37, 1–17. doi: 10.1017/S0967026201003456

Geider, R. J., Macintyre, H. L., and Kana, T. M. (1996). A dynamic model of photoadaptation in phytoplankton. *Limnol. Oceanogr.* 41, 1–15. doi: 10.4319/lo.1996.41.1.0001

Geider, R. J., Macintyre, H. L., and Kana, T. M. (1998). A dynamic regulatory model of phytoplanktonic acclimation to light, nutrients, and temperature. *Limnol. Oceanogr.* 43, 679–694. doi: 10.4319/lo.1998.43.4.0679

Geider, R. J., Mooreb, C. M., and Rossc, O. N. (2009). The role of cost-benefit analysis in models of phytoplankton growth and acclimation. *Plant Ecol. Divers.* 2, 165–178. doi: 10.1080/17550870903300949

Geider, R. J., Osborne, B. A., and Raven, J. A. (1985). Light dependence of growth and photosynthesis in *Phaeodactylum tricornutum* (Bacillariophyceae). *J. Phycol.* 21, 609–619. doi: 10.1111/j.0022-3646.1985.00609.x

Grover, J. P. (2009). Is storage an adaptation to spatial variation in resource availability? *Am. Nat.* 173, E44–E61. doi: 10.1086/595751

Grover, J. P. (2011). Resource storage and competition with spatial and temporal variation in resource availability. *Am. Nat.* 178, E124–E148. doi: 10.1086/662163

Harrison, P. J., Thompson, P. A., and Calderwood, G. S. (1990). Effects of nutrient and light limitation on the biochemical composition of phytoplankton. *J. Appl. Phycol.* 2, 45–56. doi: 10.1007/BF02179768

Hastings, W. K. (1970). Monte Carlo sampling methods using Markov chains and their applications. *Biometrika* 57, 97–109. doi: 10.1093/biomet/57.1.97

Healey, F. P. (1985). Interacting effects of light and nutrient limitation on the growth rate of *Synechococcus linearis* (Cyanophyceae). *J. Phycol.* 21, 134–146. doi: 10.1111/j.0022-3646.1985.00134.x

Ho, T., Quigg, A., Zoe, V., Milligan, A. J., Falkowski, P. G., and Morel, F. M. M. (2003). The elemental composition of some marine phytoplankton. *J. Phycol.* 39, 1145–1159. doi: 10.1111/j.0022-3646.2003.03-090.x

Inomura, K., Bragg, J., and Follows, M. J. (2017). A quantitative analysis of the direct and indirect costs of nitrogen fixation: a model based on *Azotobacter vinelandii*. *ISME J.* 11, 166–175. doi: 10.1038/ismej.2016.97

Inomura, K., Bragg, J., Riemann, L., and Follows, M. J. (2018). A quantitative model of nitrogen fixation in the presence of ammonium. *PLoS ONE* 13:e0208282. doi: 10.1371/journal.pone.0208282

Jahn, M., Vialas, V., Karlsen, J., Maddalo, G., Edfors, F., Forsström, B., et al. (2018). Growth of cyanobacteria is constrained by the abundance of light and carbon assimilation proteins. *Cell Rep.* 25, 478–486. doi: 10.1016/j.celrep.2018.09.040

Kirchhoff, H. (2014). Diffusion of molecules and macromolecules in thylakoid membranes. *Biochim. Biophys. Acta Bioenerg.* 1837, 495–502. doi: 10.1016/j.bbabio.2013.11.003

Kirchhoff, H., Haferkamp, S., Allen, J. F., Epstein, D. B., and Mullineaux, C. W. (2008). Protein diffusion and macromolecular crowding in thylakoid membranes. *Plant Physiol.* 146, 1571–1578. doi: 10.1104/pp.107.115170

Klausmeier, C. A., Litchman, E., Daufresne, T., and Levin, S. A. (2004). Optimal nitrogen-to-phosphorus stoichiometry of phytoplankton. *Nature* 429, 171–174. doi: 10.1038/nature02454

Kooijman, S. A. L. M. (2010). *Dynamic Energy Budget Theory for Metabolic Organisation, 3rd Edn*. Cambridge: Cambridge University Press. doi: 10.1017/CBO9780511805400

Laws, E., and Caperon, J. (1976). Carbon and nitrogen metabolism by *Monochrysis lutheri*: measurement of growth-rate-dependent respiration rates. *Mar. Biol.* 36, 85–97. doi: 10.1007/BF00388431

Laws, E. A., and Bannister, T. T. (1980). Nutrient- and light-limited growth of *Thalassiosira fluviatilis* in continuous culture, with implications for phytoplankton growth in the ocean. *Limnol. Oceanogr.* 25, 457–473. doi: 10.4319/lo.1980.25.3.0457

Laws, E. A., and Wong, D. C. L. (1978). Studies of carbon and nitrogen metabolism by 3 marine phytoplankton species in nitrate-limited continuous culture. *J. Phycol.* 14, 406–416. doi: 10.1111/j.1529-8817.1978.tb02460.x

Lengeler, J. W., Drews, G., and Schlegel, H. G. (1999). *Biology of the Prokaryotes.* Stuttgart: Thieme. doi: 10.1002/9781444313314

Levasseur, M., Thompson, P. A., and Harrison, P. J. (1993). Physiological acclimation of marine phytoplankton to different nitrogen sources. *J. Phycol.* 29, 587–595. doi: 10.1111/j.0022-3646.1993.00587.x

Lewis, K. M., Arntsen, A. E., Coupel, P., Joy-Warren, H., Lowry, K. E., Matsuoka, A., et al. (2019). Photoacclimation of Arctic Ocean phytoplankton to shifting light and nutrient limitation. *Limnol. Oceanogr.* 64, 284–301. doi: 10.1002/lno.11039

Lichtenthaler, H. K., Kuhn, G., Prenzel, U., and Meier, D. (1982). Chlorophyll-protein levels and degree of thylakoid stacking in radish chloroplasts from high-light, low-light and bentazon-treated plants. *Physiol. Plant.* 56, 183–188. doi: 10.1111/j.1399-3054.1982.tb00322.x

Liefer, J. D., Garg, A., Fyfe, M. H., Irwin, A. J., Benner, I., Brown, C. M., et al. (2019). The macromolecular basis of phytoplankton C:N:P under nitrogen starvation. *Front. Microbiol.* 10:763. doi: 10.3389/fmicb.2019.00763

Lin, S., Litaker, R. W., and Sunda, W. G. (2016). Phosphorus physiological ecology and molecular mechanisms in marine phytoplankton. *J. Phycol.* 52, 10–36. doi: 10.1111/jpy.12365

Litchman, E., and Klausmeier, C. A. (2008). Trait-based community ecology of phytoplankton. *Annu. Rev. Ecol. Evol. Syst.* 39, 615–639. doi: 10.1146/annurev.ecolsys.39.110707.173549

Litchman, E., Neale, P. J., and Banaszak, A. T. (2002). Increased sensitivity to ultraviolet radiation in nitrogen-limited dinoflagellates: photoprotection and repair. *Limnol. Oceanogr.* 47, 86–94. doi: 10.4319/lo.2002.47.1.0086

Litchman, E., Steiner, D., and Bossard, P. (2003). Photosynthetic and growth responses of three freshwater algae to phosphorus limitation and daylength. *Freshw. Biol.* 48, 2141–2148. doi: 10.1046/j.1365-2427.2003.01157.x

Liu, H., Bidigare, R. R., Laws, E., Landry, M. R., and Campbell, L. (1999). Cell cycle and physiological characteristics of *Synechococcus* (WH7803) in chemostat culture. *Mar. Ecol. Prog. Ser.* 189, 17–25. doi: 10.3354/meps189017

Loladze, I., and Elser, J. J. (2011). The origins of the Redfield nitrogen-to-phosphorus ratio are in a homoeostatic protein-to-rRNA ratio. *Ecol. Lett.* 14, 244–250. doi: 10.1111/j.1461-0248.2010.01577.x

McKew, B. A., Lefebvre, S. C., Achterberg, E. P., Metodieva, G., Raines, C. A., Metodiev, M. V., et al. (2013). Plasticity in the proteome of *Emiliania huxleyi* CCMP 1516 to extremes of light is highly targeted. *New Phytol.* 200, 61–73. doi: 10.1111/nph.12352

McKew, B. A., Metodieva, G., Raines, C. A., Metodiev, M. V., and Geider, R. J. (2015). Acclimation of *Emiliania huxleyi* (1516) to nutrient limitation involves precise modification of the proteome to scavenge alternative sources of N and P. *Environ. Microbiol.* 17, 4050–4062. doi: 10.1111/1462-2920.12957

Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., and Teller, A. H. (1953). Equation of state calculations by fast computing machines. *J. Chem. Phys.* 21, 1087–1092. doi: 10.1063/1.1699114

Michal, G. (1999). *Biochemical Pathways: An Atlas of Biochemistry and Molecular Biology*. Heidelberg: Wiley & Spektrum.

Mitra, A., Flynn, K. J., and Fasham, M. J. R. (2007). Accounting for grazing dynamics in nitrogen-phytoplankton-zooplankton models. *Limnol. Oceanogr.* 52, 649–661. doi: 10.4319/lo.2007.52.2.0649

Monod, J. (1949). The growth of bacterial cultures. *Ann. Rev. Mar. Sci.* 3, 371–394. doi: 10.1146/annurev.mi.03.100149.002103

Moore, J. K., Doney, S. C., and Lindsay, K. (2004). Upper ocean ecosystem dynamics and iron cycling in a global three-dimensional model. *Global Biogeochem. Cycles* 18:GB4028. doi: 10.1029/2004GB002220

Moore, L. R., and Chisholm, S. W. (1999). Photophysiology of the marine cyanobacterium *Prochlorococcus*: ecotypic differences among cultured isolates. *Limnol. Oceanogr.* 44, 628–638. doi: 10.4319/lo.1999.44.3.0628

Moreno, A. R., and Martiny, A. C. (2018). Ecological stoichiometry of ocean plankton. *Ann. Rev. Mar. Sci.* 10, 43–69. doi: 10.1146/annurev-marine-121916-063126

Neidhardt, F. C., Ingraham, J. L., and Schaechter, M. (1990). *Physiology of the Bacterial Cell: A Molecular Approach*. Sunderland, MA: Sinauer Associate, Inc.

Nicklisch, A., and Steinberg, C. E. W. (2009). RNA/protein and RNA/DNA ratios determined by flow cytometry and their relationship to growth limitation of selected planktonic algae in culture. *Eur. J. Phycol.* 44, 297–308. doi: 10.1080/09670260802578518

Omta, A. W., Talmy, D., Sher, D., Finkel, Z. V., Irwin, A. J., and Follows, M. J. (2017). Extracting phytoplankton physiological traits from batch and chemostat culture data. *Limnol. Oceanogr. Methods* 15, 453–466. doi: 10.1002/lom3.10172

Orth, J. D., Thiele, I., and Palsson, B. Ø. (2010). What is flux balance analysis? *Nat. Biotechnol.* 28, 245–248. doi: 10.1038/nbt.1614

Paasche, E. (1973). Silicon and the ecology of marine plankton diatoms. I. *Thalassiosira pseudonana* (*Cyclotella nana*) grown in a chemostat with silicate as limiting nutrient. *Mar. Biol.* 19, 117–126. doi: 10.1007/BF00353582

Pahlow, M., and Oschlies, A. (2009). Chain model of phytoplankton P, N and light colimitation. *Mar. Ecol. Prog. Ser.* 376, 69–83. doi: 10.3354/meps07748

Parsons, T. R., Takahashi, M., and Hargrave, B. (1984). *Biological Oceanographic Processes, Vol. 330*. Oxford: Pergamon Press. doi: 10.1016/B978-0-08-030765-7.50010-6

Pedersen, M. F., and Borum, J. (1996). Nutrient control of algal growth in estuarine waters. Nutrient limitation and the importance of nitrogen requirements and nitrogen storage among phytoplankton and species of macroalgae*. Mar. Ecol. Prog. Ser*. 142, 261–272. doi: 10.3354/meps142261

Perry, M. J. (1976). Phosphate utilization by an oceanic diatom in phosphorus-limited chemostat culture and in the oligotrophic waters of the central North Pacific. *Limnol. Oceanogr.* 21, 88–107. doi: 10.4319/lo.1976.21.1.0088

Platt, T., Gallegos, C. L., and Harrison, W. G. (1980). Photoinhibitioin of photosynthesis in natural assemblages of marine phytoplankton. *J. Mar. Res.* 38, 687–701.

Quigg, A., Finkel, Z. V., Irwin, A. J., Rosenthal, Y., Ho, T. Y., Reinfelder, J. R., et al. (2003). The evolutionary inheritance of elemental stoichiometry in marine phytoplankton. *Nature* 425, 291–294. doi: 10.1038/nature01953

Quigg, A., Irwin, A. J., and Finkel, Z. V (2011). Evolutionary inheritance of elemental stoichiometry in phytoplankton. *Proc. R. Soc. B* 278, 526–534. doi: 10.1098/rspb.2010.1356

Redfield, A. C. (1934). “On the proportions of organic derivations in sea water and their relation to the composition of plankton,” in *James Johnstone Memorial Volume*, ed R. J. Daniel (Liverpool: Liverpool Universtiy Press, 177–192.

Redfield, A. C. (1958). The biological control of chemical factors in the environment. *Am. Sci.* 46, 205–221.

Reimers, A.-M., Knoop, H., Bockmayr, A., and Steuer, R. (2017). Cellular trade-offs and optimal resource allocation during cyanobacterial diurnal growth. *Proc. Natl. Acad. Sci. U.S.A*. 114, E6457–E6465. doi: 10.1073/pnas.1617508114

Rhee, A. G., Gotham, I. J., and Rhee, G. (1981). The effect of environmental factors on phytoplankton growth: light and the interactions of light with nitrate limitation. *Limnol. Oceanogr.* 26, 649–659. doi: 10.4319/lo.1981.26.4.0649

Rhee, G.-Y. (1978). Effects of N:P atomic ratios and nitrate limitation on algal growth, cell compostion, and nitrate uptake. *Limnol. Oceanogr.* 23, 10–25. doi: 10.4319/lo.1978.23.1.0010

Riley, G. A. (1946). Factors controlling phytoplankton populations on Georges Bank. *J. Mar. Res.* 6, 54–73.

Rittmann, B. E., and McCarty, P. L. (2001). *Environmental Biotechnology: Principles and Applications*. New York, NY: McGraw-Hill.

Saito, M. A., Bertrand, E. M., Dutkiewicz, S., Bulygin, V. V., Moran, D. M., Monteiro, F. M., et al. (2011). Iron conservation by reduction of metalloenzyme inventories in the marine diazotroph *Crocosphaera watsonii*. *Proc. Natl. Acad. Sci. U.S.A*. 108, 2184–2189. doi: 10.1073/pnas.1006943108

Sakshaug, E., and Andersen, K. (1989). A steady state description of growth and light absorption in the marine planktonic diatom *Skeletonema costatum*. *Limnol. Oceanogr.* 34, 198–205. doi: 10.4319/lo.1989.34.1.0198

Scott, M., Gunderson, C. W., Mateescu, E. M., Zhang, Z., and Hwa, T. (2010). Interdependence of cell growth and gene expression: origins and consequences. *Science* 330, 1099–1103. doi: 10.1126/science.1192588

Shifrin, N. S., and Chisholm, S. W. (1981). Phytoplankton lipids: interspecific differences and effects of nitrate, silicate and light-dark cycles. *J. Phycol.* 17, 374–384. doi: 10.1111/j.0022-3646.1981.00374.x

Shuter, B. (1979). A model of physiological adaptation in unicellular algae. *J. Theor. Biol.* 78, 519–552. doi: 10.1016/0022-5193(79)90189-9

Smith, S. L., Pahlow, M., Merico, A., Acevedo-Trejos, E., Sasai, Y., Yoshikawa, C., et al. (2016). Flexible phytoplankton functional type (FlexPFT) model: size-scaling of traits and optimal growth. *J. Plankton Res.* 38, 977–992. doi: 10.1093/plankt/fbv038

Smith, S. L., Yamanaka, Y., Pahlow, M., and Oschlies, A. (2009). Optimal uptake kinetics: physiological acclimation explains the pattern of nitrate uptake by phytoplankton in the ocean. *Mar. Ecol. Prog. Ser.* 384, 1–12. doi: 10.3354/meps08022

Sterner, R. W., and Elser, J. J. (2002). *Ecological Stoichiometry: The Biology of Elements from Molecules to the Biosphere.* Princeton, NJ: Princeton University Press. doi: 10.1515/9781400885695

Thrane, J. E., Hessen, D. O., and Andersen, T. (2017). Plasticity in algal stoichiometry: experimental evidence of a temperature-induced shift in optimal supply N:P ratio. *Limnol. Oceanogr.* 62, 1346–1354. doi: 10.1002/lno.10500

Tozzi, S., Schofield, O., and Falkowski, P. (2004). Historical climate change and ocean turbulence as selective agents for two key phytoplankton functional groups. *Mar. Ecol. Prog. Ser.* 274, 123–132. doi: 10.3354/meps274123

Vadeboncoeur, Y., Zanden, M. J., and Vander Lodge, D. M. (2002). Putting the lake back together: reintegrating benthic pathways into lake food web models. *Bioscience* 52, 44–54. doi: 10.1641/0006-3568(2002)052[0044:PTLBTR]2.0.CO;2

Van Mooy, B. A., Fredricks, H. F., Pedler, B. E., Dyhrman, S. T., Karl, D. M., Koblízek, M., et al. (2009). Phytoplankton in the ocean use non-phosphorus lipids in response to phosphorus scarcity. *Nature* 458, 69–72. doi: 10.1038/nature07659

Werner, D., (ed.). (1977). “Silicate metabolism,” in *The Biology of the Diatoms* (Berkeley, CA: University of California Press), 110–149.

Xu, H., Paerl, H. W., Qin, B., Zhu, G., and Gao, G. (2010). Nitrogen and phosphorus inputs control phytoplankton growth in eutrophic Lake Taihu, China. *Limnol. Oceanogr.* 55, 420–432. doi: 10.4319/lo.2010.55.1.0420

Keywords: phytoplankton, elemental stoichiometry, growth rate, macromolecule, photosynthesis, protein, RNA, nutrient storage

Citation: Inomura K, Omta AW, Talmy D, Bragg J, Deutsch C and Follows MJ (2020) A Mechanistic Model of Macromolecular Allocation, Elemental Stoichiometry, and Growth Rate in Phytoplankton. *Front. Microbiol.* 11:86. doi: 10.3389/fmicb.2020.00086

Received: 10 June 2019; Accepted: 15 January 2020;

Published: 28 February 2020.

Edited by:

Pia H. Moisander, University of Massachusetts Dartmouth, United StatesReviewed by:

Ralf Steuer, Humboldt University of Berlin, GermanyJuan A. Bonachela, Rutgers, The State University of New Jersey, United States

Copyright © 2020 Inomura, Omta, Talmy, Bragg, Deutsch and Follows. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Keisuke Inomura, ag105020@gmail.com