Geographic destiny trumps taxonomy in the Roundtail Chub, Gila robusta species complex (Teleostei, Leuciscidae)

• Published in: Scientific reports 2023-09, Vol.13 (1), p.15810-15810, Article 15810
• Main Authors: Suchocki, Christopher R., Ka‘apu-Lyons, Cassie, Copus, Joshua M., Walsh, Cameron A. J., Lee, Anne M., Carter, Julie Meka, Johnson, Eric A., Etter, Paul D., Forsman, Zac H., Bowen, Brian W., Toonen, Robert J.
• Format: Article
• Language: English
• Subjects: 704/158; 704/172; Conservation; Divergence; Genetic diversity; Genetic variance; Genomes; Geography; Gila robusta; Humanities and Social Sciences; Morphology; multidisciplinary; Phylogeny; Rivers; Science; Science (multidisciplinary); Single-nucleotide polymorphism; Species; Taxonomy; Watersheds
• ISSN: 2045-2322; 2045-2322
• DOI: 10.1038/s41598-023-41719-9

The Gila robusta species complex in the lower reaches of the Colorado River includes three nominal and contested species ( G. robusta, G. intermedia, and G. nigra ) originally defined by morphological and meristic characters. In subsequent investigations, none of these characters proved diagnostic, and species assignments were based on capture location. Two recent studies applied conservation genomics to assess species boundaries and reached contrasting conclusions: an ezRAD phylogenetic study resolved 5 lineages with poor alignment to species categories and proposed a single species with multiple population partitions. In contrast, a dd-RAD coalescent study concluded that the three nominal species are well-supported evolutionarily lineages. Here we developed a draft genome (~ 1.229 Gbp) to apply genome-wide coverage (10,246 SNPs) with nearly range-wide sampling of specimens ( G. robusta N = 266, G. intermedia N = 241, and G. nigra N = 117) to resolve this debate. All three nominal species were polyphyletic, whereas 5 of 8 watersheds were monophyletic. AMOVA partitioned 23.1% of genetic variance among nominal species, 30.9% among watersheds, and the Little Colorado River was highly distinct ( F ST ranged from 0.79 to 0.88 across analyses). Likewise, DAPC identified watersheds as more distinct than species, with the Little Colorado River having 297 fixed nucleotide differences compared to zero fixed differences among the three nominal species. In every analysis, geography explains more of the observed variance than putative taxonomy, and there are no diagnostic molecular or morphological characters to justify species designation. Our analysis reconciles previous work by showing that species identities based on type location are supported by significant divergence, but natural geographic partitions show consistently greater divergence. Thus, our data confirm Gila robusta as a single polytypic species with roughly a dozen highly isolated geographic populations, providing a strong scientific basis for watershed-based future conservation.