1000 randomly subsampled posterior trees from the Bayesian species tree analyses in BEAST v2.3.0 using ribosomal ITS and three chloroplast DNA regions(trnL-F, trnH-psbA, rps16 intron), the multispecies coalescent method implemented in *BEAST, uncorrelated lognormal relaxed clock models, and two secondary calibrations under normal prior distributions.