Natural product biosynthetic potential reflects macroevolutionary diversification within a widely distributed bacterial taxon

Flavobacteriaceae spp. are key players in global biogeochemical cycling and are known for their versatile carbohydrate and peptide catabolism. However, it is currently unknown whether secondary metabolism traits underlie their broad range of occurrence across the earth's biomes. We examined 2,680 genomes to unveil an unprecedented phylogenetic signal dictating natural product biosynthesis diversification within the Flavobacteriaceae family. The distribution of secondary metabolite biosynthetic gene clusters (BGCs) across genomes usually follows macroevolutionary, genus-specific patterns. Noticeably, 88.6% of the observed BGCs were inferred to lead to the biosynthesis of likely novel natural products. We found an unanticipated, large diversity of taxon-specific BGCs encoding carotenoid and flexirubin pigments, the vast majority of which awaiting formal description. In particular, Aquimarina and Kordia spp. possessed large genomes, versatile catabolic traits, and a repertoire of BGCs possibly encoding drug-inspiring polyketides, non-ribosomal peptides, or post-translationally modified peptides. Using a machine learning approach (feature selection), we reveal that marine and non-marine Flavobacteriaceae genomes are differentially enriched in CAZymes and peptidases with distinct functionalities and molecular targets.