Growth characteristics of Synechocystis sp. PCC6803 and T. thermophila in monocultures
We designed a synthetic symbiosis system that is composed of a cyanobacterium (Synechocystis sp. PCC6803, an evolved strain created in the previous study 21) and a ciliated protozoan (T. thermophila). All cultures were illuminated for cyanobacterial photosynthesis and mixed well to simplify population dynamics, by ignoring the spatial structure. We used the synthetic medium TCM1Glc- (Materials and Methods and Table S1). This medium is obtained by removing glucose, a primary carbon source, from TCM1, where the cell populations of the cyanobacterium and the ciliate can grow independently. TCM1 was obtained by adding some components of BG11, which is the minimum medium for the cyanobacterium, to CDM15, which is the minimum medium for the ciliate, modified from a standard chemical defined medium 22. Glucose was omitted from TCM1, to enhance the nutrient dependency of the ciliate on the cyanobacterium.
We confirmed that the cyanobacterial cell population grew in monocultures using TCM1Glc-, when the initial cell concentration was higher than 106 cells/mL (Fig. 1) (the population increase at 106 cells/mL was significant; t-test, DF = 3, α = 0.05). This concentration-dependence has been revealed in a previous study 21, and has been shown to be due to the medium, in which some components (such as amino acids) are still toxic to the cyanobacterium (but necessary for the ciliate). The growth of the ciliate cell population also showed concentration-dependence in monocultures using TCM1Glc-. The cell concentration decreased when the initial cell concentration was not higher than 102.5 cells/mL. Such ciliate mortality at low concentrations has been reported previously 23. The cell concentration increased slightly when the initial cell concentration was 104 cells/mL (the increase was significant; t-test, DF = 3, α = 0.05), and we did not detect a significant increase or decrease in those intermediate ranges. The slight growth at 104 cells/mL in the medium without glucose could be because the ciliate used other substances, such as citrate, as the carbon source. As mentioned above, the growth of the cell populations of both the species was concentration-dependent (higher was better, unless near saturation) in monocultures.
Interactions In Cocultures
We investigated the effects of interactions between the two species. First, we found that the ciliate engulfed the cyanobacteria (Fig. 2a), as confirmed by the observation of cyanobacterial cells inside the ciliate cells.
Next, we investigated the effects of these interactions on cell population dynamics. When the initial ciliate cell concentration was 102.5 cells/mL [Fig. 2b-(i)], no significant difference was observed in the cyanobacterial cell concentrations between the cocultures (blue-closed circles at 3.5 d) and monocultures (blue-closed squares; t-test, DF = 25, α = 0.05). On the other hand, when the initial ciliate cell concentration was 103.5 cells/mL [Fig. 2b-(ii)], the cyanobacterial cell concentration in the cocultures (blue-closed circles at 3.5 d) was significantly lower than that in the monocultures (blue-closed square; t-test, DF = 17, α = 0.05), which suggests the effect of predation of the cyanobacterial cells by the ciliate cells. At the same time, the growth of the ciliate cell population was significantly greater in the cocultures (red-open circles with lines in Fig. 2b) than in the monocultures [red-open squares with dashed lines in Fig. 2b-(i) and without lines in Fig. 2b-(ii); one-tailed t-test, DF was 5 and 20, at initial concentrations of 103.5 cells/mL and 102.5 cells/mL, respectively, α = 0.05]. These results suggested that ciliates utilized cyanobacteria for their growth, through predation.
When the initial cell concentration of the ciliate was 102.5 cells/mL, the decrease in the population of ciliates that was observed in the monocultures [red-open squares with dashed lines in Fig. 2b-(i)] was found to be prevented in the cocultures [red-open circles with lines in Fig. 2b-(i)]. These results suggested that the cause of ciliate mortality at low cell concentrations that was observed in the monoculture was eliminated by cyanobacteria. Indeed, we confirmed that the supernatant of the cyanobacterium monoculture reduced ciliate mortality at low cell concentrations [black-open triangles with lines in Fig. 2b-(i)].
As mentioned above, we found the population decrease of the cyanobacterial cells and the population growth of the ciliate cells, which would be due to predation and the avoidance of ciliate mortality by the cyanobacterium. All these interactions were dependent on the cell concentrations in the cocultures as well as monocultures. The population growth of both the species was better at a higher cell concentration in the monoculture, but a higher cell concentration of the ciliate led to a decrease in the cyanobacterial cell concentration in the coculture. Therefore, it is necessary to quantitatively understand the cell population dynamics of the coculture, to determine the conditions for long-term co-existence.
Investigation Of Coculture Conditions Using A Mathematical Model Of Cell Population Dynamics
To find out the experimental conditions for the long-term co-existence of both species, i.e., symbiosis, we determined the range of initial cell concentrations of the cocultures, where the two species can be sustained throughout serially transferred cocultures. As mentioned above, changes in the cell populations of both the species depend on cell concentration. To quantitatively understand the cell population dynamics, we formulated a simple mathematical model for population dynamics, using a standard Monod function and concentration-dependent mortalities to match the experimental observations described above:
$$\frac{\text{d}{C}_{S}}{\text{d}\text{t}}={k}_{S}{C}_{S}-{d}_{S}\frac{{C}_{S}{C}_{T}}{{C}_{S}+{K}_{M}}$$
\(\frac{\text{d}{\text{C}}_{T}}{\text{d}\text{t}}={k}_{T}\frac{{C}_{S}{C}_{T}}{{C}_{S}+{K}_{M}}-{d}_{T}\frac{{C}_{T}}{{{K}_{I}C}_{S}+{C}_{T}}\)
where C, k, d are the cell concentration, rate constant of population growth, and rate constant of mortality, respectively, of the cyanobacterium Synechocystis sp. PCC6803 (subscript S) and the ciliate T. thermophila (subscript T). KM is the Monod constant of the predation and KI is the inhibition constant of ciliate mortality by the cyanobacterium.
For cyanobacterial cell population changes, we considered the effects of ciliate predation, in addition to independent growth. For ciliate cell population changes, we considered predation-dependent growth, self-concentration-dependent mortality, and inhibition of the ciliate mortality by cyanobacteria. All of these terms were assumed based on the experimental results shown in Fig. 2. This model is a simplified version, which includes only the terms necessary to investigate the target conditions of the coculture and ignores other effects such as the concentration-dependent growth of the cyanobacterium and cyanobacterium-independent slight growth of the ciliate. We confirmed that a more detailed model including these ignored terms could explain all the experiments carried out well, and both models gave similar results around the target conditions (Supplementary Fig. S1).
We tested the cocultures with various initial concentrations near the predicted boundary, where both cell populations increased. The values of the model constants were determined by fitting the model to the experimental results using the quasi-Newton method. Figure 3 shows the direction field calculated using the mathematical model, with the fitted constants (blue arrows) overlaid on the experimental results (red arrows). This model explains the experimental results well.
The model shows that there is no stable equilibrium point or orbit in this primitive symbiotic system, where the two species co-exist unless the system is externally controlled, such as by serial transfers. Therefore, it is necessary to find the condition in which both the species increase during serial transfers. We found that the cell populations of both the species increased after 3.5 d in cocultures, when the initial cell concentrations of the cyanobacterium and ciliate were 106.5–107 cells/mL and 102–103 cells/mL, respectively.
Below, we explain the interpretation of these experimental results using this model. The cyanobacterial population increases when predation by ciliates is not too high (Fig. 3, when the ciliate concentration is lower than the black-dashed line). The ciliate population increases when there is a sufficient concentration of the cyanobacterium for predation (Fig. 3, when the cyanobacterial concentration is higher than the black-solid line). Therefore, a balanced amount of predation is required, to allow for an increase in both the populations. When the initial cell population falls outside the range of the conditions obtained above, at least one population fails to increase, as follows: when the cyanobacterial concentration was higher than 107 cells/mL, ciliates grew rapidly, and their predation reduced the cyanobacterium steeply. When the cyanobacterial concentration is lower than 106.5 cells/mL, the ciliate population decreases, because predation is insufficient. When the ciliate concentration is higher than 103 cells/mL, the cyanobacterial population decreases because of excessive predation. When the ciliate concentration is lower than 102 cells/mL, there is no problem in this model of population dynamics, but extremely low concentrations are generally problematic in the practical scenario, because they are undetectable, and discreteness makes the experiments unstable. As described above, we obtained a range of possible conditions for the serially transferred cocultures and their interpretation from experiments and mathematical modeling.
Serial Transfers Of The Coculture
We demonstrated serial transfer of the coculture using the following procedure based on the above results. First, we started cocultures in which the initial cell concentration of the cyanobacterium was 106.5 cells/mL and that of the ciliate ranged from 101.75 to 103.5 cells/mL. Of these cocultures, we transferred the cultures where the ciliate population increased and the cyanobacterial cell concentration was higher than 106 cells/mL, after 3.5 d. The transfer was carried out by diluting the cocultures so that the cell concentration of the ciliate became the initial concentration of the next coculture. We did not control the initial cyanobacterial concentration at dilution. We set the initial ciliate cell concentration to have a range, because these experimental conditions are in the boundary region and unstable (as shown in Fig. 3). Because of the instability of the coculture, we increased the number of cultures at the time of transfer, so there were multiple lines in parallel.
Re-coculture Using Isolated Cells, After 101 Generations Of Coculture
The cells might have changed their growth characteristics in the 101 generations of coculture. We isolated and stocked post-transfer cyanobacterial and ciliate cells from the cocultures, after 41 transfers (referred to hereafter as evolved cells). We mixed these evolved cells with the original cells as re-cocultures. We set the initial ciliate concentration at 101.75 cells/mL, which made it difficult for the pair of the original cells to maintain both cell populations. We found that both cell populations were maintained in some of the four tested replicates of the re-cocultures, in which at least one population, either the cyanobacterium or the ciliate, was the evolved cell (Fig. 5). The total maintenance proportions at round 3 of these evolved cells-included re-cocultures (7/12) was significantly different from that of the original cocultures (0/4; in these experiments, note that both populations were not maintained for at least an additional of 4 more original cocultures, as shown below), as assessed using the two-proportions z-test (α = 0.05). It was still a small change that was difficult to investigate further, but we found that some significant changes associated with coculture growth occurred within the 101 generations.
Evolutionary changes can be interpreted by using mathematical models. For example, as in this case, stable growth at low ciliate concentrations depending on the cyanobacterium may have increased the KI in the model. From another point of view, in the nullcline (black solid line) of the ciliate in Fig. 3, the position when the ciliate concentration was low moved to the left. We can interpret this as an increase in the avoidance of ciliate mortality at low concentrations (black-open triangles with lines in Fig. 2b-(i) observed in monoculture). Thus, interpretation using the model contributes to verifying this mechanism. However, as mentioned above, this result is too small to verify. In addition, there is no stable equilibrium point in the novel symbiotic system constructed at present; however, an equilibrium point may appear with evolution. For example, considering that the concentration-dependence of the cyanobacterium in monoculture and the death of the ciliate at low concentrations are eliminated due to evolution, in addition to which the growth of the ciliate is saturated while assuming that the Hill coefficient is 2, as in the detailed model, there will be a condition where the equilibrium point appears at the center of the stable orbit in the mathematical model (Supplementary Figs. S1 and S3). As mentioned above, our mathematical model is useful not only for system construction, as in this case, but also for understanding future evolution.