Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Protocol
  • Published:

Cytometric fingerprinting for analyzing microbial intracommunity structure variation and identifying subcommunity function

Abstract

Functions of complex natural microbial communities are realized by single cells that contribute differently to the overall performance of a community. Usually, molecular and, more recently, deep-sequencing techniques are used for detailed but resource-consuming phylogenetic or functional analyses of microbial communities. Here we present a method for analyzing dynamic community structures that rapidly detects functional (rather than phylogenetic) coherent subcommunities by monitoring changes in cell-specific and abiotic microenvironmental parameters. The protocol involves the use of flow cytometry to analyze elastic light scattering and fluorescent cell labeling, with subsequent determination of cell gate abundance and finally the creation of a cytometric community fingerprint. Abiotic parameter analysis data are correlated with the dynamic cytometric fingerprint to obtain a time-bound functional heat map. The map facilitates the identification of activity hot spots in communities, which can be further resolved by subsequent cell sorting of key subcommunities and concurrent phylogenetic analysis (terminal restriction fragment length polymorphism, tRFLP). The cytometric fingerprint information is based on gate template settings and the functional heat maps are created using an R script. Cytometric fingerprinting and evaluation can be accomplished in 1 d, and additional subcommunity composition information can be obtained in a further 6 d.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Workflow for identifying dynamics in community structure and correlations in abundances of subsets of cells because of changes in abiotic parameter quantities.
Figure 2: Exemplary cytometric patterns after DAPI staining.
Figure 3: Gate template example.
Figure 4: Similarity analysis.
Figure 5: Gate cell number variation for ten measurements determined using cytometric fingerprinting.
Figure 6: Recognizing pitfalls in bar code interpretation.
Figure 7: Functional correlation analysis.
Figure 8: Cell sorting and tRFLP analysis of selected subcommunities.

Similar content being viewed by others

