Evolutionary processes driving spatial patterns of intraspecific genetic diversity in river ecosystems

Relationships between crayfish population genetic diversity, species richness, and abundance within impounded and unimpounded streams in Alabama, USA

Understanding the relationship between multi-scale processes driving community- and population-level diversity can guide conservation efforts. While the importance of population-level genetic diversity is widely recognized, it is not always assessed for conservation planning, and positive correlations with community-level diversity are sometimes assumed, such that only the latter is measured. We surveyed species richness and cumulative multispecies abundance of crayfishes in impounded and unimpounded streams in the southern Appalachian Mountains (Alabama, USA). We simultaneously assessed levels of population genetic diversity within two focal crayfishes (Faxonius validus and F. erichsonianus) using nuclear (nDNA; inter-simple sequence repeat (ISSR)) and mitochondrial DNA (mtDNA; mitochondrial DNA cytochrome oxidase subunit I (mtCOI)) markers. We then tested for species-genetic diversity correlations (SGDCs), species diversity-abundance correlations (i.e., more individuals hypothesis, MIH), and abundance-genetic diversity correlations (AGDCs) across sites. We also examined the relationship between each of the three different types of correlation (i.e., species richness, cumulative multispecies abundance, and population genetic diversity) and stream habitat characteristics and fragmentation. Surprisingly, based on F. validus mtDNA data, sites with the greatest multispecies abundance had the lowest genetic diversity, indicating a negative AGDC. However, no AGDC was evident from nDNA. There was no evidence of SGDCs for F. validus based on either of the two genetic data types. For F. erichsonianus, there was no evidence for SGDC or AGDC. When considering the community-level data only, there was no support for the MIH. Stream width was positively correlated with F. validus genetic diversity, but negatively correlated with multispecies abundance. Similarly, species richness was positively correlated with stream width in unimpounded streams but negatively correlated with width in impounded streams. These findings indicate that community-level diversity cannot be indiscriminately used as a proxy for population-level diversity without empirically testing this correlation on the focal group. As such, community- and population-level assessments for multiple crayfish species are needed to better understand drivers of diversity and eco-evolutionary processes which will aid in the conservation of this vulnerable taxonomic group.

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.

Highly clustered mating networks in naturally fragmented riparian tree populations

Understanding how spatial patterns of mating and gene flow respond to habitat loss and geographical isolation is a crucial aspect of forest fragmentation genetics. Naturally fragmented riparian tree populations exhibit unique characteristics that significantly influence these patterns. In this study, we investigate mating patterns, pollen-mediated gene flow, and genetic diversity in relict populations of Frangula alnus in southern Spain by testing specific hypotheses related to the riparian habitat. We employ a novel approach that combines paternity analysis, particularly suited for small and isolated populations, with complex network theory and Bayesian models to predict mating likelihood among tree pairs. Our findings reveal a prevalence of short-distance pollination, resulting in spatially driven local mating clusters with a distinct subset of trees being highly connected in the mating network. Additionally, we observe numerous pollination events over distances of hundreds of metres and considerable pollen immigration. Local neighbourhood density is the primary factor influencing within-population mating patterns and pollen dispersal; moreover, mating network properties reflect the population's size and spatial configuration. Conversely, among-population pollen dispersal is mainly determined by tree size, which influences floral display. Our results do not support a major role of directional pollen dispersal in longitudinal trends of genetic diversity. We provide evidence that long-term fragmented tree populations persist in unique environments that shape mating patterns and impose constraints to pollen-mediated gene flow. Nevertheless, even seemingly strongly isolated populations can maintain functional connectivity over extended periods, especially when animal-mediated mating networks promote genetic diversity, as in this riparian tree species.

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.

Limited congruence in phylogeographic patterns observed for riverine predacious beetles sharing distribution along the mountain rivers

Riverine predacious beetles (RPB) (Carabidae, Staphylinidae) are highly diverse and numerous elements of riverine ecosystems. Their historical and contemporary distribution and diversity are highly dependent on natural flow regimes and topography of watercourses. Despite broad knowledge of their ecology, data on population genetic diversity and connectivity are lacking. This study aimed to fill this gap in order to solve two principal hypotheses assuming (i) congruence of phylogeographic patterns observed for RPB indicating that they share a common history and the ecological adaptations to the dynamic environment, (ii) genetic structuration of populations according to river basins. The Carpathian populations of four ground beetles and three rove beetles were examined using cytochrome oxidase and arginine kinase sequencing. There are substantial differences in RPB demographic history and current genetic diversity. Star-like phylogeny of Bembidion and complex haplotype networks of Paederus/Paederidus, with some haplotypes being drainage-specific and others found in distant populations, indicate a general lack of isolation by distance. Signs of recent demographic expansion were detected for most RPB with the latest population collapse for some rove beetles. To some extent, migration of examined species has to be limited by watersheds. Observed phylogeographic patterns are essential for correctly understanding RPB meta-population functioning.

