Figure Legends:
Figure 1 Geographical distribution of Cycas bifida andC. micholitzii and distribution of their haplotypes detected from
cpDNA (a), PHYP (b), AC5 (c), PPRC (d), AAT(e) and SAMS (f). Population codes refer to Supplementary Table
S1.
Figure 2 Haplotype Networks based on the combined cpDNA data
and five low copy nuclear genes. Circle sizes are proportional to
haplotype frequency. Circle colours indicate the species origin (white =C. bifida, black = C. micholitzii ). Hypothetical
intermediates are indicated by small grey points which name after
“mv”. The steps more than one were marked near the branches.
Figure 3 Bayesian inference (K = 2, 3, 4) of microsatellite
phenotype from 11 populations of Cycas bifida and C.
micholitzii (a ).
Principal
coordinates analysis (PCoA)
(b )
and the boundaries detected using the BARRIER program based on matrices
of N ei’s (1983) unbiased genetic distance under 1000 replications
(c ).
Figure 4 Results of the identity test for Cycas bifidaand C. micholitzii . The x -axis indicates values of I and
D, and y -axis indicates number of randomizations. Null
distributions are shown by grey bars for D and black bars for I, and
arrows indicate the observed values of niche similarity (grey line for D
and black line for I).
Figure 5 Marginal distribution of posterior probabilities for
demographic parameters estimated by IM model. (a) Marginal distribution
of the effective population sizes of Cycas bifida , C.
micholitzii and their ancestral population; (b) Migration rate betweenC. bifida and C. micholitzii ; (c) The divergence time
between C. bifida and C. micholitzii .
Figure 6 Figure plot of geographical distance against genetic
distance for Cycas bifida (a), C. micholitizii (b) and the
two species (c).
Figure 7 Extended Bayesian Skyline Plot of Cycas bifida(a) and C. micholitzii (b) based on cpDNA and five low copy
nuclear genes for the effective population size fluctuation throughout
time. Black line, median estimation; area between gray lines, 95%
confidence interval.