Quantitation and Comparison of Phenotypic Heterogeneity Among Single Cells of Monoclonal Microbial Populations

Single-cell analysis reveals an active and heterotrophic microbiome in the Guaymas Basin deep subsurface with significant inorganic carbon fixation by heterotrophs

The marine subsurface is a long-term sink of atmospheric carbon dioxide with significant implications for climate on geologic timescales. Subsurface microbial cells can either enhance or reduce carbon sequestration in the subsurface, depending on their metabolic lifestyle. However, the activity of subsurface microbes is rarely measured. Here, we used nanoscale secondary ion mass spectrometry (nanoSIMS) to quantify anabolic activity in 3,203 individual cells from the thermally altered deep subsurface in the Guaymas Basin, Mexico (3-75 m below the seafloor, 0-14°C). We observed that a large majority of cells were active (83%-100%), although the rates of biomass generation were low, suggesting cellular maintenance rather than doubling. Mean single-cell activity decreased with increasing sediment depth and temperature and was most strongly correlated with porewater sulfate concentrations. Intracommunity heterogeneity in microbial activity decreased with increasing sediment depth and age. Using a dual-isotope labeling approach, we determined that all active cells analyzed were heterotrophic, deriving the majority of their cellular carbon from organic sources. However, we also detected inorganic carbon assimilation in these heterotrophic cells, likely via processes such as anaplerosis, and determined that inorganic carbon contributes at least 5% of the total biomass carbon in heterotrophs in this community. Our results demonstrate that the deep marine biosphere at Guaymas Basin is largely active and contributes to subsurface carbon cycling primarily by not only assimilating organic carbon but also fixing inorganic carbon. Heterotrophic assimilation of inorganic carbon may be a small yet significant and widespread underappreciated source of labile carbon in the global subsurface.

IMPORTANCE

The global subsurface is the largest reservoir of microbial life on the planet yet remains poorly characterized. The activity of life in this realm has implications for long-term elemental cycling, particularly of carbon, as well as how life survives in extreme environments. Here, we recovered cells from the deep subsurface of the Guaymas Basin and investigated the level and distribution of microbial activity, the physicochemical drivers of activity, and the relative significance of organic versus inorganic carbon to subsurface biomass. Using a sensitive single-cell assay, we found that the majority of cells are active, that activity is likely driven by the availability of energy, and that although heterotrophy is the dominant metabolism, both organic and inorganic carbon are used to generate biomass. Using a new approach, we quantified inorganic carbon assimilation by heterotrophs and highlighted the importance of this often-overlooked mode of carbon assimilation in the subsurface and beyond.

Stable Isotope Probing-nanoFTIR for Quantitation of Cellular Metabolism and Observation of Growth-Dependent Spectral Features

This study utilizes nanoscale Fourier transform infrared spectroscopy (nanoFTIR) to perform stable isotope probing (SIP) on individual bacteria cells cultured in the presence of 13C-labelled glucose. SIP-nanoFTIR simultaneously quantifies single-cell metabolism through infrared spectroscopy and acquires cellular morphological information via atomic force microscopy. The redshift of the amide I peak corresponds to the isotopic enrichment of newly synthesized proteins. These observations of single-cell translational activity are comparable to those of conventional methods, examining bulk cell numbers. Observing cells cultured under conditions of limited carbon, SIP- nanoFTIR is used to identify environmentally-induced changes in metabolic heterogeneity and cellular morphology. Individuals outcompeting their neighboring cells will likely play a disproportionately large role in shaping population dynamics during adverse conditions or environmental fluctuations. Additionally, SIP-nanoFTIR enables the spectroscopic differentiation of specific cellular growth phases. During cellular replication, subcellular isotope distribution becomes more homogenous, which is reflected in the spectroscopic features dependent on the extent of 13C-13C mode coupling or to specific isotopic symmetries within protein secondary structures. As SIP-nanoFTIR captures single-cell metabolism, environmentally-induced cellular processes, and subcellular isotope localization, this technique offers widespread applications across a variety of disciplines including microbial ecology, biophysics, biopharmaceuticals, medicinal science, and cancer research.

Gold-FISH enables targeted NanoSIMS analysis of plant-associated bacteria

