Prior and posterior distributions of parameters in the BPP Bayesian analysis of the nine nuclear loci.

aPriors set with Ne = 225,000 covering Watterson's [91] θ estimate (θW = 0.002, Ne = 138,889) and from branching event t = 3.31 Mya between European closely related Cobitis species [104]. Relatively fast autosomal mutation rate (μ) of 3.6×10−9 estimated in vertebrates [105] was used to transform prior expectations of θ and τ from absolute estimates of Ne and t. Both τ and θ are measured as the expected number of mutations per site.bPrior for the node age was generated from the Dirichlet distribution ([63]: equation 2).C. tae = C. taenia, C. tan = C. tanaitica, C. elo = C. elongatoides.