Rangewide molecular structuring in the Utah sucker (Catostomus ardens)

Heterochronic shift in gene expression leads to ontogenetic morphological divergence between two closely related polyploid species

Heterochrony-alteration to the rate or timing of development-is an important mechanism of trait differentiation associated with speciation. Heterochrony may explain the morphological divergence between two polyploid species, June sucker (Chasmistes liorus) and Utah sucker (Catostomus ardens). The larvae of both species have terminal mouths; however, as adults, June sucker and Utah sucker develop subterminal and ventral mouths, respectively. We document a difference in the timing of shape development and a corresponding change in the timing of gene expression, suggesting the distinctive mouth morphology in June suckers may result from paedomorphosis. Specifically, adult June suckers exhibit an intermediate mouth morphology between the larval (terminal) and ancestral (ventral) states. Endemic and sympatric Chasmistes/Catostomus pairs in two other lakes also are morphologically divergent, but genetically similar. These species pairs could have resulted from the differential expression of genes and corresponding divergence in trait development. Paedomorphosis may lead to adaptive diversification in Catostomids.

Complete mitochondrial genomes of June sucker and Utah sucker (Chasmistes liorus and Catostomus ardens)

The relationship between June sucker (Chasmistes liorus, Jordan, 1878) and Utah sucker (Catostomus ardens, Jordan & Gilbert, 1881) has been a matter of controversy since the mid 1900s. Chasmistes liorus is endemic to Utah Lake, UT and has a subterminal mouth adapted for pelagic feeding. Catostomus ardens is widely distributed throughout the Bonneville Basin and Upper Snake River Basin and has a ventral mouth adapted for benthic feeding. Chasmistes has been recognized as a separate ancient genus. Despite being morphologically distinct, no study has successfully identified residual genetic markers that separate these species. Of these studies, several have used a subset of mitochondrial genes, but no study has analyzed the complete mitochondrial genomes (mitogenomes) of these suckers (Pisces: Catostomidae). To further explore the genetic relationships between these species, we report the complete mitogenomes of Chasmistes liorus and Catostomus ardens. DNA was sequenced using an Illumina HiSeq 2500 system and mitogenomes were assembled and annotated using Geneious v. 2021.2 and MitoAnnotator, respectively. The mitogenomes of Chasmistes liorus and Catostomus ardens are both 16,623 bp and are ∼0.072% divergent. We examine the phylogenetic relationship between Chasmistes liorus and Catostomus ardens using 33 mitogenomes, representing 16 species, from Catostomidae. Our data suggest that Chasmistes liorus is sister to Catostomus ardens. Additional samples from multiple localities and/or cohorts of these species will allow us to better resolve the complicated phylogenetic relationships between these species.

Phylogeny and divergence times of suckers (Cypriniformes: Catostomidae) inferred from Bayesian total-evidence analyses of molecules, morphology, and fossils

Catostomidae ("suckers") is a diverse (76 species) and broadly distributed family of Holarctic freshwater fishes with a rich fossil record and a considerable number (∼35%) of threatened and imperiled species. We integrate DNA sequences (three mitochondrial genes, three nuclear genes), morphological data, and fossil information to infer sucker phylogenetic relationships and divergence times using Bayesian "total-evidence" methods, and then test hypotheses about the temporal diversification of the group. Our analyses resolved many nodes within subfamilies and clarified Catostominae relationships to be of the form ((Thoburniini, Moxostomatini), (Erimyzonini, Catostomini)). Patterns of subfamily relationships were incongruent, but mainly supported two placements of the Myxocyprininae; distinguishing these using Bayes factors lent strongest support to a model with Myxocyprininae sister to all remaining sucker lineages. We improved our Bayesian total-evidence dating analysis by excluding problematic characters, using a clock-partitioning scheme identified by Bayesian model selection, and employing a fossilized birth-death tree prior accommodating morphological data and fossils. The resulting chronogram showed that suckers evolved since the Late Cretaceous-Eocene, and that the Catostomini and Moxostomatini clades have accumulated species diversity since the early to mid-Miocene. These results agree with the fossil record and confirm previous hypotheses about dates for the origins of Catostomide and catostomine diversification, but reject previous molecular hypotheses about the timing of divergence of ictiobines, and between Asian-North American lineages. Overall, our findings from a synthesis of multiple data types enhance understanding of the phylogenetic relationships, taxonomic classification, and temporal diversification of suckers, while also highlighting practical methods for improving Bayesian divergence dating models by coupling phylogenetic informativeness profiling with relaxed-clock partitioning.