Genetic erosion reduces biomass temporal stability in wild fish populations

Genetic diversity sustains species adaptation. However, it may also support key ecosystems functions and services, for example biomass production, that can be altered by the worldwide loss of genetic diversity. Despite extensive experimental evidence, there have been few attempts to empirically test whether genetic diversity actually promotes biomass and biomass stability in wild populations. Here, using long-term demographic wild fish data from two large river basins in southwestern France, we demonstrate through causal modeling analyses that populations with high genetic diversity do not reach higher biomasses than populations with low genetic diversity. Nonetheless, populations with high genetic diversity have much more stable biomasses over recent decades than populations having suffered from genetic erosion, which has implications for the provision of ecosystem services and the risk of population extinction. Our results strengthen the importance of adopting prominent environmental policies to conserve this important biodiversity facet.

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.

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.

Evolution of an Amazonian Fish Is Driven by Allopatric Divergence Rather Than Ecological Divergence

Lowland central Amazonia is characterized by heterogeneous riverscapes dominated by two chemically divergent water types: black (ion-poor, rich in dissolved organic carbonate and acidic) and white (rich in nutrient and turbid) waters. Recent phylogeographic and genomic studies have associated the ecotone formed by these environments to an ecologically driven genetic divergence between fish present in both water types. With the objective of better understanding the evolutionary forces behind the central Amazonian teleostean diversification, we sampled 240 Mesonauta festivus from 12 sites on a wide area of the Amazonian basin. These sites included three confluences of black and white water environments to seek for repeated evidences of ecological divergence at the junction of these ecotones. Results obtained through our genetic assessment based on 41,268 single nucleotide polymorphism (SNPs) contrast with previous findings and support a low influence of diverging water physicochemical characteristics on the genetic structuration of M. festivus populations. Conversely, we detected patterns of isolation by downstream water current and evidence of past events of vicariance potentially linked to the Amazon River formation. Using a combination of population genetics, phylogeographic analysis and environmental association models, we decomposed the spatial and environmental genetic variances to assess which evolutionary forces shaped inter-population differences in M. festivus’ genome. Our sampling design, comprising three confluences of black and white water rivers, supports the main role of evolution by allopatry. While an ecologically driven evolution admittedly played a role in Amazonian fish diversification, we argue that neutral evolutionary processes explain most of the divergence between M. festivus populations.

Landscape connectivity and genetic structure in a mainstem and a tributary stonefly (Plecoptera) species using a novel reference genome

Understanding how environmental variation influences population genetic structure can help predict how environmental change influences population connectivity, genetic diversity, and evolutionary potential. We used riverscape genomics modelling to investigate how climatic and habitat variables relate to patterns of genetic variation in two stonefly species, one from mainstem river habitats (Sweltsa coloradensis) and one from tributaries (Sweltsa fidelis) in 40 sites in northwest Montana, USA. We produced a draft genome assembly for S. coloradensis (N50 = 0.251 Mbp, BUSCO > 95% using “insecta_ob9” reference genes). We genotyped 1930 SNPs in 372 individuals for S. coloradensis and 520 SNPs in 153 individuals for S. fidelis. We found higher genetic diversity for S. coloradensis compared to S. fidelis, but nearly identical genetic differentiation among sites within each species (both had global loci median FST = 0.000), despite differences in stream network location. For landscape genomics and testing for selection, we produced a less stringently filtered data set (3454 and 1070 SNPs for S. coloradensis and S. fidelis, respectively). Environmental variables (mean summer precipitation, slope, aspect, mean June stream temperature, land cover type) were correlated with 19 putative adaptive loci for S. coloradensis. but there was only one putative adaptive locus for S. fidelis (correlated with aspect). Interestingly, we also detected potential hybridization between multiple Sweltsa species which has never been previously detected. Studies like ours, that test for adaptive variation in multiple related species are needed to help assess landscape connectivity and the vulnerability of populations and communities to environmental change.

Landscape genetics of the tropical willow Salix humboldtiana: influence of climate, salinity, and orography in an altitudinal gradient

PREMISE

Gene flow in riparian ecosystems is influenced by landscape features such as orography, climate, and salinity. The downstream increase in genetic diversity (DIGD) hypothesis states that the unidirectionality of the watercourse causes an accumulation of genetic diversity toward downstream populations, while upstream populations are more structured and less diverse, especially in water-dispersed organisms.

METHODS