Bacteria colonize plant roots and engage in reciprocal interactions with their hosts. However, the contribution of individual taxa or groups of bacteria to plant nutrition and fitness is not well characterized due to a lack of in situ evidence of bacterial activity. To address this knowledge gap, we developed an analytical approach that combines the identification and localization of individual bacteria on root surfaces via gold-based in situ hybridization with correlative NanoSIMS imaging of incorporated stable isotopes, indicative of metabolic activity. We incubated Kosakonia strain DS-1-associated, gnotobiotically grown rice plants with 15 N-N2 gas to detect in situ N2 fixation activity. Bacterial cells along the rhizoplane showed heterogeneous patterns of 15 N enrichment, ranging from the natural isotope abundance levels up to 12.07 at% 15 N (average and median of 3.36 and 2.85 at% 15 N, respectively, n = 697 cells). The presented correlative optical and chemical imaging analysis is applicable to a broad range of studies investigating plant-microbe interactions. For example, it enables verification of the in situ metabolic activity of host-associated commercialized strains or plant growth-promoting bacteria, thereby disentangling their role in plant nutrition. Such data facilitate the design of plant-microbe combinations for improvement of crop management.

Single-cell view of deep-sea microbial activity and intracommunity heterogeneity

Microbial activity in the deep sea is cumulatively important for global elemental cycling yet is difficult to quantify and characterize due to low cell density and slow growth. Here, we investigated microbial activity off the California coast, 50-4000 m water depth, using sensitive single-cell measurements of stable-isotope uptake and nucleic acid sequencing. We observed the highest yet reported proportion of active cells in the bathypelagic (up to 78%) and calculated that deep-sea cells (200-4000 m) are responsible for up to 34% of total microbial biomass synthesis in the water column. More cells assimilated nitrogen derived from amino acids than ammonium, and at higher rates. Nitrogen was assimilated preferentially to carbon from amino acids in surface waters, while the reverse was true at depth. We introduce and apply the Gini coefficient, an established equality metric in economics, to quantify intracommunity heterogeneity in microbial anabolic activity. We found that heterogeneity increased with water depth, suggesting a minority of cells contribute disproportionately to total activity in the deep sea. This observation was supported by higher RNA/DNA ratios for low abundance taxa at depth. Intracommunity activity heterogeneity is a fundamental and rarely measured ecosystem parameter and may have implications for community function and resilience.

Photophysiological response of Symbiodiniaceae single cells to temperature stress

Photosynthetic dinoflagellates in the family Symbiodiniaceae engage in symbiosis with scleractinian corals. As coral ‘bleaching’ is partly governed by the thermal sensitivity of different Symbiodiniaceae lineages, numerous studies have investigated their temperature sensitivity. However, the systematic identification of single-cells with increased temperature resistance among these dinoflagellates has remained inaccessible, mostly due to a lack of technologies operating at the microscale. Here, we employed a unique combination of microfluidics, miniaturized temperature control, and chlorophyll fluorometry to characterize the single-cell heterogeneity among five representative species within the Symbiodiniaceae family under temperature stress. We monitored single-cell maximum quantum yields (F_v/F_m) of photosystem (PS) II under increasing temperature stress (22‒39 °C, + 1 °C every 15 min), and detected a significant F_v/F_m reduction at lineage-specific temperatures ranging from 28 °C to 34 °C alongside a 40- to 180- fold increase in intraspecific heterogeneity under elevated temperatures (>31 °C). We discovered that the initial F_v/F_m of a cell could predict the same cell’s ability to perform PSII photochemistry under moderate temperature stress (<32 °C), suggesting its use as a proxy for measuring the thermal sensitivity among Symbiodiniaceae. In combination, our study highlights the heterogeneous thermal sensitivity among photosynthetic Symbiodiniaceae and adds critical resolution to our understanding of temperature-induced coral bleaching.

Single-cell stable isotope probing in microbial ecology

Environmental and host-associated microbiomes are typically diverse assemblages of organisms performing myriad activities and engaging in a network of interactions that play out in spatially structured contexts. As the sum of these activities and interactions give rise to overall microbiome function, with important consequences for environmental processes and human health, elucidating specific microbial activities within complex communities is a pressing challenge. Single-cell stable isotope probing (SC-SIP) encompasses multiple techniques that typically utilize Raman microspectroscopy or nanoscale secondary ion mass spectrometry (NanoSIMS) to enable spatially resolved tracking of isotope tracers in cells, cellular components, and metabolites. SC-SIP techniques are uniquely suited for illuminating single-cell activities in microbial communities and for testing hypotheses about cellular functions generated for example from meta-omics datasets. Here, we illustrate the insights enabled by SC-SIP techniques by reviewing selected applications in microbiology and offer a perspective on their potential for future research.

NanoSIP: NanoSIMS Applications for Microbial Biology

