[en] Predicting the behaviour of complex microbial communities is challenging. However, this is essential for complex biotechnological processes such as those in biological wastewater treatment plants (BWWTPs), which require sustainable operation. Here we summarize 14 months of longitudinal meta-omics data from a BWWTP anaerobic tank into 17 temporal signals, explaining 91.1% of the temporal variance, and link those signals to ecological events within the community. We forecast the signals over the subsequent five years and use 21 extra samples collected at defined time intervals for testing and validation. Our forecasts are correct for six signals and hint on phenomena such as predation cycles. Using all the 17 forecasts and the environmental variables, we predict gene abundance and expression, with a coefficient of determination ≥0.87 for the subsequent three years. Our study demonstrates the ability to forecast the dynamics of open microbial ecosystems using interactions between community cycles and environmental parameters.
Research center :
Luxembourg Centre for Systems Biomedicine (LCSB): Eco-Systems Biology (Wilmes Group) Luxembourg Centre for Systems Biomedicine (LCSB): Bioinformatics Core (R. Schneider Group)
Disciplines :
Microbiology Environmental sciences & ecology
Author, co-author :
DELOGU, Francesco ; University of Luxembourg > Luxembourg Centre for Systems Biomedicine > Systems Ecology > Team Paul WILMES
Kunath, Benoit J; Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
Queirós, Pedro M; Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
HALDER, Rashi ; University of Luxembourg > Luxembourg Centre for Systems Biomedicine (LCSB) > Scientific Central Services > Sequencing Platform
LEBRUN, Laura ; University of Luxembourg > Luxembourg Centre for Systems Biomedicine > Systems Ecology > Team Paul WILMES
Pope, Phillip B ; Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway ; Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
Widder, Stefanie; Department of Medicine 1, Research Division Infection Biology, Medical University of Vienna, Vienna, Austria
MULLER, Emilie ; University of Luxembourg ; Génétique Moléculaire, Génomique, Microbiologie, UMR 7156 CNRS, Université de Strasbourg, Strasbourg, France
WILMES, Paul ; University of Luxembourg > Luxembourg Centre for Systems Biomedicine (LCSB) > Systems Ecology
External co-authors :
yes
Language :
English
Title :
Forecasting the dynamics of a complex microbial community using integrated meta-omics.
FNR - Fonds National de la Recherche FWF - Austrian Science Fund RCN - Research Council of Norway Union Européenne
Funding text :
We thank the Luxembourg National Research Fund (FNR) for supporting this work through various funding instruments. Specifically, a PRIDE doctoral training unit grant (PRIDE/15/10907093), CORE grants (CORE/17/SM/11689322), a European Union ERASysAPP grant (INTER/SYSAPP/14/05) and an ATTRACT grant (A09/03) all awarded to P.W., as well as a CORE Junior (C15/SR/10404839) grant to E.E.L.M. The project received financial support from the Integrated Biobank of Luxembourg with funds from the Luxembourg Ministry of Higher Education and Research. This work was also supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant Agreement No. 863664). The work of P.M. was funded by the ‘Plan Technologies de la Santé du Gouvernement du Grand-Duché de Luxembourg’ through the Luxembourg Centre for Systems Biomedicine (LCSB), University of Luxembourg. S.W. was supported by the Austrian Science Fund (FWF) Elise Richter V585-B31. P.B.P. is grateful for the support from the Research Council of Norway (FRIPRO programme: 250479) and the Novo Nordisk Foundation (Project No. 0054575). The authors acknowledge the ULHPC for providing and maintaining the computing resources. We duly thank Mr G. Bissen and Mr G. Di Pentima of the Syndicat Intercommunal a Vocation Ecologique (SIVEC) for access to the Schifflange wastewater treatment plant.
Martiny, J. B. H. et al. Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112 (2006). DOI: 10.1038/nrmicro1341
Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on Earth. Proc. Natl Acad. Sci. USA 115, 6506–6511 (2018). DOI: 10.1073/pnas.1711842115
Falkowski, P. G., Fenchel, T. & Delong, E. F. The microbial engines that drive Earth’s biogeochemical cycles. Science 320, 1034–1039 (2008). DOI: 10.1126/science.1153213
Larsen, P. E., Field, D. & Gilbert, J. A. Predicting bacterial community assemblages using an artificial neural network approach. Nat. Methods 9, 621–625 (2012). DOI: 10.1038/nmeth.1975
García-Jiménez, B., Muñoz, J., Cabello, S., Medina, J. & Wilkinson, M. D. Predicting microbiomes through a deep latent space. Bioinformatics 37, 1444–1451 (2021). DOI: 10.1093/bioinformatics/btaa971
O’Sullivan, J. D., Terry, J. C. D. & Rossberg, A. G. Intrinsic ecological dynamics drive biodiversity turnover in model metacommunities. Nat. Commun. 12, 3627 (2021). DOI: 10.1038/s41467-021-23769-7
Hutchins, D. A. & Fu, F. Microorganisms and ocean global change. Nat. Microbiol. 2, 17058 (2017). DOI: 10.1038/nmicrobiol.2017.58
Reed, D. C., Algar, C. K., Huber, J. A. & Dick, G. J. Gene-centric approach to integrating environmental genomics and biogeochemical models. Proc. Natl Acad. Sci. USA 111, 1879–1884 (2014). DOI: 10.1073/pnas.1313713111
Louca, S. et al. Integrating biogeochemistry with multiomic sequence information in a model oxygen minimum zone. Proc. Natl. Acad. Sci. USA 113, E5925–E5933 (2016). DOI: 10.1073/pnas.1602897113
Daims, H., Taylor, M. W. & Wagner, M. Wastewater treatment: a model system for microbial ecology. Trends Biotechnol. 24, 483–489 (2006). DOI: 10.1016/j.tibtech.2006.09.002
Dottorini, G. et al. Mass-immigration determines the assembly of activated sludge microbial communities. Proc. Natl. Acad. Sci. USA 118, e2021589118 (2021). DOI: 10.1073/pnas.2021589118
Chen, J. et al. Economic assessment of biodiesel production from wastewater sludge. Bioresour. Technol. 253, 41–48 (2018). DOI: 10.1016/j.biortech.2018.01.016
Kim, Y. K. et al. The capacity of wastewater treatment plants drives bacterial community structure and its assembly. Sci. Rep. 9, 14809 (2019). DOI: 10.1038/s41598-019-50952-0
Dueholm, M. K. D. et al. MiDAS 4: a global catalogue of full-length 16S rRNA gene sequences and taxonomy for studies of bacterial communities in wastewater treatment plants. Nat. Commun. 13, 1908 (2022). DOI: 10.1038/s41467-022-29438-7
Wade, M. J. et al. Understanding and managing uncertainty and variability for wastewater monitoring beyond the pandemic: lessons learned from the United Kingdom national COVID-19 surveillance programmes. J. Hazard. Mater. 424, 127456 (2022). DOI: 10.1016/j.jhazmat.2021.127456
Bedoya, K., Hoyos, O., Zurek, E., Cabarcas, F. & Alzate, J. F. Annual microbial community dynamics in a full-scale anaerobic sludge digester from a wastewater treatment plant in Colombia. Sci. Total Environ. 726, 138479 (2020). DOI: 10.1016/j.scitotenv.2020.138479
Frigon, D. et al. Long-term analysis of a full-scale activated sludge wastewater treatment system exhibiting seasonal biological foaming. Water Res. 40, 990–1008 (2006). DOI: 10.1016/j.watres.2005.12.015
Wang, Y. et al. Successional dynamics and alternative stable states in a saline activated sludge microbial community over 9 years. Microbiome 9, 199 (2021). DOI: 10.1186/s40168-021-01151-5
Pérez, M. V., Guerrero, L. D., Orellana, E., Figuerola, E. L. & Erijman, L. Time series genome-centric analysis unveils bacterial response to operational disturbance in activated sludge. mSystems 4, e00169-19 (2019). DOI: 10.1128/mSystems.00169-19
Sheik, A. R., Muller, E. E. L. & Wilmes, P. A hundred years of activated sludge: time for a rethink. Front. Microbiol. 5, 47 (2014). DOI: 10.3389/fmicb.2014.00047
Winkler, M. K. & Straka, L. New directions in biological nitrogen removal and recovery from wastewater. Curr. Opin. Biotechnol. 57, 50–55 (2019). DOI: 10.1016/j.copbio.2018.12.007
Alter, O., Brown, P. O. & Botstein, D. Singular value decomposition for genome-wide expression data processing and modeling. Proc. Natl Acad. Sci. USA 97, 10101–10106 (2000). DOI: 10.1073/pnas.97.18.10101
Strang, G. Linear Algebra and Its Applications. SIAM Review (Cengage Learning, 2006).
Yanai, H., Takeuchi, K. & Takane, Y. Projection Matrices, Generalized Inverse Matrices, and Singular Value Decomposition (Springer, 2011).
Hyndman, R. J. & Athanasopoulos, G. Forecasting: Principles and Practice 3rd edn (OTexts, 2021).
Roume, H. et al. A biomolecular isolation framework for eco-systems biology. ISME J. 7, 110–121 (2013). DOI: 10.1038/ismej.2012.72
Herold, M. et al. Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance. Nat. Commun. 11, 5281 (2020). DOI: 10.1038/s41467-020-19006-2
Taylor, S. J. & Letham, B. Forecasting at scale. Am. Stat. 72, 37–45 (2018). DOI: 10.1080/00031305.2017.1380080
O’Hara-Wild, M., Handyman, R. Wang, E. fable: forecasting models for tidy time series. R package https://fable.tidyverts.org/ (2023).
Roume, H. et al. Comparative integrated omics: identification of key functionalities in microbial community-wide metabolic networks. npj Biofilms Microbiomes 1, 15007 (2015). DOI: 10.1038/npjbiofilms.2015.7
Martínez Arbas, S. et al. Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics. Nat. Microbiol. 6, 123–135 (2020). DOI: 10.1038/s41564-020-00794-8
Muller, E. E. L. L. et al. Community-integrated omics links dominance of a microbial generalist to fine-tuned resource usage. Nat. Commun. 5, 5603 (2014). DOI: 10.1038/ncomms6603
Schramm, A., de Beer, D., Wagner, M. & Amann, R. Identification and activities in situ of Nitrosospira and Nitrospira spp. as dominant populations in a nitrifying fluidized bed reactor. Appl. Environ. Microbiol. 64, 3480–3485 (1998). DOI: 10.1128/AEM.64.9.3480-3485.1998
Narayanasamy, S. et al. IMP: a pipeline for reproducible reference-independent integrated metagenomic and metatranscriptomic analyses. Genome Biol. 17, 260 (2016). DOI: 10.1186/s13059-016-1116-8
Nikolenko, S. I., Korobeynikov, A. I. & Alekseyev, M. A. BayesHammer: Bayesian clustering for error correction in single-cell sequencing. BMC Genomics 14, S7 (2013). DOI: 10.1186/1471-2164-14-S1-S7
Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. MetaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834 (2017). DOI: 10.1101/gr.213959.116
Bushmanova, E., Antipov, D., Lapidus, A. & Prjibelski, A. D. rnaSPAdes: a de novo transcriptome assembler and its application to RNA-seq data. Gigascience 8, giz100 (2019). DOI: 10.1093/gigascience/giz100
West, P. T., Probst, A. J., Grigoriev, I. V., Thomas, B. C. & Banfield, J. F. Genome-reconstruction for eukaryotes from complex natural microbial communities. Genome Res. 28, 569–580 (2018). DOI: 10.1101/gr.228429.117
Krawczyk, P. S., Lipinski, L. & Dziembowski, A. PlasFlow: predicting plasmid sequences in metagenomic data using genome signatures. Nucleic Acids Res. 46, e35–e35 (2018). DOI: 10.1093/nar/gkx1321
Zhou, F. & Xu, Y. cBar: a computer program to distinguish plasmid-derived from chromosome-derived sequence fragments in metagenomics data. Bioinformatics 26, 2051–2052 (2010). DOI: 10.1093/bioinformatics/btq299
Roux, S., Enault, F., Hurwitz, B. L. & Sullivan, M. B. VirSorter: mining viral signal from microbial genomic data. PeerJ 3, e985 (2015). DOI: 10.7717/peerj.985
Ren, J. et al. Identifying viruses from metagenomic data using deep learning. Quant. Biol. 8, 64–77 (2020). DOI: 10.1007/s40484-019-0187-4
Wu, Y.-W., Tang, Y.-H., Tringe, S. G., Simmons, B. A. & Singer, S. W. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome 2, 26 (2014). DOI: 10.1186/2049-2618-2-26
Kang, D. D., Froula, J., Egan, R. & Wang, Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ 3, e1165 (2015). DOI: 10.7717/peerj.1165
Hickl, O., Queirós, P., Wilmes, P., May, P. & Heintz-Buschart, A. binny: an automated binning algorithm to recover high-quality genomes from complex metagenomic datasets. Brief. Bioinform. 23, bbac431 (2022). DOI: 10.1093/bib/bbac431
Sieber, C. M. K. et al. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat. Microbiol. 3, 836–843 (2018). DOI: 10.1038/s41564-018-0171-1
Olm, M. R., Brown, C. T., Brooks, B. & Banfield, J. F. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 11, 2864–2868 (2017). DOI: 10.1038/ismej.2017.126
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015). DOI: 10.1101/gr.186072.114
Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152 (2012). DOI: 10.1093/bioinformatics/bts565
von Meijenfeldt, F. A. B., Arkhipova, K., Cambuy, D. D., Coutinho, F. H. & Dutilh, B. E. Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT. Genome Biol. 20, 217 (2019). DOI: 10.1186/s13059-019-1817-x
Hoff, K. J. & Stanke, M. Predicting genes in single genomes with AUGUSTUS. Curr. Protoc. Bioinform. 65, e57 (2019). DOI: 10.1002/cpbi.57
Queirós, P., Delogu, F., Hickl, O., May, P. & Wilmes, P. Mantis: flexible and consensus-driven genome annotation. Gigascience 10, giab042 (2021). DOI: 10.1093/gigascience/giab042
Aramaki, T. et al. KofamKOALA: KEGG ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics 36, 2251–2252 (2020). DOI: 10.1093/bioinformatics/btz859
Selengut, J. D. et al. TIGRFAMs and genome properties: tools for the assignment of molecular function and biological process in prokaryotic genomes. Nucleic Acids Res. 35, 260–264 (2007). DOI: 10.1093/nar/gkl1043
Huerta-Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314 (2019). DOI: 10.1093/nar/gky1085
Mistry, J. et al. Pfam: the protein families database in 2021. Nucleic Acids Res. 49, D412–D419 (2021). DOI: 10.1093/nar/gkaa913
Li, W. et al. RefSeq: expanding the prokaryotic genome annotation pipeline reach with protein family model curation. Nucleic Acids Res. 49, D1020–D1028 (2021). DOI: 10.1093/nar/gkaa1105
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009). DOI: 10.1093/bioinformatics/btp324
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009). DOI: 10.1093/bioinformatics/btp352
Turro, E. et al. Haplotype and isoform specific expression estimation using multi-mapping RNA-seq reads. Genome Biol. 12, R13 (2011). DOI: 10.1186/gb-2011-12-2-r13
Turro, E., Astle, W. J. & Tavaré, S. Flexible analysis of RNA-seq data using mixed effects models. Bioinformatics 30, 180–188 (2014). DOI: 10.1093/bioinformatics/btt624
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014). DOI: 10.1093/bioinformatics/btt656
Barsnes, H. & Vaudel, M. SearchGUI: a highly adaptable common interface for proteomics search and de novo engines. J. Proteome Res. 17, 2552–2555 (2018). DOI: 10.1021/acs.jproteome.8b00175
Langella, O. et al. X!TandemPipeline: a tool to manage sequence redundancy for protein inference and phosphosite identification. J. Proteome Res. 16, 494–503 (2017). DOI: 10.1021/acs.jproteome.6b00632
Kim, S. & Pevzner, P. A. MS-GF+ makes progress towards a universal database search tool for proteomics. Nat. Commun. 5, 5277 (2014). DOI: 10.1038/ncomms6277
Eng, J. K., Jahan, T. A. & Hoopmann, M. R. Comet: an open-source MS/MS sequence database search tool. Proteomics 13, 22–24 (2013). DOI: 10.1002/pmic.201200439
Vaudel, M. et al. PeptideShaker enables reanalysis of MS-derived proteomics data sets. Nat. Biotechnol. 33, 22–24 (2015). DOI: 10.1038/nbt.3109
Schallert, K. et al. Pout2Prot: an efficient tool to create protein (sub)groups from percolator output files. J. Proteome Res. 21, 1175–1180 (2022). DOI: 10.1021/acs.jproteome.1c00685
Schiebenhoefer, H. et al. A complete and flexible workflow for metaproteomics data analysis based on MetaProteomeAnalyzer and Prophane. Nat. Protoc. 15, 3212–3239 (2020). DOI: 10.1038/s41596-020-0368-7
Van Der Boogaart, K. G. & Tolosana-Delgado, R. Compositional data analysis with ‘R’ and the package ‘compositions’. Geol. Soc. Spec. Publ. 264, 119–127 (2006). DOI: 10.1144/GSL.SP.2006.264.01.09
Jones, M. C. & Aitchison, J. The statistical analysis of compositional data. J. R. Stat. Soc. A 150, 396 (1987). DOI: 10.2307/2982045
Erb, I., Gloor, G. B. & Quinn, T. P. Editorial: Compositional data analysis and related methods applied to genomics—a first special issue from NAR Genomics and Bioinformatics. NAR Genom. Bioinform. 2, lqaa103 (2020). DOI: 10.1093/nargab/lqaa103
Royston, J. P. An extension of Shapiro and Wilk’s W test for normality to large samples. Appl. Stat. 31, 115–124 (1982). DOI: 10.2307/2347973
Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9, 559 (2008). DOI: 10.1186/1471-2105-9-559
Hothorn, T. & Zeileis, A. Diagnostic checking in regression relationships. R News 2/3, 7–10 (2011).
Sugihara, G., Ye, H., Clark, A. Deyle, E. rEDM: An R Package for Empirical Dynamic Modeling and Convergent Cross Mapping (SERDP, 2019).
Sugihara, G. Nonlinear forecasting for the classification of natural time series. Phil. Trans. R. Soc. London. A 348, 477–495 (1994). DOI: 10.1098/rsta.1994.0106
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003). DOI: 10.1101/gr.1239303
Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 36, 1925–1927 (2020). DOI: 10.1093/bioinformatics/btz848