We used chloroplast and nuclear microsatellites to characterize genetic diversity, structure, and gene flow patterns among populations of Salix humboldtiana across an elevation and salinity gradient on three rivers (Actopan, Antigua, and Blanco) in Mexico. We used optimization of resistance surface methods to determine whether genetic distances between populations are correlated with landscape features.

RESULTS

Positive F_IS values evidenced biparental inbreeding in some populations, particularly at higher elevations where lower niche availability constrains colonization and persistence. Four genetic groups were distinguished, corresponding to populations on the Actopan and Antigua rivers and upstream and downstream on the Blanco, but with high admixture between populations on the Actopan and Antigua rivers. Higher gene flow rates were found among proximate populations on the same river than among different rivers. Genetic diversity increased toward the river mouths, in support of the DIGD hypothesis, probably due to greater niche availability and larger population size. Differences among rivers in precipitation patterns and salinity, as well as geographic distance, were significant predictors of gene flow.

CONCLUSIONS

Our results depict that the DIGD and gene flow patterns in S. humboldtiana result from the complex interaction among physiography, climate, river salinity, and life-history traits of the species.

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.

Opportunities and challenges of macrogenetic studies

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

Riverscape genetics in brook lamprey: genetic diversity is less influenced by river fragmentation than by gene flow with the anadromous ecotype

Understanding the effect of human-induced landscape fragmentation on gene flow and evolutionary potential of wild populations has become a major concern. Here, we investigated the effect of riverscape fragmentation on patterns of genetic diversity in the freshwater resident European brook lamprey (Lampetra planeri) that has a low ability to pass obstacles to migration. We tested the hypotheses of (i) asymmetric gene flow following water current and (ii) an effect of gene flow with the closely related anadromous river lamprey (L. fluviatilis) ecotype on L. planeri genetic diversity. We genotyped 2472 individuals, including 225 L. fluviatilis, sampled from 81 sites upstream and downstream barriers to migration, in 29 western European rivers. Linear modelling revealed a strong positive relationship between genetic diversity and the distance from the river source, consistent with expected patterns of decreased gene flow into upstream populations. However, the presence of anthropogenic barriers had a moderate effect on spatial genetic structure. Accordingly, we found evidence for downstream-directed gene flow, supporting the hypothesis that barriers do not limit dispersal mediated by water flow. Downstream L. planeri populations in sympatry with L. fluviatilis displayed consistently higher genetic diversity. We conclude that genetic drift and slight downstream gene flow drive the genetic make-up of upstream L. planeri populations whereas gene flow between ecotypes maintains higher levels of genetic diversity in L. planeri populations sympatric with L. fluviatilis. We discuss the implications of these results for the design of conservation strategies of lamprey, and other freshwater organisms with several ecotypes, in fragmented dendritic river networks.

Range reduction of Oblong Rocksnail, Leptoxis compacta, shapes riverscape genetic patterns

Many freshwater gastropod species face extinction, including 79% of species in the family Pleuroceridae. The Oblong Rocksnail, Leptoxis compacta, is a narrow range endemic pleurocerid from the Cahaba River basin in central Alabama that has seen rapid range contraction in the last 100 years. Such a decline is expected to negatively affect genetic diversity in the species. However, precise patterns of genetic variation and gene flow across the restricted range of L. compacta are unknown. This lack of information limits our understanding of human impacts on the Cahaba River system and Pleuroceridae. Here, we show that L. compacta has likely seen a species-wide decline in genetic diversity, but remaining populations have relatively high genetic diversity. We also report a contemporary range extension compared to the last published survey. Our findings indicate that historical range contraction has resulted in the absence of common genetic patterns seen in many riverine taxa like isolation by distance as the small distribution of L. compacta allows for relatively unrestricted gene flow across its remaining range despite limited dispersal abilities. Two collection sites had higher genetic diversity than others, and broodstock sites for future captive propagation and reintroduction efforts should utilize sites identified here as having the highest genetic diversity. Broadly, our results support the hypothesis that range contraction will result in the reduction of species-wide genetic diversity, and common riverscape genetic patterns cannot be assumed to be present in species facing extinction risk.

A river runs through it: The causes, consequences, and management of intraspecific diversity in river networks