High-resolution imaging with secondary ion mass spectrometry (nanoSIMS) has become a standard method in systems biology and environmental biogeochemistry and is broadly used to decipher ecophysiological traits of environmental microorganisms, metabolic processes in plant and animal tissues, and cross-kingdom symbioses. When combined with stable isotope-labeling-an approach we refer to as nanoSIP-nanoSIMS imaging offers a distinctive means to quantify net assimilation rates and stoichiometry of individual cell-sized particles in both low- and high-complexity environments. While the majority of nanoSIP studies in environmental and microbial biology have focused on nitrogen and carbon metabolism (using ¹⁵N and ¹³C tracers), multiple advances have pushed the capabilities of this approach in the past decade. The development of a high-brightness oxygen ion source has enabled high-resolution metal analyses that are easier to perform, allowing quantification of metal distribution in cells and environmental particles. New preparation methods, tools for automated data extraction from large data sets, and analytical approaches that push the limits of sensitivity and spatial resolution have allowed for more robust characterization of populations ranging from marine archaea to fungi and viruses. NanoSIMS studies continue to be enhanced by correlation with orthogonal imaging and 'omics approaches; when linked to molecular visualization methods, such as in situ hybridization and antibody labeling, these techniques enable in situ function to be linked to microbial identity and gene expression. Here we present an updated description of the primary materials, methods, and calculations used for nanoSIP, with an emphasis on recent advances in nanoSIMS applications, key methodological steps, and potential pitfalls.

Linguistic laws in biology

Linguistic laws, the common statistical patterns of human language, have been investigated by quantitative linguists for nearly a century. Recently, biologists from a range of disciplines have started to explore the prevalence of these laws beyond language, finding patterns consistent with linguistic laws across multiple levels of biological organisation, from molecular (genomes, genes, and proteins) to organismal (animal behaviour) to ecological (populations and ecosystems). We propose a new conceptual framework for the study of linguistic laws in biology, comprising and integrating distinct levels of analysis, from description to prediction to theory building. Adopting this framework will provide critical new insights into the fundamental rules of organisation underpinning natural systems, unifying linguistic laws and core theory in biology.

Calculation and Interpretation of Substrate Assimilation Rates in Microbial Cells Based on Isotopic Composition Data Obtained by nanoSIMS

Stable isotope probing (SIP) combined with nano-scale secondary ion mass spectrometry (nanoSIMS) is a powerful approach to quantify assimilation rates of elements such as C and N into individual microbial cells. Here, we use mathematical modeling to investigate how the derived rate estimates depend on the model used to describe substrate assimilation by a cell during a SIP incubation. We show that the most commonly used model, which is based on the simplifying assumptions of linearly increasing biomass of individual cells over time and no cell division, can yield underestimated assimilation rates when compared to rates derived from a model that accounts for cell division. This difference occurs because the isotopic labeling of a dividing cell increases more rapidly over time compared to a non-dividing cell and becomes more pronounced as the labeling increases above a threshold value that depends on the cell cycle stage of the measured cell. Based on the modeling results, we present formulae for estimating assimilation rates in cells and discuss their underlying assumptions, conditions of applicability, and implications for the interpretation of intercellular variability in assimilation rates derived from nanoSIMS data, including the impacts of storage inclusion metabolism. We offer the formulae as a Matlab script to facilitate rapid data evaluation by nanoSIMS users.

Metabolic history and metabolic fitness as drivers of anabolic heterogeneity in isogenic microbial populations

Microbial populations often display different degrees of heterogeneity in their substrate assimilation, that is, anabolic heterogeneity. It has been shown that nutrient limitations are a relevant trigger for this behaviour. Here we explore the dynamics of anabolic heterogeneity under nutrient replete conditions. We applied time-resolved stable isotope probing and nanoscale secondary ion mass spectrometry to quantify substrate assimilation by individual cells of Pseudomonas putida, P. stutzeri and Thauera aromatica. Acetate and benzoate at different concentrations were used as substrates. Anabolic heterogeneity was quantified by the cumulative differentiation tendency index. We observed two major, opposing trends of anabolic heterogeneity over time. Most often, microbial populations started as highly heterogeneous, with heterogeneity decreasing by various degrees over time. The second, less frequently observed trend, saw microbial populations starting at low or very low heterogeneity, and remaining largely stable over time. We explain these trends as an interplay of metabolic history (e.g. former growth substrate or other nutrient limitations) and metabolic fitness (i.e. the fine-tuning of metabolic pathways to process a defined growth substrate). Our results offer a new viewpoint on the intra-population functional diversification often encountered in the environment, and suggests that some microbial populations may be intrinsically heterogeneous.

