Phylogeography and evolutionary history of the Panamic Clingfish Gobiesox adustus in the Tropical Eastern Pacific

[Display omitted] •Gobiesox adustus individuals is a monophyletic group.•Gobiesox adustus populations conform four genetically divergent, allopatric groups restricted to biogeographic province.•Analyses confirm the existence of four independent evolutionary units.•Habitat discontinuities, oceanographic processes and ecological characteristics seems to be causes population structure. The Panamic Clingfish Gobiesox adustus is widely distributed in the Tropical Eastern Pacific (TEP), from the central Gulf of California, Mexico to Ecuador, including the oceanic Revillagigedo Archipelago, and Isla del Coco. This cryptobenthic species is restricted to very shallow rocky-reef habitats. Here, we used one mitochondrial and three nuclear DNA markers from 155 individuals collected across the distribution range of the species in order to evaluate if geographically structured populations exist and to elucidate its evolutionary history. Phylogenetic analyses recovered a monophyletic group, with four well-supported, allopatric subgroups. Each subgroup corresponded to one of the following well-known biogeographic regions/provinces: 1) the Revillagigedo Archipelago, 2) the Cortez + Mexican provinces (Mexico), 3) the Panamic province (from El Salvador to Ecuador), and 4) Isla del Coco. A molecular-clock analysis showed a mean date for the divergence between clade I (the Revillagigedos and Cortez + Mexican provinces) and clade II (Panamic province and Isla del Coco) in the Pliocene, at ca. 5.33 Mya. Within clade I, the segregation between the Revillagigedos and Cortez + Mexican province populations was dated at ca. 1.18 Mya, during the Pleistocene. Within clade II, the segregation between samples of Isla del Coco and the Panamic province samples was dated at ca. 0.77 Mya, during the Pleistocene. The species tree, Bayesian species delimitation tests (BPP and STACEY), the ΦST, AMOVA, and the substantial genetic distances that exist between those four subgroups, indicate that they are independent evolutionary units. These cladogenetic events seem to be related to habitat discontinuities, and oceanographic and geological processes that produce barriers to gene flow for G. adustus, effects of which are enhanced by the intrinsic ecological characteristics of this species.