Skip to main page content

• HOME
• Current Issue
• Archive
• Alerts
• About ASM
• Contact us
• Tech Support
• Journals.ASM.org

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Rabouille, S. Articles by Soetaert, K. Search for Related Content PubMed PubMed Citation Articles by Rabouille, S. Articles by Soetaert, K. Agricola Articles by Rabouille, S. Articles by Soetaert, K.

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, May 2006, p. 3217-3227, Vol. 72, No. 5
0099-2240/06/$08.00+0     doi:10.1128/AEM.72.5.3217-3227.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.

Modeling the Dynamic Regulation of Nitrogen Fixation in the Cyanobacterium Trichodesmium sp.

Sophie Rabouille,^* Marc Staal, Lucas J. Stal, and Karline Soetaert

Centre for Estuarine and Marine Ecology, Netherlands Institute of Ecology, KNAW, Yerseke, The Netherlands

Received 28 December 2005/ Accepted 21 February 2006

Top Abstract Introduction Materials and Methods Results Discussion References

 
A physiological, unbalanced model is presented that explicitly describes growth of the marine cyanobacterium Trichodesmium sp. at the expense of N₂ (diazotrophy). The model involves the dynamics of intracellular reserves of carbon and nitrogen and allows the uncoupling of the metabolism of these elements. The results show the transient dynamics of N₂ fixation when combined nitrogen (NO₃^–, NH₄⁺) is available and the increased rate of N₂ fixation when combined nitrogen is insufficient to cover the demand. The daily N₂ fixation pattern that emerges from the model agrees with measurements of rates of nitrogenase activity in laboratory cultures of Trichodesmium sp. Model simulations explored the influence of irradiance levels and the length of the light period on fixation activity and cellular carbon and nitrogen stoichiometry. Changes in the cellular C/N ratio resulted from allocations of carbon to different cell compartments as demanded by the growth of the organism. The model shows that carbon availability is a simple and efficient mechanism to regulate the balance of carbon and nitrogen fixed (C/N ratio) in filaments of cells. The lowest C/N ratios were obtained when the light regime closely matched nitrogenase dynamics.

Top Abstract Introduction Materials and Methods Results Discussion References

 
In marine waters, nitrogen generally controls primary production (16). In such environments, N₂-fixing organisms have a competitive advantage over organisms that rely on the availability of combined nitrogen. In the tropical oceans, N₂-fixing cyanobacteria can be extremely abundant in surface waters and account for a considerable input of combined nitrogen into the upper mixed layer (11, 28), with a strong impact on local community production (9).

At ocean basin scales, N₂-fixing cyanobacteria affect the coupling of C-N-P cycles and contribute considerably to the net oceanic sequestration of atmospheric carbon dioxide (28). For instance, a global-scale estimate made by Lee et al. (31) of net CO₂ fixation in the absence of measurable nitrate led to the conclusion that 20 to 40% of the total new primary production in tropical and subtropical oceans could be attributed to N₂-fixing organisms. The quantitative impact of N₂-fixing organisms on nutrient cycling in the oceans and on the global carbon budget is widely acknowledged (14, 16, 26, 36, 51). Biogeochemical models have been developed to offer a dynamic view of biological and biochemical systems. At global scales, these models aim to provide estimates of the main oceanic biogeochemical fluxes, such as total, new, and regenerated production (17). Others describe the global nitrogen cycle at the ecosystem level (18, 25, 35). At smaller scales, phytoplankton growth models describe time-dependent changes in biomass or numbers as a function of one or several limiting factors, thereby offering support to hypotheses about biological and physiological processes. However, these models do not account for the physiological and regulatory mechanisms of N₂ fixation in cyanobacteria, since there is still very little knowledge about the factors that control them.

In this paper we focus on the marine nonheterocystous cyanobacterium Trichodesmium sp., identified as one of the most significant N₂-fixing organisms in oceans (6). Because nitrogenase, the enzyme complex responsible for the reduction of N₂, is very sensitive to inactivation by O₂, cyanobacteria have evolved various strategies to separate N₂ fixation from O₂-generating photosynthesis (2, 4, 21). It has been hypothesized that Trichodesmium separates N₂ fixation spatially from oxygenic photosynthesis, allowing it to fix N₂ during the day (1). In this respect, Trichodesmium has adopted a protective strategy similar to that of heterocystous cyanobacteria (4). Since C/N ratios in natural populations of Trichodesmium remain relatively constant (10), it is expected that Trichodesmium will show balanced growth.

This paper presents a model designed to assess the mechanisms that control physiological processes, in particular primary production. This model explicitly describes unbalanced growth and N₂ fixation in Trichodesmium spp. and their control by environmental factors. Different numerical simulations were performed under various conditions to test the effects of light and nutrient availability on growth and N₂ fixation under both transient and steady-state conditions. We aimed to pinpoint (i) whether light intensity and the temporal distribution of light would have an effect on the pattern and rate of N₂ fixation and (ii) the analysis of the role of carbon availability on N₂ fixation. We also analyze, in qualitative terms, the role of the nitrogen supply in (i) the interaction between nitrogen limitation and N₂ fixation and (ii) N₂ fixation dynamics. The results of these simulations are compared to measurements of nitrogenase activity carried out on exponentially growing cultures of Trichodesmium.

Top Abstract Introduction Materials and Methods Results Discussion References

 
Laboratory experiments.
Trichodesmium sp. strain IMS101 was grown at 27°C in modified YCBII medium (12) in 250-ml Erlenmeyer flasks under a 12-h light/12-h dark (L12) regime in an incubator without shaking. YCBII medium was modified by the addition of 15.9 g liter^–1 Na₂CO₃ and 1.6 x 10^–9 M Na₂SeO₃ and was devoid of combined nitrogen. The pH of the medium was 8.2. Light was provided by cool white 15-W fluorescent tubes at an incident photon irradiance of 70 µmol · m^–2 · s^–1. The cultures were swirled daily by hand in order to prevent the formation of large aggregates or wall growth. N₂ fixation was monitored by using the acetylene reduction technique. Only exponentially growing cultures were used for our experiments. Samples of 50 ml of culture were filtered on a GF/F glass fiber filter (diameter, 47 mm) under a moderate vacuum. Subsequently, the filter was placed in the incubation chamber for an online acetylene reduction assay (47). Filter incubation was preferred to liquid incubation, because it has been shown to be superior (47). The gas flow over the filter with Trichodesmium was 2 liters · h^–1 and was composed of 20% O₂, 70% N₂, and 10% C₂H₂. The incubation chamber was kept at 25°C and at the same light regime and photon flux density as in the culture. The gas flow was in line with a gas chromatograph (Chrompack-CP9000; Varian, The Netherlands) equipped with a flame ionization detector and a 500-µl sample loop for automated injections. The conditions of the gas chromatograph and other details were those described by Staal et al. (47).

