International audienceSOUL homologs constitute a heme-binding protein superfamily putatively involved in heme and tetrapyrrole metabolisms associated with a number of physiological processes. Despite their omnipresence across the tree of life and the biochemical characterization of many SOUL members, their functional role and the evolutionary events leading to such remarkable protein repertoire still remain cryptic. To explore SOUL evolution, we apply a computational phylogenetic approach, including a relevant number of SOUL homologs, to identify paralog forms and reconstruct their genealogy across the tree of life and within species. In animal lineages, multiple gene duplication or loss events and paralog functional specializations underlie SOUL evolution from the dawn of ancestral echinoderm and mollusc SOUL forms. In photosynthetic organisms, SOUL evolution is linked to the endosymbiosis events leading to plastid acquisition in eukaryotes. Derivative features, such as the F2L peptide and BH3 domain, evolved in vertebrates and provided innovative functionality to support immune response and apoptosis. The evolution of elements such as the N-terminal protein domain DUF2358, the His42 residue, or the tetrapyrrole heme-binding site is modern, and their functional implications still unresolved. This study represents the first in-depth analysis of SOUL protein evolution and provides novel insights in the understanding of their obscure physiological role