Supplementary figures. Figure S1. Graphical representation of the STRUCTURE output for the identification of unknown NGO samples for K = 2. Figure S2a. minimum Spanning Network (MSN) of C. dorsalis based on 16 composite haplotypes (collapsed haplotypes). Red circles indicate the median vector or hypothesized (ancestral) haplotypes and the length of each branch connecting sampled haplotypes indicates the genetic distance between haplotypes. b. Minimum Spanning Network (MSN) of C. callipygus based on 27 composite haplotypes (collapsed haplotypes). C. minimum Spanning Network (MSN) of P. monticola based on 25 composite haplotypes (collapsed haplotypes). Figure S3. Cytochrome c oxidase subunit I neighbour-joining bootstrap consensus phylogeny based on Tamura-Nei distances. Bootstrap values of 75% or greater are indicated at the relevant node. Figure S4a-b. Graphical representation of the STRUCTURE output for C. dorsalis for K = 2. Figure S5a-b. Graphical representation of the STRUCTURE output for C. callipygus for K = 6. Figure S6a-b. Graphical representation of the STRUCTURE output for P. monticola for K = 2. Figure S7. Environmental variables used as predictors in MAXENT and GDM analyses. Figure S8. Jackknife of regularized training gain for C. dorsalis. Fig. S9. Jackknife of regularized training gain for C. callipygus. Figure S10. Jackknife of regularized training gain for P. monticola. Figure S11. MAXENT predicted species distributions for: a) C. dorsalis b) C.callipygus and c) P. monticola. Figure S12. GDM functional response curves for C. callipygus. Figure S13. GDM functional response curves for P. monticola (PPTX 2695 kb)