N₂ fixation model. (i) Internal pools.
Depending on their structure, growth models can or cannot represent fluctuations in cellular composition and thus the adaptation of the organism to its environment. The Trichodesmium N₂ fixation model is a physiological growth model that explicitly includes reserves of carbon and nitrogen whose dynamics are unbalanced. The model describes the interaction between nitrogen and carbon metabolism and emphasizes the biotic and abiotic conditions that control N₂ fixation rates. The model is described below, and a flow diagram is given in Fig. 1. The model variables are listed in Table 1; see Table 2 for the model formulation, Table 3 for the modeled processes, and Table 4 for parameter values.

View larger version (23K): [in this window] [in a new window]

FIG. 1. Schematic representation of the four compartments in the model and their related mass fluxes. Solid arrows, C flux; dashed arrows, N flux; dotted arrows, losses due to mortality (M) and maintenance costs (R). The units (C and/or N) used to describe the content of each compartment are circled.

View this table: [in this window] [in a new window]

TABLE 1. Model variables

View this table: [in this window] [in a new window]

TABLE 2. Model formulation

View this table: [in this window] [in a new window]

TABLE 3. Modeled processes

View this table: [in this window] [in a new window]

TABLE 4. Parameter values

The cellular carbon components are decomposed as described by Lancelot and Billen (29). Cellular carbon is represented by the pool of low-molecular-weight carbon molecules (photosynthate) (C_sm), the internal carbon reserves (glycogen) (C_r), and the structural biomass (C_s). C_sm has a short turnover time, allowing fast changes of process rates. In the model this pool served as the intermediate for all metabolic processes. The C_sm consumption terms can be either (i) storage in the form of carbon reserves (C_r), (ii) incorporation in structural biomass (C_s), (iii) dissimilation, or (iv) consumption for maintenance. An exchange thus occurs between the intracellular carbon reserves (C_r) and C_sm. Similarly, internal storage of nitrogen (N_r) is included in the model in order to account for temporal variations in N₂ fixation. This nitrogen reserve groups glutamine (the early product of N₂ fixation [7]) and the pool of nitrogen reserve components (e.g., cyanophycin and phycobiliproteins). This pool represents intracellular nitrogen available for the synthesis of structural biomass (C_s) (structural carbon and nitrogen). The uptake of external combined nitrogen also feeds into the nitrogen reserve pool.

The structural biomass (C_s) pool sets the maximum rates of all processes. The uncoupling between carbon and nitrogen metabolism allows processes to be modulated (limited or inhibited) by both the carbon status and the nitrogen status of cells (C_e and N_e, respectively), which represent the availability of these two elements for metabolism. Carbon status is calculated as the relative proportion of the photosynthate (C_sm) to total carbon, excluding the minimum cellular carbon quota; nitrogen status is the relative proportion of the nitrogen reserves (N_r) to total N, excluding the minimum cellular nitrogen quota (Table 2). Hence, growth can occur independently of the periods of photosynthesis, provided that sufficient nitrogen and carbon are available to cells. The limitation functions (l_C and l_N) associated with the internal states follow a type III (sigmoid) response (Table 2). Inhibition by internal status is expressed as 1 minus sigmoid.

(ii) Formulations of photosynthesis and growth.
The carbon-specific photosynthesis rate (P) depends on the chlorophyll a (Chla) content. Light dependence [g₁(I)] is described using the exponential saturation function from reference 41, including a photoinhibition term. The actual photosynthesis rate is further limited by the extent of carbon reserves in cells (1 – l_C). The maximum rate of photosynthesis, P_m (expressed as mol of C · g of Chla^–1 · h^–1), was recalculated from Kana's measurements (27) on Trichodesmium thiebautii. The rate of carbon storage in C_r (g₃) is assumed to be proportional to the carbon status following a saturation function. The consumption rate of this reserve (g₂) is linearly related to the reserve concentration (C_r). Growth (G) is the production of structural biomass at the expense of photosynthate and nitrogen reserves. We considered a maximum theoretical specific growth rate of Trichodesmium of 0.9 day^–1. The actual growth rate results from the limitation by carbon and nitrogen internal states, following Liebig's law of the minimum. Phosphorus is not limiting in our model, and the effect of temperature is not taken into account.

(iii) Effects of nitrogen substrates and associated growth costs.
Trichodesmium is a facultative diazotroph and can take up combined inorganic nitrogen compounds (38). In the model, combined nitrogen is preferred over N₂, a condition based on the absence or lower rates of N₂ fixation when Trichodesmium is grown on NH₄⁺ or NO₃^– (24). The uptake and assimilation of combined nitrogen are described by one equation, in which the rates of consumption of NH₄⁺ (U₁) and NO₃^– (U₂) are proportional to their concentrations. The rate of consumption follows Michaelis-Menten kinetics and is modulated by the depletion of internal nitrogen reserves. As Stephens et al. (49) proposed, we simulate the potential activity of N₂ fixation as the result of a balance between an increase (g₄) and a decrease (g₈) in the potential activity rather than as a turnover of the enzyme itself. Nitrogenase activity is dependent on N₂, ATP, and reducing equivalents. The increase in the potential rate of N₂ fixation is assumed to be controlled by the nitrogen (g₆) and energy status of the cells. It is also sensitive to the concentration of NO₃^– (40). The actual rate of N₂ fixation results from carbon availability and the cellular nitrogen status as well as the potential rate of nitrogenase activity.

ATP pools are not modeled sensus stricto; the energy gained from photosynthesis is represented as the accumulation of fixed carbon. Thus, the energy for any metabolic process is provided by the dissimilation of the carbon reserves. For the model, we assume that N₂ fixation is limited by electrons originating from the mobilization of the carbon metabolites (C_sm). Respiration accounts for differences in energy costs that result from the use of different nitrogen sources. When growth cannot be sustained by the uptake of combined nitrogen, N₂ fixation is induced in order to meet nitrogen demand. The transition to N₂ fixation will moderately increase the energy cost for N assimilation, which is reflected in a higher rate of respiration. The higher energy cost of growth on NO₃^– than on NH₄⁺ is also taken into account. The minimum maintenance rate increases when the concentration of NH₄⁺ decreases.

(iv) Physical setting.
The modeled system can be run either as a well-stirred continuous culture or as a turbidostat. We emphasize that the term chemostat does not apply here, since that term usually refers to cultures in steady state controlled by a growth-limiting nutrient (often but not exclusively nitrogen) and requires a continuous light regime. In this study a light-dark regime was applied, and Trichodesmium growth was light limited and not nutrient limited.

