Evolution of prokaryotic respiratory molybdoenzymes and the frequency of their genomic co-occurrence

Molybdoenzymes are an ancient protein family found in phylogenetically and ecologically diverse prokaryotes. Under anaerobic conditions, respiratory molybdoenzymes catalyze redox reactions that transfer electrons to a variety of substrates that act as terminal electron acceptors for energy generation. Here, we used probe sequences to conduct an extensive genomic survey and phylogenetic inference for NarG, DmsA, TorA and nine other respiratory molybdoenzyme subfamilies. Our analysis demonstrates their abundance in 60% of prokaryotic phyla. In contrast to many other autonomic genetic units in prokaryotes, the major route of evolution of their predominant subfamilies is vertical gene transfer, gene duplication and divergence. Our results show the robustness of genomic co-occurrence of respiratory molybdoenzymes genes, found in the majority of studied species, for most of the enzyme subfamilies. Genomes which encode for multiple respiratory molybdoenzymes are also enriched in genes regulating replication, recombination and mobility of genetic elements. Respiratory molybdoenzymes were found in prokaryotes associated with diverse environments occupying terrestrial, aquatic, food and host-related habitats, emphasizing their essential role in adaptation of prokaryotes to changing environments. Interestingly, host-associated prokaryotes such as human pathogens more frequently carry multiple respiratory molybdoenzyme genes compared with non-host-associated prokaryotes, highlighting the importance of metabolic flexibility in host-microbiome environments.