Discriminating Bacterial Phenotypes at the Population and Single-Cell Level: A Comparison of Flow Cytometry and Raman Spectroscopy Fingerprinting

Real-time flow cytometry to assess qualitative and quantitative responses of oral pathobionts during exposure to antiseptics

Antiseptics are widely used in oral healthcare to prevent or treat oral diseases, such as gingivitis and periodontitis. However, the incidence of bacteria being tolerant to standard antiseptics has sharply increased over the last few years. This stresses the urgency for surveillance against tolerant organisms, as well as the discovery of novel antimicrobials. Traditionally, susceptibility to antimicrobials is assessed by broth micro-dilution or disk diffusion assays, both of which are time-consuming, labor-intensive, and provide limited information on the mode of action of the antimicrobials. The abovementioned limitations highlight the need for the development of new methods to monitor and further understand antimicrobial susceptibility. In this study, we used real-time flow cytometry, combined with membrane permeability staining, as a quick and sensitive technology to study the quantitative and qualitative responses of two oral pathobionts to different concentrations of chlorhexidine (CHX), cetylpyridinium chloride (CPC), or triclosan. Apart from the real-time monitoring of cell damage, we further applied a phenotypic fingerprinting method to differentiate between the bacterial subpopulations that arose due to treatment. We quantified the pathobiont damage rate of different antiseptics at different concentrations within 15 minutes of exposure and identified the conditions under which the bacteria were most susceptible. Moreover, we detected species-specific and treatment-specific phenotypic subpopulations. This proves that real-time flow cytometry can provide information on the susceptibility of different microorganisms in a short time frame while differentiating between antiseptics and thus could be a valuable tool in the discovery of novel antimicrobial compound, while at the same time deciphering their mode of action.

IMPORTANCE

With increasing evidence that microorganisms are becoming more tolerant to standard antimicrobials, faster and more accessible antimicrobial susceptibility testing methods are needed. However, traditional susceptibility assays are laborious and time-consuming. To overcome the abovementioned limitations, we introduce a novel approach to define antimicrobial susceptibility in a much shorter time frame with the use of real-time flow cytometry. Furthermore, phenotypic fingerprinting analysis can be applied on the data to study the way antiseptics affect the bacterial cell morphology over time and, thus, gain information on the mode of action of a certain compound.

Raman spectroscopic deep learning with signal aggregated representations for enhanced cell phenotype and signature identification

Feature representation is critical for data learning, particularly in learning spectroscopic data. Machine learning (ML) and deep learning (DL) models learn Raman spectra for rapid, nondestructive, and label-free cell phenotype identification, which facilitate diagnostic, therapeutic, forensic, and microbiological applications. But these are challenged by high-dimensional, unordered, and low-sample spectroscopic data. Here, we introduced novel 2D image-like dual signal and component aggregated representations by restructuring Raman spectra and principal components, which enables spectroscopic DL for enhanced cell phenotype and signature identification. New ConvNet models DSCARNets significantly outperformed the state-of-the-art (SOTA) ML and DL models on six benchmark datasets, mostly with >2% improvement over the SOTA performance of 85-97% accuracies. DSCARNets also performed well on four additional datasets against SOTA models of extremely high performances (>98%) and two datasets without a published supervised phenotype classification model. Explainable DSCARNets identified Raman signatures consistent with experimental indications.

Predicting the onset of preeclampsia by longitudinal monitoring of metabolic changes throughout pregnancy with Raman spectroscopy

Preeclampsia is a life-threatening pregnancy disorder. Current clinical assays cannot predict the onset of preeclampsia until the late 2nd trimester, which often leads to poor maternal and neonatal outcomes. Here we show that Raman spectroscopy combined with machine learning in pregnant patient plasma enables rapid, highly sensitive maternal metabolome screening that predicts preeclampsia as early as the 1st trimester with >82% accuracy. We identified 12, 15 and 17 statistically significant metabolites in the 1st, 2nd and 3rd trimesters, respectively. Metabolic pathway analysis shows multiple pathways corresponding to amino acids, fatty acids, retinol, and sugars are enriched in the preeclamptic cohort relative to a healthy pregnancy. Leveraging Pearson's correlation analysis, we show for the first time with Raman Spectroscopy that metabolites are associated with several clinical factors, including patients' body mass index, gestational age at delivery, history of preeclampsia, and severity of preeclampsia. We also show that protein quantification alone of proinflammatory cytokines and clinically relevant angiogenic markers are inadequate in identifying at-risk patients. Our findings demonstrate that Raman spectroscopy is a powerful tool that may complement current clinical assays in early diagnosis and in the prognosis of the severity of preeclampsia to ultimately enable comprehensive prenatal care for all patients.

