The architecture of river networks can drive the evolutionary dynamics of aquatic populations

Evolutionary history of the polymorphic Coatzacoalcos River endemic Thorichthys panchovillai (Cichliformes: Cichlidae)

Cichlids of the genus Thorichthys are a morphologically diverse clade of nine species occurring from Mexico south to the Motagua River in Guatemala and Honduras. Our understanding of species relationships within Thorichthys and other genera of Northern Neotropical cichlids has improved in recent years; however, phylogenetic placement of some species, as well as population-level variation, remains understudied. Thorichthys panchovillai is a polymorphic species endemic to the Coatzacoalcos River in the Atlantic slope of Mexico. The species was described from the upper reaches of this river, and since that time phenotypic variation in coloration and body shape has been observed. However, the species has never been included in a molecular phylogenetic study to understand its evolutionary relationships within the genus and better relate the observed morphological polymorphism to genetic patterns within this species. In this study we use mitochondrial data to study the phylogenetic placement of T. panchovillai in the context of all other species of Thorichthys. Results show that spatially segregated haplotypes correspond to two main morphotypes of T. panchovillai and allow for the proposal of hypotheses regarding maintenance of this polymorphism in the Coatzacoalcos River.

Balancing Inbreeding and Outbreeding Risks to Inform Translocations Throughout the Range of an Imperiled Darter

Restoring connectivity via assisted migration is a useful but currently underused approach for maintaining genetic diversity and preventing extirpations of threatened species. The use of assisted migration as a conservation strategy may be limited by the difficulty of balancing the benefits of reconnecting populations (including reduced inbreeding depression and increased adaptive capacity) with the perceived risk of outbreeding depression, which requires comprehensive knowledge of the landscape of adaptive, neutral, deleterious, and structural variation across a species' range. Using a combination of reduced-representation and whole-genome sequencing, we characterized genomic diversity and differentiation for the Arkansas Darter (Etheostoma cragini) across its range in the Midwestern US. We found strong population structure and large differences in genetic diversity and effective population sizes across drainages. The strength of genetic isolation by river distance differed among drainages, with landscape type surrounding streams and impoundments also contributing to genetic isolation. Despite low effective population sizes in some populations, there was surprisingly little evidence for recent inbreeding (based on the absence of long runs of homozygosity) or for elevated levels of deleterious variation in smaller populations. Considering neutral, adaptive, deleterious, and structural variation allowed us to identify several potential recipient populations that may benefit from translocations and potential donor sites throughout the range. Planning translocation strategies intended for restored connectivity and possible genetic rescue at earlier stages in species decline will likely increase the probability of retaining genetic diversity and population persistence over the long term while minimizing risks associated with translocation.

Spatial population genetic structure of Caquetaia kraussii (Steindachner, 1878) evidenced by species-specific microsatellite loci in the middle and low basin of the Cauca River, Colombia

The adaptative responses and divergent evolution shown in the environments habited by the Cichlidae family allow to understand different biological properties, including fish genetic diversity and structure studies. In a zone that has been historically submitted to different anthropogenic pressures, this study assessed the genetic diversity and population structure of cichlid Caquetaia kraussii, a sedentary species with parental care that has a significant ecological role for its contribution to redistribution and maintenance of sedimentologic processes in its distribution area. This study developed de novo 16 highly polymorphic species-specific microsatellite loci that allowed the estimation of the genetic diversity and differentiation in 319 individuals from natural populations in the area influenced by the Ituango hydroelectric project in the Colombian Cauca River. Caquetaia kraussii exhibits high genetic diversity levels (Ho: 0.562-0.885; He: 0.583-0.884) in relation to the average neotropical cichlids and a three group-spatial structure: two natural groups upstream and downstream the Nechí River mouth, and one group of individuals with high relatedness degree, possibly independently formed by founder effect in the dam zone. The three genetic groups show recent bottlenecks, but only the two natural groups have effective population size that suggest their long-term permanence. The information generated is relevant not only for management programs and species conservation purposes, but also for broadening the available knowledge on the factors influencing neotropical cichlids population genetics.

Riverscape community genomics: A comparative analytical approach to identify common drivers of spatial structure

Genetic differentiation among local groups of individuals, that is, genetic β-diversity, is a key component of population persistence related to connectivity and isolation. However, most genetic investigations of natural populations focus on a single species, overlooking opportunities for multispecies conservation plans to benefit entire communities in an ecosystem. We present an approach to evaluate genetic β-diversity within and among many species and demonstrate how this riverscape community genomics approach can be applied to identify common drivers of genetic structure. Our study evaluated genetic β-diversity in 31 co-distributed native stream fishes sampled from 75 sites across the White River Basin (Ozarks, USA) using SNP genotyping (ddRAD). Despite variance among species in the degree of genetic divergence, general spatial patterns were identified corresponding to river network architecture. Most species (N = 24) were partitioned into discrete subpopulations (K = 2-7). We used partial redundancy analysis to compare species-specific genetic β-diversity across four models of genetic structure: Isolation by distance (IBD), isolation by barrier (IBB), isolation by stream hierarchy (IBH), and isolation by environment (IBE). A significant proportion of intraspecific genetic variation was explained by IBH (x̄ = 62%), with the remaining models generally redundant. We found evidence for consistent spatial modularity in that gene flow is higher within rather than between hierarchical units (i.e., catchments, watersheds, basins), supporting the generalization of the stream hierarchy model. We discuss our conclusions regarding conservation and management and identify the 8-digit hydrologic unit (HUC) as the most relevant spatial scale for managing genetic diversity across riverine networks.

