Molecular systematics, species delimitation and diversification patterns of the Phyllodactylus lanei complex (Gekkota: Phyllodactylidae) in Mexico.

The description of cryptic gecko species worldwide has revealed both that many putative species are, in fact, conformed by a complex of morphologically conserved species that are genetically distinct and highly divergent, and that gecko species diversity could be underestimated. The taxonomy and species delimitation of geckos belonging to the genus Phyllodactylus is still controversial, 16 of which are distributed in Mexico and 13 are endemic. Although the large morphological variation shown by the Phyllodactylus species from Mexico has been amply documented, little is known about their genetic diversity and evolutionary relationships, and much less regarding cryptic speciation. Here, we included the most comprehensive sampling of populations and species of the Phyllodactylus lanei complex distributed in Mexico, and applied an analytical approach that included probabilistic phylogenetic analyses, jointly with species delimitation methods and Bayesian putative species validation analysis. Our results suggest the existence of 10 lineages within the complex, supporting the existence of cryptic species, and in great contrast with the current taxonomic proposal that includes only four subspecies. The most recent common ancestor (MRCA) for the P. lanei clade originated on the Early Eocene (∼54Mya), along the southern coasts of Mexico, followed by the highest diversification of the complex MRCA during the Eocene (34-56Mya). Lineages subsequently dispersed and diversified towards the northwest, and the diversification process ended with the most recent lineages inhabiting two islands on the coasts of Nayarit (Miocene; 5.5-23Mya). Our results highlight three vicariant events associated with the evolution of the lineages, two of them intimately related to the formation of the Sierra Madre del Sur and the Transmexican Volcanic Belt mountain ranges, main geographic barriers that isolated and facilitated the divergence and speciation in this group of geckos. Finally, we propose that there are 10 species in the P. lanei complex, from which four represent taxonomic changes and six are new species and require a formal description. We acknowledge that more analyses, including a detailed evaluation of morphological characters and use of more unlinked nuclear loci with enough variability, are needed to further support their taxonomic description.