flowEMMi: an automated model-based clustering tool for microbial cytometric data

Unlocking the mechanism of action: a cost-effective flow cytometry approach for accelerating antimicrobial drug development

Antimicrobial resistance is one of the greatest challenges to global health. While the development of new antimicrobials can combat resistance, low profitability reduces the number of new compounds brought to market. Elucidating the mechanism of action is crucial for developing new antimicrobials. This can become expensive as there are no universally applicable pipelines. Phenotypic heterogeneity of microbial populations resulting from antimicrobial treatment can be captured through flow cytometric fingerprinting. Since antimicrobials are classified into limited groups, the mechanism of action of known compounds can be used for predictive modeling. We demonstrate a cost-effective flow cytometry approach for determining the mechanism of action of new compounds. Cultures of Actinomyces viscosus and Fusobacterium nucleatum were treated with different antimicrobials and measured by flow cytometry. A Gaussian mixture mask was applied over the data to construct phenotypic fingerprints. Fingerprints were used to assess statistical differences between mechanism of action groups and to train random forest classifiers. Classifiers were then used to predict the mechanism of action of cephalothin. Statistical differences were found among the different mechanisms of action groups. Pairwise comparison showed statistical differences for 35 out of 45 pairs for A. viscosus and for 32 out of 45 pairs for F. nucleatum after 3.5 h of treatment. The best-performing random forest classifier yielded a Matthews correlation coefficient of 0.92 and the mechanism of action of cephalothin could be successfully predicted. These findings suggest that flow cytometry can be a cheap and fast alternative for determining the mechanism of action of new antimicrobials.IMPORTANCEIn the context of the emerging threat of antimicrobial resistance, the development of novel antimicrobials is a commonly employed strategy to combat resistance. Elucidating the mechanism of action of novel compounds is crucial in this development but can become expensive, as no universally applicable pipelines currently exist. We present a novel flow cytometry-based approach capable of determining the mechanism of action swiftly and cost-effectively. The workflow aims to accelerate drug discovery and could help facilitate a more targeted approach for antimicrobial treatment of patients.

Single-Cell Techniques in Environmental Microbiology

Environmental microbiology has been an essential part of environmental research because it provides effective solutions to most pollutants. Hence, there is an interest in investigating microorganism behavior, such as observation, identification, isolation of pollutant degraders, and interactions between microbial species. To comprehensively understand cell heterogeneity, diverse approaches at the single-cell level are demanded. Thus far, the traditional bulk biological tools such as petri dishes are technically challenging for single cells, which could mask the heterogeneity. Single-cell technologies can reveal complex and rare cell populations by detecting heterogeneity among individual cells, which offers advantages of higher resolution, higher throughput, more accurate analysis, etc. Here, we overviewed several single-cell techniques on observation, isolation, and identification from aspects of methods and applications. Microscopic observation, sequencing identification, flow cytometric identification and isolation, Raman spectroscopy-based identification and isolation, and their applications are mainly discussed. Further development on multi-technique integrations at the single-cell level may highly advance the research progress of environmental microbiology, thereby giving more indication in the environmental microbial ecology.

Tracking defined microbial communities by multicolor flow cytometry reveals tradeoffs between productivity and diversity

Cross feeding between microbes is ubiquitous, but its impact on the diversity and productivity of microbial communities is incompletely understood. A reductionist approach using simple microbial communities has the potential to detect cross feeding interactions and their impact on ecosystem properties. However, quantifying abundance of more than two microbes in a community in a high throughput fashion requires rapid, inexpensive assays. Here, we show that multicolor flow cytometry combined with a machine learning-based classifier can rapidly quantify species abundances in simple, synthetic microbial communities. Our approach measures community structure over time and detects the exchange of metabolites in a four-member community of fluorescent Bacteroides species. Notably, we quantified species abundances in co-cultures and detected evidence of cooperation in polysaccharide processing and competition for monosaccharide utilization. We also observed that co-culturing on simple sugars, but not complex sugars, reduced microbial productivity, although less productive communities maintained higher community diversity. In summary, our multicolor flow cytometric approach presents an economical, tractable model system for microbial ecology using well-studied human bacteria. It can be extended to include additional species, evaluate more complex environments, and assay response of communities to a variety of disturbances.

