Mirage: A phylogenetic mixture model to reconstruct gene-content evolutionary history using a realistic evolutionary rate model

Reconstruction of gene-content evolutionary history is an essential approach for understanding how complex biological systems have been organized. However, the existing gene-content evolutionary models cannot formulate complex and heterogeneous gene gain/loss processes, which reflect diverse evolutionary events and greatly depend on gene families. In this study, we developed Mirage (MIxture model with a Realistic evolutionary rate model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates, but reasonably limits the number of parameters to be estimated by the expectation-maximization algorithm. Simulation analysis showed that Mirage can accurately estimate complex and heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that our evolutionary model better fits genome data from various taxonomic groups than other models. Using Mirage, we revealed that gene families of metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated. The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage.