For each MCMC, we performed a run of 50 and 100 million generations (for mtDNA and nuclear genes, respectively), sampling every 1000th generation. The first 10% of the trees were discarded as burn-in for each run. We assessed whether parameter values for individual runs had reached equilibrium and convergence by visually assessing their trace plots in Tracer 1.5. Individual chains were combined using the LogCombiner v1.7.5. Finally, we analyzed combined runs using Tree Annotator v1.7.5.