“…The models of nucleotide substitution of the six DNA markers used were selected independently under the Akaike Information Criterion (
AIC ) using jModeTest v. 3.7 (
Posada 2008 ): GTR+I+Γ for
nrITS , GTR+Γ for matK , TrN+Γ for psbA-trnH , GTR+I for rbcL , TVM+Γ for rps16 and GTR+I for trnL-F . Detailed methods for
BI and
ML analyses could refer to the phylogenetic study conducted in
Yao et al (2021) , except that each of Markov Chain Monte Carlo (
MCMC ) analysis was run for 10,000,000 generations and sampling every 500 generations in
BI analysis. Number of generations in
BI analysis was sufficient, because the effective sample size (
ESS ) of all parameters were over 200 as evaluated in Tracer v. 1.6 (
Rambaut et al 2014 ), and the average standard deviations (
SD ) of split frequencies for the dataset was below 0.01.…”