Baltic Sea microbial cohorts exhibit catabolic specialization and anabolic interdependencies across environmental gradients

Abstract

Microbial communities are structured by environmental gradients and metabolic interactions, yet the genomic characteristics and metabolic functions of co-occurring populations remain underexplored. Here, we investigated co-occurring microbial cohorts across the Baltic Sea, a system characterized by strong salinity, temperature, and oxygen gradients. For this, we used a genomic catalog consisting of 701 species-representative genomes to recruit reads from 112 metagenomes and infer cohort structure, environmental distributions, and metabolic potential. We identified nine microbial cohorts that showed strong associations with environmental gradients, indicating deterministic assembly. Cohorts differed markedly in genomic traits, with the most abundant and prevalent taxa associated with smaller, streamlined genomes, while a low-oxygen cohort with larger genomes contributed disproportionately to nitrogen and sulfur transformations. Across cohorts, biosynthetic potential was unevenly distributed. Amino acid biosynthesis pathways were frequently complete, whereas B-vitamin pathways were typically incomplete and rarely encoded in full by individual genomes. Metabolites with low pathway completeness showed consistent taxonomic partitioning, with biosynthetic capabilities distributed across taxa rather than collectively encoded within cohorts. Together, these results show that Baltic Sea microbial cohorts are ecologically structured assemblages whose genomic repertoires reflect catabolic specialization and anabolic interdependencies. Our findings highlight microbial cohorts as a useful framework for linking environmental gradients, genome traits, and the organization of metabolic functions in natural microbial communities.