INFLECT: an R-package for cytometry cluster evaluation using marker modality

Background

Current methods of high-dimensional unsupervised clustering of mass cytometry data lack means to monitor and evaluate clustering results. Whether unsupervised clustering is correct is typically evaluated by agreement with dimensionality reduction techniques or based on benchmarking with manually classified cells. The ambiguity and lack of reproducibility of sequential gating has been replaced with ambiguity in interpretation of clustering results. On the other hand, spurious overclustering of data leads to loss of statistical power. We have developed INFLECT, an R-package designed to give insight in clustering results and provide an optimal number of clusters. In our approach, a mass cytometry dataset is overclustered intentionally to ensure the smallest phenotypically different subsets are captured using FlowSOM. A range of metacluster number endpoints are generated and evaluated using marker interquartile range and distribution unimodality checks. The fraction of marker distributions that pass these checks is taken as a measure of clustering success. The fraction of unimodal distributions within metaclusters is plotted against the number of generated metaclusters and reaches a plateau of diminishing returns. The inflection point at which this occurs gives an optimal point of capturing cellular heterogeneity versus statistical power.

Results

We applied INFLECT to four publically available mass cytometry datasets of different size and number of markers. The unimodality score consistently reached a plateau, with an inflection point dependent on dataset size and number of dimensions. We tested both ConsenusClusterPlus metaclustering and hierarchical clustering. While hierarchical clustering is less computationally expensive and thus faster, it achieved similar results to ConsensusClusterPlus. The four datasets consisted of labeled data and we compared INFLECT metaclustering to published results. INFLECT identified a higher optimal number of metaclusters for all datasets. We illustrated the underlying heterogeneity within labels, showing that these labels encompass distinct types of cells.

Conclusion

INFLECT addresses a knowledge gap in high-dimensional cytometry analysis, namely assessing clustering results. This is done through monitoring marker distributions for interquartile range and unimodality across a range of metacluster numbers. The inflection point is the optimal trade-off between cellular heterogeneity and statistical power, applied in this work for FlowSOM clustering on mass cytometry datasets.

Cybergenetic control of microbial community composition

The use of bacterial communities in bioproduction instead of monocultures has potential advantages including increased productivity through division of labour, ability to utilise cheaper substrates, and robustness against perturbations. A key challenge in the application of engineered bacterial communities is the ability to reliably control the composition of the community in terms of its constituent species. This is crucial to prevent faster growing species from outcompeting others with a lower relative fitness, and to ensure that all species are present at an optimal ratio during different steps in a biotechnological process. In contrast to purely biological approaches such as synthetic quorum sensing circuits or paired auxotrophies, cybergenetic control techniques - those in which computers interface with living cells-are emerging as an alternative approach with many advantages. The community composition is measured through methods such as fluorescence intensity or flow cytometry, with measured data fed real-time into a computer. A control action is computed using a variety of possible control algorithms and then applied to the system, with actuation taking the form of chemical (e.g., inducers, nutrients) or physical (e.g., optogenetic, mechanical) inputs. Subsequent changes in composition are then measured and the cycle repeated, maintaining or driving the system to a desired state. This review discusses recent and future developments in methods for implementing cybergenetic control systems, contrasts their capabilities with those of traditional biological methods of population control, and discusses future directions and outstanding challenges for the field.

Fast quantification of gut bacterial species in cocultures using flow cytometry and supervised classification

