To further gain insights into the population structure of these species, we used the BAPS software to establish the genetically differentiated groups and determine the amount of admixture between these groups. We ran the BAPS clustering model using the hierarchical manner to identify substructures inside the main clusters64. To find the optimal clustering, we ran five independent iterations with the prior upper bound of the number of clusters set to 35; clustering was performed with 4 levels in the hierarchy. The first level gave 10 clusters whereas the fourth yielded 20 clusters, for which we conducted a mixture analysis. We ran this analysis considering the minimum number of individuals for a population as one using parameters described by Castillo-Ramirez S. et al.65.