Contrasts in riverscape patterns of intraspecific genetic variation in a diverse Neotropical fish community of high conservation value

Spatial patterns of genetic variation compared across species provide information about the predictability of genetic diversity in natural populations, and areas requiring conservation measures. Due to their remarkable fish diversity, rivers in Neotropical regions are ideal systems to confront theory with observations and would benefit greatly from such approaches given their increasing vulnerability to anthropogenic pressures. We used SNP data from 18 fish species with contrasting life-history traits, co-sampled across 12 sites in the Maroni- a major river system from the Guiana Shield -, to compare patterns of intraspecific genetic variation and identify their underlying drivers. Analyses of covariance revealed a decrease in genetic diversity as distance from the river outlet increased for 5 of the 18 species, illustrating a pattern commonly observed in riverscapes for species with low-to-medium dispersal abilities. However, the mean within-site genetic diversity was lowest in the two easternmost tributaries of the Upper Maroni and around an urbanized location downstream, indicating the need to address the potential influence of local pressures in these areas, such as gold mining or fishing. Finally, the relative influence of isolation by stream distance, isolation by discontinuous river flow, and isolation by spatial heterogeneity in effective size on pairwise genetic differentiation varied across species. Species with similar dispersal and reproductive guilds did not necessarily display shared patterns of population structure. Increasing the knowledge of specific life history traits and ecological requirements of fish species in these remote areas should help further understand factors that influence their current patterns of genetic variation.

Environmental selection, rather than neutral processes, best explain regional patterns of diversity in a tropical rainforest fish

To conserve the high functional and genetic variation in hotspots such as tropical rainforests, it is essential to understand the forces driving and maintaining biodiversity. We asked to what extent environmental gradients and terrain structure affect morphological and genomic variation across the wet tropical distribution of an Australian rainbowfish, Melanotaenia splendida splendida. We used an integrative riverscape genomics and morphometrics framework to assess the influence of these factors on both putative adaptive and non-adaptive spatial divergence. We found that neutral genetic population structure was largely explainable by restricted gene flow among drainages. However, environmental associations revealed that ecological variables had a similar power to explain overall genetic variation, and greater power to explain body shape variation, than the included neutral covariables. Hydrological and thermal variables were the strongest environmental predictors and were correlated with traits previously linked to heritable habitat-associated dimorphism in rainbowfishes. In addition, climate-associated genetic variation was significantly associated with morphology, supporting heritability of shape variation. These results support the inference of evolved functional differences among localities, and the importance of hydroclimate in early stages of diversification. We expect that substantial evolutionary responses will be required in tropical rainforest endemics to mitigate local fitness losses due to changing climates.

Molecular Markers Reveal a New and Possibly Threatened Species of Cnesterodon (Poeciliidae, Cnesterodontini) from the Upper Paraná River Basin, Brazil

Cnesterodon comprises 10 valid species occurring in the major river basins of South America. Recent ichthyofaunistic studies in the Ivaí River basin, upper Paraná River system, suggested the existence of a possible new species, which was identified as Cnesterodon sp. based on morphological characters. Currently, the use of molecular tools has proved to be fundamental in aiding phylogenetics and cataloging biodiversity; therefore, in this study, we molecularly characterize a possible new species of Cnesterodon from the Ivaí River basin encoding the mitochondrial genes Cytochrome c Oxidase, subunit I (COI), and NADH dehydrogenase subunit 2 (ND2). The genetic differences found showed that this species really differs from the other Cnesterodon species, indicating that it is a distinct species, which is possibly already in serious danger of extinction since its habitat often suffers from human exploitation and its distribution is restricted to only two sites in the upper Ivaí River basin, but it has disappeared in one of them.

Contrasting population genetic responses to migration barriers in two native and an invasive freshwater fish

Habitat fragmentation impacts the distribution of genetic diversity and population genetic structure. Therefore, protecting the evolutionary potential of species, especially in the context of the current rate of human-induced environmental change, is an important goal. In riverine ecosystems, migration barriers affect the genetic structure of native species, while also influencing the spread of invasive species. In this study, we compare genetic patterns of two native and one highly invasive riverine fish species in a Belgian river basin, namely the native three-spined stickleback (Gasterosteus aculeatus) and stone loach (Barbatula barbatula), and the non-native and invasive topmouth gudgeon (Pseudorasbora parva). We aimed to characterize both natural and anthropogenic determinants of genetic diversity and population genetic connectivity. Genetic diversity was highest in topmouth gudgeon, followed by stone loach and three-spined stickleback. The correlation between downstream distance and genetic diversity, a pattern often observed in riverine systems, was only marginally significant in stone loach and three-spined stickleback, while genetic diversity strongly declined with increasing number of barriers in topmouth gudgeon. An Isolation-By-Distance pattern characterizes the population genetic structure of each species. Population differentiation was only associated with migration barriers in the invasive topmouth gudgeon, while genetic composition of all species seemed at least partially determined by the presence of migration barriers. Among the six barrier types considered (watermills, sluices, tunnels, weirs, riverbed obstructions, and others), the presence of watermills was the strongest driver of genetic structure and composition. Our results indicate that conservation and restoration actions, focusing on conserving genetic patterns, cannot be generalized across species. Moreover, measures might target either on restoring connectivity, while risking a rapid spread of the invasive topmouth gudgeon, or not restoring connectivity, while risking native species extinction in upstream populations.

