Abstract:Motivation: As the density of sampled population increases, especially as studies incorporate aspects of the spatial landscape to study evolutionary processes, efficient simulation of genetic data under the coalescent becomes a primary challenge. Beyond the computational demands, coalescence-based simulation strategies have to be reconsidered because traditional assumptions about the dynamics of coalescing lineages within local populations may be violated (e.g. more than two daughter lineages may coalesce to a parent at low population densities). Specifically, to efficiently assign n lineages to m parents, the order relation between n and m strongly affects the relevant algorithm for the coalescent simulator (e.g. only when n<2m, it is reasonable to assume that two lineages, at most, can be assigned to the same parent). Controlling the details of the simulation model as a function of n and m is then crucial to represent accurately and efficiently the assignment process, but current implementations make it difficult to switch between different types of lineage mergers at run-time or even compile-time. Results: With the described occupancy spectrum and algorithm that generates the support of the joint probability distribution of the occupancy spectrum; computation is much faster than realizing the whole assignment process under the coalescent. Using general definitions of lineage merges, which also makes the codebase reusable, we implement several variants of coalescent mergers, including an approximation where low probability spectrums are discarded. Comparison of runtimes and performance of the different C++ highly reusable coalescence mergers (binary, multiple, hybrids) are given, and we illustrate their potential utility with example applications. Availability and implementation: All components are integrated into Quetzal, an open-source C++ library for coalescence available at https://becheler.github.io/pages/quetzal.html. Supplementary information: Supplementary data are available at Bioinformatics online.

Hybrid-Lambda: simulation of multiple merger and Kingman gene genealogies in species networks and species trees

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution

Occupancy spectrum distribution: application for coalescence simulation with generic mergers

AleRax: a tool for gene and species tree co-estimation and reconciliation under a probabilistic model of gene duplication, transfer, and loss

Ultrafast learning of 4-node hybridization cycles in phylogenetic networks using algebraic invariants

Statistical Inference of Allopolyploid Species Networks in the Presence of Incomplete Lineage Sorting

Detectability of Varied Hybridization Scenarios using Genome-Scale Hybrid Detection Methods

ParGenes: a tool for massively parallel model selection and phylogenetic tree inference on thousands of genes

Anomaly zones for uniformly sampled gene trees under the gene duplication and loss model

Normalizing Kernels in the Billera-Holmes-Vogtmann Treespace

reconcILS: A gene tree-species tree reconciliation algorithm that allows for incomplete lineage sorting

ASTRAL: genome-scale coalescent-based species tree estimation

Generating Dark Matter Halo Merger Trees

A global test of hybrid ancestry from genome-scale data

Hierarchical Heuristic Species Delimitation under the Multispecies Coalescent Model with Migration

Multiple merger coalescent inference of effective population size

Statistical properties of the site-frequency spectrum associated with Lambda-coalescents

SimHyb 2: a software tool to explore and illustrate evolutionary forces in population genetics teaching and research. Application to Conservation Genetics

The Effect of Copy Number Hemiplasy on Gene Family Evolution

Reconciling Event-Labeled Gene Trees with MUL-trees and Species Networks

A Perl Package and an Alignment Tool for Phylogenetic Networks