The most accepted hypothesis has suggested that the fauna in Taiwan Island originated from South China, but some studies supported the Japan, Ryukyu Archipelago, and Taiwan Islands as a unique biogeographical district. This study examines whether the populations of freshwater snail Semisulcospira libertina in Taiwan are closer to those in Japan based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. Our study shows the populations in North Taiwan originated from Japan and the cyclic glacial caused the migrations among islands and continent repeatedly; the populations in South Taiwan might originate from South China or South Asia. Our results will not only affect the conclusions in phylogeography of freshwater species in Taiwan but also change the sampling plans in the future studies about evolutionary of freshwater species in East Asia.

We examined the genetic variation and phylogeographic relationships among 10 populations of Lunella granulata from mainland China, Penghu Archipelago, Taiwan Island, and Japan using mitochondrial COI and 16S markers. A total of 45 haplotypes were obtained in 112 specimens, and relatively high levels of haplotype diversity (h = 0.903) and low levels of nucleotide diversity (π = 0.0046) were detected. Four major phylogenetic lineage clusters were revealed and were concordant with their geographic distribution, agreeing with the haplotype network. These results suggested that geographic barrier isolating effects were occurring among the populations. This hypothesis was also supported by a significant genetic differentiation index (FST = 0.709) and by a spatial analysis of molecular variance (SAMOVA) analysis. A mismatch distribution analysis, neutrality tests and Bayesian skyline plots found a single significant population expansion. This expansion occurred on the coast of mainland China before 20–17 ka. Consequently, although the dispersal ability of the planktonic stage and the circulation of ocean currents generally promote genetic exchanges among populations, L. granulata has tended to maintain distinct genetic groups that reflect the respective geographic origins of the constituent lineages. Although the circulation of ocean currents, in principle, may still play a role in determining the genetic composition of populations, long-distance migration between regions is difficult even at the planktonic stage.

The mitochondrial DNA cytochrome c oxidase subunit I sequences from 95 specimens of Semisulcospira libertina in Taiwan were identified as two major phylogroups, exhibiting a southern and northern distribution, north of Formosa Bank and south of Miaoli Plateau. The genetic distance between these two phylogroups was 12.20%, and the distances within-phylogroups were 4.97 and 5.56%. According to a molecular clock of 1.56% per lineage per million years, the divergence time between these two major phylogroups was estimated at 4.94 million years ago (mya), with the two phylogroups forming at 3.64 and 3.75 mya, respectively. Moreover, the geological events have suggested that Taiwan Island emerged above sea level at 4-5 mya, and became its present shape at 2 mya. These results suggested that these two phylogroups might originate from two independent ancestral populations or divergent before colonizing Taiwan. Within South phylogroup, the initial colonization was hypothesized to be in Kaoping River (WT), followed by its northward. The high divergence between south- and north of WT River was influenced by the formation of the Kaoping foreland basins. Within North phylogroup, the colonization was from central sub-region through paleo-Miaoli Plateau to northern and northeastern sub-regions. This study showed that the landform changes might have shaped the genetic structure of S. libertina in concert. Apparently, two cryptic species or five different genetic stocks of S. libertina could be identified; these results are useful for the evaluation and conservation of S. libertina in Taiwan.