Genetic divergence among threespine stickleback that differ in nuptial coloration

Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.

Landscape connectivity among coastal giant salamander (Dicamptodon tenebrosus) populations shows no association with land use, fire frequency, or river drainage but exhibits genetic signatures of potential conservation concern

Determining the genetic consequences of both historical and contemporary events can clarify the effects of the environment on population connectivity and inform conservation decisions. Historical events (like glaciations) and contemporary factors (like logging) can disrupt gene flow between populations. This is especially true among species with specialized ecological requirements and low dispersal ability, like amphibians. We test for the genetic consequences of historical and contemporary disturbances in the coastal giant salamander (Dicamptodon tenebrosus) in the Pacific Northwest of the United States. We consider predictions based on the contemporary landscape (habitat connectivity, logging, forest fires, and topography), in addition to relatively ancient post-Pleistocene range expansion (following the last glacial retreat). To assess local versus larger-scale effects, we sampled 318 individuals across 23 sites, which were clustered in five sampling regions. Genetic variation was assessed using five microsatellite markers. We found evidence of (i) historical regional isolation, with decreased genetic diversity among more recently colonized northern sites, as well as (ii) high levels of inbreeding and loss of heterozygosity at local scales, despite relatively low overall population differentiation (FST) or strong evidence for population bottlenecks. Genetic diversity was not associated with contemporary disturbances (logging or fire), and there were no detectable effects on the genetic connectivity of populations based on intervening landscape features (habitat fragmentation and topography). However, lower genetic diversity in more northern regions indicates a lag in recovery of genetic diversity following post-Pleistocene expansion. Additionally, some populations had evidence of having undergone a recent genetic bottleneck or had high inbreeding (FIS) values. Lower genetic diversity in more northern sites means populations may be more vulnerable to future environmental changes, and managing for connectivity alone may not be sufficient given low mobility. Recent apparent reductions in some populations were not clearly linked to anthropogenic disturbances we examined. This suggests the type of disturbances this species is sensitive to may not be well understood.

Riverscape genomics of cichlid fishes in the lower Congo: Uncovering mechanisms of diversification in an extreme hydrological regime

Freshwater fishes are notably diverse, given that freshwater habitat represents a tiny fraction of the earth's surface, but the mechanisms generating this diversity remain poorly understood. Rivers provide excellent models to understand how freshwater diversity is generated and maintained across heterogeneous habitats. In particular, the lower Congo River (LCR) consists of a dynamic hydroscape exhibiting extraordinary aquatic biodiversity, endemicity, morphological and ecological specialization. Previous studies have suggested that the numerous high-energy rapids throughout the LCR form physical barriers to gene flow, thus facilitating diversification and speciation, generating ichthyofaunal diversity. However, this hypothesis has not been fully explored using genome-wide SNPs for fish species distributed across the LCR. Here, we examined four lamprologine cichlids endemic to the LCR that are distributed along the river without range overlap. Using genome-wide SNP data, we tested the hypotheses that high-energy rapids serve as physical barriers to gene flow that generate genetic divergence at inter- and intraspecific levels, and that gene flow occurs primarily in a downstream direction. Our results are consistent with the prediction that powerful rapids sometimes act as a barrier to gene flow but also suggest that, at certain temporal and spatial scales, they may provide multidirectional dispersal opportunities for riverine rheophilic cichlid fishes. These results highlight the complexity of diversification processes in rivers and the importance of assessing such processes across different riverscapes.

Landscape Evolution as a Diversification Driver in Freshwater Fishes

The exceptional concentration of vertebrate diversity in continental freshwaters has been termed the “freshwater fish paradox,” with > 15,000 fish species representing more than 20% of all vertebrate species compressed into tiny fractions of the Earth’s land surface area (<0.5%) or total aquatic habitat volume (<0.001%). This study asks if the fish species richness of the world’s river basins is explainable in terms of river captures using topographic metrics as proxies. The River Capture Hypothesis posits that drainage-network rearrangements have accelerated biotic diversification through their combined effects on dispersal, speciation, and extinction. Yet rates of river capture are poorly constrained at the basin scale worldwide. Here we assess correlations between fish species density (data for 14,953 obligate freshwater fish species) and basin-wide metrics of landscape evolution (data for 3,119 river basins), including: topography (elevation, average relief, slope, drainage area) and climate (average rainfall and air temperature). We assess the results in the context of both static landscapes (e.g., species-area and habitat heterogeneity relationships) and transient landscapes (e.g., river capture, tectonic activity, landscape disequilibrium). We also relax assumptions of functional neutrality of basins (tropical vs. extratropical, tectonically stable vs. active terrains). We found a disproportionate number of freshwater species in large, lowland river basins of tropical South America, Africa, and Southeast Asia, under predictable conditions of large geographic area, tropical climate, low topographic relief, and high habitat volume (i.e., high rainfall rates). However, our results show that these conditions are only necessary, but not fully sufficient, to explain the basins with the highest diversity. Basins with highest diversity are all located on tectonically stable regions, places where river capture is predicted to be most conducive to the formation of high fish species richness over evolutionary timescales. Our results are consistent with predictions of several landscape evolution models, including the River Capture Hypothesis, Mega Capture Hypothesis, and Intermediate Capture Rate Hypothesis, and support conclusions of numerical modeling studies indicating landscape transience as a mechanistic driver of net diversification in riverine and riparian organisms with widespread continental distributions.

Dispersal behaviour and riverine network connectivity shape the genetic diversity of freshwater amphipod metapopulations

