Flavin-dependent monooxygenases play a variety of key physiological roles and are also very powerful biotechnological tools. These enzymes have been classified into eight different classes (A–H) based on their sequences and biochemical features. By combining structural and sequence analysis, and phylogenetic inference, we have explored the evolutionary history of classes A, B, E, F, and G and demonstrate that their multidomain architectures reflect their phylogenetic relationships, suggesting that the main evolutionary steps in their divergence are likely to have arisen from the recruitment of different domains. Additionally, the functional divergence within in each class appears to have been the result of other mechanisms such as a complex set of single-point mutations. Our results reinforce the idea that a main constraint on the evolution of cofactor-dependent enzymes is the functional binding of the cofactor. Additionally, a remarkable feature of this family is that the sequence of the key flavin adenine dinucleotide-binding domain is split into at least two parts in all classes studied here. We propose a complex set of evolutionary events that gave rise to the origin of the different classes within this family.

Ribosome inactivating proteins are enzymes that depurinate a specific adenine residue in the alpha-sarcin-ricin loop of the large ribosomal RNA, being ricin and Shiga toxins the most renowned examples. They are widely distributed in plants and their presence has also been confirmed in a few bacterial species. According to this taxonomic distribution, the current model about the origin and evolution of RIP genes postulates that an ancestral RIP domain was originated in flowering plants, and later acquired by some bacteria via horizontal gene transfer. Here, we unequivocally detected the presence of RIP genes in fungi and metazoa. These findings, along with sequence and phylogenetic analyses, led us to propose an alternative, more parsimonious, hypothesis about the origin and evolutionary history of the RIP domain, where several paralogous RIP genes were already present before the three domains of life evolved. This model is in agreement with the current idea of the Last Universal Common Ancestor (LUCA) as a complex, genetically redundant organism. Differential loss of paralogous genes in descendants of LUCA, rather than multiple horizontal gene transfer events, could account for the complex pattern of RIP genes across extant species, as it has been observed for other genes.

The Baeyer—Villiger Monooxygenases (BVMOs) are enzymes belonging to the “Class B” of flavin monooxygenases and are capable of performing exquisite selective oxidations. These enzymes have been studied from a biotechnological perspective, but their physiological substrates and functional roles are widely unknown. Here, we investigated the origin, taxonomic distribution and evolutionary history of the BVMO genes. By using in silico approaches, 98 BVMO encoding genes were detected in the three domains of life: Archaea, Bacteria and Eukarya. We found evidence for the presence of these genes in Metazoa (Hydra vulgaris, Oikopleura dioica and Adineta vaga) and Haptophyta (Emiliania huxleyi) for the first time. Furthermore, a search for other “Class B” monooxygenases (flavoprotein monooxygenases –FMOs – and N-hydroxylating monooxygenases – NMOs) was conducted. These sequences were also found in the three domains of life. Phylogenetic analyses of all “Class B” monooxygenases revealed that NMOs and BVMOs are monophyletic, whereas FMOs form a paraphyletic group. Based on these results, we propose that BVMO genes were already present in the last universal common ancestor (LUCA) and their current taxonomic distribution is the result of differential duplication and loss of paralogous genes.