Rivers are fascinating ecosystems in which the eco-evolutionary dynamics of organisms are constrained by particular features, and biologists have developed a wealth of knowledge about freshwater biodiversity patterns. Over the last 10 years, our group used a holistic approach to contribute to this knowledge by focusing on the causes and consequences of intraspecific diversity in rivers. We conducted empirical works on temperate permanent rivers from southern France, and we broadened the scope of our findings using experiments, meta-analyses, and simulations. We demonstrated that intraspecific (genetic) diversity follows a spatial pattern (downstream increase in diversity) that is repeatable across taxa (from plants to vertebrates) and river systems. This pattern can result from interactive processes that we teased apart using appropriate simulation approaches. We further experimentally showed that intraspecific diversity matters for the functioning of river ecosystems. It indeed affects not only community dynamics, but also key ecosystem functions such as litter degradation. This means that losing intraspecific diversity in rivers can yield major ecological effects. Our work on the impact of multiple human stressors on intraspecific diversity revealed that-in the studied river systems-stocking of domestic (fish) strains strongly and consistently alters natural spatial patterns of diversity. It also highlighted the need for specific analytical tools to tease apart spurious from actual relationships in the wild. Finally, we developed original conservation strategies at the basin scale based on the systematic conservation planning framework that appeared pertinent for preserving intraspecific diversity in rivers. We identified several important research avenues that should further facilitate our understanding of patterns of local adaptation in rivers, the identification of processes sustaining intraspecific biodiversity-ecosystem function relationships, and the setting of reliable conservation plans.

Generation and application of river network analogues for use in ecology and evolution

Several key processes in freshwater ecology are governed by the connectivity inherent to dendritic river networks. These have extensively been analyzed from a geomorphological and hydrological viewpoint, yet structures classically used in ecological modeling have been poorly representative of the structure of real river basins, often failing to capture well-known scaling features of natural rivers. Pioneering work identified optimal channel networks (OCNs) as spanning trees reproducing all scaling features characteristic of natural stream networks worldwide. While OCNs have been used to create landscapes for studies on metapopulations, biodiversity, and epidemiology, their generation has not been generally accessible.Given the increasing interest in dendritic riverine networks by ecologists and evolutionary biologists, we here present a method to generate OCNs and, to facilitate its application, we provide the R-package OCNet. Owing to the stochastic process generating OCNs, multiple network replicas spanning the same surface can be built; this allows performing computational experiments whose results are irrespective of the particular shape of a single river network. The OCN construct also enables the generation of elevational gradients derived from the optimal network configuration, which can constitute three-dimensional landscapes for spatial studies in both terrestrial and freshwater realms. Moreover, the package provides functions that aggregate OCNs into an arbitrary number of nodes, calculate several descriptors of river networks, and draw relevant network features.We describe the main functionalities of the package and its integration with other R-packages commonly used in spatial ecology. Moreover, we exemplify the generation of OCNs and discuss an application to a metapopulation model for an invasive riverine species.In conclusion, OCNet provides a powerful tool to generate realistic river network analogues for various applications. It thereby allows the design of spatially realistic studies in increasingly impacted ecosystems and enhances our knowledge on spatial processes in freshwater ecology in general.

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.

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.

Disturbance reverses classic biodiversity predictions in river-like landscapes

Global analyses of biodiversity consistently reveal recurrent patterns of species distributions worldwide. However, unveiling the specific mechanisms behind those patterns remains logistically challenging, yet necessary for reliable biodiversity forecasts. Here, we combine theory and experiments to investigate the processes underlying spatial biodiversity patterns in dendritic, river-like landscapes, iconic examples of highly threatened ecosystems. We used geometric scaling properties, common to all rivers, to show that the distribution of biodiversity in these landscapes fundamentally depends on how ecological selection is modulated across space: while uniform ecological selection across the network leads to higher diversity in downstream confluences, this pattern can be inverted by disturbances when population turnover (i.e. local mortality) is higher upstream than downstream. Higher turnover in small headwater patches can slow down ecological selection, increasing local diversity in comparison to large downstream confluences. Our results show that disturbance-mediated slowing down of competitive exclusion can generate a specific transient signature in terms of biodiversity distribution when applied over a spatial gradient of disturbance, which is a common feature of many river landscapes.

Examining the role of parasites in limiting unidirectional gene flow between lake and river sticklebacks

Parasites are important selective agents with the potential to limit gene flow between host populations by shaping local host immunocompetence. We report on a contact zone between lake and river three-spined sticklebacks (Gasterosteus aculeatus) that offers the ideal biogeographic setting to explore the role of parasite-mediated selection on reproductive isolation. A waterfall acts as a natural barrier and enforces unidirectional migration from the upstream river stickleback population to the downstream river and lake populations. We assessed population genetic structure and parasite communities over four years. In a set of controlled experimental infections, we compared parasite susceptibility of upstream and downstream fish by exposing laboratory-bred upstream river and lake fish, as well as hybrids, to two common lake parasite species: a generalist trematode parasite, Diplostomum pseudospathaceum, and a host-specific cestode, Schistocephalus solidus. We found consistent genetic differentiation between upstream and downstream populations across four sampling years, even though the downstream river consisted of ~10% first-generation migrants from the upstream population as detected by parentage analysis. Fish in the upstream population had lower genetic diversity and were strikingly devoid of macroparasites. Through experimental infections, we demonstrated that upstream fish and their hybrids had higher susceptibility to parasite infections than downstream fish. Despite this, naturally sampled upstream migrants were less infected than downstream residents. Thus, migrants coming from a parasite-free environment may enjoy an initial fitness advantage, but their descendants seem likely to suffer from higher parasite loads. Our results suggest that adaptation to distinct parasite communities can influence stickleback invasion success and may represent a barrier to gene flow, even between close and connected populations.