Theory predicts that the distribution of genetic diversity in a landscape is strongly dependent on the connectivity of the metapopulation and the dispersal of individuals between patches. However, the influence of explicit spatial configurations such as dendritic landscapes on the genetic diversity of metapopulations is still understudied, and theoretical corroborations of empirical patterns are largely lacking. Here, we used microsatellite data and stochastic simulations of two metapopulations of freshwater amphipods in a 28,000 km² riverine network to study the influence of spatial connectivity and dispersal strategies on the spatial distribution of their genetic diversity. We found a significant imprint of the effects of riverine network connectivity on the local and global genetic diversity of both amphipod species. Data from 95 sites showed that allelic richness significantly increased towards more central nodes of the network. This was also seen for observed heterozygosity, yet not for expected heterozygosity. Genetic differentiation increased with instream distance. In simulation models, depending on the mutational model assumed, upstream movement probability and dispersal rate, respectively, emerged as key factors explaining the empirically observed distribution of local genetic diversity and genetic differentiation. Surprisingly, the role of site‐specific carrying capacities, for example by assuming a direct dependency of population size on local river size, was less clear cut: while our best fitting model scenario included this feature, over all simulations, scaling of carrying capacities did not increase data‐model fit. This highlights the importance of dispersal behaviour along spatial networks in shaping population genetic diversity.

A strategic sampling design revealed the local genetic structure of cold-water fluvial sculpin: a focus on groundwater-dependent water temperature heterogeneity

A key piece of information for ecosystem management is the relationship between the environment and population genetic structure. However, it is difficult to clearly quantify the effects of environmental factors on genetic differentiation because of spatial autocorrelation and analytical problems. In this study, we focused on stream ecosystems and the environmental heterogeneity caused by groundwater and constructed a sampling design in which geographic distance and environmental differences are not correlated. Using multiplexed ISSR genotyping by sequencing (MIG-seq) method, a fine-scale population genetics study was conducted in fluvial sculpin Cottus nozawae, for which summer water temperature is the determinant factor in distribution and survival. There was a clear genetic structure in the watershed. Although a significant isolation-by-distance pattern was detected in the watershed, there was no association between genetic differentiation and water temperature. Instead, asymmetric gene flow from relatively low-temperature streams to high-temperature streams was detected, indicating the importance of low-temperature streams and continuous habitats. The groundwater-focused sampling strategy yielded insightful results for conservation.

Molecular phylogeny for the Neotropical freshwater stingrays (Myliobatiformes: Potamotrygoninae) reveals limitations of traditional taxonomy

The subfamily Potamotrygoninae, the only extant clade of elasmobranchs exclusive to freshwater environments, encompasses four genera and 38 species distributed across almost every major South American river basin. Despite their importance in the ornamental fish trade, the taxonomy and evolutionary relationships within potamotrygonines have not yet been resolved. Here, we present a comprehensive molecular phylogeny for the Neotropical freshwater stingrays, based on extensive species and population sampling (35 species and > 350 individuals from drainages across South America). Our phylogeny corroborates the monophyly of the genera Paratrygon and Heliotrygon and the monophyly of the Potamotrygon + Plesiotrygon clade. Within the Potamotrygon + Plesiotrygon clade, we identify a core Potamotrygon clade characterized by short branches, low nodal support and incongruence with current species-level taxonomy. In the core Potamotrygon clade, specimens of widespread species, such as Potamotrygon motoro and Potamotrygon orbignyi, do not form monophyletic lineages; instead, specimens from these species are often closely related to those of other species from the same river basins. These patterns could be caused by inaccurate taxonomy, hybridization, incomplete lineage sorting and rapid diversification. We discuss the conservation of Neotropical freshwater stingrays from a phylogenetic perspective and suggest ways to prioritize potamotrygonid conservation efforts with respect to endemism and evolutionary distinctiveness.

Hidden or unnoticed? Multiple lines of evidence support the recognition of a new species of Pseudocorynopoma (Characidae: Corynopomini)

Species delimitation is a permanent issue in systematics. The increasing recognition of geographically isolated populations as independent lineages allowed by new methods of analysis has inflated the species-populations dilemma, which involves deciding whether to consider separate lineages as different species or structured genetic populations. This is commonly observed between fishes of adjacent river basins, with some lineages being considered allopatric sister species and others considered isolated populations or variants of the same species. Pseudocorynopoma doriae is a characid diagnosed from its single congener by the number of anal-fin rays and sexually dimorphic characters of males, including distinct fin colouration. The authors found variation in the colour pattern between isolated populations previously identified as P. doriae but no variation in scale or fin-ray counts. They analysed molecular evidence at the population level and morphological differences related to life history (e.g., colour dimorphism related to inseminating behaviour). The results provide compelling evidence for the recognition of a new species of Pseudocorynopoma despite the lack of discrete differences in meristic data. The recognition of the new species is consistent with biogeographical evidence for the long-term isolation of the respective river drainages and with differences between the ichthyofaunal communities of these rivers.

Recent and rapid anthropogenic habitat fragmentation increases extinction risk for freshwater biodiversity

Anthropogenic habitat fragmentation is often implicated as driving the current global extinction crisis, particularly in freshwater ecosystems. The genetic signal of recent population isolation can be confounded by the complex spatial arrangement of dendritic river systems. Consequently, many populations may presently be managed separately based on an incorrect assumption that they have evolved in isolation. Integrating landscape genomics data with models of connectivity that account for landscape structure, we show that the cumulative effects of multiple in-stream barriers have contributed to the recent decline of a freshwater fish from the Murray-Darling Basin, Australia. In addition, individual-based eco-evolutionary simulations further demonstrate that contemporary inferences about population isolation are consistent with the 160-year time frame since construction of in-stream barriers began in the region. Our findings suggest that the impact of very recent fragmentation may be often underestimated for freshwater biodiversity. We argue that proactive conservation measures to reconnect many riverine populations are urgently needed.

