Population Structure in the Roundtail Chub (Gila robusta Complex) of the Gila River Basin as Determined by Microsatellites: Evolutionary and Conservation Implications

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

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 (FST 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.

Taxonomic Uncertainty and the Anomaly Zone: Phylogenomics Disentangle a Rapid Radiation to Resolve Contentious Species (Gila robusta Complex) in the Colorado River

Species are indisputable units for biodiversity conservation, yet their delimitation is fraught 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 North America. Nominal species designations were originally defined according to weakly diagnostic morphological differences, but these conflicted with subsequent genetic analyses. Given this ambiguity, the complex was re-defined as a single polytypic unit, with the proposed "threatened" status under the U.S. Endangered Species Act of two elements being withdrawn. Here we re-evaluated the status of the complex by utilizing dense spatial and genomic sampling (n = 387 and >22 k loci), coupled with SNP-based coalescent and polymorphism-aware phylogenetic models. In doing so, we found that all three species were indeed supported as evolutionarily independent lineages, despite widespread phylogenetic discordance. To juxtapose this discrepancy with previous studies, we first categorized those evolutionary mechanisms driving discordance, then tested (and subsequently rejected) prior hypotheses which argued phylogenetic discord in the complex was driven by the hybrid origin of Gila nigra. The inconsistent patterns of diversity we found within G. robusta were instead associated with rapid Plio-Pleistocene drainage evolution, with subsequent divergence within the "anomaly zone" of tree space producing ambiguities that served to confound prior studies. Our results not only support the resurrection of the three species as distinct entities but also offer an empirical example of how phylogenetic discordance can be categorized within other recalcitrant taxa, particularly when variation is primarily partitioned at the species level.

Geopolitical species revisited: genomic and morphological data indicate that the roundtail chub Gila robusta species complex (Teleostei, Cyprinidae) is a single species

The Gila robusta species complex in the Lower Colorado River Basin has a complicated taxonomic history. Recent authors have separated this group into three nominal taxa, G. robusta, G. intermedia, and G. nigra, however aside from location, no reliable method of distinguishing individuals of these species currently exists. To assess relationships within this group, we examined morphology of type specimens and fresh material, and used RADseq methods to assess phylogenetic relationship among these nominal species. Maximum likelihood and Bayesian inference tree building methods reveal high concordance between tree topologies based on the mitochondrial and nuclear datasets. Coalescent SNAPP analysis resolved a similar tree topology. Neither morphological nor molecular data reveal diagnostic differences between these species as currently defined. As such, G. intermedia and G. nigra should be considered synonyms of the senior G. robusta. We hypothesize that climate driven wet and dry cycles have led to periodic isolation of population subunits and subsequent local divergence followed by reestablished connectivity and mixing. Management plans should therefore focus on retaining genetic variability and viability of geographic populations to preserve adaptability to changing climate conditions.

Can members of the south-western Gila robusta species complex be distinguished by morphological features?

The goal for this project was to re-examine key morphological characters hypothesized to differentiate Gila intermedia, Gila robusta and Gila nigra and outline methods better suited for making species designations based on morphology. Using a combination of meristic counts, morphological measurements and geometric morphometrics, morphological dissimilarities were quantified among these three putative species. Traditional meristic counts and morphological measurements (i.e. distances between landmarks) were not useful for species identification. Geometric morphometrics, however, identified differences among species, while also suggesting an effect of geographic location on morphological variation. Using canonical variate analysis for the 441 fish sampled in this study, geometric morphometrics accurately predicted true group membership 100% of the time for G. nigra, 97% of the time for G. intermedia and 91% of the time for G. robusta. These results suggest that geometric morphometric analysis is necessary to identify morphological differences among the three species. Geometric morphometric analysis used in this study can be adopted by management officials as a tool to classify unidentified individuals.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.