Introgressive Hybridization and the Evolution of Lake-Adapted Catostomid Fishes

Hybridization has been identified as a significant factor in the evolution of plants as groups of interbreeding species retain their phenotypic integrity despite gene exchange among forms. Recent studies have identified similar interactions in animals; however, the role of hybridization in the evolution of animals has been contested. Here we examine patterns of gene flow among four species of catostomid fishes from the Klamath and Rogue rivers using molecular and morphological traits. Catostomus rimiculus from the Rogue and Klamath basins represent a monophyletic group for nuclear and morphological traits; however, the Klamath form shares mtDNA lineages with other Klamath Basin species (C. snyderi, Chasmistes brevirostris, Deltistes luxatus). Within other Klamath Basin taxa, D. luxatus was largely fixed for alternate nuclear alleles relative to C. rimiculus, while Ch. brevirostris and C. snyderi exhibited a mixture of these alleles. Deltistes luxatus was the only Klamath Basin species that exhibited consistent covariation of nuclear and mitochondrial traits and was the primary source of mismatched mtDNA in Ch. brevirostris and C. snyderi, suggesting asymmetrical introgression into the latter species. In Upper Klamath Lake, D. luxatus spawning was more likely to overlap spatially and temporally with C. snyderi and Ch. brevirostris than either of those two with each other. The latter two species could not be distinguished with any molecular markers but were morphologically diagnosable in Upper Klamath Lake, where they were largely spatially and temporally segregated during spawning. We examine parallel evolution and syngameon hypotheses and conclude that observed patterns are most easily explained by introgressive hybridization among Klamath Basin catostomids.

Pluvial Drainage Patterns and Holocene Desiccation Influenced the Genetic Architecture of Relict Dace, Relictus solitarius (Teleostei: Cyprinidae)

Changing drainage patterns have played a significant role in the evolution of western North American aquatic taxa. Relict dace, Relictus solitarius, is a Great Basin endemic cyprinid with a native range that is restricted to four valleys in eastern Nevada. Relictus solitarius now occupies spring systems that are the remnants of Pleistocene-era pluvial lakes, although it may have occurred in the area for much longer. Here we use mitochondrial DNA sequence data to assess range-wide genetic diversity of R. solitarius, and to estimate divergence times to determine whether pluvial drainages played an important role in shaping intraspecific genetic diversity. Genetic diversification within R. solitarius began during the early to mid-Pleistocene, separating populations within two sets of valleys (Butte/Ruby and Goshute/Steptoe). Additional diversification in each of the two sets of valleys occurred more recently, in the mid- to late-Pleistocene. Holocene desiccation has further isolated populations, and each population sampled contains unique mtDNA haplotypes. Pluvial drainage patterns did contribute to the genetic structure observed within R. solitarius, but most of the intraspecific diversification does not appear to be associated with the Last Glacial Maximum. Holocene desiccation has also contributed to the observed genetic structure. The relict dace populations we sampled are all unique, and we recommend that future management efforts should strive to preserve as much of the genetic diversity as possible.

Fish faunal provinces of the conterminous United States of America reflect historical geography and familial composition

Although the conterminous USA has a long history of ichthyological exploration, the description of biogeographical provinces has been ad hoc. In this study we quantitatively determined fish faunal provinces and interpreted them in the context of the geological history of North America. We also evaluated influences of major river basin occupancy and contemporary environmental factors on provincial patterns. Our data set comprised 794 native fishes, which we used to generate a presence and absence matrix for U.S. Geological Survey (USGS) four-digit hydrologic units. Three nested data sets were analysed separately: primary freshwater families, continental freshwater families (including primary and secondary families) and all freshwater families (including primary, secondary and peripheral families). We used clustering analysis to delimit faunal breaks and one-way analysis of similarity (ANOSIM) to determine significance among clusters (i.e. provinces). We used an indicator-species analysis to identify species that contributed most to province delineations and a similarity-percentage (SIMPER) analysis to describe the relative influence of representatives from each category (i.e. primary, secondary, peripheral) on provincial boundaries. Lastly, we used a parsimony redundancy analysis to determine the roles of historical (i.e. major river basin) and contemporary environmental factors in shaping provinces. Analysis of the nested data sets revealed lessening provincial structure with inclusion of more families. There were 10 primary freshwater provinces, 9 continental freshwater provinces and 7 all freshwater provinces. Major basin occupancy, but not contemporary environmental factors, explained substantial variance in faunal similarities among provinces. However, provincial boundaries did not conform strictly to modern river basins, but reflected river-drainage connections of the Quaternary. Provinces represent broad-scale patterns of endemism and provide a starting point for future studies. Relative malleability of province boundaries in the continental interior highlights this region as biogeographically diverse and dynamic. Interior-core provinces of this region (Central Gulf Coastal Plains, Northern Interior) have not been recognized previously and warrant further study.

