With the advent of novel computational techniques like Markov chain Monte Carlo (MCMC) sampling, it has become feasible to integrate many of the models involved and simultaneously estimate parameters of interest. Demographic inference is a wellknown example of genealogy-based population genetics that benefited from these advances [24, 25] . Bayesian MCMC methods also enable ancestral state reconstruction while simultaneously accounting for both phylogenetic and mapping uncertainty. Although this adds much needed credibility to ancestral reconstruction [13] , phylogeographic analysis would benefit even more from fully integrating spatial, temporal and demographic inference.