F I G U R E 6 The distribution range of Turcinoemacheilus spp. in Western Asia (top). Maximum likelihood (ML) estimation of the phylogenetic relationships (bottom) based on the mitochondrial cytochrome c oxidase subunit 1 (COI) barcode region (Kimura 2-parameter (K2P) model, discrete gamma distribution for rate differences with five categories +G parameter = 0.1859). Nucleotide positions with less than 98% site coverage were eliminated, resulting in 540 analysed positions. Numbers of major nodes indicate bootstrap values of 1000 pseudoreplicates from the ML, MP, and NJ methods. "X" indicates the highest possible support from all three methods (>95). Blue bars to the right of the specimen labels indicate species boundary outcomes from the Poisson tree processes (PTP) followed by the results of the multi-rate PTP (mPTP) approach and the assemble species by automatic partitioning (ASAP) method at the rightmost.