Introduction
Magnetotactic bacteria (MTB) biomineralize magnetosomes, which are membrane-bounded magnetic iron oxide (magnetite, Fe3O4) or iron sulphide (greigite, Fe3S4) nanoparticles intracellularly (within their cells) (Bazylinski and Frankel, Reference Bazylinski and Frankel2004; Faivre and Schüler, Reference Faivre and Schüler2008; Komeili et al., Reference Komeili2012). Typically, in the cell, magnetosomes are arranged in one or more chains and/ or in preferred orientations and provide the cell with a magnetic dipole moment that enables magnetotaxis; i.e., passive orientation and then active swimming of the cell along the geomagnetic field lines of the Earth (Frankel et al., Reference Frankel, Bazylinski, Johnson and Taylor1997). This magnetic sensing mechanism is thought to aid the bacteria in locating and maintaining their optimal positions in vertical chemical concentration gradients (usually O2 gradients) that occur in their habitats. Each individual bacterial strain biomineralizes magnetic nanoparticles that have specific sizes, shapes and arrangements. By studying the magnetic behaviour of magnetosomes in appropriately chosen strains, competing effects that influence the magnetic dipole moments and magnetic microstructures of individual nanoparticles and their assemblages can be assessed and understood. The influence of nanoparticle size on magnetic domain character, the influence of shape and magnetocrystalline anisotropy on the direction of the magnetic induction within, between and surrounding the particles, the changes that are imposed on the magnetic behaviour of the particles by their interactions and the effects of structural imperfections on local magnetic properties can all be studied.
Cells of MTB have been used as natural nanoscale laboratories for the development and testing of novel magnetic imaging techniques (Dunin-Borkowski et al., Reference Dunin-Borkowski1998; Lam et al., Reference Lam2010; Le Sage et al., Reference Le Sage2013; Proksch et al., Reference Proksch1995; Staniland et al., Reference Staniland, Ward, Harrison, Van Der Laan and Telling2007) and have been proposed for applications such as genetically engineered magnonic devices (Zingsem et al., Reference Zingsem, Feggeler, Terwey, Ghaisari, Spoddig, Faivre, Meckenstock, Farle and Winklhofer2019). Off-axis electron holography (EH) is a transmission electron microscopy transmission electron microscopy (TEM) method (Kovács and Dunin-Borkowski, Reference Kovács and Dunin-Borkowski2018) that can provide images of the magnetic microstructures and quantitative measurements of the magnetic dipole moments of magnetosomes with the highest spatial resolution available to date (Pósfai and Dunin-Borkowski, Reference Pósfai and Dunin-Borkowski2009; Pósfai et al., Reference Pósfai, Kasama, Dunin-Borkowski, Nieto and Livi2013a).
In terms of both magnetism and crystal growth, one of the most intriguing aspects of the properties of magnetosomes is their well-defined, highly constrained and, in many cases, unusual crystal morphology. In general, three basic types of magnetosome habits have been described: cuboctahedral, elongated prismatic and highly elongated bullet-shaped (Pósfai et al., Reference Pósfai, Lefévre, Trubitsyn, Bazylinski and Frankel2013b). With regard to the relationship between crystallography and magnetism, there are two primary possibilities in magnetite-producing MTBs: (i) the <111> magnetocrystalline easy axis is parallel to the chain axis and (ii) a magnetocrystalline hard axis is parallel to the chain axis. These two scenarios are regarded as having cooperative and competitive anisotropy, respectively (Charilaou, Reference Charilaou2017).
EH observations of the magnetic behaviours of magnetite magnetosomes with cuboctahedral (Dunin-Borkowski et al., Reference Dunin-Borkowski1998; Dunin-Borkowski et al., Reference Dunin-Borkowski2001 and <111>-elongated prismatic (McCartney et al., Reference McCartney, Lins, Farina, Buseck and Frankel2001; Simpson et al., Reference Simpson2005) shapes in spirilla and cocci, respectively, are consistent with expectations, showing that the direction of magnetic induction within the studied magnetite particles was typically parallel to the <111> easy axis. However, since the <111> crystallographic direction coincides with both the morphological elongation axis of the magnetosome crystals and the chain axis in all known bacterial strains that produce prismatic magnetosomes, previous EH measurements have provided relatively little information about the competing effects of magnetic anisotropies associated with particle shapes and their crystallographic orientations.
Bullet-shaped magnetosome magnetite crystals provide the largest known deviations from the equilibrium morphologies of magnetite crystals. In addition, their morphological elongation axes do not necessarily coincide with the crystallographic <111> direction. Strains with magnetosome elongation directions parallel to <100> (Byrne et al., Reference Byrne2010; Li, et al., Reference Li2015; Li et al., Reference Li2010; Pósfai et al., Reference Pósfai2006), <110> (Lefèvre et al., Reference Lefèvre2011b) and <112> (Mann et al., Reference Mann, Sparks and Blakemore1987) have been described (Fig. 1), although some of these identifications remain tentative. A list of magnetosome elongation directions in different MTB strains is given in Table 1 (Pósfai et al., Reference Pósfai, Lefévre, Trubitsyn, Bazylinski and Frankel2013b).
The magnetic properties of individual bullet-shaped magnetosomes have not previously been studied in depth. EH has been used in combination with electron tomography to analyze the magnetic induction of a single chain of bullet-shaped, <111>-elongated magnetite magnetosomes in an uncultivated MTB from the River Seine in Paris (Pósfai et al., Reference Pósfai, Kasama, Dunin-Borkowski, Nieto and Livi2013a). In the latter work, the magnetic induction was found to be parallel to the elongation axes of the individual magnetosomes. For <100>-elongated magnetite magnetosomes, EH data are available from a small number of particles in an MTB that produced mostly greigite magnetosomes (Kasama et al., Reference Kasama2006). An experimental magnetic induction map has also been published for the uncultured strain MYR-1, which produces kinked but mostly <100>-elongated bullet-shaped magnetite magnetosomes (Li et al., Reference Li2015). In these studies, the magnetic induction was parallel to the elongation axis of each magnetosome instead of the [111] crystallographic easy axis. To the best of our knowledge, the magnetic behaviour of <110>-elongated bullet-shaped magnetosomes has not previously been studied.
Studies of the magnetic properties of magnetosomes not only provide fundamental insight into the magnetism of nanoparticle systems but are necessary to provide an understanding of magnetotaxis and its evolution. According to genomic analyses, MTB are widely dispersed in 16 phyla in the phylogenetic tree of the domain Bacteria (Goswami et al., Reference Goswami, He, Li, Pan, Roberts and Lin2022). However, strains that produce bullet-shaped magnetosomes have only been observed in MTB of the phyla Thermodesulfobacteriota, Desulfobacterota, Nitrospirota, Omnitrophota and Elusimicrobiota (Pósfai et al., Reference Pósfai, Lefévre, Trubitsyn, Bazylinski and Frankel2013b; Goswami et al., Reference Goswami, He, Li, Pan, Roberts and Lin2022; Uzun et al., Reference Uzun, Koziaeva, Dziuba, Alekseeva, Krutkina, Sukhacheva, Baslerov and Grouzdev2023), all of which are thought to represent ancient lineages of prokaryotes. On the assumption that magnetotaxis evolved by linear descent and not by horizontal gene transfer, the first magnetosomes were probably bullet-shaped (Lefèvre et al., Reference Lefèvre2013). By understanding the magnetic properties of bullet-shaped magnetosomes, it may, therefore, be possible to obtain information about the evolution of MTB, as well as about any evolutionary benefit/advantage from magnetotaxis that results from producing magnetite crystals with non-equilibrium morphologies.
The present study aims to analyze the magnetic properties of <100>- and <110>-elongated bullet-shaped magnetite magnetosomes and their chains. We combine quantitative magnetic imaging, electron diffraction, electron tomography and high-resolution TEM to obtain information about the structures, orientations and morphologies of chains and individual magnetosomes. We provide high-spatial-resolution and quantitative analyses of the magnetic induction and magnetization of the magnetosomes using both model-based and model-independent approaches. The results are discussed in the context of both nanoparticle magnetism and the biological functions of biomineralized magnetosomes using analytical approaches and micromagnetic modelling.
Experimental
We studied magnetosomes from the following three strains of MTB (Table 1). Desulfovibrio magneticus strain RS-1, which is affiliated to phylum Thermodesulfobacteriota, is available in pure culture and produces <100>-elongated magnetite nanoparticles (Pósfai et al., Reference Pósfai2006; Rahn-Lee et al., Reference Rahn-Lee2015; Sakaguchi et al., Reference Sakaguchi, Arakaki and Matsunaga2002). The other two strains are uncultivated and belong to the Nitrospirota phylum. Cells of the Candidatus Magnetoovum mohavensis strain LO-1, collected from sediments taken from Lake Mead, Nevada, contained <100>-elongated bullet-shaped magnetosomes that were arranged in three twisted chains in each cell (Lefèvre et al., Reference Lefèvre, Frankel, Abreu, Lins and Bazylinski2011a; Lefèvre et al., Reference Lefèvre2011b). The Candidatus Thermomagnetovibrio paiutensis strain, HSMV-1, is a moderately thermophilic vibrio that was collected from thermal springs in northern Nevada (Lefèvre et al., Reference Lefèvre2010; Lefèvre et al., Reference Lefèvre2011b) and contained <110>-elongated magnetosomes that are arranged in a single chain in each cell.
Cells of strain RS-1 were cultivated in a 6 L bioreactor (Labofors 3, Infors, Bottmingen/Switzerland) and grown in a modified version of medium 896, described in the German Collection of Microorganism and Cell Cultures GmbH (DSMZ). Magnetosome purification was performed by using a modified version of the protocol described by Mériaux et al. (Reference Mériaux, Boucher, Marty, Lalatonne, Prévéral, Motte, Lefèvre, Geffroy, Lethimonnier, Péan, Garcia, Adryanczyk-Perrier, Pignol and Ginet2015). Late exponential phase cultures (DO600nm = 1.0) were harvested by centrifugation (7500 G, 10 min, 4 °C), yielding a pellet of 7.2 g. The cells were resuspended in 40 ml of the purification buffer consisting of 20 mM HEPES (C8H18N2O4S, pH 7.5), 0.9% NaCl, 8% glycerol, 1 mM EDTA (C10H16N2O8) and a protease inhibitor cocktail. The cells were disrupted three times by the use of a French press (1000 PSI, 4 °C) and the magnetosomes were purified magnetically by placing strong magnets at the bottom of 3 × 15 ml Falcon tubes containing the suspensions. After 15 h at 4 °C, the magnetosomes had aggregated at the bottom of the tubes (against where the magnets were placed) and the supernate (soluble fraction) was carefully removed. The magnetosome pellets were pooled and resuspended in 2 ml of purification buffer without EDTA, then magnetic separation (2–3 h at 4 °C) was repeated five times in 2 ml Eppendorf tubes. The magnetosomes were finally resuspended in 200 μl of 20 mM HEPES (pH 7.5) and 8% glycerol, flash-frozen in liquid nitrogen and stored at -80 °C. The purified magnetosomes were deposited onto TEM grids. For strains LO-1 and HSMV-1, magnetically concentrated and purified cells from environmental samples were deposited directly on C-coated TEM grids and air-dried. In this way, the original arrangements of magnetosomes in chains within cells of LO-1 and HSMV-1 were mostly preserved. By contrast, different arrangements of magnetosomes from strain RS-1 on TEM grids are likely to have resulted from the self-organization of the nanoparticles and, therefore, do not represent the original spatial relationships within the cells.
Bright-field (BF) TEM images, high-angle annular dark-field (HAADF) scanning TEM (STEM) images and energy dispersive X-ray spectroscopy (EDXS) compositional maps (not shown here) were recorded using an electron probe-aberration-corrected FEI Titan ChemiSTEM instrument operated at 200 kV. For tomographic reconstruction, a series of BF and HAADF STEM images were recorded every 2º over a sample tilt range of ±64º using a single tilt Fischione analytical tomography holder in a ThermoFisher Scientific (TFS) Spectra instrument operated at 300 kV. Image alignment and visualization were performed using TFS Inspect3D and Avizo software. Off-axis EH experiments were carried out in aberration-corrected Lorentz mode in an FEI Titan 60-300 TEM, operated at 300 kV using a Fischione dual-axis tomography TEM specimen holder. The biprism voltage used for off-axis EH was typically ~150 V, resulting in a holographic interference fringe spacing of ~2.5 nm and interference fringe contrast in vacuum of ~20%. Electron holograms were recorded on a 2k × 2k CCD camera and processed using computer scripts written in the Semper image processing language (Saxton et al., Reference Saxton, Pitt and Horner1979), MATLAB and Python-based software. The magnetic contribution to the phase shift measured using EH was separated from the mean inner potential contribution to the phase by taking differences between the phases of pairs of holograms, between which the magnetization direction in the specimen was reversed in situ in the TEM by tilting the sample to ±75° and switching on the ~1.5 T magnetic field of the conventional microscope objective lens before returning to magnetic-field-free conditions and tilting the sample back to 0° for EH acquisition (Dunin-Borkowski et al., Reference Dunin-Borkowski1998). Care was taken to ensure that the magnetization direction in the sample reversed precisely and that the diffracting conditions of the crystals did not change in each pair of holograms.
Results
The morphologies of the magnetosome magnetite crystals in strains LO-1 and HSMV-1 were determined from BF TEM images of particles recorded at different specimen tilt angles, resulting in tentative morphological models (Lefèvre et al., Reference Lefèvre2011b). Although these models are based on combinations of cubic crystal forms (Fig. 1), the outlines of the particles suggest that the elongated crystals are bounded by irregular, curved surfaces parallel to their elongation directions. The ‘base’ of each crystal in RS-1 and LO-1 is, however, usually faceted with clearly-developed octahedral {111} faces, as shown in Fig. 1(b). In HSMV-1, the base of each bullet-shaped magnetosome crystal appears to be irregular and perpendicular to the long axis of the crystal, as shown in Fig. 1(c). In each strain, many of the crystals are bent. This feature is not accounted for in the models shown in Fig. 1. A two-step growth process has been suggested for similar bullet-shaped magnetite magnetosome crystals in strain MYR-1, involving the initial formation of a cuboctahedral nucleus followed by anisotropic growth along [100] (Li et al., Reference Li2015).
Figure 2 shows results taken from an analysis of the structural characteristics of magnetosome crystals from strain RS-1. The magnetosomes are either scattered or self-organized on the TEM grid in linear chains or ring-like arrangements. The TEM images revealed amorphous shell contrast around each crystal (Figs 2(a), 2(b) and 2(d)), which may be part of the membrane that remained after the magnetosomes were extracted from their cells. However, dedicated experiments would be required to confirm the unambiguous presence of a membrane. Atomic-resolution HAADF STEM images shown in Figs 2(b) and 2(c) reveal a defect-free magnetite structure. The inset digital diffractogram shown in Fig. 2(b) confirms the <100> elongation of the crystal. The degree of elongation of the RS-1 crystals was characterized by an aspect ratio defined as their length divided by their width. By measuring 40 crystals, the average aspect ratio for RS-1 was found to be 2.1 ± 0.4. Information about the three-dimensional shapes of the magnetosome crystals was obtained from tilt series of BF and HAADF STEM images using tomographic reconstruction. Figure 2(d) shows a snapshot of magnetosomes that are arranged in the form of a ring and a corresponding three-dimensional reconstruction, Fig. 2(e), which reveals the shapes and facets of the crystals.
The magnetic states and magnetic stray field distributions of the magnetite magnetosomes were studied using EH. Figures 3(a) and 3(b) show a BF TEM image and corresponding magnetic characterization of a chain of <100>-elongated magnetosome crystals from strain RS-1. The magnetic induction map was produced from the magnetic phase shift image by adding contour lines and colours. To a first approximation, each magnetosome can be seen to contain a single magnetic domain. The magnetic field lines are mostly parallel to both the elongation direction of each crystal and the chain axis. An exception is the second crystal from the lower end of the chain. Its elongation axis is misaligned with respect to the overall chain direction, resulting in curvature of the magnetic field lines within it due to the competing effects of shape anisotropy and magnetostatic interactions with its neighbours. As a result, the magnetic field lines are oriented diagonally in the last crystal. The magnetic contribution to the phase shift across one of the magnetosomes (second from the top in Fig. 3(b)) was analyzed further. Based on a line scan measurement (see Supplementary Material), the step in magnetic phase shift across this crystal was measured to be ~0.26 radians, which is approximately 75% of the phase shift that one would expect for a pure magnetite sphere of the same volume (Dunin-Borkowski et al., Reference Dunin-Borkowski, Kasama, Wei, Tripp, Hÿtch, Snoeck, Harrison and Putnis2004). It should be noted that the preparation process and shrinkage of the bacteria during drying can alter the proximity of the grains and hence their interactions when compared to hydrated bacteria in their natural environment.
In an intact MTB cell of strain LO-1, the <100>-elongated magnetosomes are typically arranged in a linear chain. A magnetic induction map and a corresponding BF TEM image of magnetosomes in a disordered chain configuration are presented in Fig. 3(c). Analysis of the measured magnetic field lines provides insight into the influence on the magnetic properties of the crystals of their morphologies and positions relative to the chain axis, which is approximately parallel to the direction of the saturating magnetic field. The first three crystals (going downwards from the top) lie partly on top of each other, effectively creating a magnetic state similar to a horseshoe magnet that results in a deviation of the magnetic induction from the chain axis. In the next four crystals, whose separation is relatively large, the magnetic induction lines approximately follow the elongation directions of the magnetosomes. The largest crystal (marked with an arrow) has an aspect ratio of 3.2 and an elongation axis that is parallel to the direction of the saturating field. The magnetic induction lines run approximately parallel to its long axis. They are most closely spaced in the middle of the crystal, increasing their separation at each end as a result of the projection of the in-plane component of the three-dimensional magnetic flux density in the electron beam direction. The remaining nine crystals lie side by side. The magnetic field lines connect many of the crystals to their immediate neighbours as a result of their stray fields and the magnetostatic interactions between them. The average aspect ratio of the magnetosome crystals from cells of strain LO-1 strain is measured to be 2.2 ± 0.6.
An HRTEM image (Fig. 4(a)) and a corresponding magnetic induction map (Fig. 4(b)) were recorded from an individual magnetosome from strain LO-1, whose morphology is similar to the model shown in Fig. 1(b). This crystal has a length of ~145 nm and a width of ~48 nm in its middle section, with an aspect ratio of ~3. It is elongated along <100>, as confirmed by a digital diffractogram (not shown) generated from the atomic-resolution TEM image. The dashed line in Fig. 4(a) marks a twin boundary in the magnetite structure (Devouard et al., Reference Devouard, Pósfai, Hua, Bazylinski, Frankel and Buseck1998). A few-nm-thick amorphous shell, which is visible on the surface of the magnetosome, may represent the remaining membrane material. Such a membrane may protect the crystal from oxidation, making the magnetosome a model system for studying the structure and magnetism of magnetite (Zhu et al., Reference Zhu, Kalirai, Hitchcock and Bazylinski2015). A magnetic induction map of the magnetosome (Fig. 4(b)) was recorded using off-axis EH after saturating it magnetically in a direction approximately 20º from its long axis. The magnetic field lines in the crystal confirm the presence of a single-domain magnetic state. They do not appear to follow any of the <111> magnetic easy axis directions, which are oriented approximately 55º from the long axis of the magnetosome. The magnetic induction lines seem to form a bundle inside the crystal and are more widely spaced near its surface and ends, suggesting the presence of a three-dimensional magnetic state. The minimum-energy magnetic configuration of an elongated nanoparticle can be described as a flower or twisted flower state close to its ends, based on finite element micromagnetic simulations (Fabian et al., Reference Fabian, Kirschner, Williams, Heider, Leibl and Huber1996; Hertel and Kronmüller, Reference Hertel and Kronmüller2002). The step in the magnetic phase shift across the middle of the magnetosome, where the magnetic field lines are most closely spaced and approximately parallel to its elongation direction, is measured to be 0.72 radians, which is slightly smaller than the value of 0.78 radians expected for a uniformly-magnetized 48-nm-diameter magnetite crystal.
Figures 5(a) and 5(b) show a BF TEM image and the corresponding magnetic field distribution of <110>-elongated magnetosomes from strain HSMV-1, which are preserved inside a dried cell. Dark contrast around the chain in the BF TEM image is associated with the remains of the dehydrated bacterial cell, whose presence is likely to have contributed to the preservation of the magnetosomes in the form of a linear chain. Some of the crystals are bent or kinked (similar to those described by Hanzlik et al. (Reference Hanzlik, Winklhofer and Petersen2002)). The average aspect ratio of the magnetosome crystals is measured to be 2.5 ± 0.4 (c.f. RS-1 = 2.1; LO-1 = 2.2). A magnetic induction map recorded from the lower part of the chain is presented in Fig. 5(b). The magnetic field is confined to be parallel to the elongation direction of each particle, which coincides closely with the chain axis. The magnetic induction does not change significantly between the magnetosomes as the gaps between adjacent crystals are small compared to their lengths. In some cases, the magnetosomes overlap in projection. The step in the magnetic contribution to the phase shift across the <110>-elongated magnetosomes in a direction perpendicular to the chain axis was measured to be approximately (1 ± 0.2) radians, which is in the same range as that predicted for a spherical 50–60 nm magnetite crystal. As the magnetic phase shift increases rapidly from 0.63 to 1.52 radians for particles with sizes of 45 to 70 nm, errors in measurements of particle size can have a significant effect on magnetic quantification.
Figures 5(c) and 5(d) show a BF TEM image and a corresponding magnetic induction map of a section of a chain that contains six <110>-elongated magnetosomes. In this arrangement, the long axes of two of the crystals are parallel to the direction of the saturating magnetic field, while the other four crystals are oriented perpendicular to the saturating magnetic field. The two chain sections become magnetized in opposite directions after the application of a saturating field (H), although within each section the magnetization directions are parallel to the long axes of the crystals. Magnetosomes with opposing magnetic polarity are a clear artifact resulting from the external magnetic field applied perpendicular to the chain axis. In earlier studies, antiparallel sub-chain polarities were observed in strain MV-1 (Kalirai et al., Reference Kalirai, Bazylinski and Hitchcock2013).
Discussion
Magnetic induction maps have been recorded from bullet-shaped magnetite crystals with <100>, <110> and <111> crystallographic directions of elongation. Results for the first two directions are presented in this study while magnetic induction maps were previously obtained from <111>-elongated magnetosomes (Pósfai et al., Reference Pósfai, Kasama, Dunin-Borkowski, Nieto and Livi2013a) in which the magnetic induction was also measured to be parallel to the particle elongation axes. Taken together, the results show that the factors that dominate magnetic properties are the elongated shapes of the crystals and magnetostatic interactions between them when they are closely spaced. Similar conclusions were drawn from results obtained using ferromagnetic resonance spectroscopy of bullet-shaped <100>-elongated magnetosome crystals from strain RS-1 (Charilaou et al., Reference Charilaou, Rahn-Lee, Kind, García-Rubio, Komeili and Gehring2015), suggesting that shape anisotropy caused by the bullet-shaped morphology compensates for the less optimal alignment of the crystallographic magnetic easy axis. Our direct measurements, recorded from chains and individual magnetosomes, provide deeper insight into the role of elongation and the persistence of the bullet-shaped magnetosome crystal morphology.
We used experimental magnetic phase shift images to record magnetic induction maps from <100>- and <110>- elongated bullet-shaped magnetosomes extracted from MTB strains RS-1, LO-1 and HSMV-1. The phase shift images can be used to provide quantitative magnetic information about individual crystals. For an individual <100>-elongated magnetosome from strain LO-1 (Fig. 4), we first used a model-independent approach based on contour integration of the magnetic phase gradient (Beleggia et al., Reference Beleggia, Kasama and Dunin-Borkowski2010) to determine the in-plane magnetic dipole moment from the phase image. The integration contours were circles of decreasing radius R2 to R1, with the smaller circle touching the magnetosome boundary, as shown in Fig. 6(a). The measurements were extrapolated to zero integration radius, as shown in Fig. 6(b), to eliminate systematic sources of error from the measurements. The measured value of the in-plane magnetic dipole moment was 4.07 × 106 μB oriented in the direction of the black arrow marked in Fig. 6(a). By estimating the volume of the magnetosome from its dimensions measured from a BF TEM image, the average saturation magnetic induction was inferred to be 0.52 ± 0.04 T, which is slightly smaller than the expected value of 0.6 T for magnetite (Dunlop and Özdemir, Reference Dunlop and Özdemir1997). We also measured the projected in-plane magnetization distribution in the magnetosome using model-based iterative reconstruction (Caron, Reference Caron2017). We previously used this approach to measure the magnetic properties of skyrmions in FeGe (Kovács et al., Reference Kovács, Caron, Savchenko, Kiselev, Shibata, Z-A, Kanazawa, Tokura, Blügel and Dunin-Borkowski2017) and precipitates in a high entropy alloy (Lan et al., Reference Lan, Kovács, Caron, Song D, Dasari, Gwalani, Chaudhary, Ramanujan, Banerjee and Dunin-Borkowski2022). With this approach, simulated magnetic phase images are calculated based on successive guesses for the projected in-plane magnetization distribution M(x,y) in the analyzed part of the specimen. It makes use of analytical solutions for the phase shifts of simple geometrical objects, with numerical discretization performed in real space to avoid artifacts generated by discretization in Fourier space. The forward simulations are used in an iterative scheme to solve the inverse problem of reconstructing the projected in-plane magnetization distribution in the sample, which can then be converted into values of in-plane magnetization based on a measurement of the sample thickness. Figure 6(c) shows the resulting projected in-plane magnetization map of the magnetosome crystal studied in Fig. 4 and Figs 6(a) and 6(b). The projected in-plane magnetization is approximately uniform, as expected for a single magnetic domain. The inset to Fig. 6(c) shows a histogram of the inferred values of magnetization, whose average is 0.54 ± 0.08 T. Uncertainties about the crystal thickness and surface layer can contribute to both statistical and systematic errors in the inferred magnetization values. The results of the different measurements are summarized in Table 2.
Slightly smaller values of magnetic phase shift and magnetization were measured from crystals in all three studied strains than expected for pure magnetite. This deviation may be explained by uncertainties in the size and volume determination of the particles, as well as by partial oxidation from magnetite to maghemite (γ-Fe2O3), the latter having a lower value of saturation magnetization (0.48 vs 0.6 T).
We now discuss the influence of magnetosome aspect ratio on the interplay between magnetocrystalline and shape anisotropy by following a similar logic to that of Körnig et al. (Reference Körnig, Winklhofer, Baumgartner, Gonzalez, Fratzl and Faivre2014). TEM imaging of the magnetosomes allows the aspect ratios of the crystals to be determined from measurements of their lengths and widths. Values obtained for <100>- and <110>-elongated crystals from strains RS-1, LO-1 and HSMV-1, respectively, are plotted in Fig. 6(d) alongside data from strain SS-5, which produces <111>-elongated magnetite crystals (Lefèvre et al., Reference Lefèvre, Viloria, Schmidt, Pósfai, Frankel and Bazylinski2012). The average crystal aspect ratio is measured at 2.1, 2.2, 2.5 and 1.2 for strains RS-1, LO-1, HSMV-1 and SS-5, respectively. This analysis assumes that the aspect ratio does not change with crystal size. The large (>2.0) aspect ratios observed here for strains RS-1, LO-1 and HSMV-1 suggest that the studied magnetosome crystals are well-developed and mature.
By considering the saturation magnetization and magnetic anisotropy constants of magnetite, the total magnetic contribution to the Helmholtz free energy density (Ftotal) can be defined as the sum of magnetocrystalline anisotropy (Fc) and demagnetizing (Fdemag) energies. By calculating the equilibrium angle of the magnetization (θ) at which Ftotal is minimized, a threshold value for the aspect ratio can be defined, above which the magnetic moments follow the elongation direction of a crystal rather than the magnetocrystalline easy axis (as also found for strain RS-1 by Körnig et al., Reference Körnig, Winklhofer, Baumgartner, Gonzalez, Fratzl and Faivre2014). Details of the analysis are provided in the Supplementary Material. Figures 6(e) and 6(f) show plots of θ vs aspect ratio for <100> and <110> crystal elongations. Our experimental phase shift and magnetization measurements (Figs 3–6) suggest that some magnetic parameters such as saturation magnetization may be reduced slightly when compared with values for ideal magnetite particles. Accordingly, blue and orange curves are plotted in Figs 6(e) and 6(f) for pure magnetite and a slightly oxidized (maghemite) structure, respectively. The threshold values of the aspect ratio for <100> and <110> elongation directions are determined to be in the range of 1.21-1.27 and 1.11-1.14, respectively. These calculations suggest that, from a magnetic perspective, a relatively small amount (~25%) of particle elongation along a magnetic hard axis is sufficient to stabilize the direction of the magnetic dipole moment in a crystal along its elongation axis, in agreement with previous modelling and experimental studies (Körnig et al., Reference Körnig, Winklhofer, Baumgartner, Gonzalez, Fratzl and Faivre2014; Charilaou et al., Reference Charilaou, Rahn-Lee, Kind, García-Rubio, Komeili and Gehring2015; Moreno et al., Reference Moreno, Poyser, Meilak, Meo, Jenkins, Lazarov, Vallejo-Fernandez, Majetich and Evans2020).
Our analysis of individual and disordered chains of closely-spaced magnetosomes shows that the chain configuration can be the most important factor affecting the direction of magnetic induction within and surrounding the magnetite crystals. Magnetostatic interactions between nanoparticles have been shown to be important in previous analyses of magnetosome chains (Dunin-Borkowski et al., Reference Dunin-Borkowski1998; Dunin-Borkowski et al., Reference Dunin-Borkowski2001; McCartney et al., Reference McCartney, Lins, Farina, Buseck and Frankel2001). The highly elongated shapes and different crystallographic elongation directions analyzed in the present study provide direct evidence that the order of importance in the determination of the magnetic induction direction in closely-spaced magnetosomes is, typically, (1) interparticle interactions, (2) shape anisotropy and (3) magnetocrystalline anisotropy.
The highly elongated shapes of the crystals studied here are likely to be beneficial for the efficiency of magnetic orientation of the bacterial cells, since the distinct particle shapes constrain the magnetic induction to be parallel to their elongation axes, which typically coincide with the chain direction. However, our results do not provide an explanation for why there are (at least) three different crystallographic elongation axes of such crystals, including <100> and <110>, which are less likely to have been selected by evolution than <111>, which lies parallel to the crystallographic magnetic easy axis of magnetite. Our magnetic induction maps show that, if the particles are in chains, then particle shape and elongation do not play the most significant role magnetically since magnetostatic interactions between adjacent magnetosomes can then have the greatest effect on the direction of the magnetic field. Therefore, the evolutionary driving force behind the development of different crystallographic directions in bullet-shaped magnetosomes is unlikely to be associated with the optimization of magnetotaxis, but may instead result from magnetosome gene variations.
The genetic background to the formation of magnetite magnetosomes has been studied extensively over the past 20 years for two model organisms – strains Magnetospirillum gryphiswaldense MSR-1 and Magnetospirillum magneticum AMB-1 – both of which belong to the Alphaproteobacteria class in the Pseudomonadota phylum and produce cuboctahedral crystals (Komeili et al., Reference Komeili2012; Lohsse et al., Reference Lohsse [Lohße]2014). In these two closely-related strains, it was shown that magnetosome formation is controlled by magnetosome membrane proteins that are encoded by so-called mam and mms genes, which cluster in a specific chromosomic region that is referred to as the ‘magnetosome gene cluster’ (MGC). Magnetosome synthesis was shown to be affected in genetically engineered mutants that lacked specific genes or operons in the MGC, with gene deletions resulting in the formation of either smaller or structurally less perfect magnetite crystals when compared to wildtype cells, or in the inability of cells to assemble magnetosomes in chains (Lohsse et al., Reference Lohsse [Lohße]2014). On the other hand, the transfer of a minimal set of magnetosome synthesis genes from MSR-1 to non-magnetotactic bacteria resulted in the formation of magnetosome chains in these organisms (Kolinko et al., Reference Kolinko2014; Dziuba et al., Reference Dziuba, Müller, Pósfai and Schüler2023). Despite this progress in understanding the functions of magnetosome proteins, almost nothing is known about genetic control over crystal morphology since strains MSR-1 and AMB-1 biomineralize cuboctahedral, equant magnetite crystals, whose morphologies approach the equilibrium shape of magnetite and are, therefore, likely to require less morphological control.
As discussed in the Introduction, bullet-shaped magnetosomes are formed in MTB that belong to the Thermodesulfobacterota, Nitrospirota, Omnitrophota and Elusimicrobiota phyla clades, which emerged before the Pseudomonadota phylum that contains all MTB that produce prismatic and cuboctahedral magnetite crystals. Bullet-shaped magnetite producers are therefore considered to be from more ancient prokaryotic lineages (Lefèvre et al., Reference Lefèvre2013). When compared to the Alphaproteobacteria, less genetic information is available from MTB from Nitrospirota (Jogler et al., Reference Jogler2011; Kolinko et al., Reference Kolinko, Jogler, Katzmann and Schüler2011; Lin et al., Reference Lin2014), to which two of the strains analyzed in the present study (LO-1 and HSMV-1) belong. In addition to mam genes, all MTB that biomineralize bullet-shaped magnetosomes have another specific group of genes in their genomes, which is referred to as mad genes. Although little is known about the function of the proteins that are encoded by mad genes, it is possible that some of these genes are involved in controlling the shapes and crystallographic orientations of anisotropic magnetite crystals. The gene mad10, found in the MGC of Desulfamplus magnetovallimortis (Descamps et al., Reference Descamps, Monteil, Menguy, Ginet, Pignol, Bazylinski and Lefèvre2017), is a cultured MTB that can biomineralize both bullet-shaped magnetite and cuboctahedral greigite and has been proposed to be involved in the elongation of anisotropic magnetite (Pohl et al., Reference Pohl, Berger, SUllan, Valverde-Tercedor, Freindl, Spiridis, Lefèvre, Menguy, Klumpp, Blank and Faivre2019). This gene is conserved in all MTB that form bullet-shaped magnetite, including the strain RS-1. The genomes of LO-1 and HSMV-1 were not available at the time of writing.
Recent studies have provided conflicting results about magnetite biomineralization in D. magneticus strain RS-1 (Baumgartner et al., Reference Baumgartner2016; Rahn-Lee et al., Reference Rahn-Lee2015). This organism possesses a MGC that differs significantly from that in the Alphaproteobacteria (Rahn-Lee et al., Reference Rahn-Lee2015). Whereas mutants of the magnetotactic Alphaproteobacteria typically produce full chains of poorly-developed magnetite crystals, mutants of strain RS-1 synthesize fully-developed but fewer magnetosomes than their wild-type counterpart. It has been suggested that magnetosomes in RS-1 are mineralized one at a time in a ‘central magnetosome factory’ and that the crystals dissociate into chains from this point of synthesis (Rahn-Lee et al., Reference Rahn-Lee2015). Such a magnetosome growth scenario could, perhaps, explain the two-step growth process of bent, bullet-shaped magnetosome magnetite crystals (Hanzlik et al., Reference Hanzlik, Winklhofer and Petersen2002), as has been suggested for strain MYR-1 (Li et al., Reference Li2015; Li et al., Reference Li2010). Whereas the nucleation and initial growth of magnetite crystals could take place in a vesicle under strict biological control, thereby producing a euhedral ‘base’ for the crystal, subsequent growth of the pointed, elongated part could result from an anisotropic flux of Fe as the particle is ejected from the ‘magnetosome factory’. A different route to magnetite mineralization in RS-1 has been suggested by Baumgartner et al. (Reference Baumgartner2016), who reported that an iron oxyhydroxide precursor, green rust, converts to magnetite in the solid state. However, neither of these crystal growth scenarios provides a satisfactory explanation for the formation of highly specific bullet-shaped crystals or for the occurrence of different crystallographic elongations of magnetite crystals in different magnetotactic strains. Future genomic analyses may provide additional information about magnetosome chain formation to resolve the genetic origins of the different crystallographic elongation directions.
Finally, we briefly speculate here on whether bullet-shaped magnetite magnetosome crystals can be used as reliable biomarkers. Magnetite crystals similar to those produced by MTB have been preserved in the geologic record and have been termed magnetofossils (Petersen et al., Reference Petersen, von Dobeneck and Vali1986; Kopp and Kirschvink, Reference Kopp and Kirschvink2008). Since bullet-shaped crystals are considered to be from MTB of ancient evolutionary lineages, as described earlier, their presence in rocks would be indicative of, and therefore a proxy for, the past presence of MTB from these groups. In addition, the presence of these crystals may also provide information about past and present environmental conditions under which they were formed (Yamazaki et al. Reference Yamazaki, Suzuki, Kouduka and Kawamura2019; Li et al. Reference Li, Menguy, Roberts, Gu, Leroy, Bourgon, Yang, Zhao, Liu, Changela and Pan2020). Bullet-shaped magnetite crystals similar to those described here have been found in several different environments, including the chimney of a hydrothermal vent (Nakano et al., Reference Nakano, Furutani, Kato, Kouduka, Yamazaki and Suzuki2023), although the interpretation of these findings is not clear. However, some trends have been observed. For example, anisotropic (e.g., bullet-shaped) magnetite crystals predominated in more reduced conditions in some marine sediments, leading to the idea that organic C fluxes control specific MTB populations and thus magnetite crystal morphologies (Yamazaki and Kawahata, Reference Yamazaki and Kawahata1998). Many or most extant MTBs are found to grow and thrive at or close to the oxic-anoxic interface (OAI) (also known as the oxic-anoxic transition zone (OATZ)), either in sediments or in water columns (Frankel et al. Reference Frankel, Bazylinski, Johnson and Taylor1997; Bazylinski and Frankel Reference Bazylinski and Frankel2004). Yamazaki et al. (Reference Yamazaki, Suzuki, Kouduka and Kawamura2019) reported the highest number of bullet-shaped magnetite magnetofossils at the OATZ and suggested that magnetofossil morphologies can be used to detect past OATZs. Before any established conclusions can be drawn about the use of specific magnetofossil magnetite morphologies as biomarkers, a key issue that needs to be addressed is a better understanding of the environmental conditions under which such crystals are preserved. The dissolution of magnetite is known to occur under reducing conditions in sediments (i.e., in suboxic (where concentrations of dissolved O2 and sulfide are very low) or anoxic zones below the Fe-redox boundary (Canfield and Berner, Reference Canfield and Berner1987; Tostevin and Poulton, Reference Tostevin and Poulton2019; Yamazaki et al., Reference Yamazaki, Abdeldayem and Ikehara2003, Reference Yamazaki, Suzuki, Kouduka and Kawamura2019). Therefore, an understanding of magnetic mineral diagenesis after deposition in sediments is essential for any interpretation of magnetofossils as biomarkers. In a controlled study, Yamazaki et al. (Reference Yamazaki, Suzuki, Kouduka and Kawamura2019) used a sediment core taken from the Japan Sea to study magnetosome magnetic crystal dissolution through the core. Partially etched crystals were common in the reductive dissolution zone. Interestingly, with progressive dissolution, the proportion of bullet-shaped crystals decreased while hexagonal prismatic forms became more dominant, probably due to differences in resistance to dissolution among crystal planes of magnetite and differences in surface area to volume ratio (Yamazaki et al., Reference Yamazaki, Suzuki, Kouduka and Kawamura2019). Thus, processes of reductive diagenesis must be understood and considered when using magnetofossil morphology as a paleoenvironmental proxy. Their use as a biomarker for the present or past presence of specific groups of MTB may be less stringent. However, because several groups of MTB biomineralize bullet-shaped magnetite crystals, a knowledge of crystallographic details (such as those presented here) may provide further help in determining which groups of MTB are present.
Future developments in the magnetic characterization of magnetite in magnetosomes using TEM can include various experimental and methodological approaches. Intact magnetosome chains in their natural environment can be studied in liquid cell TEM specimen holders (Prozorov et al., Reference Prozorov, Almeida, Kovács and Dunin-Borkowski2017). Information about local variations in magnetic moment within individual crystals and magnetic interactions between magnetosomes can be obtained quantitatively in three dimensions by using holographic tomography (Wolf et al., Reference Wolf, Schneider, Rössler, Kovács, Schmidt, Dunin-Borkowski, Büchner, Rellinghaus and Lubk2022). The strength and direction of the magnetocrystalline anisotropy in each magnetite crystal in a chain can be probed relative to contributions from shape and magnetostatic interactions as a function of temperature, both at cryogenic and at elevated temperature around the Verwey and Curie transitions. Miniature magnetizing units built inside the objective lens area of a TEM promise to allow the use of vector magnetic fields to study magnetization rotation and reversal processes in magnetosome chains in the presence of static or rapidly varying magnetic fields, thereby providing further insight into fundamental magnetic phenomena.
In summary, quantitative magnetic imaging and measurements of the magnetic dipole moment and projected in-plane magnetization of <100>- and <110>-elongated magnetite magnetosomes from strains RS-1, LO-1 and HSMV-1 were determined using off-axis EH in the TEM. Magnetic field lines in isolated bullet-shaped magnetite magnetosomes and linear chains of crystals whose elongation directions lie along the chain axes typically follow the elongation axes of the crystals and form single magnetic domain states. However, in disordered chains, the directions of the magnetic field lines can be dominated by magnetic interactions between adjacent crystals and deviate from the elongation directions of the crystals. Quantitative analysis of our experimental results suggests that the magnetization of the studied magnetosomes may be reduced slightly from that of pure magnetite. Although our study provides new information about the magnetic properties of <100>- and <110>-elongated magnetite magnetosomes, key information about the origin of such energetically unfavourable elongation directions remains unresolved.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1180/gbi.2024.3.
Acknowledgments
M.P. acknowledges support from the National Research, Development and Innovation Office of Hungary under grant no. RRF-2.3.1-21-2022-0014. C.T.L. was funded by the French National Research Agency, ANR. Z.A.L was supported by the National Natural Science Foundation of China (Grant No.11974019). The research leading to these results has received funding from the European Union’s Horizon 2020 Research and Innovation Programme (Grant No. 856538, project ‘3D MAGiC’). Support by Yanchao Dai (ThermoFisher Scientific) in acquiring the TEM image tilt series acquisition and processing is kindly acknowledged. The authors acknowledge the thorough review by Michael Winklhofer.