Biodiversity of marine microbes is safeguarded by phenotypic heterogeneity in ecological traits

Why, contrary to theoretical predictions, do marine microbe communities harbor tremendous phenotypic heterogeneity? How can so many marine microbe species competing in the same niche coexist? We discovered a unifying explanation for both phenomena by investigating a non-cooperative game that interpolates between individual-level competitions and species-level outcomes. We identified all equilibrium strategies of the game. These strategies represent the probability distribution of competitive abilities (e.g. traits) and are characterized by maximal phenotypic heterogeneity. They are also neutral towards each other in the sense that an unlimited number of species can co-exist while competing according to the equilibrium strategies. Whereas prior theory predicts that natural selection would minimize trait variation around an optimum value, here we obtained a mathematical proof that species with maximally variable traits are those that endure. This discrepancy may reflect a disparity between predictions from models developed for larger organisms in contrast to our microbe-centric model. Rigorous mathematics proves that phenotypic heterogeneity is itself a mechanistic underpinning of microbial diversity. This discovery has fundamental ramifications for microbial ecology and may represent an adaptive reservoir sheltering biodiversity in changing environmental conditions.

Novel clades of soil biphenyl degraders revealed by integrating isotope probing, multi-omics, and single-cell analyses

Most microorganisms in the biosphere remain uncultured and poorly characterized. Although the surge in genome sequences has enabled insights into the genetic and metabolic properties of uncultured microorganisms, their physiology and ecological roles cannot be determined without direct probing of their activities in natural habitats. Here we employed an experimental framework coupling genome reconstruction and activity assays to characterize the largely uncultured microorganisms responsible for aerobic biodegradation of biphenyl as a proxy for a large class of environmental pollutants, polychlorinated biphenyls. We used ¹³C-labeled biphenyl in contaminated soils and traced the flow of pollutant-derived carbon into active cells using single-cell analyses and protein–stable isotope probing. The detection of ¹³C-enriched proteins linked biphenyl biodegradation to the uncultured Alphaproteobacteria clade UBA11222, which we found to host a distinctive biphenyl dioxygenase gene widely retrieved from contaminated environments. The same approach indicated the capacity of Azoarcus species to oxidize biphenyl and suggested similar metabolic abilities for species of Rugosibacter. Biphenyl oxidation would thus represent formerly unrecognized ecological functions of both genera. The quantitative role of these microorganisms in pollutant degradation was resolved using single-cell-based uptake measurements. Our strategy advances our understanding of microbially mediated biodegradation processes and has general application potential for elucidating the ecological roles of uncultured microorganisms in their natural habitats.

Generation of Bacterial Diversity by Segregation of DNA Strands

The generation in a bacterial population of a diversity that is coherent with present and future environments is a fundamental problem. Here, we use modeling to investigate growth rate diversity. We show that the combination of (1) association of extended assemblies of macromolecules with the DNA strands and (2) the segregation of DNA strands during cell division allows cells to generate different patterns of growth rate diversity with little effect on the overall growth rate of the population and thereby constitutes an example of “order for free” on which evolution can act.

NanoSIMS sample preparation decreases isotope enrichment: magnitude, variability and implications for single-cell rates of microbial activity

The activity of individual microorganisms can be measured within environmental samples by detecting uptake of isotope-labelled substrates using nano-scale secondary ion mass spectrometry (nanoSIMS). Recent studies have demonstrated that sample preparation can decrease ¹³ C and ¹⁵ N enrichment in bacterial cells, resulting in underestimates of activity. Here, we explore this effect with a variety of preparation types, microbial lineages and isotope labels to determine its consistency and therefore potential for correction. Specifically, we investigated the impact of different protocols for fixation, nucleic acid staining and catalysed reporter deposition fluorescence in situ hybridization (CARD-FISH) on >14 500 archaeal and bacterial cells (Methanosarcina acetivorans, Sulfolobus acidocaldarius and Pseudomonas putida) enriched in ¹³ C, ¹⁵ N, ¹⁸ O, ² H and/or ³⁴ S. We found these methods decrease isotope enrichments by up to 80% - much more than previously reported - and that the effect varies by taxa, growth phase, isotope label and applied protocol. We make recommendations for how to account for this effect experimentally and analytically. We also re-evaluate published nanoSIMS datasets and revise estimated microbial turnover times in the marine subsurface and nitrogen fixation rates in pelagic unicellular cyanobacteria. When sample preparation is accounted for, cell-specific rates increase and are more consistent with modelled and bulk rates.