A bottleneck for microbial community experiments with many samples and/or replicates is the fast quantification of individual taxon abundances, which is commonly achieved through sequencing marker genes such as the 16S rRNA gene. Here, we propose a new approach for high-throughput and high-quality enumeration of human gut bacteria in a defined community, combining flow cytometry and supervised classification to identify and quantify species mixed in silico and in defined communities in vitro. We identified species in a 5-species in silico community with an F1 score of 71%. In addition, we demonstrate in vitro that our method performs equally well or better than 16S rRNA gene sequencing in two-species cocultures and agrees with 16S rRNA gene sequencing data on the most abundant species in a four-species community. We found that shape and size differences alone are insufficient to distinguish species, and that it is thus necessary to exploit the multivariate nature of flow cytometry data. Finally, we observed that variability of flow cytometry data across replicates differs between gut bacterial species. In conclusion, the performance of supervised classification of gut species in flow cytometry data is species-dependent, but is for some combinations accurate enough to serve as a faster alternative to 16S rRNA gene sequencing.

Recent advances in microbial community analysis from machine learning of multiparametric flow cytometry data

Dynamic analysis of microbial composition is crucial for understanding community functioning and detecting dysbiosis. Compositional information is mostly obtained through sequencing of taxonomic markers or whole meta-genomes, which may be productively complemented by real-time quantitative community multiparametric flow cytometry data (FCM). Patterns and clusters in FCM community data can be distinguished and compared by unsupervised machine learning. Alternatively, FCM data from preselected individual strain phenotypes can be used for supervised machine-training in order to differentiate similar cell types within communities. Both types of machine learning can quantitatively deconvolute community FCM data sets and rapidly analyse global changes in response to treatment. Procedures may further be optimized for recurrent microbiome samples to simultaneously quantify physiological and compositional states.

When FLOW-FISH met FACS: Combining multiparametric, dynamic approaches for microbial single-cell research in the total environment

In environmental microbiology, the ability to assess, in a high-throughput way, single-cells within microbial communities is key to understand their heterogeneity. Fluorescence in situ hybridization (FISH) uses fluorescently labeled oligonucleotide probes to detect, identify, and quantify single cells of specific taxonomic groups. The combination of Flow Cytometry (FLOW) with FISH (FLOW-FISH) enables high-throughput quantification of complex whole cell populations, which when associated with fluorescence-activated cell sorting (FACS) enables sorting of target microorganisms. These sorted cells may be investigated in many ways, for instance opening new avenues for cytomics at a single-cell scale. In this review, an overview of FISH and FLOW methodologies is provided, addressing conventional methods, signal amplification approaches, common fluorophores for cell physiology parameters evaluation, and model variation techniques as well. The coupling of FLOW-FISH-FACS is explored in the context of different downstream applications of sorted cells. Current and emerging applications in environmental microbiology to outline the interactions and processes of complex microbial communities within soil, water, animal microbiota, polymicrobial biofilms, and food samples, are described.

Predicting the Presence and Abundance of Bacterial Taxa in Environmental Communities through Flow Cytometric Fingerprinting

Microbiome management research and applications rely on temporally resolved measurements of community composition. Current technologies to assess community composition make use of either cultivation or sequencing of genomic material, which can become time-consuming and/or laborious in case high-throughput measurements are required. Here, using data from a shrimp hatchery as an economically relevant case study, we combined 16S rRNA gene amplicon sequencing and flow cytometry data to develop a computational workflow that allows the prediction of taxon abundances based on flow cytometry measurements. The first stage of our pipeline consists of a classifier to predict the presence or absence of the taxon of interest, with yielded an average accuracy of 88.13% ± 4.78% across the top 50 operational taxonomic units (OTUs) of our data set. In the second stage, this classifier was combined with a regression model to predict the relative abundances of the taxon of interest, which yielded an average R² of 0.35 ± 0.24 across the top 50 OTUs of our data set. Application of the models to flow cytometry time series data showed that the generated models can predict the temporal dynamics of a large fraction of the investigated taxa. Using cell sorting, we validated that the model correctly associates taxa to regions in the cytometric fingerprint, where they are detected using 16S rRNA gene amplicon sequencing. Finally, we applied the approach of our pipeline to two other data sets of microbial ecosystems. This pipeline represents an addition to the expanding toolbox for flow cytometry-based monitoring of bacterial communities and complements the current plating- and marker gene-based methods.

