Original Research ARTICLE
Temporal Patterns and Intra- and Inter-Cellular Variability in Carbon and Nitrogen Assimilation by the Unicellular Cyanobacterium Cyanothece sp. ATCC 51142
- 1Department of Earth Sciences, Utrecht University, Utrecht, Netherlands
- 2Institute of Microbiology, Czech Academy of Sciences, Centre Algatech, Třeboň, Czechia
- 3Max Planck Institute for Marine Microbiology, Bremen, Germany
- 4Sorbonne Université, CNRS, Laboratoire d’Océanographie de Villefranche, Villefranche-sur-mer, France
- 5Sorbonne Université, CNRS, Laboratoire d’Océanographie Microbienne, Banyuls-sur-mer, France
- 6Institute of Parasitology, Czech Academy of Sciences, Biology Centre, České Budějovice, Czechia
- 7Centre for Ecological Research, Balaton Limnological Institute, Tihany, Hungary
- 8Botany Department, Federal University of Santa Catarina, Campus de Trindade, Florianópolis, Brazil
- 9Department of Biology, Mount Allison University, Sackville, NB, Canada
- 10Laboratoire de Biologie des Organismes et Ecosystèmes Aquatiques, FRE 2030, Muséum National d’Histoire Naturelle, CNRS, IRD, Sorbonne Université, Université de Caen Normandie, Normandie Université, Esplanade de la Paix, France
- 11Global Change Research Institute, Czech Academy of Sciences, Brno, Czechia
- 12STU-UNIVPM Joint Algal Research Center, Marine Biology Institute, College of Sciences, Shantou University, Shantou, China
- 13NIOZ Royal Netherlands Institute for Sea Research and Utrecht University, Den Burg, Netherlands
- 14Universidade Federal de São Carlos, São Carlos, Brazil
- 15Department Experimental Limnology, Leibniz-Institute of Freshwater Ecology and Inland Fisheries, Stechlin, Germany
- 16Department of Ecology, Berlin Institute of Technology, Berlin, Germany
- 17University of Technology Sydney, Climate Change Cluster, Faculty of Science, Ultimo, NSW, Australia
- 18Department of Microbiology, Oregon State University, Corvallis, OR, United States
Unicellular nitrogen fixing cyanobacteria (UCYN) are abundant members of phytoplankton communities in a wide range of marine environments, including those with rapidly changing nitrogen (N) concentrations. We hypothesized that differences in N availability (N2 vs. combined N) would cause UCYN to shift strategies of intracellular N and C allocation. We used transmission electron microscopy and nanoscale secondary ion mass spectrometry imaging to track assimilation and intracellular allocation of 13C-labeled CO2 and 15N-labeled N2 or NO3 at different periods across a diel cycle in Cyanothece sp. ATCC 51142. We present new ideas on interpreting these imaging data, including the influences of pre-incubation cellular C and N contents and turnover rates of inclusion bodies. Within cultures growing diazotrophically, distinct subpopulations were detected that fixed N2 at night or in the morning. Additional significant within-population heterogeneity was likely caused by differences in the relative amounts of N assimilated into cyanophycin from sources external and internal to the cells. Whether growing on N2 or NO3, cells prioritized cyanophycin synthesis when N assimilation rates were highest. N assimilation in cells growing on NO3 switched from cyanophycin synthesis to protein synthesis, suggesting that once a cyanophycin quota is met, it is bypassed in favor of protein synthesis. Growth on NO3 also revealed that at night, there is a very low level of CO2 assimilation into polysaccharides simultaneous with their catabolism for protein synthesis. This study revealed multiple, detailed mechanisms underlying C and N management in Cyanothece that facilitate its success in dynamic aquatic environments.
Nitrogen fixing microorganisms (diazotrophs) are critical suppliers of bioavailable forms of nitrogen (N, e.g., ammonium) in natural ecosystems. In the surface ocean where primary production is often limited by N availability, unicellular N2-fixing cyanobacteria (UCYN) are now recognized as having key roles in biogeochemical cycles (Zehr et al., 2001; Karl et al., 2002; Montoya et al., 2004; Zehr, 2011; Wilson et al., 2017).
Because of the scarcity of cultured representatives, Cyanothece sp. ATCC 51142 (henceforth Cyanothece 51142, recently reclassified as Crocosphaera subtropica; Mareš et al., 2019) has become an experimental model of UCYN (Reddy et al., 1993; Schneegurt et al., 1994; Colón-López and Sherman, 1998; Sherman et al., 1998; Li et al., 2001). The sequenced Cyanothece 51142 genome and controlled Cyanothece culture studies have provided insights into the genetic controls underlying the temporal segregation of N2 fixation activity, regulation of photosynthesis (Stöckel et al., 2008; Toepel et al., 2008; Welsh et al., 2008), and carbon (C) metabolism (Schneegurt et al., 1994; Colón-López and Sherman, 1998; Aryal et al., 2011; Bernstein et al., 2015). Cyanothece 51142 appear to restrict N2 fixation to the night time to protect the N2 fixing enzyme complex, nitrogenase, from inactivation by molecular oxygen produced by photosynthesis during the day time (Reddy et al., 1993; Welsh et al., 2008). Intracellular storage of the newly fixed C and N facilitates temporal separation of activities. Specifically, N fixed during the night is initially stored as cyanophycin until it is used for processes such as protein and nucleic acid synthesis, whereas C fixed during the day is stored as polysaccharides until it is respired the following night to supply reducing equivalents (NADPH) and ATP needed to support N2 fixation (Schneegurt et al., 1994; Li et al., 2001; Großkopf and LaRoche, 2012; Inomura et al., 2019). C respiration during the late afternoon and night also facilitates nitrogenase activity by depleting molecular oxygen that diffuses into the cell from the environment (Dron et al., 2012; Großkopf and LaRoche, 2012; Červený et al., 2013; Inomura et al., 2019).
In addition to their ability to fix N2, Cyanothece (and other UCYN) can assimilate various forms of combined nitrogen (e.g., NH4, NO3, urea, and amino acids; Mulholland et al., 2001; Holl and Montoya, 2005; Dekaezemacker and Bonnet, 2011; Masuda et al., 2013). Combined nitrogen generally down-regulates N2 fixation because its uptake and assimilation are energetically less costly than the processes supporting nitrogenase activity (Mulholland et al., 2001; Holl and Montoya, 2005; Eichner et al., 2014). This flexibility in nitrogen metabolism is considered to be one reason why Cyanothece thrives in a variety of marine environments with widely ranging nitrogen concentrations (Rippka, 1988; Short and Zehr, 2007; Webb et al., 2009; Bonnet et al., 2011). Cyanothece 51142 appears to efficiently manage C and N storage using multiple copies of genes encoding enzymes in polysaccharide metabolism and conserved gene clusters that coordinate intersecting pathways of C and N metabolism (Welsh et al., 2008; Zhang et al., 2018). For example, the ornithine-ammonia cycle (OAC) may facilitate efficient sequestration and remobilization of N (Zhang et al., 2018). Other N-rich compounds, including arginine, can be catabolized to recycle N within the cell (Flores et al., 2019; Burnat et al., 2019). CO2 and N2 fixation activities are also at least partly regulated according to the cell cycle demands that vary across the day–night cycle (Dron et al., 2013). For example, the N demands of nucleic acid synthesis and the C and energy requirements for new cell synthesis generally restrict cell division to hours when intracellular N reserves and photosynthetic rates are high (Dron et al., 2013; Červený et al., 2013). In contrast to N2 fixation, combined N assimilation and photosynthetic C assimilation should not require temporal separation, raising questions about how N and C allocations change depending on whether cells are growing on N2 or combined N.
In this study, we characterized C and N assimilation rates over a diel cycle in Cyanothece 51142 grown under obligate diazotrophic and non-diazotrophic conditions (with added NO3). Our approach combined nanoscale secondary ion mass spectrometry (nanoSIMS) and transmission electron microscopy (TEM) to track the assimilation of 13C-labeled inorganic C and 15N-labeled N2 or NO3 into individual cells and resolve their assimilation into polysaccharides, cyanophycin, and other inclusion bodies over the day–night cycle. We observed clear differences in N allocation patterns, but also unexpected within-population heterogeneity, including wide variation in labeling of storage inclusions and day-time N2 fixation. We discuss these observations and highlight how careful evaluation of these complex nanoSIMS data revealed key metabolic mechanisms underlying C and N management in Cyanothece 51142 that facilitate its success in dynamic aquatic environments.
Materials and Methods
Bioreactor and Semi-Continuous Cultures
Cyanothece 51142 cultures were maintained in 400- or 1,000-mL flat panel photobioreactors (FMT 150, Photon Systems Instruments, Brno, Czech Republic) at 28°C and 300 μmol photons m–2 s–1 with a 14 h:10 h light:dark cycle (14L:10D), with L0 at 07:30 and D0 at 21:30. The initial and final hours of each light cycle were set to follow a sinusoidal increase or decrease in light intensity, respectively. Triplicate cultures were grown in ASP2 medium (Provasoli et al., 1957; Van Baalen, 1962) either supplemented with 17 mM NO3 (“NO3 culture”) or prepared with no inorganic N added (“N2 culture”). Cultures were bubbled with ambient air (400 mL min–1) and maintained under turbidostat mode at OD680 of ∼0.5. Another set of triplicate cultures were grown in ∼300 mL glass tubes in ASP2 medium without NO3 amendment at 28°C under 300 μmol photons m–2 s–1 (14L:10D, same light regime as above) under semi-continuous, fed-batch mode (“SC-N2 culture”). The SC-N2 cultures were cultivated and maintained to ensure a “back-up” option in case the turbidostat (N2) cultures failed during the experiment. Because both culture conditions were stable throughout the experiment, we opportunistically sampled both cultures. The SC-N2 culture exhibited a wider range of phenotypic variability compared to the N2 culture that may be akin to some more dynamic natural environments. Therefore, we chose to include results for both diazotrophic cultures. All cultures maintained constant population sizes for >21 days prior to sampling. Culture and cell characteristics for each treatment are given in Table 1. An analysis of photosynthetic electron transport and the energetic costs of N and C acquisition in the same cultures studied here are given in Rabouille et al. (In Revision).
Cell density and size distributions were determined using a Multisizer 4 Coulter Particle Counter (Beckman Coulter Inc., Brea, CA, United States). Particulate organic carbon (POC) and nitrogen (PON) were analyzed using an elemental analyzer (PerkinElmer PE2400, PerkinElmer Inc., Waltham, MA, United States) following sample collection (10 mL), centrifugation (28°C, 8,000 rpm, 7 min), and drying at 60°C. Chl a and polysaccharide contents were determined following the previously reported protocols (Zavřel et al., 2015a, b, 2018). Cyanophycin cell content was quantified by Sakaguchi reaction (Messineo, 1966), after sample concentration (30 mL) by centrifugation (28°C, 8,000 rpm, 7 min).
Stable Isotope Probing Experiments
For stable isotope probing experiments, samples were collected from one culture replicate from each condition. Incubations were performed by sub-sampling cultures into 6 mL gas-tight vials and amending the ASP2 medium with NaH13CO3 (all cultures) and either 15N2 (N2 and SC-N2 cultures) or Na15NO3 (NO3 culture). Vials were incubated under light and temperature conditions that were equivalent to culture conditions. Incubation durations were 2 h in the morning, 2.5 h during the day, and 10 h during the night, with shorter incubations in early night (3 h) and late night (5 h). Isotope labeling was calculated from the known amounts of label added to the incubation medium and measured initial concentrations of unlabeled substrate in the bioreactors (Supplementary Table 1). The 15N2 enriched stock was prepared by injecting 10 mL of 15N2 gas into 43 mL of the ASP2 medium, followed by an equilibration for >24 h. The 15N-N2 atom fraction in the incubation medium was calculated assuming that 15N2 was fully equilibrated with the stock solution. Since this may lead to an underestimation of N2 fixation rates (Mohr et al., 2010), we refrain from comparisons of N assimilation rates between N treatments (N2 vs. NO3). However, comparisons over time and among cells within each treatment are not affected since any potential underestimation would be similar in all 15N2 incubations. NO3 concentration in the incubation medium was estimated by averaging NO3 concentrations measured in the bioreactor on the respective day of the experiment. Dissolved inorganic carbon (DIC) concentration in the incubation medium was estimated by measuring DIC concentrations in the bioreactor at three time points during the day and interpolating them to the starting time points of our stable isotope incubations. Because the DIC concentration in the cultures varied depending on the time of the day, 13C-DIC atom fractions varied during our incubations, although the amounts of added NaH13CO3 were the same (Supplementary Table 1).
At the end of each isotope-labeling incubation, cells were collected and centrifuged at 2,700 rpm for 10 min at room temperature. One microliter of the pellet was mixed with 1 μL of 20% bovine serum albumin and transferred to a formvar-coated 100 mesh TEM grid. After removing the excess liquid with a filter paper, the grids were frozen in liquid ethane cooled with liquid nitrogen. Freeze-substitution was carried out in a 2% mixture of OsO4 in 100% acetone (v/v) sequentially at three temperatures: −90°C (for 96 h), −20°C (for 24 h), and 4°C (for 10 h). Temperature was increased at a rate of 5°C h–1 (from −90 to −20°C) and 3°C h–1 (from −20 to 4°C). After freeze-substitution, the samples were washed three times in acetone and infiltrated sequentially in a 2:1, 1:1, and 1:2 (v/v) mixture of acetone and low-viscosity Spurr resin (EMS) for 1 h in each step. Finally, the samples were incubated overnight in a 100% resin, transferred to embedding molds, and allowed to polymerize. Thin sections (200 nm) were cut with a diamond knife, placed on Cu-indexed TEM grids (rinsed in 30% ethanol), and contrasted for 20 min in saturated ethanolic uranyl acetate (EMS, Hatfield, United States; concentration 13 g/100 mL 50% ethanol; solution filtered before use through a 0.45 μm pore size filter). Images were taken using a JEOL 1010 TEM at 80 kV.
Nanoscale secondary ion mass spectrometry analyses were performed on two types of samples: (i) thin sections that were first imaged by TEM (as described above) and (ii) whole cells collected on polycarbonate filters. For downstream analysis of samples initially imaged by TEM, the lowest primary ion beam current (0.5 pA) was used to achieve the highest lateral resolution afforded by the instrument (∼50 nm). However, because the samples were very thin (∼200 nm), the number of imaged frames was rather low (20–50) before the cell material was sputtered away. This sometimes resulted in a poor signal-to-noise ratio (SNR) in the final secondary ion images, and thus, a low number of cells for which good quality complementary TEM and nanoSIMS images are available. Additional measurements were therefore performed on cells deposited on filters, because the imaging could be done with a stronger beam (2 pA) and over a larger area and many more frames (>200). However, the improved throughput and SNR came at the expense of a lower spatial resolution (see section “Results”). For nanoSIMS analysis of whole cells, the cells were filtered onto polycarbonate filters (2.5 cm diameter, 0.2 μm pore size, Millipore), washed three times, air-dried, and stored at room temperature. Chemical fixation was not performed thus avoiding dilution of the isotope label. Just prior to nanoSIMS analysis, filters were sputter-coated with a 10-nm gold layer, cut into small circular pieces (5 mm diameter) suitable for the nanoSIMS sample holder, and imaged with a Neoscope II JCM-6000 scanning electron microscope (JEOL, Japan) to check sample quality (cell integrity and cell density).
Nanoscale secondary ion mass spectrometry measurements were performed with the NanoSIMS 50L instrument (Cameca, France) operated at Utrecht University. Areas of interest were first pre-sputtered with Cs+-ions until secondary ion yields stabilized. Subsequently, the primary Cs+-ion beam was scanned over the sample (areas between 10 μm × 10 μm and 30 μm × 30 μm in size, dwell time of 1 ms pixel–1) while detecting secondary ions 12C–, 13C–, 16O–, 12C14N–, 12C15N–, 31P–, and 32S–. To increase the overall signal, the same area was imaged multiple times, and the resulting ion count images were aligned and accumulated.
NanoSIMS Data Processing and Quantification of Rates
Nanoscale secondary ion mass spectrometry data were processed with the Look@NanoSIMS software (Polerecky et al., 2012) to quantify 13C and 15N atom fractions, denoted as x(13C) and x(15N) (Coplen, 2011), in regions of interest (ROI’s) corresponding to cells or inclusion bodies (polysaccharide granules or cyanophycin inclusions). After drawing ROIs manually, x(13C) in the ROI was determined from the total counts of secondary ions 12C– and 13C– accumulated over the ROI pixels as x(13C) = 13C–/(12C– + 13C–). Similarly, x(15N) in the ROI was determined from the total counts of 12C15N– and 12C14N– accumulated over the ROI pixels as x(15N) = 12C14N–/(12C14N– + 12C15N–).
The C- and N-specific rates of 13C and 15N assimilation into whole cells (kC and kN, respectively) were calculated as:
while the C- and N-specific rates of 13C and 15N incorporation into polysaccharide granules (pC or pN, respectively) and cyanophycin inclusions (yC and yN, respectively) were calculated as:
In Eqs 1–4, x(13C)S and x(15N)S are atom fractions of the C and N source, respectively, and x(13C)ini and x(15N)ini are the initial atom fractions of C and N in the ROI, respectively. The isotope labeling of the C and N sources, x(13C)S and x(15N)S, was assumed to be constant during the incubation and was calculated as described above (Supplementary Table 1). x(13C)ini and x(15N)ini were determined by averaging data obtained from cells that were not exposed to the labeled substrate (i.e., control cells; x(13C)ini = 1.052 × 10–2, SD(13C)ini = 0.007 × 10–2, and x(15N)ini = 3.75 × 10–3, SD(15N)ini = 0.04 × 10–3, n = 30). A cell or an inclusion body was considered significantly enriched in 13C if the 95% confidence interval of its estimated mean 13C atom fraction did not overlap with that of the control cells, i.e., if x(13C) ± 2 × SE(13C) did not overlap with x(13C)ini ± 2 × SE(13C)ini. Here, the standard errors were calculated as SE(13C)ini = SD(13C)ini/ for the control cells, and SE(13C) = x(13C) × PE(13C) for each individual cell or inclusion body, where the relative Poisson error was calculated from the total counts of 12C– and 13C– in the cell or inclusion body as PE(13C) = [1 – x(13C)] × [1/13C– + 1/12C–]1/2 (Polerecky et al., 2012). The same approach but using the total counts of 12C14N– and 12C15N– was applied to determine significant enrichment in 15N.
Note that the C- and N-specific rates of 13C and 15N assimilation into whole cells, polysaccharide granules, and cyanophycin inclusions have units of per time (i.e., h–1 or day–1) and give the rate of 13C and 15N assimilation rates normalized to the C and N content of the ROI [i.e., mol C (mol C)–1 h–1 for kC, pC, and yC, and mol N (mol N)–1 h–1 for kN, pN, and yN]. Evaluation of the variability in assimilation rates among cells and intracellular inclusions required considering how the measured x(13C) and x(15N) in the ROI were linked to cell growth and metabolism of internal C and N pools (e.g., synthesis and degradation of polysaccharides and cyanophycin inclusions, or recycling of N from existing proteins during cyanophycin synthesis). These considerations are summarized in the Discussion section (Section “Interpreting Isotopic Enrichment Imaging Data”). The assumptions underlying the rate calculations presented above are further explored and discussed in Polerecky et al. (In Revision).
Identification of Intracellular Inclusions
Prominent intracellular inclusion bodies identified in TEM images of Cyanothece 51142 included carboxysomes, cyanophycin inclusions, polyphosphate bodies, polysaccharide granules, and thylakoid membranes (Figure 1). Some, but not all, of these inclusions could be reliably identified in nanoSIMS images when the accumulated secondary ion counts 12C14N–, 32S–, and 31P– were combined into RGB overlays. Specifically, carboxysomes, which were identified in TEM images as dark areas with a characteristic hexagonal shape (Figures 1A,C, arrows labeled “c”), had relatively higher CN– and S– but lower P– counts than the surrounding cell material. These differences caused carboxysomes to appear yellow-green in the RGB overlays (Figures 1D,F). Cyanophycin inclusions, which were identified through their oval shape and darker appearance in the TEM images, had markedly higher CN– counts, while the S– and P– counts were not different from the surrounding cell material. The relative enrichment in CN– counts gave cyanophycin inclusions an orange-to-red appearance in the RGB overlays (Figure 1, arrows labeled “cy”). Polyphosphate bodies had higher P– counts and lower CN– and S– counts relative to the surrounding cell material and appeared as bluish spots in the RGB overlays (Figure 1, arrows labeled “p”). Although polysaccharide granules could be identified in TEM images as bright oval shapes (Figure 1A, arrows labeled “ps”; see also Deschamps et al., 2008), their identification from the nanoSIMS images was not reliable. For example, most polysaccharide granules were associated with localized decreases in CN– counts (Figure 1D). However, the contrast between the polysaccharide granules and the cell matrix was low, and similar decreases in CN– counts sometimes occurred even when there was no obvious presence of polysaccharide in the TEM images. Neither could variability in S–, P–, C–, or O– counts (C– and O– data not shown) be used to distinguish polysaccharide granules in the nanoSIMS images. Similarly, thylakoid membranes were not identifiable using the nanoSIMS images although their visibility in the TEM images was often good (Figures 1A–C, arrows labeled “t”). Conversely, nucleoids were clearly observed as the violet-colored regions in the RGB overlays due to P counts being markedly higher than the surrounding cell material and CN– and S– counts that were similar to the surrounding cell material (Figures 1D–F, arrows labeled “n”); however, direct nucleoid identification in the TEM images was not possible.
Figure 1. Correlative microscopy of thin sections of Cyanothece 51142 cells. Shown are examples of TEM images (A–C) and the corresponding nanoSIMS images (D–F). The dashed polygon in panels (D–F) shows the boundary of the TEM image aligned within the NanoSIMS image. Shown are images from samples collected at 09:30 following a morning incubation (07:30–09:30) from the N2 culture (A,D) and NO3 culture (B,E), and from a sample collected at 07:00 following a night-time incubation (21:45–07:00) from the NO3 culture (C,F). Arrows in the images point to sub-cellular structures identified as carboxysomes (c), cyanophycin inclusions (cy), polyphosphate bodies (p), polysaccharide granules (ps), thylakoid membranes (t), and nucleoid (n). NanoSIMS images are shown as RGB overlays of secondary ion counts 12C14N– (red), 32S– (green), and 31P– (blue). Note that the hues are not comparable among the images because, for each overlay, the contrast for the three color channels was modified so as to enhance the overall visibility of the intracellular variability. In addition to P–, the polyphosphate bodies had markedly increased O– ion counts (data not shown).
Carbon and Nitrogen Assimilation Rates and Allocation Patterns
Daytime C fixation in Cyanothece 51142 grown under diazotrophic conditions in turbidostat mode (“N2 culture”) was observed in all but one of the 126 cells imaged (Figures 2A,B). The C-specific rates of 13C assimilation in whole cells, kC, were highest in the morning and declined on average by about 80% in the afternoon (Table 2). 13C enrichment was highest in polysaccharide granules and lower and diffusely spread throughout the cell matrix (Figures 3A,B). The C-specific rates of 13C assimilation in polysaccharide granules, pC, varied (CV≈32%), with 61% of the variance explained by differences among cells and 39% of the variance explained by differences within cells (Supplementary Figure 1A). Moreover, the relative area of the cell sections covered by polysaccharide granules varied among cells (range: 0.08–0.28, CV≈43%) and was significantly positively correlated with the kC values (R = 0.68, p = 0.002; Supplementary Figure 2B).
Figure 2. Element-specific rates of C and N assimilation in Cyanothece 51142. Data points in panels (A–J) show element-specific assimilation rates of 13C [kC, in mol C (mol C)–1 h–1] and 15N [kN, in mol N (mol N)–1 h–1] for individual cells from three parallel cultures (N2, diazotrophic culture grown in turbidostat mode; NO3, culture grown on nitrate in turbidostat mode; SC-N2, diazotrophic culture grown in semi-continuous batch mode). Cells were incubated during time intervals indicated above the graphs. L0 was 07:30 and D0 was 21:30. Data points are colored according to the amounts of C and N assimilated (see color scheme). The number of measured cells in each subpopulation is given in the legend. Subpopulations labeled P1–P9 are described in further detail in text. Note different scaling of the x and y axes among the different panels.
Figure 3. Images of isotopic enrichment in Cyanothece 51142. Shown are overlays of 13C (red) and 15N (green) enrichments measured in cells from three parallel cultures (N2, diazotrophic culture grown in turbidostat mode; NO3, culture grown on nitrate in turbidostat mode; SC-N2, diazotrophic culture grown in semi-continuous batch mode). Cells were incubated during time intervals indicated above the images. Images shown in panels (A,B,E–G) were obtained from thin cell sections analyzed by TEM [panels (A,E,G) correspond to panels (A–C) in Figure 1], whereas panels (C,D,H–J) show images of cells collected on polycarbonate membrane filters. Additional images of filtered cells are shown in Supplementary Figure 2. In each image, the intensity (“brightness”) of the red and green color scales linearly with the 13C and 15N enrichment, respectively, with black indicating no enrichment. Note, however, that because the scaling for the red and green colors was optimized independently for each image to enhance the visibility of the intracellular heterogeneity, the intensities of the red and green colors are not comparable among the images. In panels (A,C,E–G,I,J), examples of cyanophycin inclusions (cy) and carboxysomes (c) are marked with dashed-line and solid-line arrows, respectively. In panels (A,B,E–G), areas and spots with pronounced 13C enrichment (red) correspond to polysaccharide granules. Black areas in panels (F,G) correspond to voids (v) due to artifacts associated with the preparation of the thin cell sections.
The majority of cells in the N2 culture fixed N2 during the night (Figure 2C), with the exception of one cell (out of 104 imaged) that fixed N2 in the morning (green circle in Figure 2A; “asynchronous diazotrophic” cell in Figure 3A). During the early night, cells grouped into two clear subpopulations exhibiting different activities: P1 showed significant N2 fixation and accounted for 83% of cells, while P2 showed no significant N2 fixation and accounted for 17% of cells (compare green and red circles in Figure 2C). P1 and P2 showed low but significant C fixation during the early night incubation (Figure 2C). The kC values in P1 and P2 did not differ at night [ANOVA, F(1,77) = 1.64, p = 0.204] and were about 1–1.5% of morning kC values. The kC and kN values did not correlate in subpopulation P1 (R = 0.043, p = 0.73), and kC:kN ranged from 0.008 to 0.035 (mean = 0.016, SD = 0.008, CV≈50%) among the cells. During the late night, only about 3% of cells fixed N2 (green pluses in Figure 2C), and the average kN was about 1% of the average kN in the early night (Table 2). Significant C fixation was observed in 5% of cells, and these cells again grouped into subpopulations depending on whether they also fixed N2 (3%) or not (2%) (compare green and red pluses in Figure 2C). The remaining 95% of cells showed no detectable C nor N2 fixation during the late night (black pluses in Figure 2C).
The C and N2 fixation patterns observed separately during the day and night were reflected in cells incubated with isotopes over the full night + day cycle (Figures 2D, 3D). The majority of cells (∼87%) fixed both C and N2 (P3, orange circles in Figure 2D), whereas the remaining 13% of cells had kC values similar to P3 but showed no significant N2 fixation (P4, red circles in Figure 2D). Average kC and kN values were similar for the majority of cells incubated for the full night + day cycle, but these values were only about 5% of the peak kC and kN values observed during the morning and early night, respectively (Table 2). These dramatic shifts in metabolism over the day caused kC in the morning to overestimate specific growth rate (μ = 0.22 day–1) by 7.7-fold. Theoretically, the average night + day kC value should estimate μ, but was only 0.09 day–1. Similarly, kN at night overestimated μ by up to 6-fold, and the average night + day kN was only 0.08 day–1. The discrepancies between μ and the average kC and kN from the night + day (∼20 h) incubations were likely caused by differences in incubation conditions, including gas flow and medium exchange in the turbidostat that were not possible in the isotope labeling incubation. Finally, kC:kN varied widely among cells (range 0.6–3.3, mean = 1.24, SD = 0.78, CV = 63%).
Cyanophycin inclusions showed the greatest 13C and 15N enrichment compared to other inclusions within individual cells from the early night incubation (Figures 3C, 4A–C). Similar to the data clustering observed for the whole cells, individual cyanophycin inclusions grouped into two clear subpopulations, one with significant 15N enrichment (P1) and one with no significant 15N enrichment (P2) (compare green and red circles in Figure 4B). There was no significant correlation between the C- and N-specific 13C and 15N assimilation rates in cyanophycin granules, yC and yN, in P1 (R = 0.10, p = 0.10), and yC:yN were highly variable between cyanophycin granules among cells (range: 0–0.09, mean = 0.015, SD = 0.016, CV = 107%) and within individual cells (Figure 4A). In one cell where we could clearly resolve all relevant intracellular structures, we observed low but significant 15N enrichment in the carboxysomes and a slightly greater 15N enrichment in the nucleoid in addition to the strong 15N enrichment in the cyanophycin inclusions (Figure 4C).
Figure 4. Within-cell heterogeneity of N assimilation in N2-fixing Cyanothece 51142. Shown are data for cells from the N2 culture (diazotrophic culture grown in turbidostat mode) incubated during early night [21:45–00:45; (A–C)] and from the SC-N2 (diazotrophic culture grown in semi-continuous batch mode) incubated during morning [7:30–9:30; (D)] and afternoon [14:45–17:15; (E,F)]. L0 was 07:30 and D0 was 21:30. Images show an overlay of the 13C (red) and 15N (green) atom fractions (A), the 15N atom fraction (C), and the 13C atom fraction (F). The corresponding TEM images of selected cells are shown in the inset. Dashed-line arrows point to cyanophycin inclusions (cy), while solid-line arrows point to carboxysomes (c), voids due to artifacts associated with the preparation of the cell sections (v), and areas where the atom fractions could not be quantified due to low secondary ion counts (lc). In panel A, numbers indicate the ratios of C and N-specific rates of 13C and 15N assimilation in individual cyanophycin granules, yC:yN. In panel C, the 15N atom fractions in the marked inclusion bodies range between 0.022 and 0.035 for cyanophycin (cy), 0.0047 and 0.0052 for carboxysomes (c), and 0.0082 and 0.0083 for nucleoid (n), and are all significantly greater than in control cells [x(15N)ini = 0.00375]. (B,D,E) Scatter plots of yC vs. yN in cyanophycin granules. Values depicted with open and filled symbols were derived from images obtained from thin cell sections and whole cells deposited on a filter, respectively. Data points are colored as in Figure 2. Correlation coefficients (R) and the corresponding p-values are also shown for green symbols in panel (B) (P1, corresponding to values from the P1 subpopulation shown in Figure 2C) and orange symbols in panel (D) (P6, corresponding to values from the P6 subpopulation shown in Figure 2H). Red data points correspond to cyanophycin inclusions with no significant 15N enrichment.
Cyanothece 51142 grown under non-diazotrophic conditions in turbidostat mode (with NO3 added; “NO3 culture”) showed daily patterns of C fixation that were similar to the N2 culture. Morning kC was 2.7-fold greater than μ of 0.31 day–1. In the afternoon, kC decreased on average by 50% and reached about 0.5% of the morning values during the night time (Figures 2E–G and Table 2). Cells in the NO3 culture always assimilated newly fixed C into polysaccharide granules (Figures 3E–G), whereas in the N2 culture the highest 13C enrichment was observed in polysaccharide granules during the day but in cyanophycin inclusions during the night (compare Figures 3A–C).
The daily patterns of N assimilation and intracellular allocation were more complex in the NO3 culture than in the N2 culture. Values of kN were highest in the morning (10.3-fold greater than μ) and decreased by about 60% in the afternoon and by 90% during the night (Figures 2E–G and Table 2). In the morning, 15N was often accumulated in what appear to be newly synthesized cyanophycin inclusions (Figure 3E and Supplementary Figure 3D, arrows 1) or was added to existing cyanophycin inclusions as manifested by a 15N-rich “shell” surrounding a 15N-poor core (Figure 3E, arrow 2). In many cells, 15N enrichment was clearly present in the cell matrix but not in existing cyanophycin inclusions (Supplementary Figure 3D, arrows 3). In the afternoon, 15N enrichment was distributed relatively homogeneously within the cell matrix, and 15N enrichment in cyanophycin inclusions was sometimes greater but more often lower than in the cell matrix (Figures 3E,F). In the morning, kC and kN were significantly correlated (R = 0.43, p = 0.045), whereas no significant correlation was observed in the afternoon (R = −0.02, p = 0.90). For both morning and afternoon incubations, kC:kN varied among individual cells from 0.15 to 0.5 (mean = 0.31, SD = 0.08, CV = 26%).
Night-time NO3 assimilation resulted in 15N being homogeneously enriched within the cell matrix, but carboxysomes were notably more enriched in 15N than the cell matrix, and cyanophycin inclusions showed no 15N enrichment (Figure 3G). This intracellular N allocation pattern was observed in all cells (Supplementary Figure 3F) despite the large intercellular variability in kN values. In contrast to the N2 culture, kC and kN in the NO3 culture were strongly correlated during the night (R = 0.94, p < 10–4; Figure 2G) even though the newly assimilated 13C and 15N were allocated into different cell compartments (Figure 3G). Night-time kC:kN values varied only slightly in the NO3 culture among individual cells (range 0.01–0.03, mean = 0.015, SD = 0.005, CV = 33%).
In the majority of Cyanothece 51142 cells grown under diazotrophic conditions in semi-continuous, fed-batch mode (“SC-N2 culture”), temporal patterns in kC and kN values as well as in C and N allocation were similar to the N2 culture grown in turbidostat mode. In these cells, the highest kC and kN values were observed during the morning and night, respectively (Figures 2H–J and Table 2). As in the N2 culture, the highest 13C and 15N enrichments were detected in the polysaccharide granules and cyanophycin inclusions, respectively (Figures 3H–J). In the afternoon, there was no measurable N2 fixation, and kC values decreased by about 80% compared to the morning values (Figure 2I).
There were, however, two notable differences in the behaviors of the N2 (turbidostat) and SC-N2 cultures. Firstly, in the SC-N2 culture there was a large subpopulation of cells (∼46%) that fixed N2 during the morning (subpopulation P6 in Figure 2H). In these asynchronous diazotrophic cells, 13C and 15N enrichments were concentrated in cyanophycin inclusions (Figure 3H), and kN and kC as well as yC and yN were significantly negatively correlated (R = −0.81, p < 10–4, Figure 2H; R = −0.36, p < 10–5, Figure 4D, respectively). Additionally, the average kN measured in these morning N2-fixing cells was 6.4-fold higher than the average kN during the night time. Consequently, the average kN for the SC-N2 culture in the morning was about 3.1-fold higher than at night (Table 2). The second notable difference was that cells from the SC-N2 culture incubated in the afternoon had some cyanophycin inclusions that were significantly more enriched in 13C compared to the polysaccharide granules and cell matrix (Figures 3I, 4F).
In all cultures and incubations, kC and kN were markedly heterogeneous among individual cells (Figure 2). Intercellular heterogeneity in day-time kC values was similar with CV = 16–22% (Figure 5) across the turbidostat cultures (N2 and NO3 cultures). Heterogeneity in night-time kN values was also similar between the turbidostat cultures but was about 3-fold greater (CV = 53–57%) than heterogeneity in day-time kC values. Heterogeneity in day-time kN in NO3 cultures decreased from the morning to the afternoon (CV declined from ∼30 to 15%). Moreover, heterogeneity in kN in the N2 culture during the night time was higher than in the NO3 culture during the morning.
Figure 5. Intercellular variability of C and N assimilation rates in Cyanothece 51142. Shown are the coefficients of variation (CV) of C-specific (A) and N-specific (B) rates rates of 13C and 15N assimilation for whole cells, kC and kN, derived from data shown in Figure 2.
The semi-continuous diazotrophic batch culture (SC-N2) showed considerably greater intercellular heterogeneity than the turbidostat cultures (N2 and NO3), especially for rates of N2 fixation (Figure 5). In the SC-N2 culture, the differentiation of cells into subpopulations in the morning was reflected in the high heterogeneity in kC (CV = 85%) and kN (CV = 155%). In the afternoon, heterogeneity in kC in the SC-N2 culture decreased to a level similar to the N2 culture (CV = 20%).
Interpreting Isotopic Enrichment Imaging Data
The data yielded from nanoSIMS analyses hold valuable information about metabolic strategies used by cells across time and space. However, knowledge about cells’ activities that influence their isotopic composition at the end of an SIP incubation is critical to properly interpret nanoSIMS data. Carbon and nitrogen assimilation in Cyanothece mainly occurred during short and intensive periods either in the few hours after dawn or during the night. These periods of rapid assimilation activities fueled the majority of the cells’ C and N needs for growth and were followed by long periods of very low assimilation rates. Our results also show a wide range of C- and N-specific rates of 13C and 15N assimilation within and between whole cells and among inclusions. Here, we critically evaluate the factors that can lead to variations in these measured rates and discuss several important, and to our knowledge previously unrecognized, considerations for using the spatially resolved 13C and 15N enrichment data obtained by nanoSIMS to infer rates of substrate assimilation. We first focus on principles of labeling as applied to subcellular structures, such as inclusion bodies, and then discuss these principles in the context of whole cell assimilation rates. A more comprehensive and mathematical analysis of these considerations can be found in Polerecky et al. (In Revision).
13C or 15N enrichment depends upon the amount of labeled C or N added to a structure during an incubation relative to the amount of unlabeled C or N present prior to the incubation. Any structure that is newly synthesized during an incubation will have 13C and 15N enrichments that match those of the enriched C and N sources. However, C or N that is added to an existing structure during the incubation will cause the average 13C and 15N enrichment measured in the structure to be lower than in the C and N sources. The deviation between structure enrichment and source enrichment will decrease with incubation time as a function of the rate of biosynthesis but increase with the initial C and N content of the structure. Consequently, variation in the initial C and N content of structures will lead to apparent differences in 13C and 15N enrichments among structures present in an incubation even though the rates of biosynthesis and accumulation of 13C and 15N may have been the same.
13C or 15N enrichment of a structure will also vary if the C and N used for its biosynthesis are derived from unlabeled sources of C and N, including the turnover of cellular macromolecules, in addition to the labeled sources external to the cell. One way to detect the relative importance of internal macromolecular recycling is to quantify the ratio of C- and N-specific rates of 13C and 15N incorporation into a structure (e.g., yC:yN for a cyanophycin granule). Because the C:N ratio of many compounds comprising cell structures is well defined (e.g., cyanophycin has a C:N ratio of two), their biosynthesis will preserve the 13C:15N signature of the C and N sources (except for the minute deviations linked to kinetic isotope fractionation). Consequently, if only the external, labeled pools of C and N are utilized for biosynthesis, the ratio of the C- and N-specific rates of 13C and 15N assimilation must be equal to 1. Any departure of this ratio from 1 implies that some of the C or N in the structure originated from a different source (i.e., with a different 13C:15N signature than that of the externally supplied sources).
Similar reasoning is needed when analyzing the isotopic enrichment of whole cells. Average 13C or 15N enrichment of a cell depends upon the amount of labeled C or N taken up during an incubation relative to the amount of unlabeled C or N present in the cell prior to the incubation. Consequently, intercellular variability in the amounts of C and N storage compounds will lead to apparent differences in the cellular 13C and 15N enrichments among cells even if the rates of 13C and 15N assimilation into cells were same. The influence of varying storage compound content can be revealed by analyzing the ratio between the C- and N-specific rates of 13C and 15N assimilation into whole cells, kC:kN. For an individual cell, this ratio will be 1 (or very close to 1, if the subtle effects of kinetic isotope fractionation are included) provided the cell is in balanced growth, and the externally supplied labeled pools of C and N were the only sources of C and N assimilated by the cell. Any deviation from 1 indicates that (1) the cell assimilated an additional, unlabeled external source; (2) a storage product was preferentially synthesized over another (e.g., polysaccharides over cyanophycin); or (3) the cell recycled internal, unlabeled stores of C or N.
The foregoing analysis of enrichment sources, biosynthesis, and initial C and N content highlights that kC:kN can reveal the presence, synthesis, or mobilization of intracellular C and N stores. With these ideas in mind, we evaluate the isotope enrichment results for Cyanothece 51142 cultures grown with different N sources across a day–night cycle to understand their C and N assimilation processes and allocation strategies.
Roles of Internal C and N Recycling in Cyanophycin Synthesis
The majority of C used for cyanophycin synthesis at night in cells growing diazotrophically originated from recycling of existing C compounds within the cells. Cyanophycin synthesis involved some 13C (external C source) assimilation, but the kC:kN and yC:yN values were much lower than 1 (Figures 2C, 4A,B). Thus, the bulk of the CO2 incorporated into cyanophycin was likely derived from polysaccharide catabolism needed to simultaneously provide energy and reducing power (ATP and NADPH) for N2-fixation. Nevertheless, cyanophycin synthesis was detected via their enrichment in 13C rather than 15N in about 17% of cells (Figure 3C). 13C is assimilated via two CO2 fixation steps leading to synthesis of the non-ribosomal peptide, cyanophycin, which is comprised of aspartate and arginine (Flores et al., 2019). Specifically, CO2 is incorporated via (i) pyruvate carboxylase or phosphoenolpyruvate carboxykinase yielding oxaloacetate which is transaminated by glutamate to form aspartate and (ii) carbamoyl phosphate synthase together with ornithine transcarbamoylase operating to generate arginine (Zhang et al., 2018).
The differences in yC:yN values between and within N2-fixing cells in the same culture (Figures 3C, 4A,B) were also caused by variations in the relative contributions to cyanophycin synthesis of 15N assimilated during the incubation and unlabeled N assimilated prior to the incubation. Unlabeled N may originate from efficient recycling of N in polyamines, including degradation of arginine via the recently described AgrE/PutA pathway (Burnat et al., 2019; Lee and Rhee, 2020). We hypothesize that the enrichment patterns observed in the N2 cultures at night were caused by a variable fraction of cyanophycin-N that originated from active N2 fixation (P1) or from protein degradation (P2). These findings suggest that diazotrophy demands internal N redistribution at night in all cells (as indicated by the similar 13C enrichment in cyanophycin inclusions), but distinct subpopulations emerge depending on their rates of N2-fixation.
Some cyanophycin inclusions in cells from the SC-N2 culture showed pronounced 13C but no 15N labeling during the afternoon incubation (Figures 3I, 4E,F). We speculate that these cells are part of the asynchronous diazotrophic morning subpopulation that fixed N2 into cyanophycin (P6, Figure 3H) but that later synthesized cyanophycin using internal (and unlabeled) N sources. This idea is supported by the observation that about half of the population fixed N2 in the morning and exhibited pronounced 13C enrichment in cyanophycin granules (Figures 3H,I). These observations suggest that cyanophycin synthesis in Cyanothece can occur throughout the entire light period, with cyanophycin-N derived either from N2-fixation or internal N (e.g., via protein degradation). Our high-resolution imaging shows that cyanophycin effectively collects, stores, and redistributes N to facilitate ongoing protein synthesis and catabolism.
Cyanophycin Biosynthesis Is Prioritized, but N Can Flow Directly to Protein During Growth on NO3
Whether growing diazotrophically or with NO3, Cyanothece prioritized cyanophycin synthesis when the rates of N assimilation were at their highest. When N assimilation rates were lower, most cells growing on NO3 assimilated N into the cell matrix and carboxysomes, suggesting that N was used directly for protein synthesis without prior storage in cyanophycin. Night-time NO3 uptake into carboxysomes indicates that this new N was used immediately for the synthesis of RubisCO to maintain its content throughout the diel cycle (Nassoury et al., 2001). How cyanophycin synthesis is prioritized over protein synthesis is not clear, but our data suggest that once the cell has acquired sufficient N storage into cyanophycin, this storage step can be bypassed in favor of direct incorporation into proteins.
Night-time assimilation of N in the NO3 culture was accompanied by small but detectable assimilation of CO2, indicative of pyruvate carboxylase activity during the night. Surprisingly, this new C was directed into polysaccharides (Figure 3G). Typically, the pyruvate carboxylase reaction is considered important to ensure availability of oxaloacetate for citrate synthase to initiate the TCA cycle that produces NADH and amino acid precursors. However, in NO3-grown cells, existing polysaccharides appear to have supplied all of the C for protein synthesis (because no 13C was incorporated into the carboxysomes coincident with 15N; see above). These results suggest that the amphibolic nature of the glycolytic/gluconeogenic pathway is directional with respect to the flow of newly fixed C into polysaccharides: newly fixed C flows through gluconeogenesis into polysaccharides simultaneous to glycolytic catabolism of “old C” stored in polysaccharides for use in protein synthesis. Glycolytic and gluconeogenesis pathways are strictly controlled so that they cannot be both highly active at the same time, which would create a futile cycle. The highly sensitive detection of labeled C and N afforded by the stable isotope probing and nanoSIMS technologies combined with TEM allowed us to view these unexpected cell activities that occur at very low levels but that support the careful modulation of C and N storage and re-mobilization in Cyanothece.
Intercellular Heterogeneity in C and N Metabolism
Within-population heterogeneity in 13C and 15N enrichments (kC and kN), such as that observed in N2 and SC-N2 treatments, has been reported in previous nanoSIMS-based studies (Foster et al., 2013; Mohr et al., 2013; Masuda et al., 2020), although the causes remain poorly understood (Ackermann, 2015). Intercellular heterogeneity has been attributed to stochastic gene expression or state switching in fast growing bacteria and yeast (Elowitz et al., 2002; Blake et al., 2003; Raj and van Oudenaarden, 2008; Raser and O’Shea, 2013; Sanchez et al., 2013; Damodaran et al., 2015). Cyanothece metabolism is strongly regulated by circadian rhythms, and thus, the intercellular heterogeneity observed in our study is at least partly associated with the regulation of C and N fixation determined by the light period and cell cycle (Caudron and Barral, 2013; Bach and Taucher, 2019).
We find that the metabolism of internal C and N storage compounds is another mechanism contributing to cell-to-cell heterogeneity in isotopic enrichment. Dual-label stable isotope probing combined with sub-cellular resolution imaging enabled us to identify internal recycling of N during cyanophycin synthesis, which led to variation in kN values during night-time N2 fixation (CV≈55%; Figure 2C) and to the wide ranges of yN (CV≈72%; Figure 4B) and yC:yN (Figures 4A,B) values.
The intercellular heterogeneity in kC values during morning C fixation could be caused by the variable polysaccharide content among cells (see Section “Interpreting Isotopic Enrichment Imaging Data”). The limited cellular volume probed by the nanoSIMS measurement also likely contributes to an apparent population heterogeneity. These alternatives are supported by the variable content and uneven distribution of polysaccharide granules within cells (Figures 1, 3A,B,E–G) and the fact that a large fraction of the variability in kC values (∼50%) was explained by the areal coverage by polysaccharide granules (Supplementary Figure 1B).
Differences in the turnover rates of storage inclusions with different C:N contents (cyanophycin vs. polysaccharides) may underlie the three-fold mismatch between the cell-to-cell variation in kN and kC. Other cyanobacteria and pico-eukaryotes have exhibited similar differences in kN and kC (Berthelot et al., 2019; Masuda et al., 2020). One explanation is that day-time acquisition of C reserves was insufficient to fuel N-fixation and other night-time metabolisms (Dron et al., 2013). However, pre-dawn cells were never completely depleted of polysaccharide granules (data not shown), making it unlikely that C-reserves limited N2 fixation. The large difference in cellular inclusion content also suggests that their subcellular metabolism influences kN and kC. Cyanophycin comprised ∼3% of cellular N in the N2 culture, whereas polysaccharides comprised ∼30% of cellular C (Table 1). Moreover, the early night yN values were considerably higher than the morning pC values, suggesting that the turnover rate of cyanophycin is considerably faster than polysaccharide turnover at times of highest N and C assimilation, respectively. Thus, it appears that during diazotrophy, cells retain large pools of C storage with slow turnover rates and small pools of N storage with high turnover rates to manage their C and N demands. This strategy could result in a greater range of enrichment in cyanophycin inclusions compared to the larger and less dynamic pool of polysaccharides. A high turnover rate of cyanophycin also helps explain why variation in kN values in NO3-grown cells was higher in the morning, when the cells assimilated N into cyanophycin, and lower in the afternoon, when the cyanophycin pool was bypassed. Differences in the turnover rates of N-rich proteins might also explain the large variation in kN in NO3-grown cells during the night when N assimilation again bypassed cyanophycin.
Differences in the timing of N2 fixation revealed a surprising amount of within-population cell-to-cell heterogeneity in diazotrophic cultures. While the majority of cells fixed N2 at night as expected (Mitsui et al., 1986; Gallon, 1992; Tuit et al., 2004; Wilson et al., 2017), subpopulations in both the turbidostat-grown and semicontinuous batch cultures fixed N2 in the morning. Asynchronous diazotrophy has been suggested to occur when the amount of N2 fixed at night is insufficient to support growth in the following day (Dron et al., 2013; Rabouille et al., 2014; Rabouille and Claquin, 2016). In Cyanothece 51142, asynchronous diazotrophy coincided with the diel maxima in population-level C fixation (Figure 2), but single-cell analysis revealed that morning N2 fixation was limited to cells whose C fixation rates, and thus presumably intracellular O2 concentrations, were low compared to the rest of the population. This behavior may be associated with prolonged deactivation of PSII through the early morning hours (Rabouille and Claquin, 2016). How these activities are regulated is not yet known.
Semi-continuous, fed-batch cultures are exposed to a wider range of nutrient and light concentrations compared to turbidostat cultures. These variations could result in a greater range of cell physiologies within a population. Together with previous reports of N2 fixation in UCYN occurring during a subjective dark phase under continuous light (Colón-López and Sherman, 1998; Pennebaker et al., 2010; Dron et al., 2013), our findings suggest that the timing of N2 fixation is not only regulated by the circadian rhythm or light/dark cycle but also by the cell’s ability to balance N and light energy demands. The greater range of heterogeneity within the SC-N2 culture compared to the N2 culture is also consistent with the idea that cell-to-cell metabolic heterogeneity facilitates rapid population adjustment to environmental changes (Ackermann, 2015; Schreiber et al., 2016) such as those present in coastal environments (Rippka, 1988).
Dual labeling combined with nanoSIMS imaging enabled a much richer and more complex view of cell activities than previously observed using measurements of bulk activities. Specifically, we observed significant cell-to-cell variation, which we attribute to differences in (1) the degree to which internal storage compounds are used as sources of C and N for cyanophycin synthesis, (2) the turnover rates of different storage pools, (3) the range of environmental conditions experienced by a population over a day–night cycle, and (4) the timing of N2 fixation. The intercellular heterogeneity potentially reflects adaptive mechanisms that allow Cyanothece to thrive in dynamic environments.
Additional details of C and N metabolism were also elucidated by evaluation of 13C and 15N labeling patterns across the day–night cycle. Cyanophycin synthesis is a highly effective N-scavenging pathway that assimilates N from protein degradation as well as external sources (NO3 or N2). Whether growing on N2 or NO3, cells prioritize cyanophycin synthesis when N assimilation rates are highest. In NO3-growing cells, N assimilation switches from cyanophycin synthesis to RubisCO synthesis, suggesting that there is a cyanophycin requirement that, once met, can be bypassed in favor of protein synthesis. In NO3-grown cells, night-time CO2 was assimilated into polysaccharides simultaneous with catabolism of polysaccharides used for protein synthesis, suggesting that one way these cells control C is to maintain a directional flow of new carbon entering the cell: CO2 → gluconeogenesis → polysaccharides → glycolysis → protein.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
OP organized the experimental part of this study conducted during the 10th Group for Aquatic Productivity (GAP) workshop in August 2017. LP, TM, ME, and KH designed the study. TM, ME, and LP performed the SIP experiment. MK and LP performed the nanoSIMS analysis. MV performed the TEM analysis. LP, KH, ME, and TM drafted the manuscript, and all authors provided input during writing of the manuscript. All authors contributed to sampling and data interpretation.
The NanoSIMS facility at Utrecht University was financed through a large infrastructure grant by the Netherlands Organization for Scientific Research (NWO; grant no. 175.010.2009.011). Electron microscopy was supported by the Czech Ministry of Education, Youth and Sport through the BioImaging Research Infrastructure project LM2015062. GB was supported by the Hungarian Academic Foundation (OTKA, grant no. K128950). TZ and JČ were supported by the Ministry of Education, Youth and Sports of the Czech Republic (OP RDE grant number CZ.02.1.01/0.0/0.0/16–026/0008413) and by the Czech Science Foundation (GA CR, grant number 18–24397S). SR acknowledges support from the LEFE-INSU funding program. TM, OP, and ME were supported by the Czech Science Foundation (GACR 20-17627S and 20-02827Y). SS was supported by the ILES (Illuminating Lake Ecosystem) project funded through the Leibniz Competition (SAW-2015-IGB-1).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work results from the 10th Group for Aquatic Photosynthesis (GAP) workshop held in Třeboň, Czech Republic, in August 2017 and subsequent discussions. With this contribution, we wish to honor the late MG and JK for their collegiality, friendship, and extensive contributions to the field of phytoplankton ecophysiology.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.620915/full#supplementary-material
Aryal, U. K., Stöckel, J., Krovvidi, R. K., Gritsenko, M. A., Monroe, M. E., Moore, R. J., et al. (2011). Dynamic proteomic profiling of a unicellular cyanobacterium Cyanothece ATCC 51142 across light-dark diurnal cycles. BMC Syst. Biol. 5:194. doi: 10.1186/1752-0509-5-194
Bach, L. T., and Taucher, J. (2019). CO2 effects on diatoms: a synthesis of more than a decade of ocean acidification experiments with natural communities. Ocean Sci. 15, 1159–1175. doi: 10.5194/os-15-1159-2019
Bernstein, H. C., Charania, M. A., McClure, R. S., Sadler, N. C., Melnicki, M. R., Hill, E. A., et al. (2015). Multi-omic dynamics associate oxygenic photosynthesis with nitrogenase-mediated H2 production in Cyanothece sp. ATCC 51142. Sci. Rep. 5:16004. doi: 10.1038/srep16004
Berthelot, H., Duhamel, S., L’Helguen, S., Maguer, J.-F., Wang, S., Cetinić, I., et al. (2019). NanoSIMS single cell analyses reveal the contrasting nitrogen sources for small phytoplankton. ISME J. 13, 651–662. doi: 10.1038/s41396-018-0285-288
Bonnet, S., Grosso, O., and Moutin, T. (2011). Planktonic dinitrogen fixation along a longitudinal gradient across the Mediterranean Sea during the stratified period (BOUM cruise). Biogeosciences 8, 2257–2267. doi: 10.5194/bg-8-2257-2011
Burnat, M., Picossi, S., Valladares, A., Herrero, A., and Flores, E. (2019). Catabolic pathway of arginine in anabaena involves a novel bifunctional enzyme that produces proline from arginine. Mol. Microbiol. 111, 883–897. doi: 10.1111/mmi.14203
Červený, J., Sinetova, M. A., Valledor, L., Sherman, L. A., and Nedbal, L. (2013). Ultradian metabolic rhythm in the diazotrophic cyanobacterium Cyanothece sp. ATCC 51142. Proc. Natl. Acad. Sci. U S A. 110, 13210–13215. doi: 10.1073/pnas.1301171110
Colón-López, M. S., and Sherman, L. A. (1998). Transcriptional and translational regulation of photosystem I and II genes in light-dark- and continuous-light-grown cultures of the unicellular cyanobacterium Cyanothece sp. strain ATCC 51142. J. Bacteriol. 180, 519–526. doi: 10.1128/jb.180.3.519-526.1998
Damodaran, S. P., Eberhard, S., Boitard, L., Rodriguez, J. G., Wang, Y., Bremond, N., et al. (2015). A millifluidic study of cell-to-cell heterogeneity in growth-rate and cell-division capability in populations of isogenic cells of Chlamydomonas reinhardtii. PLoS One 10:e0118987. doi: 10.1371/journal.pone.0118987
Dekaezemacker, J., and Bonnet, S. (2011). Sensitivity of N 2 fixation to combined nitrogen forms ( and ) in two strains of the marine diazotroph Crocosphaera watsonii (Cyanobacteria). Mar. Ecol. Prog. Ser. 438, 33–46. doi: 10.3354/meps09297
Deschamps, P., Colleoni, C., Nakamura, Y., Suzuki, E., Putaux, J.-L., Buleon, A., et al. (2008). Metabolic symbiosis and the birth of the plant kingdom. Mol. Biol. Evol. 25, 536–548. doi: 10.1093/molbev/msn053
Dron, A., Rabouille, S., Claquin, P., Chang, P., Raimbault, V., Talec, A., et al. (2012). Light:dark (12:12 h) quantification of carbohydrate fluxes in Crocosphaera watsonii. Aquat. Microb. Ecol. 68, 43–55. doi: 10.3354/ame01600
Dron, A., Rabouille, S., Claquin, P., Talec, A., Raimbault, V., and Sciandra, A. (2013). Photoperiod length paces the temporal orchestration of cell cycle and carbon-nitrogen metabolism in Crocosphaera watsonii. Environ. Microbiol. 15, 3292–3304. doi: 10.1111/1462-2920.12163
Eichner, M., Kranz, S. A., and Rost, B. (2014). Combined effects of different CO2 levels and N sources on the diazotrophic cyanobacterium Trichodesmium. Physiol. Plant. 152, 316–330. doi: 10.1111/ppl.12172
Foster, R. A., Sztejrenszus, S., and Kuypers, M. M. M. (2013). Measuring carbon and N2 fixation in field populations of colonial and free-living unicellular cyanobacteria using nanometer-scale secondary ion mass spectrometry1. J. Phycol. 49, 502–516. doi: 10.1111/jpy.12057
Großkopf, T., and LaRoche, J. (2012). Direct and indirect costs of dinitrogen fixation in Crocosphaera watsonii WH8501 and possible implications for the nitrogen cycle. Front. Microbiol. 3:236. doi: 10.3389/fmicb.2012.00236
Holl, C. M., and Montoya, J. P. (2005). Interactions between nitrate uptake and nitrogen fixation in continuous cultures of the marine diazotroph Trichodesmium (Cyanobacteria). J. Phycol. 41, 1178–1183. doi: 10.1111/j.1529-8817.2005.00146.x
Inomura, K., Deutsch, C., Wilson, S. T., Masuda, T., Lawrenz, E., Bučinská, L., et al. (2019). Quantifying oxygen management and temperature and light dependencies of nitrogen fixation by Crocosphaera watsonii. mSphere 4:e00531-19. doi: 10.1128/msphere.00531-519
Karl, D., Michaels, A., Bergman, B., Capone, D., Carpenter, E., Letelier, R., et al. (2002). “Dinitrogen fixation in the world’s oceans,” in The Nitrogen Cycle at Regional to Global Scales, eds E. W. Boyer and R. W. Howarth (Dordrecht: Springer), doi: 10.1007/978-94-017-3405-9_2
Lee, H., and Rhee, S. (2020). Structural and mutational analyses of the bifunctional arginine dihydrolase and ornithine cyclodeaminase AgrE from the cyanobacterium Anabaena. J. Biol. Chem. 295, 5751–5760. doi: 10.1074/jbc.RA120.012768
Li, H., Sherman, D. M., Bao, S., and Sherman, L. A. (2001). Pattern of cyanophycin accumulation in nitrogen-fixing and non-nitrogen-fixing cyanobacteria. Arch. Microbiol. 176, 9–18. doi: 10.1007/s002030100281
Mareš, J., Johansen, J. R., Hauer, T., Zima, J., Ventura, S., Cuzman, O., et al. (2019). Taxonomic resolution of the genus Cyanothece (Chroococcales, Cyanobacteria), with a treatment on Gloeothece and three new genera, Crocosphaera, Rippkaea, and Zehria. J. Phycol. 55, 578–610. doi: 10.1111/jpy.12853
Masuda, T., Furuya, K., Kodama, T., Takeda, S., and Harrison, P. J. (2013). Ammonium uptake and dinitrogen fixation by the unicellular nanocyanobacterium Crocosphaera watsonii in nitrogen-limited continuous cultures. Limnol. Oceanogr. 58, 2029–2036. doi: 10.4319/lo.2013.58.6.2029
Masuda, T., Inomura, K., Takahata, N., Shiozaki, T., Sano, Y., Deutsch, C., et al. (2020). Heterogeneous nitrogen fixation rates confer energetic advantage and expanded ecological niche of unicellular diazotroph populations. Commun. Biol. 3:172. doi: 10.1038/s42003-020-0894-894
Messineo, L. (1966). Modification of the Sakaguchi reaction: spectrophotometric determination of arginine in proteins without previous hydrolysis. Arch. Biochem. Biophys. 117, 534–540. doi: 10.1016/0003-9861(66)90094-90094
Mitsui, A., Kumazawa, S., Takahashi, A., Ikemoto, H., Cao, S., and Arai, T. (1986). Strategy by which nitrogen-fixing unicellular cyanobacteria grow photoautotrophically. Nature 323, 720–722. doi: 10.1038/323720a0
Mohr, W., Intermaggio, M. P., and LaRoche, J. (2010). Diel rhythm of nitrogen and carbon metabolism in the unicellular, diazotrophic cyanobacterium Crocosphaera watsonii WH8501. Environ. Microbiol. 12, 412–421. doi: 10.1111/j.1462-2920.2009.02078.x
Mohr, W., Vagner, T., Kuypers, M. M. M., Ackermann, M., and LaRoche, J. (2013). Resolution of conflicting signals at the single-cell level in the regulation of cyanobacterial photosynthesis and nitrogen fixation. PLoS One 8:e66060. doi: 10.1371/journal.pone.0066060
Montoya, J. P., Holl, C. M., Zehr, J. P., Hansen, A., Villareal, T. A., and Capone, D. G. (2004). High rates of N2 fixation by unicellular diazotrophs in the oligotrophic Pacific Ocean. Nature 430, 1027–1031. doi: 10.1038/nature02824
Mulholland, M. R., Ohki, K., and Capone, D. G. (2001). Nutrient controls on nitrogen uptake and metabolism by natural populations and cultures of Trichodesmium (Cyanobacteria). J. Phycol. 37, 1001–1009. doi: 10.1046/j.1529-8817.2001.00080.x
Nassoury, N., Fritz, L., and Morse, D. (2001). Circadian changes in ribulose-1,5-bisphosphate carboxylase/oxygenase distribution inside individual chloroplasts can account for the rhythm in dinoflagellate carbon fixation. Plant Cell 13, 923–934. doi: 10.2307/3871349
Pennebaker, K., Mackey, K. R. M., Smith, R. M., Williams, S. B., and Zehr, J. P. (2010). Diel cycling of DNA staining and nifH gene regulation in the unicellular cyanobacterium Crocosphaera watsonii strain WH 8501 (Cyanophyta). Environ. Microbiol. 12, 1001–1010. doi: 10.1111/j.1462-2920.2010.02144.x
Polerecky, L., Adam, B., Milucka, J., Musat, N., Vagner, T., and Kuypers, M. M. M. (2012). Look@NanoSIMS - a tool for the analysis of nanoSIMS data in environmental microbiology. Environ. Microbiol. 14, 1009–1023. doi: 10.1111/j.1462-2920.2011.02681.x
Polerecky, L., Eichner, M., Masuda, T., Zavřel, T., Rabouille, S., Campbell, D. A., et al. (In Revision). Calculation and Interpretation of Substrate Assimilation Rates in Microbial Cells Based on Isotopic Composition Data Obtained by Nanosims.
Rabouille, S., and Claquin, P. (2016). Photosystem-II shutdown evolved with Nitrogen fixation in the unicellular diazotroph Crocosphaera watsonii. Environ. Microbiol. 18, 477–485. doi: 10.1111/1462-2920.13157
Rabouille, S., Van de Waal, D. B., Matthijs, H. C. P., and Huisman, J. (2014). Nitrogen fixation and respiratory electron transport in the cyanobacterium Cyanothece under different light/dark cycles. FEMS Microbiol. Ecol. 87, 630–638. doi: 10.1111/1574-6941.12251
Reddy, K. J., Haskell, J. B., Sherman, D. M., and Sherman, L. A. (1993). Unicellular, aerobic nitrogen-fixing cyanobacteria of the genus Cyanothece. J. Bacteriol. 175, 1284–1292. doi: 10.1128/jb.175.5.1284-1292.1993
Schneegurt, M. A., Sherman, D. M., Nayar, S., and Sherman, L. A. (1994). Oscillating behavior of carbohydrate granule formation and dinitrogen fixation in the cyanobacterium Cyanothece sp. strain ATCC 51142. J. Bacteriol. 176, 1586–1597. doi: 10.1128/jb.176.6.1586-1597.1994
Schreiber, F., Littmann, S., Lavik, G., Escrig, S., Meibom, A., Kuypers, M. M. M., et al. (2016). Phenotypic heterogeneity driven by nutrient limitation promotes growth in fluctuating environments. Nat. Microbiol. 1:16055. doi: 10.1038/nmicrobiol.2016.55
Sherman, L. A., Meunier, P., and Colón-López, M. S. (1998). Diurnal rhythms in metabolism: a day in the life of a unicellular, diazotrophic cyanobacterium. Photosynth. Res. 58, 25–42. doi: 10.1023/A:1006137605802
Stöckel, J., Welsh, E. A., Liberton, M., Kunnvakkam, R., Aurora, R., and Pakrasi, H. B. (2008). Global transcriptomic analysis of Cyanothece 51142 reveals robust diurnal oscillation of central metabolic processes. Proc. Natl. Acad. Sci. U S A. 105, 6156–6161. doi: 10.1073/pnas.0711068105
Toepel, J., Welsh, E., Summerfield, T. C., Pakrasi, H. B., and 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. doi: 10.1128/JB.00206-208
Webb, E. A., Ehrenreich, I. M., Brown, S. L., Valois, F. W., and Waterbury, J. B. (2009). Phenotypic and genotypic characterization of multiple strains of the diazotrophic cyanobacterium, Crocosphaera watsonii, isolated from the open ocean. Environ. Microbiol. 11, 338–348. doi: 10.1111/j.1462-2920.2008.01771.x
Welsh, E. A., Liberton, M., Stöckel, J., Loh, T., Elvitigala, T., Wang, C., et al. (2008). The genome of Cyanothece 51142, a unicellular diazotrophic cyanobacterium important in the marine nitrogen cycle. Proc. Natl. Acad. Sci. U S A. 105, 15094–15099. doi: 10.1073/pnas.0805418105
Wilson, S. T., Aylward, F. O., Ribalet, F., Barone, B., Casey, J. R., Connell, P. E., et al. (2017). Coordinated regulation of growth, activity and transcription in natural populations of the unicellular nitrogen-fixing cyanobacterium Crocosphaera. Nat. Microbiol. 2:17118. doi: 10.1038/nmicrobiol.2017.118
Zavřel, T., Sinetova, M. A., Búzová, D., Literáková, P., and Červený, J. (2015a). Characterization of a model cyanobacterium Synechocystis sp. PCC 6803 autotrophic growth in a flat-panel photobioreactor. Eng. Life Sci. 15, 122–132. doi: 10.1002/elsc.201300165
Zehr, J. P., Waterbury, J. B., Turner, P. J., Montoya, J. P., Omoregie, E., Steward, G. F., et al. (2001). Unicellular cyanobacteria fix N2 in the subtropical north Pacific Ocean. Nature 412, 635–638. doi: 10.1038/35088063
Keywords: Crocosphaera subtropica (former Cyanothece sp. ATCC 51142), Cyanothece, photosynthesis, carbon fixation, nitrogen fixation, nanoSIMS, TEM
Citation: Polerecky L, Masuda T, Eichner M, Rabouille S, Vancová M, Kienhuis MVM, Bernát G, Bonomi-Barufi J, Campbell DA, Claquin P, Červený J, Giordano M, Kotabová E, Kromkamp J, Lombardi AT, Lukeš M, Prášil O, Stephan S, Suggett D, Zavřel T and Halsey KH (2021) Temporal Patterns and Intra- and Inter-Cellular Variability in Carbon and Nitrogen Assimilation by the Unicellular Cyanobacterium Cyanothece sp. ATCC 51142. Front. Microbiol. 12:620915. doi: 10.3389/fmicb.2021.620915
Received: 24 October 2020; Accepted: 11 January 2021;
Published: 04 February 2021.
Edited by:Petra M. Visser, University of Amsterdam, Netherlands
Reviewed by:Anne W. Thompson, Portland State University, United States
Deng Liu, Washington University in St. Louis, United States
Enrique Flores, Consejo Superior de Investigaciones Científicas (CSIC), Spain
Copyright © 2021 Polerecky, Masuda, Eichner, Rabouille, Vancová, Kienhuis, Bernát, Bonomi-Barufi, Campbell, Claquin, Červený, Giordano, Kotabová, Kromkamp, Lombardi, Lukeš, Prášil, Stephan, Suggett, Zavřel and Halsey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.