Lineage isolation in the face of active gene flow in the coastal plant wild radish is reinforced by differentiated vernalisation responses

Table 2 Maximum-likelihood estimates (MLEs) and the 95 % highest posterior density (HPD) intervals of demographic parameters from pairwise IMa2 multi-locus analyses

	θ _S	θ _N	θ _A	m _N>S	m _S>N	2N_Nm_N	2N_Sm_S	t	T (approx.year)
L-model
MLE	0.094	0.082	1.582	51	11	3.188	1.351	0.027	17859
HPD95Lo	0.030	0.026	0.914	0	0	0	0	0.006	4173
HPD95Hi	0.246	0.200	2.642	201	165	7.162	7.064	0.071	47237
M-model 1
MLE	0.094	0.078	1.574	45	15	3.710	2.828	0.026	17922
HPD95Lo	0.034	0.026	0.914	0	0	0	0.073	0.006	3865
HPD95Hi	0.242	0.202	2.634	201	157	7.237	6.365	0.072	50251
M-model 2
MLE	0.094	0.086	1.582	37	11	3.581	3.140	0.027	18846
HPD95Lo	0.030	0.026	0.930	0	0	0.095	0	0.007	4623
HPD95Hi	0.250	0.202	2.634	205	167	6.920	7.090	0.072	50850
M-model 3
MLE	0.094	0.082	1.590	57	11	3.013	1.351	0.029	20269
HPD95Lo	0.034	0.030	0.898	0	0	0.163	0	0.007	4623
HPD95Hi	0.246	0.206	2.654	199	171	7.313	7.625	0.070	49428

θ _S θ _N and θ _A indicated the effective population size of the southern region, the northern region and the ancestral populations, respectively. The migration rate is indicated as backwards in time, e.g. migration parameter (m _N>S) indicates migration from northern to southern population in coalescent time (i.e. from population southern to northern forward in time). 2N _N m _N population migration rate from southern to northern population, 2N _S m _S population migration rate from to northern to southern population. The divergence time (T) was scaled by geometric means of substation rates per locus (1.498 × 10⁻⁶)

ISSN: 2730-7182