An integrative insight into the diversity, distribution, and biogeography of the freshwater endemic Glade of the Ponticola syrman group (Teleostei: Gobiidae) in the Caucasus biodiversity hotspot

Freshwater habitats of the Caucasus biodiversity hotspot represent a center of endemism for the gobiid genus Ponticola IIjin, 1927. Hitherto, large-scale molecular studies, owing to restricted taxon and geographical sampling, have failed to give an elaborate picture of diversity and evolutionary history of these species. Here, to contribute to filling this gap, we assessed taxonomic diversity, phylogeography and evolutionary history for the south Caspian populations of Ponticola presently classified as P. iranicus and P. patimari, using an integrative taxonomic approach comprising an entire geographic range sampling, and analyses of mitochondrial DNA haplotypes, the head lateral line system, otolith shape, and meristic and morphometric variation. All freshwater samples of the P. syrman group belong to a monophyletic Glade with two main subclades: a small subclade confined to the upper Sefidroud sub-basin including the type locality of P. iranicus and a large subclade with three geographically constrained haplogroups (Hg1, Hg2, and Hg3), comprising the rest of the distribution. Hg1 showed an eastern distribution including the type locality of P. patimari, while Hg2 and Hg3 are sister groups with central and western-central distributions, respectively. The freshwater Glade diverged from P. syrman during the Tyurkyanian low stand (similar to 150 m b.s.l. lasting similar to 0.1 Myr), while the divergence of P. iranicus and P. patimari and radiations within P. patimari took place during the Bakunian high stand (up to 50 m a.s.l. lasting similar to 378-480 kya). Species delimitation analyses indicated two distinct species, corresponding to each main subclade. Although the otolith shape and lateral line analyses did not reflect with phylogeographic pattern, PCA and DFA plots of meristic and morphometric data showed a clear separation of the two major subclades corresponding to P. iranicus and P. patimari, suggesting the presence of significant morphological variation meriting formal taxonomic recognition. Overall, our findings (i) reveal the presence of two freshwater endemic species in the P. syrman group, and pending further investigation, hypothesize the presence of a third cryptic species; (ii) revise and document a narrow distributional range and low diversity for P. iranicus, in contrast to a wider distributional range and high diversity for P. patimari; (iii) suggest that the climatic oscillations of the Pleistocene were associated with the cladogenesis within the P. syrman group; and (iv) allowed for the recognition of conservation units and proposition of management measures.