The well-known vicariance and dispersal models dominate in understanding the allopatric pattern for related species and presume the simultaneous occurrence of speciation and biogeographic events. However, the formation of allopatry may postdate the species divergence. We examined this hypothesis using DNA sequence data from 3 chloroplast fragments and 5 nuclear loci of Dipelta floribunda and D. yunnanensis, two shrub species with the circum Sichuan Basin distribution, combining the climatic niche modeling approach. The best-fit model supported by the approximate Bayesian computation (ABC) analysis indicated that, D. floribunda and D. yunnanensis diverged during the mid-Pleistocene period, consistent with the largest glacial period in the Qinghai-Tibet Plateau (QTP). The historically inter-specific gene flow was identified but seemed to have ceased after the last interglacial period (LIG), when the range of D. floribunda moved northward from the south of the Sichuan Basin. Further, populations of D. floribunda had expanded obviously in the north of the Sichuan Basin after the last glacial maximum (LGM). Relatively, the range of D. yunnanensis expanded before the LGM, reduced during the post-LGM especially in the north of the Sichuan Basin, reflecting the asynchronous responses of related species to the contemporary climate changes. Our results suggested that complex topography should be considered in understanding the distributional patterns even for closely related species and their demographic responses.