Watershed boundaries and geographic isolation: patterns of diversification in cutthroat trout from western North America

Population genetics of the African snakehead fish Parachanna obscura along West Africa's water networks: Implications for sustainable management and conservation

An essential factor for aquatic conservation is genetic diversity or population divergence, which in natural populations reflects the interplay between geographical isolation with restricted gene flow and local evolution of populations. The long geological history of Africa may induce stronger among-population divergence and lower within-population divergence in fish populations of African watersheds. As an example, we studied population structure of the African snakehead fish Parachanna obscura. Our study aimed: (1) to develop a set of highly polymorphic microsatellite markers suitable for the analysis of genetic diversity among P. obscura and (2) to study the genetic diversity and structure of P. obscura populations from the West Africa region and mostly from Côte d'Ivoire, with respect to the effects of climate region and geographical distance on the genetic differentiation. A total of 259 specimens from 15 locations of P. obscura were collected over the West Africa region reflecting a high variability of pairwise geographical distances and variability of hydrological connectivity of the area. We developed a set of 21 polymorphic microsatellite markers for studying population genetics of the fish. The results showed relatively low intragenetic diversity for all the 15 locations, certainly attributable to confinement of fish in segregated catchments, resulting in limited gene flow. We also found evidence for local adaptation processes, suggested by five out of 21 microsatellite loci being under putative selection, according to BAYESCAN analysis. In turn, there was strong genetic differentiation (F _ST > 0.5) among fish from most locations, reflecting the allopatric development in watersheds without hydraulic connectivity. Neighbor-joining dendrogram, Principal Coordinate Analysis, and analysis of ancestral groups by STRUCTURE suggested that the 15 locations can be divided into three clusters, mainly matching the dominant climate zones and the segregation of the watersheds in the region. The distinct genetic structure of the fish from the 15 locations obtained in this study suggests that conservation and sustainable management actions of this fish resource should avoid genetic mixing of potentially locally adapted populations.

Estimating the genetic diversity of Pacific salmon and trout using multigene eDNA metabarcoding

Genetic diversity underpins species conservation and management goals, and ultimately determines a species' ability to adapt. Using freshwater environmental DNA (eDNA) samples, we examined mitochondrial genetic diversity using multigene metabarcode sequence data from four Oncorhynchus species across 16 sites in Oregon and northern California. Our multigene metabarcode panel included targets commonly used in population genetic NADH dehydrogenase 2 (ND2), phylogenetic cytochrome c oxidase subunit 1 (COI) and eDNA (12S ribosomal DNA) screening. The ND2 locus showed the greatest within-species haplotype diversity for all species, followed by COI and then 12S rDNA for all species except Oncorhynchus kisutch. Sequences recovered for O. clarkii clarkii were either identical to, or one mutation different from, previously characterized haplotypes (95.3% and 4.5% of reads, respectively). The greatest diversity in O. c. clarkii was among coastal watersheds, and subsets of this diversity were shared with fish in inland watersheds. However, coastal streams and the Umpqua River watershed appear to harbour unique haplotypes. Sequences from O. mykiss revealed a disjunction between the Willamette watershed and southern watersheds suggesting divergent histories. We also identified similarities between populations in the northern Deschutes and southern Klamath watersheds, consistent with previously hypothesized connections between the two via inland basins. Oncorhynchus kisutch was only identified in coastal streams and the Klamath River watershed, with most diversity concentrated in the coastal Coquille watershed. Oncorhynchus tshawytscha was only observed at one site, but contained multiple haplotypes at each locus. The characterization of genetic diversity at multiple loci expands the knowledge gained from eDNA sampling and provides crucial information for conservation actions and genetic management.

Genetic differentiation among populations of the blackfin goodea Goodea atripinnis (Cyprinodontiformes: Goodeidae): implications for its evolutionary history

