Fig. 4
Demographical history of T. fargesii var. fargesii and T. fargesii var. yunnanensis was simulated by (A) coalescent simulation (fastsimcoal2) and (B) isolation-with-migration (IM) model. Effective population size (Ne) for current and historical populations were depicted by the ratio of the width of the gray bars (A) or blank boxes (B), and the effective population sizes were showed in parentheses next to each population name; divergence time (Ma, million years ago) between the two varieties was presented on a horizontal dotted lines (A) or line (B); migration rate (m or 2Nm) were showed on a solid arrow, and the arrow indicated the direction of migration. The 95% highest posterior density intervals were showed with the double-sided arrows and blank boxes in grey for effective population size and divergence time in IM model (B). The detailed posterior interval for each demographic parameter can be found in Additional file 1: Table S11, S14. Asterisk (*) presented the migration rate at significant level (P < 0.05)