Combining population genomics and forward simulations to investigate stocking impacts: A case study of Muskellunge (Esox masquinongy) from the St. Lawrence River basin

Understanding the genetic and evolutionary impacts of stocking on wild fish populations has long been of interest as negative consequences such as reduced fitness and loss of genetic diversity are commonly reported outcomes. In an attempt to sustain a fishery, managers implemented nearly five decades of extensive stocking of over a million Muskellunge (Esox masquinongy), a native species in the Lower St. Lawrence River (Québec, Canada). We investigated the effect of this stocking on population genetic structure and allelic diversity in the St. Lawrence River in addition to tributaries and several stocked inland lakes. Using genotype by sequencing, we genotyped 643 individuals representing 22 locations and combined this information with forward simulations to investigate the genetic consequences of long-term stocking. Individuals native to the St. Lawrence watershed were genetically differentiated from stocking sources and tributaries, and inland lakes were naturally differentiated from the main river. Empirical data and simulations within the St. Lawrence River revealed weak stocking effects on admixture patterns. Our data suggest that the genetic structure associated with stocked fish was diluted into its relatively large effective population size. This interpretation is also consistent with a hypothesis that selection against introgression was in operation and relatively efficient within the large St. Lawrence River system. In contrast, smaller populations from adjacent tributaries and lakes displayed greater stocking-related admixture that resulted in comparatively higher heterozygosity than the St. Lawrence. Finally, individuals from inland lakes that were established by stocking maintained a close affinity with their source populations. This study illustrated a benefit of combining extensive genomic data with forward simulations for improved inference regarding population-level genetic effects of long-term stocking, and its relevance for fishery management decision making.

Genetic changes caused by restocking and hydroelectric dams in demographically bottlenecked brown trout in a transnational subarctic riverine system

Habitat discontinuity, anthropogenic disturbance, and overharvesting have led to population fragmentation and decline worldwide. Preservation of remaining natural genetic diversity is crucial to avoid continued genetic erosion. Brown trout (Salmo trutta L.) is an ideal model species for studying anthropogenic influences on genetic integrity, as it has experienced significant genetic alterations throughout its natural distribution range due to habitat fragmentation, overexploitation, translocations, and stocking. The Pasvik River is a subarctic riverine system shared between Norway, Russia, and Finland, subdivided by seven hydroelectric power dams that destroyed about 70% of natural spawning and nursing areas. Stocking is applied in certain river parts to support the natural brown trout population. Adjacent river segments with different management strategies (stocked vs. not stocked) facilitated the simultaneous assessment of genetic impacts of dams and stocking based on analyses of 16 short tandem repeat loci. Dams were expected to increase genetic differentiation between and reduce genetic diversity within river sections. Contrastingly, stocking was predicted to promote genetic homogenization and diversity, but also potentially lead to loss of private alleles and to genetic erosion. Our results showed comparatively low heterozygosity and clear genetic differentiation between adjacent sections in nonstocked river parts, indicating that dams prevent migration and contribute to genetic isolation and loss of genetic diversity. Furthermore, genetic differentiation was low and heterozygosity relatively high across stocked sections. However, in stocked river sections, we found signatures of recent bottlenecks and reductions in private alleles, indicating that only a subset of individuals contributes to reproduction, potentially leading to divergence away from the natural genetic state. Taken together, these results indicate that stocking counteracts the negative fragmentation effects of dams, but also that stocking practices should be planned carefully in order to ensure long-term preservation of natural genetic diversity and integrity in brown trout and other species in regulated river systems.

A spatial stream-network approach assists in managing the remnant genetic diversity of riparian forests

Quantifying the genetic diversity of riparian trees is essential to understand their chances to survive hydroclimatic alterations and to maintain their role as foundation species modulating fluvial ecosystem processes. However, the application of suitable models that account for the specific dendritic structure of hydrographic networks is still incipient in the literature. We investigate the roles of ecological and spatial factors in driving the genetic diversity of Salix salviifolia, an Iberian endemic riparian tree, across the species latitudinal range. We applied spatial stream-network models that aptly integrate dendritic features (topology, directionality) to quantify the impacts of multiple scale factors in determining genetic diversity. Based on the drift hypothesis, we expect that genetic diversity accumulates downstream in riparian ecosystems, but life history traits (e.g. dispersal patterns) and abiotic or anthropogenic factors (e.g. drought events or hydrological alteration) might alter expected patterns. Hydrological factors explained the downstream accumulation of genetic diversity at the intermediate scale that was likely mediated by hydrochory. The models also suggested upstream gene flow within basins that likely occurred through anemophilous and entomophilous pollen and seed dispersal. Higher thermicity and summer drought were related to higher population inbreeding and individual homozygosity, respectively, suggesting that increased aridity might disrupt the connectivity and mating patterns among and within riparian populations.