Channel network structure determines genetic connectivity of landward-seaward Avicennia marina populations in a tropical bay

Mangrove ecosystems along the East African coast are often characterized by a disjunct zonation pattern of seaward and landward Avicennia marina trees. This disjunct zonation may be maintained through different positions in the tidal frame, yielding different dispersal settings. The spatial configuration of the landscape and coastal processes such as tides and waves is expected to largely influence the extent of propagule transport and subsequent regeneration. We hypothesized that landward sites would keep a stronger genetic structure over short distances in comparison with enhanced gene flow among regularly flooded seaward fringes. We tested this hypothesis from densely vegetated A. marina transects of a well-documented mangrove system (Gazi Bay, Kenya) and estimated local gene flow and kinship-based fine-scale genetic structure. Ten polymorphic microsatellite markers in 457 A. marina trees revealed no overall significant difference in levels of allele or gene diversities between sites that differ in hydrological proximity. Genetic structure and connectivity of A. marina populations however indicated an overall effect of geographic distance and revealed a pronounced distinction between channels and topographic setting. Migration models allowed to infer gene flow directionality among channels, and indicated a bidirectional steppingstone between seaward and nearest located landward stands. Admixed gene pools without any fine-scale structure were found within the wider and more exposed Kidogoweni channel, suggesting open systems. Elevated kinship values and structure over 5 to 20 m distance were only detected in two distant landward and seaward transects near the mouth of the Mkurumuji River, indicating local retention and establishment. Overall, our findings show that patterns of A. marina connectivity are explained by hydrological proximity, channel network structure, and hydrokinetic energy, rather than just their positioning as disjunct landward or seaward zones.

River capture or ancestral polymorphism: an empirical genetic test in a freshwater fish using approximate Bayesian computation

A headwater or river capture is a phenomenon commonly invoked to explain the absence of reciprocal monophyly of genetic lineages among isolated hydrographic basins in freshwater fish. Under the assumption of river capture, a secondary contact between populations previously isolated in different basins explains the observed genetic pattern. However, the absence of reciprocal monophyly could also arise under population isolation through the retention of ancestral of polymorphisms. Here, we applied an approximate Bayesian computation (ABC) framework for estimating the relative probability of scenarios with and without secondary contact. We used Cnesterodon decemmaculatus as a study model because of the multiple possible cases of river capture and the demographic parameters estimated in a previous mitochondrial DNA study that are useful for simulating scenarios to test both hypotheses using the ABC framework. Our results showed that, in general, mitochondrial DNA is useful for distinguishing between these alternative demographic scenarios with reasonable confidence, but in extreme cases (e.g. recent divergence or large population size) there is no power to discriminate between scenarios. Testing hypotheses of drainage rearrangement under a statistically rigorous framework is fundamental for understanding the evolution of freshwater fish fauna as a complement to, or in the absence of, geological evidence.

Riverscape properties contribute to the origin and structure of a hybrid zone in a Neotropical freshwater fish

Understanding the structure of hybrid zones provides valuable insights about species boundaries and speciation, such as the evolution of barriers to gene flow and the strength of selection. In river networks, studying evolutionary processes in hybrid zones can be especially challenging, given the influence of past and current river properties along with biological species-specific traits. Here, we suggest that a natural hybrid zone between two divergent lineages of the sexually dimorphic Neotropical fish Nematocharax venustus was probably established by secondary contact as a result of a river capture event between the Contas and Pardo river basins. This putative river capture is supported by hydrogeological evidence of elbows of capture, wind gaps and geological faults. The morphological (colour pattern) and genetic (mtDNA and RADseq) variation reveal a clinal transition between parental lineages along the main river, with predominance of F2 hybrids at the centre of the hybrid zone, absence of early generation backcrosses and different levels of hybridization in the tributaries. We highlight that different sources of information are crucial for understanding how the riverscape spatial history influences the connectivity between and within river systems and, consequently, the dynamics of gene flow between freshwater lineages/species.

Is Niagara Falls a barrier to gene flow in riverine fishes? A test using genome-wide SNP data from seven native species

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post-glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water-flow drive the downstream accumulation of genetic diversity, whereas post-glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls - a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls' role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome-wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post-glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species' genetic diversity and highlight the need to better understand post-glacial dispersal pathways.

Quaternary climate changes as speciation drivers in the Amazon floodplains

The role of climate as a speciation driver in the Amazon has long been discussed. Phylogeographic studies have failed to recover synchronous demographic responses across taxa, although recent evidence supports the interaction between rivers and climate in promoting speciation. Most studies, however, are biased toward upland forest organisms, while other habitats are poorly explored and could hold valuable information about major historical processes. We conducted a comparative phylogenomic analysis of floodplain forest birds to explore the effects of historical environmental changes and current connectivity on population differentiation. Our findings support a similar demographic history among species complexes, indicating that the central portion of the Amazon River basin is a suture zone for taxa isolated across the main Amazonian sub-basins. Our results also suggest that changes in the fluvial landscape induced by climate variation during the Mid- and Late Pleistocene drove population isolation, leading to diversification with subsequent secondary contact.

