Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.