Central Mexico is characterized by a complex topography that is the result of historic and contemporary tectonic and climatic factors. These events have influenced the evolutionary history of numerous freshwater fishes in the region. Nonetheless, recent studies have shown that life-history traits and ecological characteristics of species may influence dispersal capabilities and the degree of genetic connectivity. Goodea (Cyprinodontiformes: Goodeidae) is one of the most widely distributed and environmentally tolerant genera of goodeids. In this study, the authors analysed variation in the mitochondrial cytochrome b gene to evaluate the phylogeographic relationships, genetic structure, genetic diversity and demographic history of Goodea from across its distribution range. They found low genetic differentiation and identified shared haplotypes among several regions. Geographic segregation was found in samples southwest and northeast of the Lower Lerma region, with some internal isolated groups showing phylogeographic differentiation and unique haplotypes. The AMOVA best explained genetic structure when grouped by haplogroups rather than when grouped by recognized biogeographic regions. Several regions showed null genetic diversity, raising the possibility of dispersal mediated by humans. Finally, Bayesian Skyline Plot analysis showed a population expansion for the Southwest haplogroup, except for the Armería population and sub-group II of the Northeast haplogroup. All this suggests a recent colonization of Goodea atripinnis throughout some of the biogeographic regions currently inhabited by this species.

Morphological and genetic concordance of cutthroat trout (Oncorhynchus clarkii) diversification from western North America

The cutthroat trout (Oncorhynchus clarkii (Richardson, 1836)) is one of the most widely distributed species of freshwater fish in western North America. Occupying a diverse range of habitats, they exhibit significant phenotypic variability that is often recognized by intraspecific taxonomy. Recent molecular phylogenies have described phylogenetic diversification across cutthroat trout populations, but no study has provided a range-wide morphological comparison of taxonomic divisions. In this study, we used linear- and geometric-based morphometrics to determine if phylogenetic and subspecies divisions correspond to morphological variation in cutthroat trout, using replicate populations from throughout the geographic range of the species. Our data indicate significant morphological divergence of intraspecific categories in some, but not all, cutthroat trout subspecies. We also compare morphological distance measures with distance measures of mtDNA sequence divergence. DNA sequence divergence was positively correlated with morphological distance measures, indicating that morphologically more similar subspecies have lower sequence divergence in comparison to morphologically distant subspecies. Given these results, integrating both approaches to describing intraspecific variation may be necessary for developing a comprehensive conservation plan in wide-ranging species.

Molecular phylogeny and biogeography of Triplophysa stone loaches in the Central Chinese Mountains

The geological structure and history of mountains often cause highly complex patterns of evolution. The origin and diversity of the genus Triplophysa is considered to be closely associated with the rapid and persistent rise of the Qinghai-Tibet Plateau (QTP). Species of the genus are widely distributed in the QTP and adjacent regions, and are strongly adapted to its extreme environments. However, the Central Chinese Mountains (CCM), which are far apart from the QTP, also have many Triplophysa populations and yet their taxonomy, phylogeny and origin are unknown. We used three mitochondrial genes of 266 individuals from 29 sites in the CCM to shed light on the evolution of these Triplophysa populations. Phylogenetic analyses and species delimitation revealed seven distinct molecular operational taxonomic units (mOTUs) in the CCM, each with low genetic diversity. Molecular clock and biogeographic analyses suggested that vicariance events resulting from the rapid uplifting of the Qinling and Taihang Mountains initialized diversification in the Pliocene. During the early Pleistocene, dispersal events in the central parts of the Yellow River drainage altered genetic diversity and led to the currently observed biogeographic patterns. Our results suggest the presence of at least six species of Triplophysa in the CCM and more are likely to be discovered in the future.

Population dynamics of threatened Lahontan cutthroat trout in Summit Lake, Nevada

Summit Lake, Nevada (USA) is the last high-desert terminal lake to have a native self-sustaining population of threatened Lahontan cutthroat trout (Oncorhynchus clarkii henshawi). From spring 2015 to fall 2017, we quantified adult abundance and survival and the total annual spawning run. Abundance and survival were estimated with mark-recapture using PIT tags, and the annual spawning run was estimated with PIT tag detections and counts of spawners. Adult abundance fluctuated from 830 (95% CI 559–1248) to 1085 (95% CI 747–1614), with no overall temporal trend, as a decrease in male abundance was generally offset by an equal increase in female abundance. Estimated mean adult survival was 0.51 (95% CI 0.44–0.58). The spawning run increased from 645 (2015) to 868 (2016), but then decreased slightly to 824 (2017, mean = 789 ± 118). Female spawners increased in 2016 but decreased slightly in 2017, whereas male spawners decreased each year. In addition, the proportion of adults that spawned each year increased overall. Our study suggests that the adult population remained stable although most of the study period included the recent, severe regional drought in the western United States (2012–2016).