Single-cell phenotypes revealed as a key biomarker in bacterial-fungal interactions: a case study of Staphylococcus and Malassezia

IMPORTANCE

Evaluating bacterial-fungal interactions is important for understanding ecological functions in a natural habitat. Many studies have defined bacterial-fungal interactions according to changes in growth rates when co-cultivated. However, the current literature lacks detailed studies on phenotypic changes in single cells associated with transcriptomic profiles to understand the bacterial-fungal interactions. In our study, we measured the single-cell phenotypes of bacteria co-cultivated with fungi using Raman spectroscopy with its transcriptomic profiles and determined the consequence of these interactions in detail. This rapid and reliable phenotyping approach has the potential to provide new insights regarding bacterial-fungal interactions.

Development of an Automated Online Flow Cytometry Method to Quantify Cell Density and Fingerprint Bacterial Communities

Cell density is an important factor in all microbiome research, where interactions are of interest. It is also the most important parameter for the operation and control of most biotechnological processes. In the past, cell density determination was often performed offline and manually, resulting in a delay between sampling and immediate data processing, preventing quick action. While there are now some online methods for rapid and automated cell density determination, they are unable to distinguish between the different cell types in bacterial communities. To address this gap, an online automated flow cytometry procedure is proposed for real-time high-resolution analysis of bacterial communities. On the one hand, it allows for the online automated calculation of cell concentrations and, on the other, for the differentiation between different cell subsets of a bacterial community. To achieve this, the OC-300 automation device (onCyt Microbiology, Zürich, Switzerland) was coupled with the flow cytometer CytoFLEX (Beckman Coulter, Brea, USA). The OC-300 performs the automatic sampling, dilution, fixation and 4',6-diamidino-2-phenylindole (DAPI) staining of a bacterial sample before sending it to the CytoFLEX for measurement. It is demonstrated that this method can reproducibly measure both cell density and fingerprint-like patterns of bacterial communities, generating suitable data for powerful automated data analysis and interpretation pipelines. In particular, the automated, high-resolution partitioning of clustered data into cell subsets opens up the possibility of correlation analysis to identify the operational or abiotic/biotic causes of community disturbances or state changes, which can influence the interaction potential of organisms in microbiomes or even affect the performance of individual organisms.

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.

Phenotypic convergence of bacterial adaption to sub-lethal antibiotic treatment

Microorganisms can adapt quickly to changes in their environment, leading to various phenotypes. The dynamic for phenotypic plasticity caused by environmental variations has not yet been fully investigated. In this study, we analyzed the time-series of phenotypic changes in Staphylococcus cells during adaptive process to antibiotics stresses using flow cytometry and Raman spectroscopy. The nine antibiotics with four different mode of actions were treated in bacterial cells at a sub-lethal concentration to give adaptable stress. Although the growth rate initially varied depending on the type of antibiotic, most samples reached the maximum growth comparable to the control through the short-term adaptation after 24 h. The phenotypic diversity, which showed remarkable changes depending on antibiotic treatment, converged identical to the control over time. In addition, the phenotype with cellular biomolecules converted into a bacterial cell that enhance tolerance to antibiotic stress with increases in cytochrome and lipid. Our findings demonstrated that the convergence into the phenotypes that enhance antibiotic tolerance in a short period when treated with sub-lethal concentrations, and highlight the feasibility of phenotypic approaches in the advanced antibiotic treatment.

High-Pressure Microfluidics for Ultra-Fast Microbial Phenotyping

Here, we present a novel methodology based on high-pressure microfluidics to rapidly perform temperature-based phenotyping of microbial strains from deep-sea environments. The main advantage concerns the multiple on-chip temperature conditions that can be achieved in a single experiment at pressures representative of the deep-sea, overcoming the conventional limitations of large-scale batch metal reactors to conduct fast screening investigations. We monitored the growth of the model strain Thermococcus barophilus over 40 temperature and pressure conditions, without any decompression, in only 1 week, whereas it takes weeks or months with conventional approaches. The results are later compared with data from the literature. An additional example is also shown for a hydrogenotrophic methanogen strain (Methanothermococcus thermolithotrophicus), demonstrating the robustness of the methodology. These microfluidic tools can be used in laboratories to accelerate characterizations of new isolated species, changing the widely accepted paradigm that high-pressure microbiology experiments are time-consuming.

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.

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.

