Display omitted •The wild-ginger, series Sakawanum shows clinal floral variation.•Evolutionary history of four taxa was inferred by molecular markers.•Two admixture events may have occurred between morphologically distinct taxa.•Geographic range shift to glacial refugia could have allowed secondary contact.•A complicated demographic history may have shaped intertaxonomic clinal variation. Clinal variation is a major pattern of observed phenotypic diversity and identifying underlying demographic processes is a necessary step to understand the establishment of clinal variation. The wild ginger series Sakawanum (genus Asarum) comprises four taxa, which exhibit intertaxonomic clinal variation in calyx lobe length across two continental islands isolated by a sea strait. To test alternative hypotheses of the evolutionary history and to determine the implications for the formation of clinal variation, we conducted approximate Bayesian computation (ABC) analysis and ecological niche modeling (ENM). ABC analysis indicated that the scenario assuming multiple admixture events was strongly supported. This scenario assumed two admixture events occurred between morphologically distinct taxa, likely leading to the generation of intermediate taxa. One of the admixture events was estimated to have occurred during the last glacial maximum (LGM), during which the taxa were estimated to have formed a common refugia in southern areas by ENM analysis. Although four taxa are currently distributed allopatrically on different islands and trans-oceanic dispersal appears unlikely, the formation of a land bridge and the geographic range shift to refugia would have allowed secondary contact between previously isolated taxa. This study suggests that clinal variation can be shaped by demographic history including multiple admixtures due to climatic oscillations.