2.3. Demographic analyses

The expected mismatch distribution analysis under pure demographic growth (Rogers & Harpending, 1992), as well as neutrality tests [Fu's Fs (Fu, 1997) and Tajima's D (Tajima, 1989)], was performed for each marker to test the deviation from demographic equilibrium. These analyses were employed in Arlequin 3.5.2.2 (Excoffier & Lischer, 2010). We also performed extended Bayesian skyline plot (EBSP) in BEAST2 (Bouckaert et al., 2014) using trnS‐trnG and PHYC datasets, assuming the models for nucleotide evolution and molecular clock identified for each partition. The substitution rate was available only for plastid sequences. Depending on the analysis, the MCMC runs were carried for 20 to 80 million generations sampling every 2,000 steps, with a 15% burn‐in. The results of EBSP were analyzed in TRACER 1.6 (available from http://beast.bio.ed.ac.uk/Tracer).