Abstract
The Japanese archipelago as a continental island of the Eurasia continent and harboring high levels of plant species diversity provides an ideal geographical setting for investigating vicariant allopatric speciation due to the sea-level fluctuations associated with climatic oscillations during the Quaternary. In this study, three chloroplast DNA regions and 14 nuclear loci were sequenced for 31 individuals from three populations of Torreya nucifera var. nucifera and 52 individuals from three populations of T. nucifera var. radicans. Population genetic analyses (Network, STRUCTURE and phylogeny) revealed that the genetic boundaries of the two varieties are distinct, with high genetic differentiation (F(ST)) of 0.9619 in chloroplast DNA and 0.6543 in nuclear loci. The relatively ancient divergence times between the two varieties were estimated to 3.03 Ma by DIYABC and 1.77 Ma by IMa2 when dated back to the late Pliocene and the early Pleistocene, respectively. The extremely weak gene flow (2Nm = 0.1) between the two varieties was detected by IMa2, which might be caused by their population expansion since the early Pleistocene (~2.0 Ma) inferred in the Bayesian skyline plots and DIYABC. Niche modeling showed that the two varieties had significant ecological differentiation (p < 0.001) since the Last Interglacial even earlier. These results demonstrate that vicariant allopatric speciation due to sea-level fluctuations may be a common mode of speciation in the Japanese archipelago. This finding provides insights into the understanding of species diversification in the Japanese Archipelago and even East Asian flora under climatic oscillations during the Quaternary.