Disclosure statement:Authors have nothing to disclose.Abstract (242/250) Species are an indisputable unit for biodiversity conservation, yet their delimitation is fraught 2 with both conceptual and methodological difficulties. A classic example is the taxonomic controversy surrounding the Gila robusta complex in the lower Colorado River of southwestern 4 North America. Three species were originally defined according to subtle morphological differences with weak diagnostic power, with traditional genetic analyses failing to support the 6 nominal species. Consequently, the complex was re-defined as a single polytypic unit, with the proposed 'threatened' status of two putative species withdrawn at the federal level. Here, we 8 utilized dense spatial and genomic sampling (N=387 and >22k loci) to re-evaluate the status of the complex, based on SNP-based coalescent and polymorphism-aware phylogenetic models. In 10 doing so, all three species were supported as evolutionarily independent lineages, despite widespread phylogenetic discordance. To understand this discrepancy with past studies, we 12 categorized evolutionary mechanisms driving discordance. We tested (and subsequently rejected) prior hypotheses suggesting that phylogenetic discord in the complex was hybridization-driven. 14 Instead, we found the G. robusta complex to have diverged within the 'anomaly zone' of tree space and, as such, have accumulated inconsistent patterns of diversity which have confounded 16 prior studies. After extending these analyses with phylogeographic modeling, we propose that this is reflective of a rapid radiation promoted by Plio-Pleistocene tectonism. Our results not only 18 support resurrection of the three species as distinct entities, but also offer an empirical example of how phylogenetic discordance can be categorized in other recalcitrant taxa. 20 between G. robusta from the upper and lower basins of the Colorado River ecosystem (Chafin et 130 al. 2019). One species, Gila elegans, is extirpated from the wild, thus samples were provided by the 132 Southwestern Native Aquatic Resources and Recovery Center (Dexter, NM). The genus Ptychocheilus served to root the Gila clade within the broader context of western leuciscids 134 (Schönhuth et al. 2012, 2014, 2018). 8 136 Reduced-representation sequencing Genomic DNA was extracted using either PureGene® or DNeasy® kits (Qiagen Inc.) and 138 quantified via fluorometer (Qubit™; Thermo-Fisher Scientific). Library preparations followed the published ddRAD protocol (Peterson et al. 2012). Restriction enzyme and size-selection 140 ranges were first screened using an in silico procedure (Chafin et al. 2018), with the target fragment sizes further optimized by quantifying digests for 15 representative samples on an 142 Agilent 2200 TapeStation. Final library preparations were double-digested using a high-fidelity PstI (5'-CTGCAG-3') and MspI (5'-CCGG-3') following manufacturer's protocols (New 144 England Biosciences). Digests were purified using bead purification (Ampure XP; Beckman-Coul...