Integrating metagenome-scale metabolic modelling and metabolomics to identify biochemical interactions in Microcystis phycospheres

Favoured by global changes, freshwater cyanobacterial harmful blooms generate major ecological, economical and public health challenges. Microcystis, one of the most widespread cyanobacterial genera, grows within a phycosphere where specialised interactions with its microbiome occur, and are suspected to influence bloom appearance and its potential toxicity. Using a combination of metagenomic, metabolomic and metabolic modelling, we characterised the phycospheres of twelve Microcystis strains isolated from a French pond. The distribution of metabolic reactions within Microcystis was consistent with their genospecies, whereas the metabolic landscape at the community level diverged from cyanobacterial phylogeny indicating functional decoupling between cyanobacteria and their associated microbiomes. Phycosphere-associated bacteria substantially expand the metabolic repertoire of the system, while maintaining functional redundancy within and across communities. On the other hand, metabolomic profiles were largely driven by cyanobacterial metabolic outputs. Metabolic modelling, together with the identification of toxic specialised metabolites produced by specific biosynthetic gene clusters, further highlighted differences in metabolic potential among phycospheres. Together, these findings deepen the understanding of Microcystis phycosphere functioning, demonstrate the value of multi-omics systems biology approaches, and underscore the ecological relevance of interspecies and inter-phycosphere metabolic interactions as a structuring process in bloom-associated microbiomes.