Figure S1 and Tables S1-S2

Figure S1. 99-sequence phylogeny (topology identical to Fig. 3) depicting ancestral state reconstructions for the caste determination phenotype (see Methods).

Table S1. Significance (P-values) from 2024 bootstrap permutations for all pairwise FST comparisons among the P. barbatus subgroups as calculated in Arlequin (see Methods).

Table S2. Significance (P-values) from 2024 bootstrap permutations for all pairwise FST comparisons among the P. rugosus subgroups as calculated in Arlequin (see Methods).

Appendix S1. Summary of discrepancies and proposed amendments regarding the distribution of P. barbatus and P. rugosus morphospecies.

Appendix S2. Detailed analysis of P. rugosus and P. barbatus phylogeography.