Global analyses of biosynthetic gene clusters in phytobiomes reveal strong phylogenetic conservation of terpenes and aryl polyenes

There are gaps in our understandings on how did the evolutionary relationships among members of the phytobiomes shape their ability to produce tremendously complex specialized metabolites under the influence of plant host. To determine these relationships, we investigated the phylogenetic conservation of biosynthetic gene clusters (BGCs) on a global collection of 4,519 high-quality and nonredundant (out of 12,181) bacterial isolates and metagenome-assembled genomes from 47 different plant hosts and soil, by adopting three independent phylogenomic approaches (D-test, Pagel’s λ, and consenTRAIT). We report that the BGCs are phylogenetically conserved to varying strengths and depths in their different classes. We show that the ability to produce specialized metabolites qualifies as a complex trait, and the depth of conservation is equivalent to ecologically relevant complex microbial traits. Interestingly, terpene and aryl polyene BGCs had the strongest phylogenetic conservation in the phytobiomes, but not in the soil microbiomes. Furthermore, we showed that terpenes are largely uncharacterized in phytobiomes and pinpointed specific clades that harbor potentially novel terpenes. Taken together, this study sheds light on the evolution of specialized metabolites’ biosynthesis potential in phytobiomes under the influence of plant hosts and presents strategies to rationally guide the discovery of potentially novel classes of metabolites. IMPORTANCE: This study expands our understandings of the biosynthetic potential of phytobiomes by using such worldwide and extensive collection of microbiomes from plants and soil. Apart from providing such vital resource for the plant microbiome researchers, this study provides fundamental insights into the evolution of biosynthetic gene clusters (BGCs) in phytobiomes under the influence of plant host. Specifically, we report that the strength of phylogenetic conservation in microbiomes varies for different classes of BGCs and is influenced as a result of plant host association. Furthermore, our results indicate that biosynthetic potential of specialized metabolites is deeply conserved equivalent to other complex and ecologically relevant microbial traits. Finally, for the most conserved class of specialized metabolites (terpenes), we identified clades harboring potentially novel class of molecules. Future studies could focus on plant–microbe coevolution and interactions through specialized metabolites building upon these findings.