Tempo and mode of allopatric divergence in the weakly electric fish Sternopygus dariensis in the Isthmus of Panama

Spatial isolation is one of the main drivers of allopatric speciation, but the extent to which spatially-segregated populations accumulate genetic differences relevant to speciation is not always clear. We used data from ultraconserved elements (UCEs) and whole mitochondrial genomes (i.e., mitogenomes) to explore genetic variation among allopatric populations of the weakly electric fish Sternopygus dariensis across the Isthmus of Panama. We found strong genetic divergence between eastern and western populations of S. dariensis. Over 77% of the UCE loci examined were differentially fixed between populations, and these loci appear to be distributed across the species' genome. Population divergence occurred within the last 1.1 million years, perhaps due to global glaciation oscillations during the Pleistocene. Our results are consistent with a pattern of genetic differentiation under strict geographic isolation, and suggest the presence of incipient allopatric species within S. dariensis. Genetic divergence in S. dariensis likely occurred in situ, long after the closure of the Isthmus of Panama. Our study highlights the contribution of spatial isolation and vicariance to promoting rapid diversification in Neotropical freshwater fishes. The study of spatially-segregated populations within the Isthmus of Panama could reveal how genetic differences accumulate as allopatric speciation proceeds.

Connections and containers: Using genetic data to understand how watershed evolution and human activities influence cutthroat trout biogeography

Species with large geographic distributions often exhibit complex patterns of diversity that can be further complicated by human activities. Cutthroat trout (Oncorhynchus clarkii) are one of the most widely distributed freshwater fish species in western North America exhibiting substantial phenotypic and genetic variability; however, fish stocking practices have translocated populations outside of their native range and may have obscured intraspecific boundaries. This study focuses on cutthroat trout populations representing three distinct evolutionary clades that are found intermixed within a contact zone between the Bonneville and upper Snake River watersheds in the western United States. We used mitochondrial and microsatellite genetic data, as well as historical stocking records, to evaluate whether populations of cutthroat trout in the contact zone are native or are introduced. We found significant genetic differentiation and fine-scale genetic population structure that was organized primarily by watershed boundaries. While we detected increased genetic diversity in some areas in close proximity to the greatest number of stocking events, the highly organized population structure both within and between areas of the contact zone indicates that the populations are native to the watersheds. Intermixing of distinct evolutionary lineages of cutthroat trout appears to be the result of historical connections between paleodrainages. Our analyses provide a context for understanding how genetic data can be used to assess the status of populations as native or introduced.

Patterns of hybridization among cutthroat trout and rainbow trout in northern Rocky Mountain streams

Introgressive hybridization between native and introduced species is a growing conservation concern. For native cutthroat trout and introduced rainbow trout in western North America, this process is thought to lead to the formation of hybrid swarms and the loss of monophyletic evolutionary lineages. Previous studies of this phenomenon, however, indicated that hybrid swarms were rare except when native and introduced forms of cutthroat trout co‐occurred. We used a panel of 86 diagnostic, single nucleotide polymorphisms to evaluate the genetic composition of 3865 fish captured in 188 locations on 129 streams distributed across western Montana and northern Idaho. Although introgression was common and only 37% of the sites were occupied solely by parental westslope cutthroat trout, levels of hybridization were generally low. Of the 188 sites sampled, 73% contained ≤5% rainbow trout alleles and 58% had ≤1% rainbow trout alleles. Overall, 72% of specimens were nonadmixed westslope cutthroat trout, and an additional 3.5% were nonadmixed rainbow trout. Samples from seven sites met our criteria for hybrid swarms, that is, an absence of nonadmixed individuals and a random distribution of alleles within the sample; most (6/7) were associated with introgression by Yellowstone cutthroat trout. In streams with multiple sites, upstream locations exhibited less introgression than downstream locations. We conclude that although the widespread introduction of nonnative trout within the historical range of westslope cutthroat trout has increased the incidence of introgression, sites containing nonadmixed populations of this taxon are common and broadly distributed.

