Unraveling the diversification and systematic puzzle of the highly polymorphic Psammobates tentorius (Bell, 1828) complex (Reptilia: Testudinidae) through phylogenetic analyses and species delimitation approaches

The high level of phenotypic diversity in southern African tent tortoises (Psammobates tentorius complex) has for decades prevented systematists from developing a stable taxonomy for the group. Here, we used a comprehensive DNA sequence dataset (mtDNA: Cytb, ND4, ND4 adjacent tRNA-His, and tRNA-Ser, 12S, 16S; and nDNA: PRLR gene) of 455 specimens, and the latest phylogenetic and species delimitation analytical procedures, to unravel the long-standing P. tentorius complex systematic puzzle. Our results for mtDNA and nDNA were incongruent, with the poorly supported nDNA phylogeny differentiating the three recognized subspecies, and showing potential hybridization in some regions. In contrast, the concatenated mtDNA phylogeny identified seven operational taxonomic units, with strong support. Clades 1, 4, 5, and 7 corresponded to tortoises identified as P. t. tentorius, clade 3 to P. t. trimeni, and clades 2 and 6 to P. t. verroxii. Our analyses showed conflicting topologies for the placement of C6 (P. t. verroxii north of the Orange River), with stronger support for it being sister to C2 + C3 than to the other clades. Clades 1, 2, and 6 had significantly higher genetic diversity than clades 3, 4, 5, and 7, perhaps because these clades inhabit substantially larger areas. The potential for future cladogenic radiations seems high in C1 and C6, particularly in C6 for which the within-clade diversification level was highest. Further research involving microsatellite DNA, phylogeographic evaluations, and morphological variation among clades is crucial for understanding the adaptive radiation of the P. tentorius complex and for modifying their taxonomy.