Abstract High-throughput sequencing projects generate genome-scale sequence data for species-level phylogenies1–3. However, state-of-the-art Bayesian methods for inferring timetrees are computationally limited to small datasets and