Reserve network planning for fishes in the middle and lower Yangtze River basin by systematic conservation approaches

Although China has established more than 600 wetland nature reserves, conservation gaps still exist for many species, especially for freshwater fishes. Underlying this problem is the fact that top-level planning is missing in the construction of nature reserves. To promote the development of nature reserves for fishes, this study took the middle and lower reaches of the Yangtze River basin (MLYRB) as an example to carry out top-level reserve network planning for fishes using approaches of systematic conservation planning. Typical fish species living in freshwater habitats were defined and considered in the planning. Based on sample data collected from large quantities of literatures, continuous distribution patterns of 142 fishes were obtained with species distribution modeling and subsequent processing, and the distributions of another eleven species were artificially designated. With the distribution pattern of species, Marxan was used to carry out conservation planning. To obtain ideal solutions with representativeness, persistence, and efficiency, parameters were set with careful consideration regarding existing wetland reserves, human disturbances, hydrological connectivity, and representation targets of species. Marxan produced the selection frequency of planning units (PUs) and a best solution. Selection frequency indicates the relative protection importance of a PU. The best solution is a representative of ideal fish reserve networks. Both of the PUs with high selection frequency and those in the best solution have low proportions included in existing wetland nature reserves, suggesting that there are significant conservation gaps for fish species in MLYRB. The best solution could serve as a reference for establishing a fish reserve network in the MLYRB. There is great flexibility for replacing selected PUs in the solution, and such flexibility facilitates the implementation of the solution in reality in case of unexpected obstacles. Further, we suggested adopting a freshwater management framework in the implementation of such solution.

Climate oscillations, glacial refugia, and dispersal ability: factors influencing the genetic structure of the least salmonfly, Pteronarcella badia (Plecoptera), in Western North America

Background

Phylogeographic studies of aquatic insects provide valuable insights into mechanisms that shape the genetic structure of communities, yet studies that include broad geographic areas are uncommon for this group. We conducted a broad scale phylogeographic analysis of the least salmonfly Pteronarcella badia (Plecoptera) across western North America. We tested hypotheses related to mode of dispersal and the influence of historic climate oscillations on population genetic structure. In order to generate a larger mitochondrial data set, we used 454 sequencing to reconstruct the complete mitochondrial genome in the early stages of the project.

Results

Our analysis revealed high levels of population structure with several deeply divergent clades present across the sample area. Evidence from five mitochondrial genes and one nuclear locus identified a potentially cryptic lineage in the Pacific Northwest. Gene flow estimates and geographic clade distributions suggest that overland flight during the winged adult stage is an important dispersal mechanism for this taxon. We found evidence of multiple glacial refugia across the species distribution and signs of secondary contact within and among major clades.

Conclusions

This study provides a basis for future studies of aquatic insect phylogeography at the inter-basin scale in western North America. Our findings add to an understanding of the role of historical climate isolations in shaping assemblages of aquatic insects in this region. We identified several geographic areas that may have historical importance for other aquatic organisms with similar distributions and dispersal strategies as P. badia. This work adds to the ever-growing list of studies that highlight the potential of next-generation DNA sequencing in a phylogenetic context to improve molecular data sets from understudied groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0553-4) contains supplementary material, which is available to authorized users.

The Dual Challenges of Generality and Specificity When Developing Environmental DNA Markers for Species and Subspecies of Oncorhynchus

Environmental DNA (eDNA) sampling is a powerful tool for detecting invasive and native aquatic species. Often, species of conservation interest co-occur with other, closely related taxa. Here, we developed qPCR (quantitative PCR) markers which distinguish westslope cutthroat trout (Oncorhynchus clarkii lewsi), Yellowstone cutthroat trout (O. clarkii bouvieri), and rainbow trout (O. mykiss), which are of conservation interest both as native species and as invasive species across each other’s native ranges. We found that local polymorphisms within westslope cutthroat trout and rainbow trout posed a challenge to designing assays that are generally applicable across the range of these widely-distributed species. Further, poorly-resolved taxonomies of Yellowstone cutthroat trout and Bonneville cutthroat trout (O. c. utah) prevented design of an assay that distinguishes these recognized taxa. The issues of intraspecific polymorphism and unresolved taxonomy for eDNA assay design addressed in this study are likely to be general problems for closely-related taxa. Prior to field application, we recommend that future studies sample populations and test assays more broadly than has been typical of published eDNA assays to date.