Riverscape genetic variation, migration patterns, and morphological variation of the threatened Round Rocksnail, Leptoxis ampla

Within riverine systems, headwater populations are hypothesized to harbour higher amounts of genetic distinctiveness than populations in the main stem of a river and display increased genetic diversity in large, downstream habitats. However, these hypotheses were mostly developed with insects and fish, and they have not been tested on many invertebrate lineages. Pleuroceridae gastropods are of particular ecological importance to rivers of eastern North America, sometimes comprising over 90% of macroinvertebrate biomass. Yet, virtually nothing is known of pleurocerid landscape genetics, including whether genetic diversity follows predictions made by hypotheses developed on more mobile species. Moreover, the commonly repeated hypothesis that intraspecific morphological variation in gastropods results from ecophenotypic plasticity has not been well tested on pleurocerids. Using 2bRAD-seq to discover single nucleotide polymorphisms, we show that the threatened, Cahaba River endemic pleurocerid, Leptoxis ampla, has limited gene flow among populations and that migration is downstream-biased, conflicting with previous hypotheses. Both tributary and main stem populations harbour unique genomic profiles, and genetic diversity was highest in downstream populations. Furthermore, L. ampla shell morphology was more correlated with genetic differences among individuals and populations than habitat characteristics. We anticipate similar genetic and demographic patterns to be seen in other pleurocerids, and hypotheses about gene flow and population demographics that were based on more mobile taxa often, but not always, apply to freshwater gastropods. From a conservation standpoint, genetic structure of L. ampla populations suggests distinctive genetic diversity is lost with localized extirpation, a phenomenon common across the range of Pleuroceridae.

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.

The Robust Classification Model Based on Combinatorial Features

Analyzing the disease data from the view of combinatorial features may better characterize the disease phenotype. In this study, a novel method is proposed to construct feature combinations and a classification model (CFC-CM) by mining key feature relationships. CFC-CM iteratively tests for differences in the feature relationship between different groups. To do this, it uses a modified $k$k-top-scoring pair (M-$k$k-TSP) algorithm and then selects the most discriminative feature pairs in the current feature set to infer the combinatorial features and build the classification model. Compared with support vector machines, random forests, least absolute shrinkage and selection operator, elastic net, and M-$k$k-TSP, the superior performance of CFC-CM on nine public gene expression datasets validates its potential for more precise identification of complex diseases. Subsequently, CFC-CM was applied to two metabolomics datasets, it obtained accuracy rates of $88.73\pm 2.06\%$88.73±2.06% and $79.11\pm 2.70\%$79.11±2.70% in distinguishing between hepatocellular carcinoma and hepatic cirrhosis groups and between acute kidney injury (AKI) and non-AKI samples, results superior to those of the other five methods. In summary, the better results of CFC-CM show that in contrast to molecules and combinations constituted by just two features, the combinations inferred by appropriate number of features could better identify the complex diseases.

Founder effects drive the genetic structure of passively dispersed aquatic invertebrates

Populations of passively dispersed organisms in continental aquatic habitats typically show high levels of neutral genetic differentiation despite their high dispersal capabilities. Several evolutionary factors, including founder events, local adaptation, and life cycle features such as high population growth rates and the presence of propagule banks, have been proposed to be responsible for this paradox. Here, we have modeled the colonization process to assess the impact of migration rate, population growth rate, population size, local adaptation and life-cycle features on the population genetic structure in these organisms. Our simulations show that the strongest effect on population structure are persistent founder effects, resulting from the interaction of a few population founders, high population growth rates, large population sizes and the presence of diapausing egg banks. In contrast, the role of local adaptation, genetic hitchhiking and migration is limited to small populations in these organisms. Our results indicate that local adaptation could have different impact on genetic structure in different groups of zooplankters.

Multilocus phylogeny of the zebra mussel family Dreissenidae (Mollusca: Bivalvia) reveals a fourth Neotropical genus sister to all other genera

