Molecular dating chronogram using cytochrome b with a HKY + I + G model (ML and BI topology presents similar topology).

Statistical support higher than 0.90 posterior probability (PP) and 70 bootstraps is identified next to the nodes divided in two sides, the left side represents PP values and the right side the bootstrap values; black represent 1 PP and 100 bootstrap, very dark gray represent 0,99 PP and 91–99 bootstrap, dark gray represent values between 0.96–0.98 PP and 80–90 bootstrap, and light gray represents values between 0.95–0.96 PP and 70–80 bootstrap. The letters B, C, D, E, F, G, H, J and I identify clades formed in P. guyannensis. The letters K, L, M, N, O, P, Q, R, S, T, U, V and X identify clades formed in P. cuvieri. Each terminal is identified for sample followed by country of origin. The numbers in the parentheses identified the locality in the Fig 1. We also identified the state of origin in the case of Brazil. Clusters recovered by BAPS are identified after the terminals for P. guyannensis (G1, G2, G3, G4) and for P. cuvieri (C1, C2, C3, C4). Inferred spliting dates are shown on the nodes.