Citizen-Science Data Reveal Global Diversity Patterns in the Gecko Genus Hemidactylus

Understanding how species richness is distributed across space is a central goal of macroecology and biogeography. The gecko genus Hemidactylus ranks among the most speciose and broadly distributed of all squamate genera, with species spanning tropical and subtropical environments on most continents. Yet the large-scale structure of its diversity has rarely been examined at a global scale. Here we assembled georeferenced occurrence data for the genus from the Global Biodiversity Information Facility (GBIF Access data is 14 March 2026) and, following cleaning and quality control, retained 143,858 records covering 157 currently recognized species (approximately 79% of the genus) to characterize broad-scale biodiversity patterns. Our analyses addressed the latitudinal diversity gradient, the distribution of species range sizes, continental richness, spatial clustering of diversity, and regional beta diversity. Richness declined steadily from the tropics toward the poles, producing a pronounced latitudinal diversity gradient (Pearson r = −0.71, p < 0.001). A generalized additive model captured a strong non-linear association between latitude and richness (deviance explained = 98.1%), consistent with tropical regions acting as principal centers of diversity. Range-size distributions were markedly right-skewed, with the majority of species restricted to comparatively small areas; however, these estimates co-varied strongly with per-species sampling effort (r = 0.90) and warrant cautious interpretation. We found only a weak, marginally significant association between range size and latitudinal midpoint, offering little support for Rapoport’s rule in the genus. At the continental scale, Hemidactylus diversity was overwhelmingly concentrated in Afro-Asian tropical regions (Asia, 93 species; Africa, 66), whereas Europe and Oceania supported few species (four each). Spatial analyses indicated significant autocorrelation in richness (Moran’s I = 0.56, p < 0.001) and well-defined diversity hotspots in East Africa, the Indian subcontinent, and Southeast Asia, and regional beta diversity was high (Jaccard dissimilarity up to 0.97), reflecting pronounced faunal turnover among continents.