Registering the evolutionary history in individual-based models of speciation

Understanding the emergence of biodiversity patterns in nature is a central problem in biology. Theoretical models of speciation have addressed this question in the macroecological scale, but little has been done to connect microevolutionary processes with macroevolutionary patterns. Knowledge of the evolutionary history allows the study of patterns underlying the processes being modeled, revealing their signatures and the role of speciation and extinction in shaping macroevolutionary patterns. In this paper we introduce two algorithms to record the evolutionary history of populations and species in individual-based models of speciation, from which genealogies and phylogenies can be constructed. The first algorithm relies on saving ancestor-descendant relationships, generating a matrix that contains the times to the most recent common ancestor between all pairs of individuals at every generation (the Most Recent Common Ancestor Time matrix, MRCAT). The second algorithm directly records all speciation and extinction events throughout the evolutionary process, generating a matrix with the true phylogeny of species (the Sequential Speciation and Extinction Events, SSEE). We illustrate the use of these algorithms in a spatially explicit individual-based model of speciation. We compare the trees generated via MRCAT and SSEE algorithms with trees inferred by methods that use only genetic distance between individuals of extant species, commonly used in empirical studies and applied here to simulated genetic data. Comparisons between trees are performed with metrics describing the overall topology, branch length distribution and imbalance degree. We observe that both MRCAT and distance-based trees differ from the true phylogeny, with the first being closer to the true tree than the second.