References

  1. Staley, J.T. & Konopka, A. Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and tertestrial habitats. Ann. Rev. Microbiol. 39, 321–346 (1985).

    Article  CAS  Google Scholar 

  2. Clardy, J. Discovery of new compounds in nature. Proc. Am. Philos. Soc. 151, 201–210 (2007).

    Google Scholar 

  3. Manz, W., Wagner, M., Amann, R. & Schleifer, K.-H. In situ characteriszation of the microbial consortia active in two wastewater treatment plants. Water Res. 28, 1715–1723 (1994).

    Article  CAS  Google Scholar 

  4. Hammes, F. et al. Flow-cytometric total bacterial cell counts as a descriptive microbiological parameter for drinking water treatment processes. Water Res. 42, 269–277 (2008).

    Article  CAS  Google Scholar 

  5. Müller, S. & Nebe-von-Caron, G. Functional single-cell analyses: flow cytometry and cell sorting of microbial populations and communities. FEMS Microbiol. Rev. 34, 554–587 (2010).

    Article  Google Scholar 

  6. Wang, Y., Hammes, F., De Roy, K., Verstraete, W. & Boon, N. Past, present and future applications of flow cytometry in aquatic microbiology. Trends Biotechnol. 28, 416–424 (2010).

    Article  CAS  Google Scholar 

  7. Günther, S. et al. Correlation of community dynamics and process parameters as a tool for the prediction of the stability of wastewater treatment. Environ. Sci. Technol. 46, 84–92 (2012).

    Article  Google Scholar 

  8. Müller, S. Modes of cytometric bacterial DNA pattern: a tool for pursuing growth. Cell Prolif. 40, 621–639 (2007).

    Article  Google Scholar 

  9. Müller, S., Harms, H. & Bley, T. Origin and analysis of microbial population heterogeneity in bioprocesses. Curr. Opin. Biotechnol. 21, 100–113 (2010).

    Article  Google Scholar 

  10. Wiacek, C., Müller, S. & Benndorf, D. A cytomic approach reveals population heterogeneity of Cupriavidus necator in response to harmful phenol concentrations. Proteomics 6, 5983–5994 (2006).

    Article  CAS  Google Scholar 

  11. Müller, S., Hübschmann, T., Kleinsteuber, S. & Vogt, C. High resolution single cell analytics to follow microbial community dynamics in anaerobic ecosystems. Methods 57, 338–349 (2012).

    Article  Google Scholar 

  12. Harnisch, F. et al. Revealing the electrochemically driven selection in natural community derived microbial biofilms using flow-cytometry. Energy Environ. Sci. 4, 1265–1267 (2011).

    Article  CAS  Google Scholar 

  13. Patil, S.A. et al. Electroactive mixed culture derived biofilms in microbial bioelectrochemical systems: the role of pH on biofilm formation, performance and composition. Bioresour. Technol. 102, 9683–9690 (2011).

    Article  CAS  Google Scholar 

  14. Hamady, M. & Knight, R. Microbial community profiling for human microbiome projects: tools, techniques, and challenges. Genome Res. 19, 1141–1152 (2009).

    Article  CAS  Google Scholar 

  15. Martinez-Garcia, M. et al. High-throughput single-cell sequencing identifies photoheterotrophs and chemoautotrophs in freshwater bacterioplankton. ISME J. 6, 113–123 (2012).

    Article  CAS  Google Scholar 

  16. Hartmann, M. & Widmer, F. Reliability for detecting composition and changes of microbial communities by T-RFLP genetic profiling. FEMS Microbiol. Ecol. 63, 249–260 (2007).

    Article  Google Scholar 

  17. Engelbrektson, A. et al. Experimental factors affecting PCR-based estimates of microbial species richness and evenness. ISME J. 4, 642–647 (2010).

    Article  CAS  Google Scholar 

  18. Eagon, R.G. Pseudomonas natrigens, a marine bacterium with a generation time of less than 10 minutes. J. Bacteriol. 83, 736–737 (1962).

    Article  CAS  Google Scholar 

  19. Ulitzur, S. Vibrio parahaemolyticus and Vibrio alginolyticus: short generation-time marine bacteria. Microb. Ecol. 1, 127–135 (1974).

    Article  CAS  Google Scholar 

  20. Liu, B. & Pop, M. MetaPath: identifying differentially abundant metabolic pathways in metagenomic datasets. BMC Proc. 5, S9 (2011).

    Article  Google Scholar 

  21. Kuczynski, J. et al. Experimental and analytical tools for studying the human microbiome. Nat. Rev. Genet. 13, 47–58 (2012).

    Article  CAS  Google Scholar 

  22. Teeling, H. et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611 (2012).

    Article  CAS  Google Scholar 

  23. Wallner, G., Erhart, R. & Amann, R. Flow cytometric analysis of activated sludge with rRNA-targeted probes. Appl. Environ. Microbiol. 61, 1859–1866 (1995).

    Article  CAS  Google Scholar 

  24. Günther, S. et al. Fixation procedures for flow cytometric analysis of environmental bacteria. J. Microbiol. Methods 75, 127–134 (2008).

    Article  Google Scholar 

  25. Bombach, P. et al. Resolution of natural microbial community dynamics by community fingerprinting, flow cytometry, and trend interpretation analysis. In Adv. in Biochem. Eng. Biotechnol. 124, 151–181 Springer, (2011).

    Google Scholar 

  26. R Development Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2011).

  27. Vogt, C., Lösche, A., Kleinsteuber, S. & Müller, S. Population profiles of a stable, commensalistic bacterial culture grown with toluene under sulphate-reducing conditions. Cytometry A 66A, 91–102 (2005).

    Article  CAS  Google Scholar 

  28. Günther, S. et al. Dynamics of polyphosphate-accumulating bacteria in wastewater treatment plant microbial communities detected via DAPI (4′,6′-diamidino-2-phenylindole) and tetracycline labeling. Appl. Environ. Microbiol. 75, 2111–2121 (2009).

    Article  Google Scholar 

  29. Bolboaca, S.D. & Jäntschi, L. Pearson versus Spearman, Kendall's Tau correlation analysis on structure-activity relationships of biologic active compounds. Leonardo J. Sci. 9, 179–200 (2006).

    Google Scholar 

  30. Ramette, A. Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 62, 142–160 (2007).

    Article  CAS  Google Scholar 

  31. Oksanen, J. et al. Vegan: community ecology package. R package version 1.17-10. http://CRAN.R-project.org/package=vegan (2011).