Global determinants of freshwater and marine fish genetic diversity

Genetic diversity is estimated to be declining faster than species diversity under escalating threats, but its spatial distribution remains poorly documented at the global scale. Theory predicts that similar processes should foster congruent spatial patterns of genetic and species diversity, but empirical studies are scarce. Using a mined database of 50,588 georeferenced mitochondrial DNA barcode sequences (COI) for 3,815 marine and 1,611 freshwater fish species respectively, we examined the correlation between genetic diversity and species diversity and their global distributions in relation to climate and geography. Genetic diversity showed a clear spatial organisation, but a weak association with species diversity for both marine and freshwater species. We found a predominantly positive relationship between genetic diversity and sea surface temperature for marine species. Genetic diversity of freshwater species varied primarily across the regional basins and was negatively correlated with average river slope. The detection of genetic diversity patterns suggests that conservation measures should consider mismatching spatial signals across multiple facets of biodiversity.

Biogeographic patterns of genetic diversity are poorly documented, especially for fish species. Here the authors show that (mitochondrial) genetic diversity has global spatial organization patterns with different environmental drivers for marine and freshwater fishes, where genetic diversity is only partly congruent with species richness.

Small fish, large river: Surprisingly minimal genetic structure in a dispersal-limited, habitat specialist fish

Genetic connectivity is expected to be lower in species with limited dispersal ability and a high degree of habitat specialization (intrinsic factors). Also, gene flow is predicted to be limited by habitat conditions such as physical barriers and geographic distance (extrinsic factors). We investigated the effects of distance, intervening pools, and rapids on gene flow in a species, the Tuxedo Darter (Etheostoma lemniscatum), a habitat specialist that is presumed to be dispersal-limited. We predicted that the interplay between these intrinsic and extrinsic factors would limit dispersal and lead to genetic structure even at the small spatial scale of the species range (a 38.6 km river reach). The simple linear distribution of E. lemniscatum allowed for an ideal test of how these factors acted on gene flow and allowed us to test expectations (e.g., isolation-by-distance) of linearly distributed species. Using 20 microsatellites from 163 individuals collected from 18 habitat patches, we observed low levels of genetic structure that were related to geographic distance and rapids, though these factors were not barriers to gene flow. Pools separating habitat patches did not contribute to any observed genetic structure. Overall, E. lemniscatum maintains gene flow across its range and is comprised of a single population. Due to the linear distribution of the species, a stepping-stone model of dispersal best explains the maintenance of gene flow across its small range. In general, our observation of higher-than-expected connectivity likely stems from an adaptation to disperse due to temporally unstable and patchy habitat.

Hybridization drives genetic erosion in sympatric desert fishes of western North America

Many species have evolved or currently coexist in sympatry due to differential adaptation in a heterogeneous environment. However, anthropogenic habitat modifications can either disrupt reproductive barriers or obscure environmental conditions which underlie fitness gradients. In this study, we evaluated the potential for an anthropogenically-mediated shift in reproductive boundaries that separate two historically sympatric fish species (Gila cypha and G. robusta) endemic to the Colorado River Basin using ddRAD sequencing of 368 individuals. We first examined the integrity of reproductive isolation while in sympatry and allopatry, then characterized hybrid ancestries using genealogical assignment tests. We tested for localized erosion of reproductive isolation by comparing site-wise genomic clines against global patterns and identified a breakdown in the drainage-wide pattern of selection against interspecific heterozygotes. This, in turn, allowed for the formation of a hybrid swarm in one tributary, and asymmetric introgression where species co-occur. We also detected a weak but significant relationship between genetic purity and degree of consumptive water removal, suggesting a role for anthropogenic habitat modifications in undermining species boundaries or expanding historically limited introgression. In addition, results from basin-wide genomic clines suggested that hybrids and parental forms are adaptively nonequivalent. If so, then a failure to manage for hybridization will exacerbate the long-term extinction risk in parental populations. These results reinforce the role of anthropogenic habitat modification in promoting interspecific introgression in sympatric species by relaxing divergent selection. This, in turn, underscores a broader role for hybridization in decreasing global biodiversity within rapidly deteriorating environments.

To remain or leave: Dispersal variation and its genetic consequences in benthic freshwater invertebrates

Variation in dispersal capacity may influence population genetic variation and relatedness of freshwater animals thus demonstrating how life-history traits influence patterns and processes that in turn influence biodiversity. The majority of studies have focused on the consequences of dispersal variation in taxa inhabiting riverine systems whose dendritic nature and upstream/downstream gradients facilitate characterizing populations along networks. We undertook extensive, large-scale investigations of the impacts of hydrological connectivity on population genetic variation in two freshwater bryozoan species whose dispersive propagules (statoblasts) are either attached to surfaces (Fredericella sultana) or are released as buoyant stages (Cristatella mucedo) and that live primarily in either lotic (F. sultana) or lentic environments (C. mucedo). Describing population genetic structure in multiple sites characterized by varying degrees of hydrological connectivity within each of three (or four) UK regions enabled us to test the following hypotheses: (1) genetic diversity and gene flow will be more influenced by hydrological connectivity in populations of C. mucedo (because F. sultana dispersal stages are retained); (2) populations of F. sultana will be characterized by greater genetic divergence than those of C. mucedo (reflecting their relative dispersal capacities); and (3) genetic variation will be greatest in F. sultana (reflecting a propensity for genetic divergence as a result of its low dispersal potential). We found that hydrological connectivity enhanced genetic diversity and gene flow among C. mucedo populations but not in F. sultana while higher overall measures of clonal diversity and greater genetic divergence characterized populations of F. sultana. We suggest that genetic divergence over time within F. sultana populations reflects a general constraint of releasing propagules that might eventually be swept to sea when taxa inhabit running waters. In contrast, taxa that primarily inhabit lakes and ponds may colonize across hydrologically connected regions, establishing genetically related populations. Our study contributes more nuanced views about drivers of population genetic structures in passively dispersing freshwater invertebrates as outlined by the Monopolization Hypothesis (Acta Oecologica, 23, 2002, 121) by highlighting how a range of demographic and evolutionary processes reflect life-history attributes of benthic colonial invertebrates (bryozoans) and cyclically parthenogenetic zooplankton. In addition, growing evidence that genetic divergence may commonly characterize populations of diverse groups of riverine taxa suggests that organisms inhabiting lotic systems may be particularly challenged by environmental change. Such change may predispose riverine populations to extinction as a result of genetic divergence combined with limited dispersal and gene flow.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.1tm8705.