Dreissenidae is one of the most economically and ecologically important families of freshwater and estuarine mollusks. Fourteen extant species and three genera are currently recognized: Congeria contains three species from karst caves along the eastern Adriatic coast and one from the Orinoco River of Venezuela, Dreissena contains six species native to Eastern European rivers and estuaries, and Mytilopsis contains three species from the Gulf of Mexico, Caribbean, and northwestern coast of South America and one from the Tocantins River of Brazil. Previous molecular phylogenetic studies have examined all species except those from South American rivers, and found each genus to be monophyletic with Congeria and Mytilopsis forming a clade sister to Dreissena. We present the first multilocus phylogeny of Dreissenidae inclusive of South American riverine species. Bayesian and maximum likelihood analyses of a 3085 bp alignment consisting of mitochondrial (COI and 16S) and nuclear (18S and 28S) gene regions found Neotropical species to be consistently and strongly supported as sister to all other dreissenids, although incomplete sequencing of the single Orinoco specimen obscured Neotropical monophyly. Our intergeneric relationships are inconsistent with an extensive fossil record suggesting that dreissenids originated in Europe approximately 30 My before dispersing to the Western Hemisphere. Fossil-calibrated analyses indicated that Neotropical dreissenids diverged from European lineages in the mid to late Eocene (∼39.3 Ma), and Brazilian and Guiana shield populations diversified during the Oligocene to Miocene. We erect the new genus Rheodreissena for all Neotropical freshwater dreissenids and present haplotype data indicative of at least three species. Widespread anthropogenic alteration of the middle Xingu River and lower Amazon threatens the persistence of these endemic, poorly studied mussels and may facilitate introduction beyond their native range.

Comparative genetic isolation patterns for multiple headwater fishes in three geographic regions

Headwater-resident fishes may be prone to a high rate of isolation and a pronounced hierarchical genetic structure due to a combination of niche preference and dendritic effects of river networks. Genetic isolation patterns were compared using microsatellites in six headwater fishes, Fundulus olivaceus, Semotilus atromaculatus, Erimyzon claviformis, Etheostoma artesiae, Etheostoma whipplei and Etheostoma parvipinne, in three geographic regions that included drainages of small, medium and large sizes in the southern United States. All species showed hierarchical nesting of genetic populations and there were clear and mostly consistent differences between species and regions that were identified through summary statistics derived from two independent analyses. For species comparisons, a high isolation grouping (increased number of isolated genetic clusters or sites within regions) and a low-isolation grouping (decreased number of clusters or sites) were identified. Species group placement was related to niche breadth along the river continuum and presumed abundance and variability of preferred microhabitats, with increased headwater specialization among species being associated with placement in the high-isolation grouping. There was a weakly significant positive effect of drainage size on the number of isolated clusters or sites across all species. Regional patterns were shared in two species, with the region containing the smallest drainages having lower rates of isolation in both datasets. This study shows the effects of regional and species characteristics on genetic isolation for headwater species, which are especially prone to isolation due to spatial, dendritic effects of river networks.

Inferring Causalities in Landscape Genetics: An Extension of Wright's Causal Modeling to Distance Matrices

Identifying landscape features that affect functional connectivity among populations is a major challenge in fundamental and applied sciences. Landscape genetics combines landscape and genetic data to address this issue, with the main objective of disentangling direct and indirect relationships among an intricate set of variables. Causal modeling has strong potential to address the complex nature of landscape genetic data sets. However, this statistical approach was not initially developed to address the pairwise distance matrices commonly used in landscape genetics. Here, we aimed to extend the applicability of two causal modeling methods-that is, maximum-likelihood path analysis and the directional separation test-by developing statistical approaches aimed at handling distance matrices and improving functional connectivity inference. Using simulations, we showed that these approaches greatly improved the robustness of the absolute (using a frequentist approach) and relative (using an information-theoretic approach) fits of the tested models. We used an empirical data set combining genetic information on a freshwater fish species (Gobio occitaniae) and detailed landscape descriptors to demonstrate the usefulness of causal modeling to identify functional connectivity in wild populations. Specifically, we demonstrated how direct and indirect relationships involving altitude, temperature, and oxygen concentration influenced within- and between-population genetic diversity of G. occitaniae.

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.

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation‐by‐distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation‐by‐distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.

Genetic diversity and gene flow decline with elevation in montane mayflies