Highly variable reproductive isolation among pairs of Catostomus species

Hybridization between diverged taxa tests the strength of reproductive isolation and can therefore reveal mechanisms of reproductive isolation. However, it remains unclear how consistent reproductive isolation is across species' ranges and to what extent reproductive isolation might remain polymorphic as species diverge. To address these questions, we compared outcomes of hybridization across species pairs of Catostomus fishes in three rivers in the Upper Colorado River basin, where an introduced species, C. commersoni, hybridizes with at least two native species, C. discobolus and C. latipinnis. We observed substantial heterogeneity in outcomes of hybridization, both between species pairs and across geographically separate rivers within each species pair. We also observed hybridization of additional related species with our focal species, suggesting that reproductive isolation in this group involves interactions of multiple evolutionary and ecological factors. These findings suggest that a better understanding of the determinants of variation in reproductive isolation is needed and that studies of reproductive isolation in hybrids should consider how the dynamics and mechanisms of reproductive isolation vary over ecological space and over evolutionary time. Our results also have implications for the conservation and management of native catostomids in the Colorado River basin. Heterogeneity in outcomes of hybridization suggests that the threat posed by hybridization and genetic introgression to the persistence of native species probably varies with extent of reproductive isolation, both across rivers and across species pairs.

Investigating the effects of Pleistocene events on genetic divergence within Richardsonius balteatus, a widely distributed western North American minnow

Background

Biogeographers seek to understand the influences of global climate shifts and geologic changes to the landscape on the ecology and evolution of organisms. Across both longer and shorter timeframes, the western North American landscape has experienced dynamic transformations related to various geologic processes and climatic oscillations, including events as recently as the Last Glacial Maximum (LGM; ~20 Ka) that have impacted the evolution of the North American biota. Redside shiner is a cyprinid species that is widely distributed throughout western North America. The species’ native range includes several well-documented Pleistocene refugia. Here we use mitochondrial DNA sequence data to assess phylogeography, and to test two biogeographic hypotheses regarding post-glacial colonization by redside shiner: 1) Redside shiner entered the Bonneville Basin at the time of the Bonneville Flood (Late Pleistocene; 14.5 Ka), and 2) redside shiner colonized British Columbia post-glacially from a single refugium in the Upper Columbia River drainage.

Results

Genetic diversification in redside shiner began in the mid to late Pleistocene, but was not associated with LGM. Different clades of redside shiner were distributed in multiple glacial age refugia, and each clade retains a signature of population expansion, with clades having secondary contact in some areas.

Conclusions

Divergence times between redside shiner populations in the Bonneville Basin and the Upper Snake/Columbia River drainage precedes the Bonneville Flood, thus it is unlikely that redside shiner invaded the Bonneville Basin during this flooding event. All but one British Columbia population of redside shiner are associated with the Upper Columbia River drainage with the lone exception being a population near the coast, suggesting that the province as a whole was colonized from multiple refugia, but the inland British Columbia redside shiner populations are affiliated with a refugium in the Upper Columbia River drainage.

Stream hierarchy defines riverscape genetics of a North American desert fish

