Comparative Phylogeography of Three Marine Species with Different PLD Modes Reveals Two Genetic Breaks Across the Southern Caribbean Sea

Supplementary Materials: Figure S1. Bar graphs of the population ancestry coefficient from the STRUCTURE program indicating the probable number of populations (K = 2 - 6) for each species: (a) Acanthemblemaria rivasi, (b) Cittarium pica, and (c) Nerita tessellata; Figure S2.Correlogram plots of correlation coefficients and averages of estimated Ф_ST at each linear geographic distance class mark (right) for (a) Acanthemblemaria rivasi, (b) Cittarium pica, and (c) Nerita tessellata; Figure S3. Dendrograms of Ф_ST values constructed with the UPGMA method for A. rivasi, C. pica, and N. tessellata and spatial autocorrelation analysis between the Ф_ST matrices relating the localities across the southern and southwestern Caribbean, Colombia sector. (a) Analysis between A. rivasi and C. pica; (b) A. rivasi and N. tessellata; (c) C. pica and N. tessellata. FC = topological concordance factor; W = Kendall's statistic for topological (Wt) and genetic (Wg) congruence; rm = Mantel correlation coefficient; p = significance level. The black, blue, and red colors of the site names indicate high, medium, and low concordance among their locations in the dendrograms; Table S1. The genetic diversity indices calculated for Acanthemblemaria rivasi, Cittarium pica, and Nerita tessellata in five localities distributed across the southern Caribbean Sea, Colombia sector. Num: number of alleles, Eff_num: effective number of alleles, Ho: observed heterozygosity, Hs: heterozygosity within populations, and Gis: inbreeding coefficient.