ABSTRACT
The cyanobacterium Trichodesmium is an important contributor of new nitrogen (N) to the surface ocean, but its strategies for protecting the nitrogenase enzyme from inhibition by oxygen (O2) remain poorly understood. We present a dynamic physiological model to evaluate hypothesized conditions that would allow Trichodesmium to carry out its two conflicting metabolic processes of N2 fixation and photosynthesis. First, the model indicates that managing cellular O2 to permit N2 fixation requires high rates of respiratory O2 consumption. The energetic cost amounts to ∼80% of daily C fixation, comparable to the observed diminution of the growth rate of Trichodesmium relative to other phytoplankton. Second, by forming a trichome of connected cells, Trichodesmium can segregate N2 fixation from photosynthesis. The transfer of stored C to N-fixing cells fuels the respiratory O2 consumption that protects nitrogenase, while the reciprocal transfer of newly fixed N to C-fixing cells supports cellular growth. Third, despite Trichodesmium lacking the structural barrier found in heterocystous species, the model predicts low diffusivity of cell membranes, a function that may be explained by the presence of Gram-negative membrane, production of extracellular polysaccharide substances (EPS), and “buffer cells” that intervene between N2-fixing and photosynthetic cells. Our results suggest that all three factors—respiratory protection, trichome formation, and diffusion barriers—represent essential strategies that, despite their energetic costs, facilitate the growth of Trichodesmium in the oligotrophic aerobic ocean and permit it to be a major source of new reactive nitrogen.
IMPORTANCE Trichodesmium is a major nitrogen-fixing cyanobacterium and exerts a significant influence on the oceanic nitrogen cycle. It is also a widely used model organism in laboratory studies. Since the nitrogen-fixing enzyme nitrogenase is extremely sensitive to oxygen, how these surface-dwelling plankton manage the two conflicting processes of nitrogen fixation and photosynthesis has been a long-standing question. In this study, we developed a simple model of metabolic fluxes of Trichodesmium capturing observed daily cycles of photosynthesis, nitrogen fixation, and boundary layer oxygen concentrations. The model suggests that forming a chain of cells for spatially segregating nitrogen fixation and photosynthesis is essential but not sufficient. It also requires a barrier against oxygen diffusion and high rates of oxygen scavenging by respiration. Finally, the model indicates that the life span of intracellular oxygen is extremely short, thus enabling cells to instantly create a low-oxygen environment upon deactivation of photosynthesis.
INTRODUCTION
Biological dinitrogen (N2) fixation provides bioavailable nitrogen (N) to the marine biosphere, supporting up to half of net community production in otherwise nutrient-depleted environments (1). The process of N2 fixation by the enzyme nitrogenase requires large amounts of energy and electrons (2–4). Furthermore, the nitrogenase enzyme contains metal cofactors that are irreversibly disabled in the presence of even trace levels of oxygen (O2) (5, 6). The mechanisms by which cells of a few micrometers in size maintain an active nitrogenase enzyme in an O2-rich environment are diverse and not fully understood (7). Some nitrogen fixers form a thick glycolipid layer of specialized cells (heterocysts) that prevent O2 diffusion into the N2-fixing cells (8). Other nitrogen fixers maintain high respiration rates to counteract the passive O2 diffusion (9–12). The metabolic strategies that enable N2 fixation to function in an oxygenated environment occur at the expense of other physiological activities, including growth. The growth rate handicap of diazotrophs is considered a key ecological trade-off (13–17) with important implications for the global N cycle (18, 19).
In the marine environment, a major contributor to N2 fixation is the photosynthetic diazotroph Trichodesmium, mainly observed in oligotrophic tropical and subtropical oceans (20–23). Although this species forms trichomes, they do not contain heterocysts to protect from O2 invasion (24, 25). Moreover, they are observed to fix N2 during the day, when photosynthetic production of O2 is also occurring (26, 27). Some studies show Trichodesmium respiration rates exceeding those of non-N2-fixing cyanobacteria (10, 26), despite a positive net O2 evolution rate during the daytime (25). Despite over a century of research on Trichodesmium, there is no unequivocal explanation for how N2 fixation occurs when the cells are photosynthetically active and O2 should be at its highest levels.
In response to this physiological enigma, it has been hypothesized that N2 fixation and photosynthesis are temporally and/or spatially segregated (28, 29). Spatial segregation is a highly debated strategy, as some previous work revealed the presence of nitrogenase in almost all cells (29, 30), while other reports showed nitrogenase occurred in about 10% to 20% of the cells (24, 31, 32). In support of temporal segregation, it has been shown that the rates of N2 fixation and respiration increase, while the rate of photosynthesis decreases, during the middle of the light period (26). Whether this temporal segregation is sufficient for photosynthesis and N2 fixation to occur simultaneously remains unclear (25). A recent approach to the Trichodesmium paradox has been to track 13C and 15N uptake at the cellular level using high-resolution secondary ion mass spectrometry (33, 34). However, even with near-hourly resolution measurements, it was not possible to determine spatial segregation along the trichome because the redistribution of newly fixed N occurs on a time scale of minutes (33).
Compiling previous studies reveals a common general feature of diurnal cycles in Trichodesmium physiology (Fig. 1). Rates of photosynthesis increase at sunrise and peak in early morning (Fig. 1A). The rate decreases during midday, increases slightly again toward evening, and decreases to nearly zero at night. The rate of N2 fixation, on the other hand, reaches its maximum value during the midday and its minimum (approximately zero) during the night (Fig. 1B). Similarly to N2 fixation, respiration rates peak during the midday, yielding a dip in near-cell O2 within a colony. Whereas some studies show reduced O2 within a colony of Trichodesmium (28), recent work shows that photosynthesis causes the boundary layer to have O2 levels that are 20% to 30% higher than those seen in the ambient water, where O2 is nearly saturated (34). Therefore, Trichodesmium needs to manage O2 fluxes not only directly from the photosynthetic cells (P cells) but also from the boundary layer environment within a colony. The question remains as to how these observed trends relate to the temporal coexistence of photosynthesis and N2 fixation.
In this study, we take a fresh approach to investigation of the Trichodesmium paradox by modeling the physiology of Trichodesmium (cell flux model of Trichodesmium [CFM-Tricho]) over a diurnal cycle to evaluate the hypothesized spatial and temporal strategies of the cyanobacteria for maintaining N2 fixation and photosynthesis. It has been suggested that colony formation by Trichodesmium plays an important role in creating a low-O2 environment to facilitate N2 fixation (28, 35), though some evidence contradicts assertions of such a role (34, 36). Recent studies show that the majority of Trichodesmium exists as filaments (37), which have higher rates of N2 fixation per cell than colony-maintained Trichodesmium cells (38). We focused the model on a single trichome due to its simple morphology. However, the considerations are equally valid for colonies, as the model resolves the near-cell environment (referred to here as the “boundary layer”) where O2 concentrations are influenced by cellular metabolism. The model simulates cellular resource allocation by combining a model of the cellular reserves of C and N (39) with a representation of O2 management critical to nitrogenase activity (12). The primary mechanisms of O2 protection include high respiration rates (respiratory protection [including dark respiration and light-dependent respiration]), segregation of N2-fixing cells from photosynthetic cells (trichome formation), and low diffusivity between cells (diffusion barriers). The diurnal variation of metabolism can be explained by fluctuations in the relative abundances of photosynthetic cells and nonphotosynthetic cells. The quantitative model requirements are evaluated against our current knowledge of Trichodesmium and other diazotrophs.
RESULTS
Simulating cellular differentiation.The model resolves two types of cells: photosynthetic cells (P cells) and nonphotosynthetic cells (N cells) (Fig. 2). The fractions of P and N cells are represented by fP and fN, respectively. The P cells fix carbon and make it available for growth, storage, and respiration. The N cells use stored C obtained from P cells for O2 consumption and N2 fixation. The proportion of cells carrying out each metabolic function determines the rates within and fluxes from the trichome as a whole. We are interested in the rates of N2 fixation (FNfix; here, “F” indicates fluxes), which we assume depend on O2 and the standing stocks of stored N and C, as represented in equation 1:
Schematic depiction of molecular pools and fluxes in the model. The model consists of the following four boxes: photosynthetic cells (P), nonphotosynthetic cells (N), boundary layer (B), and environment (E). In each cell type, the model solves for the concentration of stored C and N and the concentration of O2. The model assumes that the supply of CO2 and N2 does not limit the rates of photosynthesis and N2 fixation, the assumption made in most ecological models. Each flux symbol corresponds to those used in the model equations (see Materials and Methods).
Here,
To resolve the variation in O2 concentrations outside the cells due to photosynthesis and respiration, the O2 balance in the boundary layer environment is included. The model normalizes fluxes and molecules to the volume of the entire trichome. Therefore, the model can represent any number of cells with a certain proportion of N2-fixing cells. The model assumes that the supply of CO2 and N2 does not limit the rate of photosynthesis and N2 fixation, a common assumption of most diazotroph models. Instead, photosynthesis is a function of light and stored C, and N2 fixation is limited by stored C and N (equation 1). Further details are provided in Materials and Methods and Text S1 in the supplemental material.
TEXT S1
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
We run the model under a 12-h/12-h light/dark cycle, with fP prescribed as a step function in time (Fig. 3A). These transitions occur smoothly in nature and in experimental observations (26, 33, 34). Because our goal is to elucidate mechanisms, rather than to simulate precise details of particular experiments, we choose abrupt transitions that can be clearly discerned in the model output. As a check on the broad applicability of our diurnal forcing, we compared the trend of averaged Fv/Fm (photosynthetic quantum yield of photosystem II) to previous observations (26) (Fig. 3B). Fv/Fm is an indication on the efficiency of light use, which tends to be lower when cells are actively fixing nitrogen (26, 40). Also, it indicates the activity of photosystem II and thus the production of O2 (26, 41). The general trend in observations is captured by Fv/Fm = 0.5 and 0.1 for P cells and N cells, respectively. Similar values were observed during photosynthetic and nonphotosynthetic periods in Crocosphaera (40).
Modeled diurnal cycle of (A) the fraction of photosynthetic cells (fP) and (B) quantum yield (Fv/Fm). In panel B, model results (blue curves) are compared with observed data (red circles) (26). Light period, 0 to 12 h; dark period, 12 to 24 h (indicated by gray shading).
Mechanisms of O2 management.The maintenance of low intracellular O2 levels to permit N2 fixation can be achieved through several potential mechanisms, including the following: (i) trichome formation, (ii) respiratory protection, and (iii) diffusion barriers. All of these factors need to be included in order to reproduce the observed diurnal variation in metabolic rates (Fig. 1). Here we describe model results that evaluate the importance of each mechanism.
Trichome formation.Due to the production of O2, N2 fixation cannot occur in the photosynthetic cells. Thus, trichome formation with differently functioning cells is a key factor for temporal coexistence of N2 fixation and photosynthesis. We reproduce the observed daily cycle of photosynthesis and N2 fixation (Fig. 4) by simulating the daily cycle of the fractions of P cells and N cells (Fig. 3A).
Rates of photosynthesis and N2 fixation for a simulated trichome. (A) C-based rates of photosynthesis. Data represent C fixation rates normalized from results reported in references 33 (Data1) and 26 (Data2 and Data3). Data3 results represent conversion from gross O2 evolution data. (B) N-based rates of nitrogen fixation. Data represent N2 fixation rates normalized from results reported in references 33 (Data1) and 26 (Data2 and Data3). For Data1, only daytime data are plotted. The normalized data are plotted using magnitudes similar to those of the model outputs. Model fluxes are normalized to biomass carbon levels. Light period, 0 to 12 h; dark period, 12 to 24 h (indicated by gray shading).
The observed cycle of photosynthesis with midmorning and midafternoon peaks (26, 33) (Fig. 1) is reproduced by the model. This is achieved by an increase in levels of N2-fixing cells (decreased fP) during the middle of the day (from h 3 to h 9). In the early morning, photosynthesis increases rapidly because of increased light levels, followed by a sharp decrease during the middle of the light period (Fig. 4A) due to decreasing levels of photosynthetic cells. After the colony shifts toward greater numbers of photosynthetic cells in the evening (h 9), the photosynthesis rate increases again but at lower rates than in early morning, a difference that is consistent with observations (Fig. 4A). In the model, this difference is due to C storage reaching the maximum capacity of finite cell volumes at the end of the light period. The levels of N cells increase at the expense of P cells, and the rate of N2 fixation rises. This temporal physiological shift has been experimentally observed when a decrease in the rate of photosynthesis coincided with an increase in the rate of N2 fixation (Fig. 4B) (26, 33). The model captures this trend with increased fN levels (and thus decreased fP levels), but does not do so if fN and fP levels are held constant (see Fig. S1 in the supplemental material), suggesting that diurnal shifts in metabolic function of cells are important.
FIG S1
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
The rate of N2 fixation peaks during the middle of the day due to increased levels of nonphotosynthetic cells and accumulated C. During the period from ∼h 0 to ∼h 3, C stores accumulates due to high rates of photosynthesis, leading to a gradual increase in the rate of N2 fixation (Fig. 4B). When the trichomes shift toward greater numbers of N2-fixing cells (fP levels having decreased from 0.8 to 0.55) later in the morning (h 3), the rate of N2 fixation almost doubles. A smaller subsequent decrease is due to loss of C storage and simultaneously N storage getting closer to its maximum, both of which act to reduce FNfix rates (equation 1). In the evening (h 9), the rate is decreased by almost half due to the decreased presence of nonphotosynthetic cells before finally dropping to zero at the onset of the dark period.
Respiratory protection.Light harvesting is essential for providing the organic C and electrons needed to maintain low O2 concentrations via respiratory activity. Over 80% of C is used for respiratory protection during the light period, except for storage accumulation (Fig. 5). In comparison, the consumption of C to supply energy and electrons for N2 fixation was less than 10%. The fraction of C corresponding to respiratory protection was above 70% during the dark period (Fig. 5). This large quantity of C used for respiratory protection explains the previously observed differences between apparent levels of C fixation and O2 evolution. Under optimal conditions, the mean O2 production can reach rates of 26 mg O2 (mg Chl-a)−1 h−1, while the net C fixation rate is 4.5 mg C (mg Chl-a)−1 h−1 (10). The difference in respiration between day and night corresponds to light-dependent respiration levels; the predicted ratio of light-dependent respiration to dark respiration (Fig. 5) is close to what has been previously observed (42). The model also explains the unusually high basal rates of dark respiration observed in Trichodesmium: 0.18 μmol O2 (μg Chl-a)−1 h−1 (10, 42). Estimated rates of basal dark respiration for Skeletonema costatum and Pavlova lutheri are ∼5% to ∼6% of those for Trichodesmium (43, 44). This high respiration rate indicates that Trichodesmium maintains low intracellular O2 levels even during the dark period, possibly to maintain and/or synthesize nitrogenase.
Diurnal allocation of C fluxes to modeled cellular functions. (A) C fluxes in moles of C per mole of C biomass per day. FCsto, C storage production; FRP, respiratory protection; FResN2, respiration for providing energy for N2 fixation; FNfix, carbon consumption for providing electrons for N2 fixation; FResBio (=
The use of C/electron for respiratory protection comes at the expense of cellular growth and therefore may also explain the low growth rate of Trichodesmium. While non-nitrogen fixers have nutrient replete growth rate (μmax) of over 1 day−1 (45), that for Trichodesmium is about 0.1 to 0.5 day−1 (10, 45). Under conditions of nutrient repletion, the μmax rate for the cell can be described as follows:
Diffusion barriers.Even in the absence of photosynthesis, diffusion of O2 from the ambient environment may result in O2 permeating the cell, disabling the active site of the nitrogenase enzyme. Published observations of O2 concentrations within the colony of Trichodesmium reveal strong diurnal variations in O2 concentrations in the boundary layer (Fig. 6A, red circles) (34), with maxima in the early morning and evening, echoing variations in photosynthesis rate. Although O2 is consistently supersaturated in the boundary layer during the day, its levels become slightly lower during the night. These observations provide a key constraint on the model representation of O2 diffusion and the importance of the Trichodesmium strategy to protect N2 fixation by minimizing diffusive O2 fluxes (Fig. 6).
Model-data comparison of O2 concentrations. (A) Modeled time series of O2 with low diffusion coefficient based on the predicted low diffusivity; see “Model diffusivity” in panel B. (B) Modeled daytime average O2 concentrations for various levels of diffusivity of cell membrane layers relative to that of water at 25°C. The key in the center of the figure applies to both panel A and panel B. “P cell,” “N cell,” “Boundary,” and “Environ.” represent model results corresponding to boxes representing photosynthetic cells, nonphotosynthetic cells, the boundary layer, and the environment. Points in panel A represent time series data corresponding to the boundary layer (±300 μm from the center of the colony [34]); compare to “Boundary”. The horizontal dashed line in panel B indicates the median of the daytime data presented in panel A. Light period, 0 to 12 h; dark period, 12 to 24 h (indicated by gray shading in panel A).
The observed trends in levels of O2 in the boundary layer can be captured when we apply low diffusivity for cell membrane layers (lower than 10−3 of that in water) (here referred to as “model diffusivity”) (Fig. 6A and B). First, we predict consistently higher levels of O2 in the boundary layer than in the environment, as observed (34) (Fig. 6A). Despite the respiratory protection, net O2 production rates are still positive, increasing the levels of boundary layer O2. Second, [O2]B reached two peaks: each peak during the early and later parts of the day. This high O2 concentration is due to high rates of photosynthesis (Fig. 4A). The peak value during early light period is slightly higher (∼400 μM) than that in the evening, consistent with the observations (34). The O2 levels decrease during the middle of the day, due to the decreased fraction of photosynthetic cells. However, the concentration (∼300 μM) is still higher than the O2 concentration in the environment [O2]E, also consistent with the observations. The higher value occurs based on the balance between respiration rates and photosynthesis; if the number of photosynthetic cells decreases, [O2]B can become lower than [O2]E. During the dark period, since there is no photosynthesis, the model predicts lower [O2]B (< 200 μM) than [O2]E, as previously observed (34), due to respiration.
Intracellular and boundary layer O2 concentrations are highly sensitive to the levels of diffusivity (Fig. 6B), supporting the importance of strong diffusion barriers. Decreasing diffusivity would increase the passive uptake of O2 by N cells, requiring a higher amount of C. For diffusivity levels exceeding twice that of our default model diffusivity, intracellular O2 cannot be maintained at the minimum level. Also, when the diffusivity is higher than three times the default value, the boundary layer O2 concentration becomes similar to that in the environment, failing to reproduce higher concentrations in the boundary layer environment. If we assume diffusivity of water for the cellular membrane (47), all the boxes have similar averaged O2 levels (Fig. 6B), and N2 fixation cannot be maintained. The predicted low diffusivity is qualitatively consistent with the results of a recent study (38).
DISCUSSION
Potential explanations for low diffusivity.One explanation for the low level of model diffusivity is the low diffusivity of the bacterial membranes. Trichodesmium is a Gram-negative bacterium whose cell envelope has multiple layers (48), with inner and outer lipid membranes separated by a periplasm containing peptidoglycan. The outer lipid membrane is connected to lipopolysaccharide (LPS), creating a capsule. Possibly due to the presence of these layers, the O2 diffusivity of bacterial membranes is predicted to be ∼2 to ∼3 orders of magnitude lower than that of water (12, 49, 50). The diffusivity of cells with simple lipid bilayers is generally within the same order of magnitude as that of water (50). Thus, it is likely that the layers of LPS or peptidoglycan provide strong barriers against O2 in Trichodesmium.
Although the Gram-negative membrane might be sufficient, another potential mechanism for the diffusion barrier may be production of extracellular polymeric substances (EPS). Recent studies suggest involvement of EPS in protecting heterotrophic nitrogen fixers of the species Azotobacter vinelandii (12, 51, 52). This has been supported by a laboratory study where respiration decreases with EPS production (53). Trichodesmium produces EPS especially under nutrient-limited conditions (54–56). Here, given the low model diffusivity, another purpose of EPS might be the management of O2 from the environment, which accounts for over 80% of O2 input in model simulations (see Fig. S2 in the supplemental material). The remaining ∼10% to ∼20% of O2 input is directly transported from neighboring photosynthetic cells, and this intercellular flux must also be minimized.
FIG S2
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
To reproduce the observed fluctuations of O2 concentration in the boundary layer, the model required 45% of nonphotosynthetic cells during the middle of the light period; fractions below 45% result in insufficient diel fluctuations of O2. This value sits between two observational results; some studies show almost all the cells having nitrogenase (30, 57), whereas other laboratory studies show that only ∼10% to ∼20% of cells actually contain nitrogenase (24, 31, 32). In the latter case, nitrogenase may be confined to a subset of nonphotosynthetic cells.
Given the predicted low diffusivity, certain cells may play a role in reducing O2 diffusion from photosynthetic cells to N2-fixing cells. On the basis of this hypothesized function, we refer to them as “buffer cells.” Because cell membranes have significantly slower diffusion and higher viscosity than water, having more membranes between P cells and N cells containing nitrogenase would decrease O2 transport between them. Because EPS can prevent O2 diffusion from the environment only, such buffer cells may represent an alternative to the glycolipid layers seen in heterocysts and thus may be essential for the coexistence of photosynthesis and N2 fixation (Fig. 7).
Model-based hypothesis of how Trichodesmium manages O2 during the light period in comparison to heterocystous cyanobacteria. (A) General understanding of how heterocystous cyanobacteria manage O2. (B) Proposed O2 management by Trichodesmium. The key at the right side of the figure applies to both panel A and panel B. “EPS” represents extracellular polymeric substances (Extra. Poly. Sub.). The model predicts that N2-fixing cells and buffer cells in Trichodesmium would have high respiration for O2 scavenging (respiratory protection). The Gram-negative membranes between cells are doubled. The locations of these cells may be switching in the time scale of minutes, but an anoxic environment can be created instantly (i.e., within seconds) after such a switch due to the high rate of respiration relative to the intracellular O2 concentration.
The double layer of Gram-negative membranes between two adjacent cells could effectively reduce O2 transport. However, such cells and membranes could also slow the transport of carbon or other reduced molecules. To circumvent such problems, we hypothesize that cells may selectively transfer molecules. Such a mechanism seems to exist in heterocystous cyanobacteria, which rapidly transfer sucrose through cell junctions (58). It is also possible that sites for N2 fixation occur where C storage is sufficiently accumulated. This would explain why the lowest fractions of N2-fixing cells are observed just before dawn, when cellular C storage is most depleted (24, 31).
The residence time of O2 is significantly shorter than time scales of metabolic switching.Fluorescence kinetics data show various states of photosynthetic activity throughout a trichome (41). However, until now, it has not been resolved whether cellular O2 management strategies can keep up with changes in photosynthetic activity, which occur over time scales of minutes (25). Since our model focuses on the averaged state of a trichome based on the fraction of different metabolisms (photosynthetic or nitrogen fixing), it does not resolve such dynamic locational shifts of photosynthetic activity. It does show, however, that the residence time of O2 is extremely short, on the order of 1 s. The predicted daytime rate of respiration is ∼2 mol O2 mol C−1 day−1, equivalent to 0.424 mol O2 m−3 s−1 with a cellular C concentration of 1.83 × 104 mol m−3 (59). For a typical concentration of O2 in tropical surface water (∼0.2 mol m−3), the time scale of O2 turnover is ∼0.5 s. Even if we apply our predicted O2 concentration in photosynthetic cells (∼0.6 mol m−3), the time scale is only ∼1.4 s. This simple calculation indicates the potential for photosynthetic cells to quickly switch to nitrogen-fixing cells as long as nitrogenase can be preserved within the photosynthetic cells. The capability of preserving nitrogenase may be supported by evidence indicating that nitrogenase is found in most of the cells (30, 57). However, how nitrogenase is preserved during photosynthesis needs to be further investigated. Potential mechanisms include “conformational protection,” in which the activity of nitrogenase is rapidly and reversibly switched off in response to high O2 concentrations (2, 60, 61), thereby limiting the time-consuming processes of resynthesizing or repairing nitrogenases.
Broader context: how can photosynthesis and N2 fixation occur simultaneously in general?Nitrogenase is highly sensitive to O2 (5, 6, 9); therefore, photosynthetic O2 production is detrimental to N2 fixation. However, some major nitrogen fixers and their symbioses enable these conflicting processes to occur simultaneously; e.g., Trichodesmium, diatom diazotroph association (DDA) (62, 63), and UCYN-A and phytoplankton symbiosis (64). Here, on the basis of a simple model of Trichodesmium, we discuss how these processes can coexist over short distances.
First, N2 fixation and photosynthesis should not occur in the same cell simultaneously. Nitrogenase exists in the cytoplasm, while photosynthesis occurs on the intracytoplasmic membranes. Thus, photosynthesis would likely release O2 into cytoplasm, damaging nitrogenase. Also, we predict that rates of respiration exceeding that of photosynthesis would be required for consuming this O2. This would not be sustainable in cells that lack metabolic specialization. On the other hand, if these processes occur in different cells, the amount of O2 that needs to be managed decreases significantly. Thus, it would make sense that such functional separation has been selected through evolution.
Second, diffusivity control of O2 and high respiration rates must occur together, unless diffusivity is perfectly controlled. The model uses low diffusivity through membranes, but it still requires high respiration rates. High respiration has been reported not only in Trichodesmium, but also Crocosphaera, unicellular photoautotrophic nitrogen fixers (11). Organisms with heterocysts may not require such high levels of respiration with sophisticated O2 barriers. It is likely that decreasing the diffusion can affect the uptake/transfer of other important nutrients/metabolites. Thus, for such strong barriers, selective transfer of organic C substrates to fuel respiration must be involved.
Conclusions.With a mechanistic model of Trichodesmium, we investigate how these photosynthetic cyanobacteria manage to fix nitrogen while also fixing carbon. The model reproduces observed patterns of diel cycling in Trichodesmium physiology and indicates that these conflicting processes must occur in different cells and also need to be separated by a barrier to O2 diffusion. These results support the hypothesis that a Gram-negative membrane could represent an essential diffusion barrier against O2. EPS and buffer cells might additionally decrease O2 diffusion. The remaining O2 must be consumed through respiratory protection for N2 fixation to occur, and such respiration may explain the observed low growth rate of Trichodesmium. The model indicates that the residence time of O2 within trichomes is a few seconds, significantly shorter than the time scale of switching on and off photosynthesis. Our model and the results provide a theoretical basis for investigating how these two conflicting processes occur in one of the major sources of new N to the oligotrophic surface ocean.
MATERIALS AND METHODS
Model equations.In this section, we explain the assumptions and equations used in the model. This model resolves C, N, and O2 budgets (Fig. 2) within cells, as well as the O2 budget in the immediate vicinity (boundary layer) of those cells. The budgets consist of biochemical rates within each type of cell and exchange between cells of each type and between cells and the environment. At any given time, these rates depend on the fraction of the trichome that is assigned to each metabolic function. We denote fP and fN as the fractions of photosynthetic and nonphotosynthetic cells, respectively. Of the nonphotosynthetic cells (fN), a constant fraction (fNITROGE) is assumed to contain nitrogenase and thus to be able to fix N2. The detailed nomenclature (with units) is provided in Table S1 in the supplemental material. To determine the rate of N2 fixation (equation 1), we need to specify N and C storage and cellular O2 in the N2-fixing cells. Since these quantities are themselves governed by an allocation of resources that varies with environmental conditions, we model these variables with a set of differential equations. The difference between C production (photosynthesis) and consumption (respiration, N2 fixation, and biomass production) leads to the change in C storage (CSto). Similarly, the change in N storage (NSto) is governed by the balance of N2 fixation and biomass production as represented in the following equations:
TABLE S1
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
CSto and NSto are governed by 2 common rates (FNfix and FBio), whereas CSto has additional input from C fixation and respiration that need to be modeled. The CSto and NSto values are computed with a finite-difference method and therefore cause temporal variations in C:N values.
We assume that the rate of biomass production (FBio) is a minimum function of available storage resources (CSto and NSto), calculated as follows:
The rate of C fixation (FCfix) is proportional to the level of chlorophyll available for photosynthesis:
Here C fixation saturates with light and slows as C storage approaches its maximum capacity. We assume that the level of chlorophyll per P cell is constant and proportional to the fraction of P cells fP (= 1 – fN) since the model runs at constant light during the light period as follows:
The total respiration rate (FRes) is the sum of the rates in P cells (
We assume that biosynthesis occurs mostly in P cells, and the respiration for supporting the biosynthesis assumed to be proportional with a constant ratio of YRes:Bio, leading to the following equation:
To compute
Here, Vi represents the volume of box i (i = P, N, or B, representing boxes for photosynthetic cells, nonphotosynthetic cells, or boundary layer, respectively); V = VP+VN, [O2]i represents the O2 concentration of box i (i = P, N, or B);
Since the time scale of O2 equilibration is much smaller than that of N and C, we assume a pseudo-steady state of O2. This assumption largely applies since the cellular concentrations of C and N are on the order of 104 and 103 (mol m−3), respectively, while that of O2 is generally below 1 (mol m−3) despite the fact that the magnitudes of the fluxes are similar. The pseudo-steady-state assumption gives a single solution for
Parameterization.The model requires specification of several parameters which are obtained from previous studies (Table S2). These include elemental stoichiometry of biomass and boundary conditions such as light intensities and O2 concentrations in the environment. The remaining parameters are selected to reproduce the observed metabolic rates of Trichodesmium as compiled from multiple published studies (10, 27, 29, 42, 65–67) (Table S3). In general, the key metabolic processes follow a well-defined order, with N2 fixation being the slowest (∼0.006 to ∼0.146 mol N mol biomass C−1 day−1), O2 production being the fastest (∼1.10 to ∼288 mol O2 mol biomass C−1 day−1), and net C fixation being intermediate between the other two (∼0.16 to ∼2.57 mol C mol biomass C−1 day−1). The remaining parameters have been tuned to bring averaged model output values within these ranges.
TABLE S2
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S3
Copyright © 2019 Inomura et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Model availability.The model in this study was written in Python 3 and is freely available at https://zenodo.org/record/1245128.
ACKNOWLEDGMENTS
We thank Michael J. Follows, Jonathan P. Zehr, and Meri Eichner for useful discussion.
This study was supported by the Simons Foundation (Life Sciences-Simons Postdoctoral Fellowships in Marine Microbial Ecology, award 544338, K.I.), the Simons Collaboration on Ocean Processes and Ecology (SCOPE; grant 329108 to David M. Karl), National Science Foundation (OCE-1756524 to S.T.W.), and the Gordon and Betty Moore Foundation (GBMF grant 3775, C.D.).
FOOTNOTES
- Received March 22, 2019.
- Accepted July 15, 2019.
- Copyright © 2019 Inomura et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