A palette of fluorophores that are differentially accumulated by wild-type and mutant strains of Escherichia coli: surrogate ligands for profiling bacterial membrane transporters

Our previous work demonstrated that two commonly used fluorescent dyes that were accumulated by wild-type Escherichia coli MG1655 were differentially transported in single-gene knockout strains, and also that they might be used as surrogates in flow cytometric transporter assays. We summarize the desirable properties of such stains, and here survey 143 candidate dyes. We eventually triage them (on the basis of signal, accumulation levels and cost) to a palette of 39 commercially available and affordable fluorophores that are accumulated significantly by wild-type cells of the 'Keio' strain BW25113, as measured flow cytometrically. Cheminformatic analyses indicate both their similarities and their (much more considerable) structural differences. We describe the effects of pH and of the efflux pump inhibitor chlorpromazine on the accumulation of the dyes. Even the 'wild-type' MG1655 and BW25113 strains can differ significantly in their ability to take up such dyes. We illustrate the highly differential uptake of our dyes into strains with particular lesions in, or overexpressed levels of, three particular transporters or transporter components (yhjV, yihN and tolC). The relatively small collection of dyes described offers a rapid, inexpensive, convenient and informative approach to the assessment of microbial physiology and phenotyping of membrane transporter function.

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.

Raman Spectroscopy-Based Measurements of Single-Cell Phenotypic Diversity in Microbial Populations

Microbial cells that live in the same community can exist in different physiological and morphological states that change as a function of spatiotemporal variations in environmental conditions. This phenomenon is commonly known as phenotypic heterogeneity and/or diversity. Measuring this plethora of cellular expressions is needed to better understand and manage microbial processes. However, most tools to study phenotypic diversity only average the behavior of the sampled community. In this work, we present a way to quantify the phenotypic diversity of microbial samples by inferring the (bio)molecular profile of its constituent cells using Raman spectroscopy. We demonstrate how this tool can be used to quantify the phenotypic diversity that arises after the exposure of microbes to stress. Raman spectroscopy holds potential for the detection of stressed cells in bioproduction.

Microbial cells experience physiological changes due to environmental change, such as pH and temperature, the release of bactericidal agents, or nutrient limitation. This has been shown to affect community assembly and physiological processes (e.g., stress tolerance, virulence, or cellular metabolic activity). Metabolic stress is typically quantified by measuring community phenotypic properties such as biomass growth, reactive oxygen species, or cell permeability. However, bulk community measurements do not take into account single-cell phenotypic diversity, which is important for a better understanding and the subsequent management of microbial populations. Raman spectroscopy is a nondestructive alternative that provides detailed information on the biochemical makeup of each individual cell. Here, we introduce a method for describing single-cell phenotypic diversity using the Hill diversity framework of Raman spectra. Using the biomolecular profile of individual cells, we obtained a metric to compare cellular states and used it to study stress-induced changes. First, in two Escherichia coli populations either treated with ethanol or nontreated and then in two Saccharomyces cerevisiae subpopulations with either high or low expression of a stress reporter. In both cases, we were able to quantify single-cell phenotypic diversity and to discriminate metabolically stressed cells using a clustering algorithm. We also described how the lipid, protein, and nucleic acid compositions changed after the exposure to the stressor using information from the Raman spectra. Our results show that Raman spectroscopy delivers the necessary resolution to quantify phenotypic diversity within individual cells and that this information can be used to study stress-driven metabolic diversity in microbial populations.

IMPORTANCE Microbial cells that live in the same community can exist in different physiological and morphological states that change as a function of spatiotemporal variations in environmental conditions. This phenomenon is commonly known as phenotypic heterogeneity and/or diversity. Measuring this plethora of cellular expressions is needed to better understand and manage microbial processes. However, most tools to study phenotypic diversity only average the behavior of the sampled community. In this work, we present a way to quantify the phenotypic diversity of microbial samples by inferring the (bio)molecular profile of its constituent cells using Raman spectroscopy. We demonstrate how this tool can be used to quantify the phenotypic diversity that arises after the exposure of microbes to stress. Raman spectroscopy holds potential for the detection of stressed cells in bioproduction.