Biomass generates an extinction of light (e_I). Hence, the higher the biomass, the lower the resultant average irradiance () within the vessel becomes, lowering the specific photosynthesis rate. Biomass is considered to be well mixed and as a result homogeneously distributed. Growth will be light limited at the equilibrium state. Dilution removes part of the biomass and prevents culture senescence, since light is never completely extinguished. When the model reaches equilibrium, the daily average of all the specific rates becomes constant and the daily net specific growth rate equals the dilution rate. To qualify this equilibrium, we do not use here the term "steady state," which refers to a perfect constancy of variables in time. In the model, the light regime induces a periodicity in the metabolism, and at equilibrium, variables are still fluctuating within 24 h but their daily average is constant from one day to another. Such an equilibrium is called "dynamic equilibrium." In this model, the dynamic equilibrium is determined by light conditions and by the dilution rate. When the model is run as a turbidostat, the dilution rate becomes a state variable that fluctuates until the average irradiance () in the culture reaches a chosen value, I_eq (Table 2). I_eq characterizes the light regime reached in the turbidostat at the equilibrium.

(v) Implementation.
The model is formulated as a system of ordinary differential equations. It is implemented in the FORTRAN language, in the FEMME simulation environment (44), running on Windows. The model is available for download from the FEMME website (http://www.nioo.knaw.nl/cemo/femme/).

Simulations.
The model was run as a continuous culture and as a turbidostat. In the continuous-culture runs, the daily average of the net specific growth rate equaled the dilution rate when the model reached the dynamic equilibrium. The turbidostat allowed better control of the light intensity in the culture. In each turbidostat simulation, the average desired irradiance in the culture at equilibrium (I_eq) was set as a parameter. The dilution rate was not constant in these experiments but changed at each time step, resulting in a constant daily average of the biomass concentration that created the light conditions in the vessel equal to I_eq.

(i) Transient behavior experiments.
The first two simulations consisted of a well-mixed continuous culture in which the initial concentration of combined nitrogen in the culture vessel was 10 µM and the initial total biomass was 6.8 mmol of C · m^–3. The medium supply did not contain combined nitrogen. The culture was exposed to a constant incident photon irradiance of 120 µmol · m^–2 · s^–1 during the light period.

In the L12 simulation, the on-off light rhythm applied was 12 h of light and 12 h of darkness (hence, the daily average was 60 µmol · m^–2 · s^–1). The model was first run as a continuous culture with a dilution rate set to 0.0085 h^–1. Then, under the same conditions, the model was run as a turbidostat, with an I_eq of 12.5 µmol · m^–2 · s^–1 during the light period.

In the L16 simulation, the length of the light period was changed to 16 h of light and 8 h of darkness. The same incident photon irradiance was applied (120 µmol · m^–2 · s^–1 during the 16 h of light), and the model was run as a turbidostat, with an I_eq of 9.375 µmol · m^–2 · s^–1 during the light period. The total average daily light dose experienced by the turbidostat culture at the dynamic equilibrium was then the same as that in the turbidostat culture in the L12 simulation.

(ii) Equilibrium state analysis.
One hundred simulations were performed in order to investigate the influence of the average light intensity in the culture vessel on the stoichiometry of cellular carbon and nitrogen. The model was run as a turbidostat, and I_eq was uniformly varied between predefined ranges.

Series 1 consisted of 100 simulations similar to the L12 experiments (i.e., under a 12-h light/12-h dark cycle) with I_eq values ranging from 5 to 100 µmol · m^–2 · s^–1. The incident photon irradiance applied during the light phase was again constant in all runs (120 µmol · m^–2 · s^–1); the total daily light dose remained the same.

Series 2 consisted of 100 simulations to analyze the coinfluence of the average light intensity and the length of the light period on cell metabolism. Both I_eq (5 to 100 µmol · m^–2 · s^–1) and the length of the light period (8 h to 18 h) were changed during these simulations, while the incident photon irradiance applied during the light phase was again constant in all runs (120 µmol · m^–2 · s^–1).

Top Abstract Introduction Materials and Methods Results Discussion References

 
Figure 2 displays experimental data; Fig. 3 shows simulation results with the model run as a continuous culture; and Fig. 4 to 6 show results of the simulation with the model run as a turbidostat.

View larger version (15K): [in this window] [in a new window]

FIG. 2. Measurement of nitrogenase activity in Trichodesmium sp. strain IMS101 over 24 h.

View larger version (34K): [in this window] [in a new window]

FIG. 3. Simulation of the successive use of different nitrogen sources in a continuous culture of Trichodesmium sp. strain IMS101 grown under a regime consisting of 12 h of light and 12 h of darkness. Shown are the concentration of DIN in the medium, the DIN uptake rate, and the N2 fixation rate. Horizontal black bars on the x axis show the dark period.

View larger version (42K): [in this window] [in a new window]

FIG. 4. L12 (left panels) and L16 (right panels) turbidostat experiments. Shown are daily dynamics when the model has reached the dynamic equilibrium, for simulations run with a regime consisting of 12 h of light and 12 h of darkness (L12) or 16 h of light and 8 h of darkness (L16). (a) Availability of carbon reserves, and synthesis and decay of nitrogenase activity. (b) Nitrogenase activity, rate of N2 fixation, and availability of nitrogen reserves. (c) Resultant average irradiance (Ieq) in the vessel and C/N ratio. Horizontal black bars represent the duration of the dark period.

View larger version (17K): [in this window] [in a new window]

FIG. 6. Combined influences of the average irradiance (Ieq) in the culture and the length of the light period at the dynamic equilibrium (series 2). Daily averages of variables at the equilibrium state are plotted as a function of the resultant daily average irradiance in the vessel. Lines represent changes in the daily average C/N ratio for a given light period: 8 h (L8), 12 h (L12), or 16 h (L16). Note that the x axis gives the irradiance (Ieq) averaged over 24 h.

Dynamics of nitrogenase activity.
Figure 2 depicts the daily pattern of nitrogenase activity in a culture of Trichodesmium sp. strain IMS101, measured by the online acetylene reduction technique. Nitrogenase activity was observed essentially during the light period of the light-dark cycle. The rate of N₂ fixation increased throughout the first half of the light period and reached its maximum 4 to 5 h after the onset of the light period. During the second half of the light period, activity decreased steadily, and activity quickly ceased when it became dark.

Transient dynamics of nitrogen concentration, rate of N₂ fixation, and consumption of nitrogen in a simulated experiment (L12 experiment).
The model behaved qualitatively in the same way whether it was run in the light-limited continuous-culture mode or in the turbidostat mode. This was true for the transient phase as well as for the equilibrium phase (data not shown). This was because the system was limited by light in both modes. The only difference was the resultant light intensity in the culture vessel. The results of the L12 continuous culture run are depicted in Fig. 3. The dilution experiment exhibited a progressive shift from growth at the expense of combined nitrogen to diazotrophic growth when the external combined nitrogen levels fell below the semisaturation concentration for nitrogen uptake. Nitrogen uptake and N₂ fixation occurred during the transition period. When available, combined nitrogen was preferred. Subsequently, when uptake decreased, N₂ fixation became progressively more important (Fig. 3). During this transient phase, both dissolved inorganic nitrogen (DIN) uptake and N₂ fixation were strongly stimulated when the light was switched on. A gradual and slight decrease in DIN uptake was then observed during the light phase (due to inhibition by the nitrogen reserves [N_r]), while N₂ fixation continued to increase due to the increasing potential activity of the nitrogenase. After the light was switched off, DIN uptake, which is dependent on energy from the catabolism of carbon reserves, dropped to half the value obtained during the light period, while N₂ fixation dropped to zero because nitrogenase broke down during the night. Several days were required for the population to deplete the DIN from the environment, because the inoculum was small.

Influence of the photoperiod on the daily pattern (L12 versus L16 run).
In the turbidostat as well as in the continuous-culture experiment, the model displayed a daily pattern of N₂ fixation with increasing activity during the first half of the light period and a decrease during the second half. The daily behavior at equilibrium is depicted comparatively for the two light regimes in Fig. 4. Nitrogenase synthesis was activated at the beginning of the day by the depletion of the nitrogen reserves (Fig. 4a and b). When energy became available, the rate of N₂ fixation increased, and so did the nitrogen reserves (Fig. 4b). Eventually, the cellular nitrogen status caused a decay of active nitrogenase in the middle of the light period, thus decreasing the N₂ fixation potential (Fig. 4a). N₂ fixation continued during the rest of the light period but decreased (i) because the "synthesis" of nitrogenase decreased as the internal nitrogen reserves became replete during the course of the day and (ii) due to the decay of nitrogenase. The decay of nitrogenase is proportional to the nitrogenase activity as well as to the shortage of available carbohydrates in cells (Fig. 4a). When the light was switched off, carbon (C_sm) was no longer supplied by photosynthesis and net consumption led to a decrease of the available C_sm reserves.

The resultant average irradiance in the vessel chosen in L16 was such that the total daily average light experienced by colonies was the same as that in the L12 turbidostat. Due to the lower average irradiance (I_eq) experienced by colonies in the L16 culture, a lower carbon-specific maximum rate of C_s synthesis was observed. However, photosynthesis occurred for a longer period, so that the daily net C_s synthesis rate was higher in L16. As a result, L16 showed higher total intracellular carbon and nitrogen contents (in structural biomass as well as in reserves) than L12 (data not shown).

Figure 4 shows that N₂ fixation activity occurred as long as the culture was exposed to light. Biomass-specific nitrogenase activity still showed a maximum approximately halfway into the light period, but this maximum value was lower than that in the L12 experiment. The daily N₂ fixation rate was higher in L16 than in L12. Cells in L16 accumulated more carbon and nitrogen than cells in L12. However, proportionally, more carbon was stored in C_r and C_sm reserves in L16, while more was allocated to N_r and C_s in L12. As a result, the C/N ratio observed in L16 was somewhat higher than that in L12.

Influence of the available irradiance on cell metabolism (series 1).
In the series 1 runs, the average irradiance (I_eq) in the culture at the equilibrium was set and the dilution rate progressively changed until the amount of biomass gave the desired value I_eq. Beyond a critical dilution rate, which equals the maximum growth rate, washout occurs. An unbalanced growth rate was achieved by changes in the cellular carbon and nitrogen contents, which were expressed as the relative proportions of the different molecules (Table 2). Thus, we define the relative proportion of structural carbon, pC_s, as its content relative to total intracellular carbon: pC_s = [C_s]/C_tot (Table 2). The relative content of the carbon reserve, pC_r, was calculated in the same way. The relative proportion of stored nitrogen, pN_r, is the content of N_r, expressed as carbon, versus total carbon content. In Fig. 5 the daily average of different variables at equilibrium is plotted as a function of the daily average photon flux density (). Because at dynamic equilibrium the biomass dilution equaled the net growth rate, a higher imposed irradiance provoked a higher growth rate and therefore higher dilution rates at equilibrium (Fig. 5a). With increasing irradiances, more carbon was fixed and stored in C_sm and C_r. Hence, the proportion of carbon storage (pC_sm and pC_r) increased, while the proportion of carbon in structural biomass (C_s) decreased (Fig. 5b and c). N₂ fixation is regulated by the availability of C_sm and therefore increased with light, as did the proportion of nitrogen stored (N_r) (Fig. 5c). Although growth decreased the cellular carbon content (because of losses due to carbon dissimilation), carbon assimilation increased faster than nitrogen assimilation. As a result, the C/N ratio increased with the average irradiance (Fig. 5a).

View larger version (23K): [in this window] [in a new window]

FIG. 5. Overview of 100 simulations similar to the L12 experiments (series 1), except that the dilution rate changes. The daily averages of variables at the equilibrium state are plotted as a function of the daily average irradiance in the vessel. (a) Average dilution rate and C/N ratio. (b) pCs and pCr. (c) pCsm and pNr. Rates are given per volume of medium.

Coinfluence of the average irradiance in the culture and the duration of the light period (series 2).
Figure 6 displays the fluctuations of the C/N ratio for different light periods. Three lines were drawn that represent, respectively, L8, L12, and L16. These curves show that higher C/N ratios resulted from the extra carbon assimilated during longer day lengths. This phenomenon is attributed to the saturation in the photosynthesis-light relationship. In this experiment, the total daily amount of light in the culture vessel (I_eq x day length) was equal in all runs. Hence, when the light period was short, light intensity was high. This resulted in a maximum rate of photosynthesis, but not all energy was used. In contrast, when the light period was long, incident irradiance was lower and photosynthesis was submaximal. Last, the relative increase in the C/N ratio with irradiance is slightly higher at long day lengths (Fig. 6). For instance, the C/N ratio at 34.8 µmol · m^–2 · s^–1 is 113.8% of the value at 3.7 µmol · m^–2 · s^–1 in L8, while it is 119% of that value in L16.

Top Abstract Introduction Materials and Methods Results Discussion References

 
Little is known about the regulatory mechanisms of N₂ fixation in Trichodesmium. Nitrogenase concentrations in populations are rarely quantified, and even when they are, the active and inactive forms of this enzyme complex are not easily distinguished from each other (12, 22, 54). Furthermore, the number of N₂-fixing cells in a trichome of Trichodesmium is mostly unknown. Therefore, it seems almost impossible to quantify the nitrogenase content or its potential activity in individual nitrogen-fixing cells. Due to these uncertainties, nitrogenase content is not a suitable state variable for the model. Therefore, we chose to model nitrogenase activity (potential and actual rates of N₂ fixation) rather than the enzyme content. N₂ fixation is an energy-costly process and therefore is not preferred when combined nitrogen is sufficiently available. As emphasized by Gallon (21), use of N₂ instead of DIN entails some extra costs, associated with (i) the synthesis of nitrogenase and (ii) the maintenance of anaerobic or microaerobic conditions to protect nitrogenase from inactivation by O₂. High respiration rates have been observed in Trichodesmium populations during periods of N₂ fixation (4, 8, 27, 42), and it has been suggested that additional O₂-consuming reactions are enhanced when N₂ fixation and photosynthesis occur simultaneously (2). In addition, combined nitrogen availability in the environment seems to elicit a down-regulation of the synthesis and activity of nitrogenase (40). A representation of intracellular reserves of both carbon and nitrogen is essential in order to accurately describe the dynamics of N₂ fixation and its regulation by the environment. To date, different physiologically based models already describe the uncoupling of carbon and nitrogen assimilation in the phytoplankton. In such unbalanced models, internal reserves consisting of either energy reserves (carbon) or intracellular nitrogen are included as an intermediate between photosynthesis and growth (15, 29, 50). These formulations have proved their accuracy in models related to cultures (23) as well as to natural environments (19, 30, 45, 52). Stephens et al. (49) devised a model that simulates the interaction between the different nitrogen sources in the nonheterocystous diazotrophic cyanobacterium Gloeothece. In this model, uncoupling of carbon and nitrogen assimilation was achieved through the description of the intracellular concentration of a nitrogen reserve (glutamine) and included the description of nitrogenase activity. Thus, in the Gloeothece model, the changes in the C/N ratio were driven by fluctuations of the glutamine pool (nitrogen reserve), while the accumulation of storage carbon (and hence energy reserves) was not explicitly included. Therefore, the representation of the interactions led to a complex set of equations. In order to be able to more easily interpret the analysis of the model behavior and to consider N₂ fixation from a mechanistic point of view, we have chosen to use a simpler formulation. We included decomposition of cellular carbon components as proposed by Lancelot and Billen (29) in which an intermediary pool of low-molecular-weight carbon molecules (C_sm), the product of photosynthesis, is described in addition to the internal reserves of carbon. The C_sm pool is an important component of the model, since its rapid turnover allows fast fluctuations of the cellular processes. Added to that description, an internal storage of nitrogen reserves was included to account for N₂ fixation. Hence, the simulated mechanisms include a down-regulation of active nitrogenase by available combined nitrogen in the environment, while the actual fixation rate is further regulated by the internal nitrogen and carbon status. Carbon is allocated differently according to fluctuations in the environment. This model has a flexible structure, resulting in a dynamic response of cells to external forcing. Since both photosynthesis and N₂ fixation depend on light and on the internal status of cells, it is possible to analyze whether differences in the light patterns can result in different carbon acquisition and growth.

Diurnal dynamics of N₂ fixation in Trichodesmium.
Because Trichodesmium fixes N₂ during the day, the energy required for diazotrophy is provided mainly by photosynthesis. Such a relation has been clearly pointed out for heterocystous cyanobacteria (46) and for Trichodesmium (26, 48). Trichodesmium showed a rapid decline in the N₂ fixation rate when photosynthesis ceased, suggesting that an immediate shortage of energy occurred after darkness (Fig. 2). Hence, our results confirm the dependence of Trichodesmium sp. strain IMS101 on ATP provided by photosynthesis, highlighting the importance of the different substrates required for N₂ fixation (N₂, ATP, and reducing equivalents). In the model, the lowering of the available C_sm reserves in cells at the end of the day resulted in an increased decay of active nitrogenase during the dark (Fig. 4b). The actual N₂ fixation rate drops to zero at the end of the day while nitrogenase is still active, indicating that activity is limited by the shortage of carbohydrates in the form of C_sm. N₂ fixation, driven by the "small-metabolite" pool (C_sm) in the model, shows a fast response to light. The role of this compartment is to take into account fast exchanges of carbon to represent, for instance, the steep decline in N₂ fixation when the light is turned off. However, this decline is less steep than that measured in the laboratory culture. Furthermore, C_sm content originates both from photosynthesis and from storage/catabolism. When cells are exposed to the dark, C_sm is fuelled not by photosynthesis but only by catabolism of cellular carbohydrates. Provided that the C_sm reserve is available to meet the energy demand, the model also shows that diazotrophy in Trichodesmium can occur during part of the night (Fig. 4b). These results disagree with observations by Mulholland and Bernhardt (37), who did not observe any night N₂ fixation in cultures of Trichodesmium sp. strain IMS101. On the other hand, our results are in accordance with the behavior of cultured populations of Trichodesmium sp. strain GBRTRLI101 (20) and with studies of natural Trichodesmium populations from the Pacific (32), supporting the idea that nitrogenase activity is not strictly restricted to the light period in Trichodesmium.

Although photosynthesis and N₂ fixation both occur during the day, unbalanced growth leads to a fluctuation of the C/N ratio in cells at a diel scale (Fig. 4c). The maximum hourly N₂ fixation rate [expressed in mol of N · (g of Chla)^–1 · h^–1] and carbon fixation rate [expressed in mol of C · (g of Chla)^–1 · h^–1] calculated by the model were, respectively, 0.043 and 0.24 in L12 and 0.036 and 0.19 in L16 cultures. These values are in the range reported for different studies on Trichodesmium, listed by Mulholland and Bernhardt (37). The daily averages of the molar C/N₂ fixation ratio were 22 in the L12 culture (where the daily average of the dilution rate is 0.28 day^–1) and 18.7 in the L16 culture (where the daily average of the dilution rate is 0.29 day^–1), consistent with results from acetylene reduction assays reported from continuous-culture experiments with Trichodesmium (37).

Influence of the light rhythm.
Fluctuations of the C/N ratio as a function of light indicate that both the irradiance and the day length affect N₂ fixation and cellular elemental composition in Trichodesmium (Fig. 6). L12 and L16 simulations received the same total daily light dose (I_eq x day length) but experienced different light regimes. When the day length increased, nitrogenase was active in cells for a longer period during the day, and although it was progressively down-regulated by the increasing nitrogen reserves, activity proceeded toward the end of the light period because carbon was still available to fuel N₂ fixation. Light energy is more efficiently used by photosynthesis at long day lengths: more carbon is fixed per cell and C_s production is higher at the daily scale (and so is pC_r), but cells proportionately allocate less nitrogen to their structural biomass. Hence, at shorter day lengths, photosynthesis was better synchronized to nitrogenase requirements and resulted in a lower C/N ratio. In contrast, with longer light periods, an imbalance occurred between photosynthesis and biomass production, and cells redirected the extra carbon fixed into storage. The higher diurnal accumulation of carbon in reserves in L16 also indicates that much less energy was used by the N₂ fixation process at the beginning and end of the period of activity. Because nitrogenase activity is not constant during the day, our model simulations emphasize the importance of energy availability when nitrogenase approaches its maximum potential activity, i.e., mainly before midday. The temporal light patterns may thus be part of the factors that influence the proliferation of Trichodesmium populations. The model suggests that cell C/N stoichiometry is tightly regulated, and the range of C/N ratios calculated by the model is in agreement with field observations of natural Trichodesmium populations (3).

Regulation in the modeled system versus laboratory culture.
With the model, we could generate a daily N₂ fixation pattern qualitatively similar to the daily patterns measured by us in a laboratory culture of Trichodesmium (Fig. 2) and patterns known from natural populations and other cultures (39).

A special feature of Trichodesmium is that it fixes N₂ during the day, concomitantly with photosynthesis. N₂ fixation rates are related to photosynthetic activity on a short time scale. Incident irradiance was kept constant during the light period in both the laboratory culture and the model simulations. The observed diel pattern of N₂ fixation is not correlated with incident irradiance, suggesting that the process of N₂ fixation is not linearly related to light. Rather, the daily pattern had modulations, indicating the occurrence of internal regulatory mechanisms that result from cell history. In the laboratory culture, the increase in the rate of N₂ fixation observed in the morning can be due either to (i) an increase in the fixation potential within N₂-fixing cells, (ii) an increase in the number of cells that actually fix N₂, or (iii) an increase in availability of the substrate (carbon) limiting the N₂ fixation rate. The first two possibilities agree with the modeled increase in nitrogenase activity (whether in one cell or in several cells), whereas the third is the modulation of the N₂ fixation rate by carbon availability. In the same way, declining rates in the afternoon follow deactivation of the enzyme but could also reflect a decrease in numbers of N₂-fixing cells.

Because N₂ fixation in the model depends only on physiological factors, enzymatic activities are not confined to the light period. Rather, we assume that diazotrophic activity is determined by nitrogen limitation and cell C/N stoichiometry. This assumption is based on two main ideas: (i) the synthesis and decay of the enzyme are activated and modulated by the internal nutrient and energy states of the cells and (ii) the potential fixation rate is enhanced by nitrogen starvation but limited by electrons originating from the small carbon metabolites or ATP produced by respiration or photosynthesis. As a result, the onset of enzyme activity starts before the onset of the light, and the activity varies along the light period although the applied irradiance is constant. Daily fluctuations calculated by the model proved to be in good agreement with the dynamics of N₂ fixation observed in Trichodesmium populations (3), as well as with our own measurements (Fig. 2). Hitherto, models available in the literature have displayed patterns of nitrogenase activity that were most often directly related to light intensity (26, 33, 34). However the synthesis and activity of nitrogenase are probably under the control of a circadian clock (12, 43). In Synechococcus sp. strain Rf-1, the cyclic nitrogenase activity resulted from de novo synthesis of nitrogenase during the dark (13). In Trichodesmium, part of the enzymatic complex is destroyed during the afternoon and synthesized de novo toward the end of the night and the beginning of the day through genetic regulatory mechanisms (5, 53, 54). This transcriptional regulation cannot be described explicitly within our model, because it occurs at different scales than the processes we have focused on. However, regulation is still exerted by the metabolism in the model, which integrates processes at lower levels. In particular, the gene transcriptional level is implicitly involved in the calculation of the potential activity of nitrogenase through the feedback regulation exerted by the internal nitrogen status and a presumed regulatory mechanism based on the C_sm and N_r stoichiometry (Fig. 4a and 5a). The model is then able to reproduce the daily patterns of the metabolism that result from these coupled influences of the environment and the internal status of cells.

Finally, the model shows that carbon availability is a simple and efficient mechanism for regulating the balance of carbon and nitrogen fixed (C/N ratio) in cells with a spatial separation of N₂ fixation. We conclude that N₂ fixation is ultimately determined by substrate availability and that this direct control prevents high fluctuations on a daily basis.

Conclusions.
The dynamics of N₂ fixation are complex because they involve the influence of external as well as internal factors. The mathematical model presented here aimed to provide a scientific reflection tool conceptualizing the ecological determinants of the N₂ fixation process in Trichodesmium. The combined representation of intracellular reserves of carbon and nitrogen and the description of the activity of nitrogenase are essential parts of the model formulation. The modeling is unconventional in that it precedes observations that could be used to test, validate, or falsify its assumptions. However, where available, the model simulation results agree qualitatively with observations. The model generates a correct timing of the onset of enzyme activity synthesis. Moreover, it also proved to represent the co-occurrence and transition between nitrogen uptake and N₂ fixation, and we were able to reproduce a daily pattern of nitrogenase activity as observed in nature. This model is suitable for exploring the influence of environmental factors on the metabolism of this important cyanobacterium and the conditions promoting N₂ fixation. In particular, we emphasize the strong relation between N₂ fixation and carbon availability. The C/N ratio appeared to be dependent on the light regime. Fluctuations in cell C/N stoichiometry resulted from changes in carbon uptake, assimilation into the different cellular compartments, and energy consumption, according to growth requirements. With the simple mechanism of carbon-limited N₂ fixation, the C/N ratio was balanced over days in the model. This modeling exercise thus emphasizes the need for acquisition of data on nitrogenase activity under various environmental conditions, as it may be substrate (carbon) limited. Substrate limitation can directly control C/N ratios and could thus precede the delayed genetic regulatory mechanisms based on the cellular C/N status.

 
We are very grateful to J. Zehr for his contructive comments on the manuscript.

This work was funded by a grant from the European Commission, contract MEIF-CT-2003-500516.

 
* Corresponding author. Present address: Ocean Sciences Department, University of California—Santa Cruz, 1156 High Street, Santa Cruz, CA 95064. Phone: (831) 459-5152. Fax: (831) 459-4882. E-mail: srabouille{at}pmc.ucsc.edu.

Top Abstract Introduction Materials and Methods Results Discussion References

 

1
1. Bergman, B., and E. J. Carpenter. 1991. Nitrogenase confined to randomly distributed trichomes in the marine cyanobacterium Trichodesmium thiebautii. J. Phycol. 27:158-165.[CrossRef]
2
2. Bergman, B., J. R. Gallon, A. N. Rai, and L. J. Stal. 1997. N₂ fixation by non-heterocystous cyanobacteria. FEMS Microbiol. Rev. 19:139-185.[CrossRef]
3
3. Berman-Frank, I., J. T. Cullen, Y. Shaked, R. M. Sherrell, and P. G. Falkowski. 2001. Iron availability, cellular iron quotas, and nitrogen fixation in Trichodesmium. Limnol. Oceanogr. 46:1249-1260.
4
4. Berman-Frank, I., P. Lundgren, Y.-B. Chen, H. Kupper, Z. Kolber, B. Bergman, and P. Falkowski. 2001. Segregation of nitrogen fixation and oxygenic photosynthesis in the marine cyanobacterium Trichodesmium. Science 294:1534-1537.[Abstract/Free Full Text]
5
5. Capone, D. G., J. M. O'Neil, J. Zehr, and E. J. Carpenter. 1990. Basis for diel variation in nitrogenase activity in the marine planktonic cyanobacterium Trichodesmium thiebautii. Appl. Environ. Microbiol. 56:3532-3536.[Abstract/Free Full Text]
6
6. Capone, D. G., J. P. Zehr, H. W. Paerl, B. Bergman, and E. J. Carpenter. 1997. Trichodesmium, a globally significant marine cyanobacterium. Science 276:1221-1229.[Abstract/Free Full Text]
7
7. Carpenter, E. J., B. Bergman, R. Dawson, P. J. A. Siddiqui, E. Soderback, and D. G. Capone. 1992. Glutamine-synthetase and nitrogen cycling in colonies of the marine diazotrophic cyanobacteria Trichodesmium spp. Appl. Environ. Microbiol. 58:3122-3129.[Abstract/Free Full Text]
8
8. Carpenter, E. J., and T. Roenneberg. 1995. The marine planktonic cyanobacteria Trichodesmium spp.—photosynthetic rate measurements in the SW Atlantic Ocean. Mar. Ecol. Prog. Ser. 118:267-273.
9
9. Carpenter, E. J., and K. Romans. 1991. Major role of the cyanobacterium Trichodesmium in nutrient cycling in the North Atlantic Ocean. Science 254:1356-1358.[Abstract/Free Full Text]
10
10. Carpenter, E. J., A. Subramaniam, and D. G. Capone. 2004. Biomass and primary productivity of the cyanobacterium Trichodesmium spp. in the tropical N Atlantic Ocean. Deep-Sea Res. Pt. I 51:173-203.
11
11. Cavender-Bares, K. K., D. M. Karl, and S. W. Chisholm. 2001. Nutrient gradients in the western North Atlantic Ocean: relationship to microbial community structure and comparison to patterns in the Pacific Ocean. Deep-Sea Res. Pt. I 48:2373-2395.
12
12. Chen, Y. B., J. P. Zehr, and M. Mellon. 1996. Growth and nitrogen fixation of the diazotrophic filamentous nonheterocystous cyanobacterium Trichodesmium sp. IMS 101 in defined media: evidence for a circadian rhythm. J. Phycol. 32:916-923.[CrossRef]
13
13. Chow, T. J., and F. R. Tabita. 1994. Reciprocal light-dark transcriptional control of Nif and Rbc expression and light-dependent posttranslational control of nitrogenase activity in Synechococcus sp. strain Rf-1. J. Bacteriol. 176:6281-6285.[Abstract/Free Full Text]
14
14. Coles, V. J., R. R. Hood, M. Pascual, and D. G. Capone. 2004. Modeling the impact of Trichodesmium and nitrogen fixation in the Atlantic Ocean. J. Geophys. Res. Oceans 109:C06007. doi:10.1029/2002JC001754.
15
15. Droop, M. R. 1973. Some thoughts on nutrient limitation in algae. J. Phycol. 9:264-272.[CrossRef]
16
16. Falkowski, P. G. 1997. Evolution of the nitrogen cycle and its influence on the biological sequestration of CO₂ in the ocean. Nature 387:272-275.
17
17. Faugeras, B., M. Levy, L. Memery, J. Verron, J. Blum, and I. Charpentier. 2003. Can biogeochemical fluxes be recovered from nitrate and chlorophyll data? A case study assimilating data in the Northwestern Mediterranean Sea at the JGOFS-DYFAMED station. J. Mar. Syst. 40:99-125.[CrossRef]
18
18. Fennel, K., Y. H. Spitz, R. M. Letelier, M. R. Abbott, and D. M. Karl. 2002. A deterministic model for N₂ fixation at stn. ALOHA in the subtropical North Pacific Ocean. Deep-Sea Res. Pt. II 49:149-174.
19
19. Flynn, K. J., and M. J. R. Fasham. 2003. Operation of light-dark cycles within simple ecosystem models of primary production and the consequences of using phytoplankton models with different abilities to assimilate N in darkness. J. Plankt. Res. 25:83-92.
20
20. Fu, F. X., and P. R. F. Bell. 2003. Factors affecting N₂ fixation by the cyanobacterium Trichodesmium sp. GBR-TRLI101. FEMS Microbiol. Ecol. 45:203-209.[CrossRef]
21
21. Gallon, J. R. 2001. N₂ fixation in phototrophs: adaptation to a specialized way of life. Plant Soil 230:39-48.[CrossRef]
22
22. Gallon, J. R., J. J. Cheng, L. J. Dougherty, V. A. Gallon, H. Hilz, D. M. Pederson, H. M. Richards, S. Ruggeberg, and C. J. Smith. 2000. A novel covalent modification of nitrogenase in a cyanobacterium. FEBS Lett. 468:231-233.[CrossRef][Medline]
23
23. Geider, R. J., H. L. MacIntyre, and T. M. Kana. 1998. A dynamic regulatory model of phytoplanktonic acclimation to light, nutrients, and temperature. Limnol. Oceanogr. 43:679-694.
24
24. Holl, C. M., and J. P. Montoya. 2005. Interactions between nitrate uptake and nitrogen fixation in continuous cultures of the marine diazotroph Trichodesmium (Cyanobacteria). J. Phycol. 41:1178-1183.[CrossRef]
25
25. Hood, R. R., N. R. Bates, D. G. Capone, and D. B. Olson. 2001. Modeling the effect of nitrogen fixation on carbon and nitrogen fluxes at BATS. Deep-Sea Res. Pt. II 48:1609-1648.
26
26. Hood, R. R., A. Subramaniam, L. R. May, E. J. Carpenter, and D. G. Capone. 2002. Remote estimation of nitrogen fixation by Trichodesmium. Deep-Sea Res. Pt. II 49:123-147.
27
27. Kana, T. M. 1993. Rapid oxygen cycling in Trichodesmium thiebautii. Limnol. Oceanogr. 38:18-24.
28
28. Karl, D., A. Michaels, B. Bergman, D. Capone, E. Carpenter, R. Letelier, F. Lipschultz, H. Paerl, D. Sigman, and L. Stal. 2002. Dinitrogen fixation in the world's oceans. Biogeochemistry 57/58:47-98.[CrossRef]
29
29. Lancelot, C., and G. Billen. 1985. Carbon-nitrogen relationships in nutrient metabolism of coastal marine ecosystems. Adv. Aquat. Microbiol. 3:263-321.
30
30. Lancelot, C., G. Billen, C. Veth, S. Becquevort, and S. Mathot. 1991. Modeling carbon cycling through phytoplankton and microbes in the Scotia-Weddell Sea area during sea ice retreat. Mar. Chem. 35:305-324.
31
31. Lee, K., D. M. Karl, R. Wanninkhof, and J.-Z. Zhang. 2002. Global estimates of net carbon production in the nitrate-depleted tropical and subtropical oceans. Geophys. Res. Lett. 29:1907. doi:10.1029/2001GL014198.[CrossRef]
32
32. Letelier, R. M., and D. M. Karl. 1996. Role of Trichodesmium spp. in the productivity of the subtropical North Pacific Ocean. Mar. Ecol. Prog. Ser. 133:263-273.
33
33. Levine, S. N., and W. M. Lewis, Jr. 1984. Diel variation of nitrogen fixation in Lake Valencia, Venezuela. Limnol. Oceanogr. 29:887-893.
34
34. Levine, S. N., and W. M. Lewis, Jr. 1987. A numerical model of nitrogen fixation and its application to Lake Valencia. Freshwat. Biol. 17:265-274.[CrossRef]
35
35. Moore, J. K., S. C. Doney, D. M. Glover, and I. Y. Fung. 2002. Iron cycling and nutrient-limitation patterns in surface waters of the World Ocean. Deep-Sea Res. Pt. II 49:463-507.
36
36. Moore, J. K., S. C. Doney, J. A. Kleypas, D. M. Glover, and I. Y. Fung. 2002. An intermediate complexity marine ecosystem model for the global domain. Deep-Sea Res. Pt. II 49:403-462.
37
37. Mulholland, M. R., and P. W. Bernhardt. 2005. The effect of growth rate, phosphorus concentration, and temperature on N₂ fixation, carbon fixation, and nitrogen release in continuous cultures of Trichodesmium IMS101. Limnol. Oceanogr. 50:839-849.
38
38. Mulholland, M. R., and D. G. Capone. 1999. Nitrogen fixation, uptake and metabolism in natural and cultured populations of Trichodesmium spp. Mar. Ecol. Prog. Ser. 188:33-49.
39
39. Mulholland, M. R., and D. G. Capone. 2000. The nitrogen physiology of the marine N₂-fixing cyanobacteria Trichodesmium spp. Trends Plant Sci. 5:148-153.[CrossRef][Medline]
40
40. Mulholland, M. R., K. Ohki, and D. G. Capone. 2001. Nutrient controls on nitrogen uptake and metabolism by natural populations and cultures of Trichodesmium (Cyanobacteria). J. Phycol. 37:1001-1009.[CrossRef]
41
41. Platt, T., C. L. Gallegos, and W. G. Harrison. 1980. Photoinhibition of photosynthesis in natural assemblages of marine phytoplankton. J. Mar. Res. 38:687-701.
42
42. Roenneberg, T., and E. J. Carpenter. 1993. Daily rhythm of O₂-evolution in the cyanobacterium Trichodesmium thiebautii under natural and constant conditions. Mar. Biol. 117:693-697.[CrossRef]
43
43. Schneegurt, M. A., D. L. Tucker, J. K. Ondr, D. M. Sherman, and L. A. Sherman. 2000. Metabolic rhythms of a diazotrophic cyanobacterium, Cyanothece sp. strain ATCC 51142, heterotrophically grown in continuous dark. J. Phycol. 36:107-117.[CrossRef]
44
44. Soetaert, K., V. deClippele, and P. Herman. 2002. FEMME, a flexible environment for mathematically modelling the environment. Ecol. Model. 151:177-193.[CrossRef]
45
45. Soetaert, K., P. M. J. Herman, J. J. Middelburg, C. Heip, C. L. Smith, P. Tett, and K. Wild-Allen. 2001. Numerical modelling of the shelf break ecosystem: reproducing benthic and pelagic measurements. Deep-Sea Res. Pt. II 48:3141-3177.
46
46. Staal, M., S. T. L. Hekkert, P. Herman, and L. J. Stal. 2002. Comparison of models describing light dependence of N₂ fixation in heterocystous cyanobacteria. Appl. Environ. Microbiol. 68:4679-4683.[Abstract/Free Full Text]
47
47. Staal, M., S. T. Lintel-Hekkert, F. Harren, and L. Stal. 2001. Nitrogenase activity in cyanobacteria measured by the acetylene reduction assay: a comparison between batch incubation and on-line monitoring. Environ. Microbiol. 3:343-351.[CrossRef][Medline]
48
48. Staal, M., F. J. R. Meysman, and L. J. Stal. 2003. Temperature excludes N₂-fixing heterocystous cyanobacteria in the tropical oceans. Nature 425:504-507.[CrossRef][Medline]
49
49. Stephens, N., K. J. Flynn, and J. R. Gallon. 2003. Interrelationships between the pathways of inorganic nitrogen assimilation in the cyanobacterium Gloeothece can be described using a mechanistic mathematical model. New Phytol. 160:545-555.[CrossRef]
50
50. Tett, P. 1998. Parameterising a microplankton model. Napier University, Edinburgh, United Kingdom.
51
51. Tyrrell, T. 1999. The relative influences of nitrogen and phosphorus on oceanic primary production. Nature 400:525-531.[CrossRef]
52
52. Van den Meersche, K., J. J. Middelburg, K. Soetaert, P. van Rijswijk, H. T. S. Boschker, and C. H. R. Heip. 2004. Carbon-nitrogen coupling and algal-bacterial interactions during an experimental bloom: modeling a C-13 tracer experiment. Limnol. Oceanogr. 49:862-878.
53
53. Wyman, M., J. P. Zehr, and D. G. Capone. 1996. Temporal variability in nitrogenase gene expression in natural populations of the marine cyanobacterium Trichodesmium thiebautii. Appl. Environ. Microbiol. 62:1073-1075.[Abstract]
54
54. Zehr, J. P., M. Wyman, V. Miller, L. Duguay, and D. G. Capone. 1993. Modification of the Fe protein of nitrogenase in natural populations of Trichodesmium thiebautii. Appl. Environ. Microbiol. 59:669-676.[Abstract/Free Full Text]

---

Applied and Environmental Microbiology, May 2006, p. 3217-3227, Vol. 72, No. 5
0099-2240/06/$08.00+0     doi:10.1128/AEM.72.5.3217-3227.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.

This article has been cited by other articles:

• Finzi-Hart, J. A., Pett-Ridge, J., Weber, P. K., Popa, R., Fallon, S. J., Gunderson, T., Hutcheon, I. D., Nealson, K. H., Capone, D. G. (2009). Fixation and fate of C and N in the cyanobacterium Trichodesmium using nanometer-scale secondary ion mass spectrometry. Proc. Natl. Acad. Sci. USA 106: 6345-6350 [Abstract] [Full Text]  
• Toepel, J., Welsh, E., Summerfield, T. C., Pakrasi, H. B., Sherman, L. A. (2008). Differential Transcriptional Analysis of the Cyanobacterium Cyanothece sp. Strain ATCC 51142 during Light-Dark and Continuous-Light Growth. J. Bacteriol. 190: 3904-3913 [Abstract] [Full Text]  

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Rabouille, S. Articles by Soetaert, K. Search for Related Content PubMed PubMed Citation Articles by Rabouille, S. Articles by Soetaert, K. Agricola Articles by Rabouille, S. Articles by Soetaert, K.