IMPORTANCE Monitoring of microbial community composition is crucial for both microbiome management research and applications. Existing technologies, such as plating and amplicon sequencing, can become laborious and expensive when high-throughput measurements are required. In recent years, flow cytometry-based measurements of community diversity have been shown to correlate well with those derived from 16S rRNA gene amplicon sequencing in several aquatic ecosystems, suggesting that there is a link between the taxonomic community composition and phenotypic properties as derived through flow cytometry. Here, we further integrated 16S rRNA gene amplicon sequencing and flow cytometry survey data in order to construct models that enable the prediction of both the presence and the abundances of individual bacterial taxa in mixed communities using flow cytometric fingerprinting. The developed pipeline holds great potential to be integrated into routine monitoring schemes and early warning systems for biotechnological applications.

Assessing Biodegradability of Chemical Compounds from Microbial Community Growth Using Flow Cytometry

Compound biodegradability tests with natural microbial communities form an important keystone in the ecological assessment of chemicals. However, biodegradability tests are frequently limited by a singular focus either on the chemical and potential transformation products or on the individual microbial species degrading the compound. Here, we investigated a methodology to simultaneously analyze community compositional changes and biomass growth on dosed test compound from flow cytometry (FCM) data coupled to machine-learned cell type recognition. We quantified the growth of freshwater microbial communities on a range of carbon dosages of three readily biodegradable reference compounds, phenol, 1-octanol, and benzoate, in comparison to three fragrances, methyl jasmonate, myrcene, and musk xylene (as a nonbiodegradable control). Compound mass balances with between 0.1 to 10 mg C · liter-1 phenol or 1-octanol, inferred from cell numbers, parent compound analysis, and CO2 evolution, as well as use of 14C-labeled compounds, showed between 6 and 25% mg C · mg C-1 substrate incorporation into biomass within 2 to 4 days and 25 to 45% released as CO2 In contrast, similar dosage of methyl jasmonate and myrcene supported slower (4 to 10 days) and less (2.6 to 6.6% mg C · mg C-1 with 4.9 to 22% CO2) community growth. Community compositions inferred from machine-learned cell type recognition and 16S rRNA amplicon sequencing showed substrate- and concentration-dependent changes, with visible enrichment of microbial subgroups already at 0.1 mg C · liter-1 phenol and 1-octanol. In general, community compositions were similar at the start and after the stationary phase of the microbial growth, except at the highest used substrate concentrations of 100 to 1,000 mg C · liter-1 Flow cytometry cell counting coupled to deconvolution of communities into subgroups is thus suitable to infer biodegradability of organic chemicals, permitting biomass balances and near-real-time assessment of relevant subgroup changes.IMPORTANCE The manifold effects of potentially toxic compounds on microbial communities are often difficult to discern. Some compounds may be transformed or completely degraded by few or multiple strains in the community, whereas others may present inhibitory effects. In this study, we benchmark a new method based on machine-learned microbial cell recognition to rapidly follow dynamic changes in aquatic communities. We further determine productive biodegradation upon dosing of a number of well-known readily biodegradable tester compounds at a variety of concentrations. Microbial community growth was quantified using flow cytometry, and the multiple cell parameters measured were used in parallel to deconvolute the community on the basis of similarity to previously standardized cell types. Biodegradation was further confirmed by chemical analysis, showing how distinct changes in specific populations correlate to degradation. The method holds great promise for near-real-time community composition changes and deduction of compound biodegradation in natural microbial communities.

PhenoGMM: Gaussian Mixture Modeling of Cytometry Data Quantifies Changes in Microbial Community Structure