Montane environments around the globe are biodiversity 'hotspots' and important reservoirs of genetic diversity. Montane species are also typically more vulnerable to environmental change than their low-elevation counterparts due to restricted ranges and dispersal limitations. Here we focus on two abundant congeneric mayflies (Baetis bicaudatus and B. tricaudatus) from montane streams over an elevation gradient spanning 1400 m. Using single-nucleotide polymorphism genotypes, we measured population diversity and vulnerability in these two species by: (i) describing genetic diversity and population structure across elevation gradients to identify mechanisms underlying diversification; (ii) performing spatially explicit landscape analyses to identify environmental drivers of differentiation; and (iii) identifying outlier loci hypothesized to underlie adaptive divergence. Differences in the extent of population structure in these species were evident depending upon their position along the elevation gradient. Heterozygosity, effective population sizes and gene flow all declined with increasing elevation, resulting in substantial population structure in the higher elevation species (B. bicaudatus). At lower elevations, populations of both species are more genetically similar, indicating ongoing gene flow. Isolation by distance was detected at lower elevations only, whereas landscape barriers better predicted genetic distance at higher elevations. At higher elevations, dispersal was restricted due to landscape effects, resulting in greater population isolation. Our results demonstrate differentiation over small spatial scales along an elevation gradient, and highlight the importance of preserving genetic diversity in more isolated high-elevation populations.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

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.

A portrait of a sucker using landscape genetics: how colonization and life history undermine the idealized dendritic metapopulation

Dendritic metapopulations have been attributed unique properties by in silico studies, including an elevated genetic diversity relative to a panmictic population of equal total size. These predictions have not been rigorously tested in nature, nor has there been full consideration of the interacting effects among contemporary landscape features, colonization history and life history traits of the target species. We tested for the effects of dendritic structure as well as the relative importance of life history, environmental barriers and historical colonization on the neutral genetic structure of a longnose sucker (Catostomus catostomus) metapopulation in the Kogaluk watershed of northern Labrador, Canada. Samples were collected from eight lakes, genotyped with 17 microsatellites, and aged using opercula. Lakes varied in differentiation, historical and contemporary connectivity, and life history traits. Isolation by distance was detected only by removing two highly genetically differentiated lakes, suggesting a lack of migration-drift equilibrium and the lingering influence of historical factors on genetic structure. Bayesian analyses supported colonization via the Kogaluk's headwaters. The historical concentration of genetic diversity in headwaters inferred by this result was supported by high historical and contemporary effective sizes of the headwater lake, T-Bone. Alternatively, reduced allelic richness in headwaters confirmed the dendritic structure's influence on gene flow, but this did not translate to an elevated metapopulation effective size. A lack of equilibrium and upstream migration may have dampened the effects of dendritic structure. We suggest that interacting historical and contemporary factors prevent the achievement of the idealized traits of a dendritic metapopulation in nature.

Applying gene flow science to environmental policy needs: a boundary work perspective

One application of gene flow science is the policy arena. In this article, we describe two examples in which the topic of gene flow has entered into the U.S. national environmental policymaking process: regulation of genetically engineered crops and clarification of the jurisdictional scope of the Clean Water Act. We summarize both current scientific understanding and the legal context within which gene flow science has relevance. We also discuss the process by which scientific knowledge has been synthesized and communicated to decision-makers in these two contexts utilizing the concept of 'boundary work'. Boundary organizations, the work they engage in to bridge the worlds of science, policy, and practice, and the boundary objects they produce to translate scientific knowledge existed in both examples. However, the specific activities and attributes of the objects produced varied based on the needs of the decision-makers. We close with suggestions for how scientists can contribute to or engage in boundary work with policymakers.

Spatial patterns of genetic diversity, community composition and occurrence of native and non-native amphipods in naturally replicated tributary streams

Background

Worldwide, natural communities are invaded by non-native species, with potentially devastating effects on the native communities. A large part of past research aimed at finding traits and characteristics of the invading species or the invaded community explaining observed invasions. Only recently, the focus shifted on the spatial patterns during invasions per se. Empirical data, however, are limited, as invasions are often unique incidences of a complex spatio-temporal process. In order to identify generalities of invasion patterns, we studied 13 naturally replicated tributary streams draining into Lake Constance, and studied the occurrence of native and non-native amphipods along linear transects from the stream outlets to the upstream headwater reaches.

Results

We found repeated spatial patterns of community composition and the occurrence of native and non-native amphipod species across two different years. Specifically, occurrence as well as abundance of two non-native amphipod species decreased from the stream outlets at the lake site towards upstream headwater reaches. Populations of the most common native amphipod species were largest at the uppermost headwater reaches. All populations of this native species, however, showed significant signals of recent genetic bottlenecks, irrespective of the stream position and occurrence of non-native species. Contrary to our expectations, this native species also showed no longitudinal genetic differentiation within individual tributaries as postulated for headwater versus outlet populations.

Conclusions

Our results indicate that invasions of river-systems may overall follow predictable patterns on the level of spatial distributions and community composition. However, effects of invading organisms on the genetic diversity and genetic structure of native populations observed at larger scales may not necessarily be directly reflected at the scale of smaller tributaries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-016-0079-7) contains supplementary material, which is available to authorized users.