Cryptic Lineages and a Population Dammed to Incipient Extinction? Insights into the Genetic Structure of a Mekong River Catfish

An understanding of the genetic composition of populations across management boundaries is vital to developing successful strategies for sustaining biodiversity and food resources. This is especially important in ecosystems where habitat fragmentation has altered baseline patterns of gene flow, dividing natural populations into smaller subpopulations and increasing potential loss of genetic variation through genetic drift. River systems can be highly fragmented by dams built for flow regulation and hydropower. We used reduced-representation sequencing to examine genomic patterns in an exploited catfish, Hemibagrus spilopterus, in a hotspot of biodiversity and hydropower development—the Mekong River basin. Our results revealed the presence of 2 highly divergent coexisting genetic lineages which may be cryptic species. Within the lineage with the greatest sample sizes, pairwise FST values, principal component analysis, and a STRUCTURE analysis all suggest that long-distance migration is not common across the Lower Mekong Basin, even in areas where flood-pulse hydrology has limited genetic divergence. In tributaries, effective population size estimates were at least an order of magnitude lower than in the Mekong mainstream indicating these populations may be more vulnerable to perturbations such as human-induced fragmentation. Fish isolated upstream of several dams in one tributary exhibited particularly low genetic diversity, high amounts of relatedness, and a level of inbreeding (GIS = 0.51) that has been associated with inbreeding depression in other outcrossing species. Our results highlight the importance of assessing genetic structure and diversity in riverine fisheries populations across proposed dam development sites for the preservation of these critically important resources.

Systematic conservation planning for intraspecific genetic diversity

Intraspecific diversity informs the demographic and evolutionary histories of populations, and should be a main conservation target. Although approaches exist for identifying relevant biological conservation units, attempts to identify priority conservation areas for intraspecific diversity are scarce, especially within a multi-specific framework. We used neutral molecular data on six European freshwater fish species (Squalius cephalus, Phoxinus phoxinus, Barbatula barbatula, Gobio occitaniae, Leuciscus burdigalensis and Parachondrostoma toxostoma) sampled at the riverscape scale (i.e. the Garonne-Dordogne river basin, France) to determine hot- and coldspots of genetic diversity, and to identify priority conservation areas using a systematic conservation planning approach. We demonstrate that systematic conservation planning is efficient for identifying priority areas representing a predefined part of the total genetic diversity of a whole landscape. With the exception of private allelic richness (PA), classical genetic diversity indices (allelic richness, genetic uniqueness) were poor predictors for identifying priority areas. Moreover, we identified weak surrogacies among conservation solutions found for each species, implying that conservation solutions are highly species-specific. Nonetheless, we showed that priority areas identified using intraspecific genetic data from multiple species provide more effective conservation solutions than areas identified for single species or on the basis of traditional taxonomic criteria.

Biological Diversity in Headwater Streams

Headwaters, the sources of all stream networks, provide habitats that are unique from other freshwater environments and are used by a specialised subset of aquatic species. The features of headwaters that provide special habitats include predator-free or competitor-free spaces; specific resources (particularly detrital based); and moderate variations in flows, temperature and discharge. Headwaters provide key habitats for all or some life stages for a large number of species across just about all freshwater phyla and divisions. Some features of headwaters, including isolation and small population sizes, have allowed for the evolutionary radiation of many groups of organisms within and beyond those habitats. As small and easily engineered physical spaces, headwaters are easily degraded by streambank development, ditching and even burial. Headwater streams are among the most sensitive of freshwater ecosystems due to their intimate linkage with their catchments and how easily they are impacted. As a unique ecosystem with many specialist species, headwater streams deserve better stewardship.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Genetic diversity in fishes is influenced by habitat type and life-history variation

Populations of fishes are increasingly threatened by over-exploitation, pollution, habitat destruction, and climate change. In order to better understand the factors that can explain the amount of genetic diversity in wild populations of fishes, we collected estimates of genetic diversity (mean heterozygosity and mean rarefied number of alleles per locus) along with habitat associations, conservation status, and life-history information for 463 fish species. We ran a series of phylogenetic generalized least squares models to determine which factors influence genetic diversity in fishes after accounting for shared evolutionary history among related taxa. We found that marine fishes had significantly higher genetic diversity than freshwater fishes with marine fishes averaging 11.3 more alleles per locus than their freshwater counterparts. However, contrary to our expectations, genetic diversity was not found to be lower in threatened versus not-threatened fishes. Finally, we found that both age at maturity and fecundity were negatively related to genetic variation in both marine and freshwater fishes. Our results demonstrate that both life-history characteristics and habitat play a role in shaping patterns of genetic diversity in fishes and should be considered when prioritizing species for conservation.