Fragmentation and thermal risks from climate change interact to affect persistence of native trout in the Colorado River basin

Impending changes in climate will interact with other stressors to threaten aquatic ecosystems and their biota. Native Colorado River cutthroat trout (CRCT; Oncorhynchus clarkii pleuriticus) are now relegated to 309 isolated high-elevation (>1700 m) headwater stream fragments in the Upper Colorado River Basin, owing to past nonnative trout invasions and habitat loss. Predicted changes in climate (i.e., temperature and precipitation) and resulting changes in stochastic physical disturbances (i.e., wildfire, debris flow, and channel drying and freezing) could further threaten the remaining CRCT populations. We developed an empirical model to predict stream temperatures at the fragment scale from downscaled climate projections along with geomorphic and landscape variables. We coupled these spatially explicit predictions of stream temperature with a Bayesian Network (BN) model that integrates stochastic risks from fragmentation to project persistence of CRCT populations across the upper Colorado River basin to 2040 and 2080. Overall, none of the populations are at risk from acute mortality resulting from high temperatures during the warmest summer period. In contrast, only 37% of populations have a ≥90% chance of persistence for 70 years (similar to the typical benchmark for conservation), primarily owing to fragmentation. Populations in short stream fragments <7 km long, and those at the lowest elevations, are at the highest risk of extirpation. Therefore, interactions of stochastic disturbances with fragmentation are projected to be greater threats than warming for CRCT populations. The reason for this paradox is that past nonnative trout invasions and habitat loss have restricted most CRCT populations to high-elevation stream fragments that are buffered from the potential consequences of warming, but at risk of extirpation from stochastic events. The greatest conservation need is for management to increase fragment lengths to forestall these risks.

Single nucleotide polymorphism discovery in cutthroat trout subspecies using genome reduction, barcoding, and 454 pyro-sequencing

Background

Salmonids are popular sport fishes, and as such have been subjected to widespread stocking throughout western North America. Historically, stocking was done with little regard for genetic variation among populations and has resulted in genetic mixing among species and subspecies in many areas, thus putting the genetic integrity of native salmonid populations at risk and creating a need to assess the genetic constitution of native salmonid populations. Cutthroat trout is a salmonid species with pronounced geographic structure (there are 10 extant subspecies) and a recent history of hybridization with introduced rainbow trout in many populations. Genetic admixture has also occurred among cutthroat trout subspecies in areas where introductions have brought two or more subspecies into contact. Consequently, management agencies have increased their efforts to evaluate the genetic composition of cutthroat trout populations to identify populations that remain uncompromised and manage them accordingly, but additional genetic markers are needed to do so effectively. Here we used genome reduction, MID-barcoding, and 454-pyrosequencing to discover single nucleotide polymorphisms that differentiate cutthroat trout subspecies and can be used as a rapid, cost-effective method to characterize the genetic composition of cutthroat trout populations.

Results

Thirty cutthroat and six rainbow trout individuals were subjected to genome reduction and next-generation sequencing. A total of 1,499,670 reads averaging 379 base pairs in length were generated by 454-pyrosequencing, resulting in 569,060,077 total base pairs sequenced. A total of 43,558 putative SNPs were identified, and of those, 125 SNP primers were developed that successfully amplified 96 cutthroat trout and rainbow trout individuals. These SNP loci were able to differentiate most cutthroat trout subspecies using distance methods and Structure analyses.

Conclusions

Genomic and bioinformatic protocols were successfully implemented to identify 125 nuclear SNPs that are capable of differentiating most subspecies of cutthroat trout from one another. The ability to use this suite of SNPs to identify individuals of unknown genetic background to subspecies can be a valuable tool for management agencies in their efforts to evaluate the genetic structure of cutthroat trout populations prior to constructing and implementing conservation plans.