Microorganisms are vital components in various ecosystems on Earth. In order to investigate the microbial diversity, researchers have largely relied on the analysis of 16S rRNA gene sequences from DNA.

Microbial flow cytometry can rapidly characterize the status of microbial communities. Upon measurement, large amounts of quantitative single-cell data are generated, which need to be analyzed appropriately. Cytometric fingerprinting approaches are often used for this purpose. Traditional approaches either require a manual annotation of regions of interest, do not fully consider the multivariate characteristics of the data, or result in many community-describing variables. To address these shortcomings, we propose an automated model-based fingerprinting approach based on Gaussian mixture models, which we call PhenoGMM. The method successfully quantifies changes in microbial community structure based on flow cytometry data, which can be expressed in terms of cytometric diversity. We evaluate the performance of PhenoGMM using data sets from both synthetic and natural ecosystems and compare the method with a generic binning fingerprinting approach. PhenoGMM supports the rapid and quantitative screening of microbial community structure and dynamics.

IMPORTANCE Microorganisms are vital components in various ecosystems on Earth. In order to investigate the microbial diversity, researchers have largely relied on the analysis of 16S rRNA gene sequences from DNA. Flow cytometry has been proposed as an alternative technology to characterize microbial community diversity and dynamics. The technology enables a fast measurement of optical properties of individual cells. So-called fingerprinting techniques are needed in order to describe microbial community diversity and dynamics based on flow cytometry data. In this work, we propose a more advanced fingerprinting strategy based on Gaussian mixture models. We evaluated our workflow on data sets from both synthetic and natural ecosystems, illustrating its general applicability for the analysis of microbial flow cytometry data. PhenoGMM supports a rapid and quantitative analysis of microbial community structure using flow cytometry.

Computational Analysis of Microbial Flow Cytometry Data

Flow cytometry is an important technology for the study of microbial communities. It grants the ability to rapidly generate phenotypic single-cell data that are both quantitative, multivariate and of high temporal resolution. The complexity and amount of data necessitate an objective and streamlined data processing workflow that extends beyond commercial instrument software. No full overview of the necessary steps regarding the computational analysis of microbial flow cytometry data currently exists. In this review, we provide an overview of the full data analysis pipeline, ranging from measurement to data interpretation, tailored toward studies in microbial ecology. At every step, we highlight computational methods that are potentially useful, for which we provide a short nontechnical description. We place this overview in the context of a number of open challenges to the field and offer further motivation for the use of standardized flow cytometry in microbial ecology research.

Bacterial mock communities as standards for reproducible cytometric microbiome analysis

Flow cytometry has recently established itself as a tool to track short-term dynamics in microbial community assembly and link those dynamics with ecological parameters. However, instrumental configurations of commercial cytometers and variability introduced through differential handling of the cells and instruments frequently cause data set variability at the single-cell level. This is especially pronounced with microorganisms, which are in the lower range of optical resolution. Although alignment beads are valuable to generally minimize instrumental noise and align overall machine settings, an artificial microbial cytometric mock community (mCMC) is mandatory for validating lab workflows and enabling comparison of data between experiments, thus representing a necessary reference standard for the reproducible cytometric characterization of microbial communities, especially in long-term studies. In this study, the mock community consisted of two Gram-positive and two Gram-negative bacterial strains, which can be assembled with respective subsets of cells, including spores, in any selected ratio or concentration. The preparation of the four strains takes a maximum of 5 d, and the stains are storable with either PFA/ethanol fixation at -20 °C or drying at 4 °C for at least 6 months. Starting from this stock, an mCMC can be assembled within 1 h. Fluorescence staining methods are presented and representatively applied with two high-resolution cell sorters and three benchtop flow cytometers. Benchmarked data sets allow the use of bioinformatic evaluation procedures to decode community behavior or convey qualified cell sorting decisions for subsequent high-resolution sequencing or proteomic routines.