Geomorphology Influencing the Diversification of Fish in Small-Order Rivers of Neighboring Basins

The current analysis investigates whether the uplift of the Serra da Esperança and the Ponta Grossa Arch in the Serra Geral resulted in ichthyofaunistic changes in adjacent basins. For this, we describe the phylogeographic structure among populations of Trichomycterus collected in hydrographic basins in southern Brazil by using partial nucleotide sequences of the mitochondrial gene Cytochrome C Oxidase subunit I. Analyses revealed that the nomenclature Trichomycterus davisi fails to contain the whole genetic diversity range found in the collected specimens and indicates at least six genetic lineages in Trichomycterus. Diagnostic morphological characteristics not associated to T. davisi could be identified in some specimens from the Iguaçu Piquiri haplogroup, indicating the occurrence of species Trichomycterus stawiarski. The lack of morphological differences among the other clades clearly suggests a cryptic species case. Molecular analyses revealed at least five new species besides T. davisi in the hydrographic basins and support the interpretation that genetic structure in T. davisi species complex is explained by tectonic events intrinsic to the areas of influence of Serra da Esperança and the Ponta Grossa Arch which occurred around 1.7 My.

Loss of dendritic connectivity in southern California's urban riverscape facilitates decline of an endemic freshwater fish

Life history adaptations and spatial configuration of metapopulation networks allow certain species to persist in extreme fluctuating environments, yet long-term stability within these systems relies on the maintenance of linkage habitat. Degradation of such linkages in urban riverscapes can disrupt this dynamic in aquatic species, leading to increased extinction debt in local populations experiencing environment-related demographic flux. We used microsatellites and mtDNA to examine the effects of collapsed network structure in the endemic Santa Ana sucker Catostomus santaanae of southern California, a threatened species affected by natural flood-drought cycles, "boom-and-bust" demography, hybridization and presumed artificial transplantation. Our results show a predominance of drift-mediated processes in shaping population structure and that reverse mechanisms for counterbalancing the genetic effects of these phenomena have dissipated with the collapse of dendritic connectivity. We use approximate Bayesian models to support two cases of artificial transplantation and provide evidence that one of the invaded systems better represents the historic processes that maintained genetic variation within watersheds than any remaining drainages where C. santaanae is considered native. We further show that a stable dry gap in the northern range is preventing genetic dilution of pure C. santaanae persisting upstream of a hybrid assemblage involving a non-native sucker and that local accumulation of genetic variation in the same drainage is influenced by position within the network. This work has important implications for declining species that have historically relied on dendritic metapopulation networks to maintain source-sink dynamics in phasic environments, but no longer possess this capacity in urban-converted landscapes.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Genomic signatures of paleodrainages in a freshwater fish along the southeastern coast of Brazil: genetic structure reflects past riverine properties

Past shifts in connectivity in riverine environments (for example, sea-level changes) and the properties of current drainages can act as drivers of genetic structure and demographic processes in riverine population of fishes. However, it is unclear whether the same river properties that structure variation on recent timescales will also leave similar genomic signatures that reflect paleodrainage properties. By characterizing genetic structure in a freshwater fish species (Hollandichthys multifasciatus) from a system of basins along the Atlantic coast of Brazil we test for the effects of paleodrainages caused by sea-level changes during the Pleistocene. Given that the paleodrainage properties differ along the Brazilian coast, we also evaluate whether estimated genetic diversity within paleodrainages can be explained by past riverine properties (i.e., area and number of rivers in a paleodrainage). Our results demonstrate that genetic structure between populations is not just highly concordant with paleodrainages, but that differences in the genetic diversity among paleodrainages correspond to the joint effect of differences in the area encompassed by, and the number of rivers, within a paleodrainage. Our findings extend the influence of current riverine properties on genetic diversity to those associated with past paleodrainage properties. We discuss how these findings may explain the inconsistent support for paleodrainages in structuring divergence from different global regions and the importance of taking into account past conditions for understanding the high species diversity of freshwater fish that we currently observe in the world, and especially in the Neotropics.

River network architecture, genetic effective size and distributional patterns predict differences in genetic structure across species in a dryland stream fish community

Dendritic ecological network (DEN) architecture can be a strong predictor of spatial genetic patterns in theoretical and simulation studies. Yet, interspecific differences in dispersal capabilities and distribution within the network may equally affect species' genetic structuring. We characterized patterns of genetic variation from up to ten microsatellite loci for nine numerically dominant members of the upper Gila River fish community, New Mexico, USA. Using comparative landscape genetics, we evaluated the role of network architecture for structuring populations within species (pairwise F_ST ) while explicitly accounting for intraspecific demographic influences on effective population size (N_e ). Five species exhibited patterns of connectivity and/or genetic diversity gradients that were predicted by network structure. These species were generally considered to be small-bodied or habitat specialists. Spatial variation of N_e was a strong predictor of pairwise F_ST for two species, suggesting patterns of connectivity may also be influenced by genetic drift independent of network properties. Finally, two study species exhibited genetic patterns that were unexplained by network properties and appeared to be related to nonequilibrium processes. Properties of DENs shape community-wide genetic structure but effects are modified by intrinsic traits and nonequilibrium processes. Further theoretical development of the DEN framework should account for such cases.