Global climate change is apparent within the Arctic and the south-western deserts of North America, with record drought in the latter reflected within 640,000 km(2) of the Colorado River Basin. To discern the manner by which natural and anthropogenic drivers have compressed Basin-wide fish biodiversity, and to establish a baseline for future climate effects, the Stream Hierarchy Model (SHM) was employed to juxtapose fluvial topography against molecular diversities of 1092 Bluehead Sucker (Catostomus discobolus). MtDNA revealed three geomorphically defined evolutionarily significant units (ESUs): Bonneville Basin, upper Little Colorado River and the remaining Colorado River Basin. Microsatellite analyses (16 loci) reinforced distinctiveness of the Bonneville Basin and upper Little Colorado River, but subdivided the Colorado River Basin into seven management units (MUs). One represents a cline of three admixed gene pools comprising the mainstem and its lower-gradient tributaries. Six others are not only distinct genetically but also demographically (i.e. migrants/generation <9.7%). Two of these (i.e. Grand Canyon and Canyon de Chelly) are defined by geomorphology, two others (i.e. Fremont-Muddy and San Raphael rivers) are isolated by sharp declivities as they drop precipitously from the west slope into the mainstem Colorado/Green rivers, another represents an isolated impoundment (i.e. Ringdahl Reservoir), while the last corresponds to a recognized subspecies (i.e. Zuni River, NM). Historical legacies of endemic fishes (ESUs) and their evolutionary potential (MUs) are clearly represented in our data, yet their arbiter will be the unrelenting natural and anthropogenic water depletions that will precipitate yet another conservation conflict within this unique but arid region.

Y chromosome phylogeny for cutthroat trout (Oncorhynchus clarkii) subspecies is generally concordant with those of other markers

Sequence divergence was evaluated in the non-recombining, male-specific OmyY1 region of the Y chromosome among the subspecies of cutthroat trout (Oncorhynchus clarkii) in the western United States. This evaluation identified subspecies-discriminating OmyY1-haplotypes within a ∼1200bp region of the OmyY1 locus and localized the region to the end of the Y chromosome by FISH analysis. OmyY1 sequences were aligned and used to reconstruct a phylogeny of the cutthroat trout subspecies and related species via maximum-parsimony and Bayesian analyses. In the Y-haplotype phylogeny, clade distributions generally corresponded to the geographic distributions of the recognized subspecies. This phylogeny generally corresponded to a mitochondrial tree obtained for these subspecies in a previous study. Both support a clade of trout vs. Pacific salmon, of rainbow trout, and of a Yellowstone cutthroat group within the cutthroat trout. In our OmyY1 tree, however, the cutthroat "clade", although present topologically, was not statistically significant. Some key differences were found between trees obtained from the paternally-inherited OmyY1 vs. maternally-inherited mitochondrial haplotypes in cutthroat trout compared to rainbow trout. Other findings are: The trout OmyY1 region evolves between 3 and 13 times slower than the trout mitochondrial regions that have been studied. The Lahontan cutthroat trout had a fixed OmyY1 sequence throughout ten separate populations, suggesting this subspecies underwent a severe population bottleneck prior to its current dispersal throughout the Great Basin during the pluvial phase of the last ice age. The Yellowstone group is the most derived among the cutthroat trout and consists of the Yellowstone, Bonneville, Colorado, Rio Grande and greenback subspecies. Identification of subspecies and sex with this Y-chromosome marker may prove useful in conservation efforts.

Genetic structuring in the freshwater mussel Anodonta corresponds with major hydrologic basins in the western United States

Freshwater mussels (unionids) are increasingly recognized as important providers of ecosystem services, yet are among the most endangered fauna in the world. Because unionids are generally sessile and require specific fish hosts for development and dispersal, they are particularly vulnerable to habitat degradation. Surprisingly, little is known about the distribution of genetic diversity in freshwater mussels and this gap has a negative impact on taxonomy, monitoring, conservation and ecological research in these species. Here, we focus on western North American Anodonta, one of only three genera known to exist in this broad landscape and which contains three highly divergent lineages. We describe phylogeographical subdivision in the most widespread and diverse of these lineages, which includes Anodonta californiensis and Anodonta nuttalliana and occurs from Canada to Mexico. Using mitochondrial and nuclear data, we found that genetic structuring within this clade is inconsistent with morphologically based species designations, but instead follows patterns of vicariance among major hydrogeologic basins. Furthermore, there was a strong tendency for population diversity within drainage systems to increase downstream, implying greater habitat or host fish availability in this direction. Microsatellite results indicated that sampling locations were all genetically distinct, even at short distances. Many of our sample populations showed evidence of a recent demographic bottleneck, although this effect seemed to be very local and not drainage or basin-specific. This study provides a foundation for the establishment of appropriate management units and future research on adaptive differentiation and host fish relationships.

