Multicellular magnetotactic bacterial consortia are metabolically differentiated and not clonal

Abstract

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing eight new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nano-scale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal non-canonical amino acid tagging (BONCAT) we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

Fig. 1.

Morphology and structure of MMB. (A) Cartoon depicting the morphology and internal organization of a MMB consortium. At the center of each MMB consortium lies an acellular space that is surrounded by a single layer of cells. Each cell harbors magnetosome organelles (black polygons aligned along cytoskeleton-like filaments), compartments for carbon or energy storage (gray circles), as well as other, currently unidentified structures. Scale bar ca. 1 μm. (B) Scanning electron microscopy (SEM) image of two MMB magnetically enriched from LSSM, possibly undergoing division. Scale bar, 1 μm. (C) Backscatter electron microscopy image of magnetosome chains within MMB cells (arrow). Magnetosome minerals appear to have 4–8 visible facets and are approximately 30–60 nm in diameter. Scale bar, 300 nm. Contrast and brightness of image (C) was increased for better visualization.

Fig. 2.

Genomic and phylogenetic analysis of all publicly available MMB MAGs and the 22 SCMs generated in this study. (A) Maximum-likelihood tree, inferred with FastTree, using a concatenated set of six conserved COGs (Table S3) present in all entries. Ultrafast bootstrap support values and selected genome statistics are listed. The color codes for the SCM Groups remain the same throughout all figures. (B) Average full length 16S rRNA gene identity and (C) average genome nucleotide identity heat maps of the eight newly identified MMB species compared to two available MMB reference genomes (Ca. M. multicellularis and Ca. Magnetomorum sp. HK-1). For a phylogenetic tree of all publicly available MMB 16S rRNA gene sequences, see Fig. S2. For an exhaustive sequence identity analyses of 16S rRNA and whole genomes of MMB see Figs. S3–5.

Fig. 3.

Clonality analysis of individual MMB consortia. (A) Individual reads were mapped to the same genome bin for each of the 22 SCMs. This analysis revealed that the genomes of cells within MMB consortia have a higher single nucleotide polymorphism rate (SNP expressed as Variations per kb) as compared to a clonal Pseudomonas sp. control (p < 7.3 × 10⁻⁹, n = 10, 30, 60, and 100 Pseudomonas cells) and other environmental cells (p < 2.4 × 10⁻⁶, e.g. “Other”). (B) The three sample categories showed no statistically significant difference in terms of their ratio of non-synonymous to synonymous substitutions (dN/dS). Values near 0 indicate that substitutions are neutral and there is no positive selection of the protein-coding genes in which the SNPs reside. The color of each SCM corresponds to the color identifying each unique species in Fig. 2.

Fig. 4.

Metabolic potential of the eight MMB species in LSSM. Arrows without circles indicate presence of the respective enzyme or pathway in all bins. Circles indicate complete presence (black), partial presence (gray), or missing (white) genes in each species. A full list of genes used to construct this figure can be found in Table S5.

Fig. 5.

NanoSIMS analysis of the cellular ¹³C-content of MMB consortia after in situ incubation with isotopically light or heavy carbon sources, specifically 1,2-¹³C₂-acetate, ¹³C-bicarbonate, 1,2-¹³C₂-propionate, or 1,2-¹³C_2-succinate, for 24 hours. The kill control contained magnetically enriched MMB that had been fixed in 4% paraformaldehyde prior to ¹³C-acetate addition. The negative control was sediment containing MMB without substrate addition. The dotted line shows the natural abundance of ¹³C. For further description of boxplots, see SI Appendix Text. Inset images show representative NanoSIMS hue saturated images (HSI) for each ¹³C-labeled substrate analyzed. Color scales in HSI images are 1.1% - 5% atom percent ¹³C. Scale bars are 5 μm. Fig. S1CD show the incubation setup. For a comparison of the anabolic activity of MMB groups 1, 3, and 4 see Fig. S11. Fig. S12 provides an example for correlative microscopy analysis of MMB. SI Materials and Methods detail the calculation of atom percent. For ROIs, see Fig. S13.

Fig. 6.

NanoSIMS analysis of MMB consortia incubated with 1,2-¹³C₂-acetate and ²H₂O. Hotspots within individual consortia were auto-segmented in ImageJ and the isotope ratios of hotspots compared to the value for the whole consortium and negative controls. The ¹³C and ²H hotspots showed significantly higher isotopic enrichment when compared to the values for the respective whole consortium (p <1.3×10⁻³ and <2.2×10⁻⁸, respectively), indicating they are metabolically differentiated. For further description of boxplots, see SI Appendix Text. Inset images show NanoSIMS HSI of the same MMB consortium analyzed using mass ratio ¹³C¹²C/¹²C₂ and ²H/¹H, revealing cell-to-cell differentiation. The HSI are scaled to show the atom percent of the respective isotope. For an example of the correlative microscopy workflow used to study MMB see Fig. S12. For ROIs, see Fig. S14.

Fig. 7.

Heterogeneity in anabolic activity within individual MMB consortia as revealed by BONCAT. (A) The averaged intensity profile across the diameter of 57 rotationally averaged BONCAT-labeled MMB with standard deviation shown in gray. Relative fluorescence intensity (RFI) and diameter of each MMB was scaled as a ratio (0 to 1) to account for differences in fluorescence intensity between consortia and size of consortia. The boxplots show the averaged RFI for each quarter section of the radius with a pairwise statistical difference of p < 1.0×10⁻¹⁰. For further description of boxplots, see SI Appendix Text. (B) Gray scale confocal microscopy image of a BONCAT labeled MMB showing proteins that had been synthesized over a 24-hour period. (C) Image of the MMB shown in (B) that has been rotationally averaged prior to quantification in Eman2. The red dotted line shows each quarter analyzed for the boxplots shown in (A). For raw and rotationally averaged images of all 57 MMB, see Fig. S16.