INTRODUCTION
Almost all freshwater and marine environments host a broad spectrum of special kinds of bacteria, called magnetotactic bacteria (MTB) (
1), that synthesize membrane-bound, single-domain-sized (35- to 120-nm) magnetite and/or greigite crystals called magnetosomes (
2–4). Chains of magnetosomes that contain permanent magnetic moments are fixed within the cell, which enable MTB to swim along magnetic field lines (
5,
6). MTB host a wide variety of cell morphotypes, including cocci, vibrios, spirilla, ovoids, rods, and multicellular prokaryotes (
7).
MTB generally live within the oxic-anoxic transition zone, yet they are also found in anoxic environments (
8). Upon death, they become entrained in the sediment column, where eventually, through lithification, the magnetosomes become important recorders of the geomagnetic field (
9,
10). Magnetosome concentrations in sediments also serve as paleoclimate proxies (
11–14). As they contain 2 to 4% (dry weight) Fe (
15,
16), MTB could influence iron cycles in nature by acting to sequester Fe in sediments (
17). Several studies have addressed the spatiotemporal variation of MTB along vertical profiles (
16,
18,
19), yet much less is known about how their horizontal distribution varies in space and time. Specifically, we question how well MTB populations measured in a given place (a spatial niche) represent a particular environment through time; e.g., do the populations remain fixed in a limited area, move randomly from place to place, or some combination of the two, and over what time scales? We therefore carried out relatively dense spatial and temporal sampling of decameter-sized environments and quantified natural MTB populations over a period of 198 days. This enabled us to map their distribution in space and to observe how the population evolved over time.
DISCUSSION
There are few studies on the horizontal distribution of MTB in sediment microcosms over time. We found a mostly nonuniform and constantly varying horizontal distribution of MTB, consistent with data reported previously by Jogler et al. (
18), who also found uneven population distributions of MTB rods in seven habitats. Simmons et al. (
20) found sharp changes in the abundances of MTB cocci and greigite (iron sulfide)-producing MTB over small length scales in a natural pond.
MTB rods and MTB cocci were initially concentrated in the northern part of the aquaria. After aquarium A was rotated 180° on day 44, MTB rods did not flourish in the “new” northern part of the aquarium or die off toward the south. Because the sediment was mixed into the aquaria in the magnetic field, and because MTB swimming velocities are higher in water than in sediment (
21), the fact that MTB rods were concentrated in the north could indicate that pouring the sediment into the aquaria in the ambient field might have biased the initial results. On the other hand, the finding that MTB rods died off in aquarium A and not in aquarium B, whereas MTB cocci and MTB spirilla were less affected by the rotation, suggests that the change in the magnetic environment influenced the MTB rod community. Mao (
22) found that the abundance of MTB rods dropped after decreasing the ambient magnetic field to near null, whereas the MTB rods recovered to normal concentrations after the geomagnetic field was restored. Lin et al. (
23) analyzed more than 900 MTB 16S rRNA gene sequences from 25 locations around the world based on the UniFrac and Sørensen indices. They found that the geomagnetic field strength influenced MTB activity and diversity to the same extent as did salinity, sulfate, temperature, and E
h. The fact that only MTB rods showed a response to the 180° rotation of the aquarium may also be due to the fact that MTB rods have a magnetic moment that is larger than that of MTB cocci by a factor of 10 (
24). Correspondingly, a higher sensitivity of MTB rods to changes in the magnetic field could be expected.
To test whether morphotype populations were correlated with O
2 concentrations, we measured nine O
2 profiles in each aquarium on day 120 (see Fig. S1 in the supplemental material). Oxygen concentrations in water ranged from 100 to 150 μmol/liter; O
2 disappeared at a depth of 3 mm in the sediment. Pearson correlation analysis revealed no significant relationship between bacterial abundances and O
2 concentrations at a depth of 1 mm. A similar analysis based on the O
2 concentration in water just above the sediment also yielded no significant degree of correlation. MTB populations were uncorrelated with local O
2 conditions in the aquarium environment in the plane parallel to the surface. Flies et al. (
16) found that MTB distributions in the vertical dimension were restricted to a narrow sediment layer overlapping or closely below the maximum oxygen penetration depth. Different species showed various preferences within vertical gradients, but over 60% of MTB were detected within the suboxic zone, which begins ∼1 to 2 mm below the sediment-water interface in our study. Our sampling integrates a volume of
h · π ·
r2 (where
r is 2.5 mm and
h is 10 mm) spanning the oxic-anoxic transition zone where most bacteria should live. The idea that the numbers of MTB do not correlate with O
2 concentrations might not be surprising considering that O
2 restricts MTB to live in a specific region in the vertical direction but does not control the absolute numbers of MTB within that region. Our experiments suggest that the absolute numbers of MTB could vary considerably if one measured them in the horizontal plane at a given depth or constant horizon of chemical activity.
MTB cocci initially grew in the corners of the aquaria and ended up being more concentrated in the center of the aquaria toward the end of the experiment. The finding that the total number of MTB increased over time suggests that the changes in the horizontal distribution do not stem from migration but rather are due to changes in cell division (growth) and death rates, which then trended toward relatively equal numbers at the latter stages of the experiment in both aquaria. This is consistent with the observations of Mao et al. (
21), who found that MTB move via a slightly biased random walk in the external magnetic field in the sediment. While the swimming velocity of MTB cocci in water can reach 112 μm/s (
25), their swimming velocity in sediment is unknown.
Another way to quantify changes in MTB populations is to examine the relative population over time (
Fig. 3). Viewing them in this way, one sees that MTB coccus and MTB spirillum populations were generally anticorrelated at the beginning stages of the experiment. Making correlations with data from other studies needs to be done with caution, as there are multiple species of cocci and spirilla (
3,
7), and each species might respond differently to a given microenvironment.
The horizontal distributions of MTB coccus and MTB spirillum population centers were in continuous flux during incubation, and the spatial positioning of these centers was consistently anticorrelated (
Fig. 2). Initially, MTB rods and MTB spirilla predominated, but MTB cocci became the dominant group after day 37, which has been found in natural habitats (
16,
19,
26,
27). Moreover, when MTB cocci had similar spatial distributions between adjacent days of counting, MTB spirilla were concentrated in different places between those same adjacent days of counting and vice versa. Pearson correlation analysis suggests that MTB cocci and MTB spirilla were significantly anticorrelated on day 198 in aquarium B (
r = −0.523), with a lower probability of anticorrelation in aquarium A (
r = −0.301). The anticorrelation between MTB cocci and MTB spirilla leads us to conclude that these two morphotypes occupy spatially distinct niches whose horizontal positioning undergoes constant flux.
We calculated the relative growth rate (
k) of MTB by using the formula
k = (
Nn+1 −
Nn)/[
Nn × (
tn+1 −
tn)], where
Nn+1 and
Nn are the amounts of MTB at time
n + 1 (
tn+1) and time
n (
tn), respectively (
Fig. 4). MTB cocci had a higher average growth rate between days 7 and 67, whereas MTB spirilla and MTB rods in aquarium B had higher growth rates from day 120. The maximum growth rates at all 72 sites in both aquaria for MTB cocci, MTB spirilla, and MTB rods were 10, 3, and 1 day
−1, respectively. Average
k values were 0.24, 0.11, and 0.02 day
−1, respectively. The minimum growth rate was −0.03 day
−1. In comparison, the generation times for cultured magnetotactic spirilla are between 6 and 26 h (
28), which correspond to
k values of 15 and 1 day
−1 and hence are significantly higher than our average growth rates but on the same order as the maximum growth rates. Nutrient availability and/or predator-prey relationships between natural sediment and culture environments likely account for this difference.
The two aquaria experienced identical environmental conditions and parallel initial conditions. The spatial distributions of MTB cocci and MTB spirilla were fairly similar in each aquarium, independent of the rotation at day 44. Changes in the distribution of MTB cannot be related to temperature, light, or O
2 concentration, as was also concluded by Jogler et al. (
18). Lin and Pan (
26) speculated that the variation in nitrite-oxidizing and ammonia-oxidizing bacteria could change the concentration of nitrate in the sediment, which may have had an effect on MTB communities. Different areas of the sediment may contain different predators, like phages or eukaryotic grazers (
18), which could account for the uneven distribution of MTB. On the other hand, some workers found that the abundance of spherical, mulberry-like, magnetotactic multicellular prokaryotes correlated with the concentration of organic matter (
29–31). This could explain our results if a local depletion of organic matter led to population deceases, whereas population booms represent the exploitation of new sources of organic matter, although other variables such as nitrate (
26), etc., could play equally significant roles.
Sobrinho et al. (
32) found that iron and bioavailable sulfur concentrations regulate the magnetotactic multicellular prokaryote density in the Araruama Lagoon in Brazil. Other workers found that salinity, nitrate concentrations, and sulfate concentrations correlate with the amounts of MTB (
16,
19,
26,
33). We observed worms and small aquatic plants living in the sediment. Bioturbation by invertebrate burrowing as well as the roots of aquatic plants can increase the heterogeneity of microenvironments. In our particular case, the aquaria likely had an uneven distribution of nutrients.
The method of counting MTB has some limitations, as only highly motile MTB swam out of the sediment into the water drop, and only live cells were counted. Nevertheless, the relative changes that we observed should be robust, as the same method was applied across all samples and time points. The anticorrelation of MTB coccus and MTB spirillum population centers cannot be attributed solely to migration or swimming speed. Rather, changes in the growth and death rates of MTB must be contributing factors, with the population centers undergoing constant flux throughout the time period. These results suggest that many MTB lie dormant and then wake up at different times, either in a stochastic fashion (
34) or in response to newly available nutrients such as iron (
35). This appears to vary on the centimeter scale, given the anticorrelation of the fluxes in growth and death rates of the two morphotypes. Our study suggests that quantifying the presence and abundance of MTB in nature should be carried out with high-resolution sampling in space and time. The highly dynamic niche partitioning shown here implies that results based on environmental samples taken from a single point in space (e.g., for metagenomics or 16S rRNA gene sequencing) from benthic ecosystems should be interpreted with caution.