Phylogenetic relationships of the western North American cyprinid genus Richardsonius, with an overview of phylogeographic structure

Diversification of many North American taxa, including freshwater fishes, has been heavily influenced by the effects of complex geological and climatic events throughout the Cenozoic that have significantly altered the landscape. Here, we employ an array of phylogenetic analyses using a multiple gene tree approach to address several questions regarding the phylogenetic relationships of the North American cyprinid genus Richardsonius and two other closely related genera, Clinostomus and Iotichthys. We also use divergence time estimates generated using fossil calibrations to qualitatively assess the phylogeographic implications of evolution within the group. Mitochondrial and nuclear DNA sequences show a sister relationship between Iotichthys and Richardsonius, with Clinostomus being sister to an Iotichthys-Richardsonius clade, hence the currently recognized sister relationship between Clinostomus and Richardsonius is not supported. These genera appear to be monophyletic lineages, and sister species within genera appear to be reciprocally monophyletic. The two species within the genus Richardsonius both exhibit phylogeographic structure that is worthy of further investigation. Divergence time estimates between genera and species are Miocene or Pliocene in age, and divergence between phylogroups within species occurred in the late Pliocene to Pleistocene. These splits coincide with documented geological and climatic events.

Morphological and genetic structuring in the Utah Lake sucker complex

Population decline in the federally endangered June sucker (Chasmistes liorus), a lakesucker unique to Utah Lake, Utah, has been attributed in part to hybridization with the more widespread Utah sucker (Catostomus ardens). As a group, suckers in Utah Lake exhibit considerable external morphological variation. Meristic and morphological ambiguities, presumably the result of hybridization, create a continuum of intermediate forms between Chasmistes and Catostomus extremes and prevent definitive identification to species. Here we describe and evaluate the morphological and genetic variation in suckers in Utah Lake by comparing a morphological analysis with amplified fragment length polymorphism and microsatellite analyses. Suckers were morphologically differentiated using mouth characters associated with different feeding strategies: planktivory (June sucker) and benthivory (Utah sucker). Although we found no genetic evidence for a deep divergence between June and Utah morphs, significant, but slight population structuring accompanied the substantial morphological variation. Bayesian model-based genetic clustering analyses detected two sucker populations in Utah Lake; however, these clusters were not strongly concordant with morphological groupings or between marker systems. The suckers in Utah Lake present an interesting dilemma regarding conservation: should one conserve (breed and stock) a subset of the morphotypic variation in the Utah Lake sucker complex, focusing on the endangered June sucker morphotype, or should one conserve both June sucker and Utah sucker morphotypes in this complex, possibly maximizing evolutionary potential? We explore this question in the context of current genetic and morphological variation in the Utah Lake sucker complex as well as historical information on this complex and other lakesuckers.

An introduced and a native vertebrate hybridize to form a genetic bridge to a second native species

The genetic impacts of hybridization between native and introduced species are of considerable conservation concern, while the possibility of reticulate evolution affects our basic understanding of how species arise and shapes how we use genetic data to understand evolutionary diversification. By using mitochondrial NADH dehydrogenase subunit 2 (ND2) sequences and 467 amplified fragment-length polymorphism nuclear DNA markers, we show that the introduced white sucker (Catostomus commersoni) has hybridized with two species native to the Colorado River Basin--the flannelmouth sucker (Catostomus latipinnis) and the bluehead sucker (Catostomus discobolus). Hybrids between the flannelmouth sucker and white sucker have facilitated introgression between the two native species, previously isolated by reproductive barriers, such that individuals exist with contributions from all three genomes. Most hybrids had the mitochondrial haplotype of the introduced white sucker, emphasizing its pivotal role in this three-way hybridization. Our findings highlight how introduced species can threaten the genetic integrity of not only one species but also multiple previously reproductively isolated species. Furthermore, this complex three-way reticulate (as opposed to strictly bifurcating) evolution suggests that seeking examples in other vertebrate systems might be productive. Although the present study involved an introduced species, similar patterns of hybridization could result from natural processes, including stream capture or geological formations (e.g., the Bering land bridge).