Download references

Acknowledgements

We thank C. Süring, A. Schwarzer, D. Schlosser and G. Weichert for technical assistance. We thank I. Fetzer for helpful discussions regarding R. This work is integrated in the internal research and development program of the UFZ and the CITE (Chemicals In The Environment) program.

Author information

Authors and Affiliations

Authors

Contributions

C.K., S.G. and A.F.D. performed the experiments with contributions of T.H.; C.K. and S.G. developed the bar code and correlation tools; C.K., S.G., T.H. and S.M. designed and performed the research, and interpreted the data; C.K. and S.M. wrote the manuscript.

Corresponding author

Correspondence to Susann Müller.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Supplementary information

Supplementary Figure 1

Effects of duration of DAPI staining. The duration of the staining period has a profound impact on cytometric community analyses. The use of low DAPI concentrations and appropriate incubation times allows for very specific and reproducible DNA staining. The increasing abundance of cells with enhanced DAPI fluorescence after 30 and 150 minutes is shown. The pattern generated after 150 min incubation remained stable for several hours. (PDF 2480 kb)

Supplementary Table 1

Method reproducibility. The error introduced during sample preparation and measurement was tested. Two samples from a natural community were taken, treated identically and stained for DNA analysis with DAPI. For each sample three cytometric measurements were performed. A gate template with four representative gates covering the most dominant sub-communities was applied to all measurements and the number of DAPI stained cells determined (values given in % of all stained cells). It was found that the methodical variation of the cytometric measurement itself introduces a variation between 0.82 and 3.26 %. The sample preparation bias causes only variation up to 2.23 %. From these results it can be stated that a variation found above 3.26 % was due to a real variation in the community structure. Community variations below this value were covered by methodical variations. (PDF 262 kb)

Supplementary Table 2

Original data for batches 1 to 10 of SBR cultivation of a microbial community obtained from Ethiopian tannery wastewater. The table summarizes all cytometric and abiotic data that were used within the ANTICIPATED RESULTS section. (PDF 271 kb)

Supplementary Methods

Experimental setup of the SBR. The general setup of the SBR is presented including details about the determination of abiotic parameters and the tanning agents. (PDF 377 kb)

Supplementary Data 1

barcode.r This R-script can be used to perform step 12 of the procedure. It creates a heat map that displays the normalized cell abundance values for each gate using a color code with 21 gradations (blue, white, red). A box plot visualizes the actual cell number per gate as percentage of all measured cells. (TXT 2 kb)

Supplementary Data 2

nmds.r This R-script can be used to perform step 13 of the procedure. It visualizes the distance between samples based on their cytometrically measured cell abundance information using nonmetric multidimensional scaling (NMDS). (TXT 1 kb)

Supplementary Data 3

correlation.r This R-script can be used to perform step 15 of the procedure. It determines relationships between biotic and abiotic parameters based on Spearman's rank-order correlation coefficient and visualizes them as heat map with positive correlations in red, negative correlations in blue and neutral correlations in white. (TXT 2 kb)

Rights and permissions

Reprints and permissions

About this article

Cite this article

Koch, C., Günther, S., Desta, A. et al. Cytometric fingerprinting for analyzing microbial intracommunity structure variation and identifying subcommunity function. Nat Protoc 8, 190–202 (2013). https://doi.org/10.1038/nprot.2012.149

Download citation

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nprot.2012.149

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing Microbiology

Sign up for the Nature Briefing: Microbiology newsletter — what matters in microbiology research, free to your inbox weekly.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing: Microbiology