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

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Does life history mediate discharge as a driver of multi-decadal changes in populations of freshwater fish?

Understanding how and why the size of populations varies is critical knowledge for conservation and management. While considerable work has explored how different demographic parameters affect population growth, less is known the drivers of variability in these parameters. Long-term time series tracking population size that are coupled with empirical data to examine the relative importance of different drivers are rare, especially in freshwater systems. Even rarer are studies that collect this information concurrently from multiple species with contrasting life history strategies in the same system to assess whether population size and the relative importance of drivers also vary. We studied changes in the abundance and size structure of four native freshwater fish species in the Murray River, southeastern Australia, over a continuous 19-yr period. Two species with traits typical of "equilibrium" species (Murray cod Maccullochella peelii and trout cod Maccullochella macquariensis) and two with traits of "periodic" species (golden perch Macquaria ambigua and silver perch Bidyanus bidyanus) were sampled annually and capture-mark-recapture modeling was used to ask (1) how did population size change during this period, (2) how were changes in population size related to variability in hydrology, and (3) how were changes in population size driven by different processes (local recruitment or migration events)? Populations of all four species varied throughout the study, and our results are consistent with the notion that local recruitment is an important driver of this variability for Murray cod and trout cod, whereas immigration is more important for the two other species. Increases in spring river discharge strongly influenced these responses for trout cod and golden perch. Our study provides fundamental insights into the population dynamics of these valued species, and how management strategies might differ based on their life histories. Management should focus on allowing connectivity for golden and silver perch, and on promoting local scale recruitment and survival for Murray cod and trout cod. More generally, our study highlights the importance of understanding the processes underpinning population persistence, how these processes may vary for different species, and ultimately how this knowledge can inform targeted management actions.

Fish out of water: Genomic insights into persistence of rainbowfish populations in the desert

How populations of aquatic fauna persist in extreme desert environments is an enigma. Individuals often breed and disperse during favorable conditions. Theory predicts that adaptive capacity should be low in small populations, such as in desert fishes. We integrated satellite-derived surface water data and population genomic diversity from 20,294 single-nucleotide polymorphisms across 344 individuals to understand metapopulation persistence of the desert rainbowfish (Melanotaenia splendida tatei) in central Australia. Desert rainbowfish showed very small effective population sizes, especially at peripheral populations, and low connectivity between river catchments. Yet, there was no evidence of population-level inbreeding and a signal of possible adaptive divergence associated with aridity was detected. Candidate genes for local adaptation included functions related to environmental cues and stressful conditions. Eco-evolutionary modeling showed that positive selection in refugial subpopulations combined with connectivity during flood periods can enable retention of adaptive diversity. Our study suggests that adaptive variation can be maintained in small populations and integrate with neutral metapopulation processes to allow persistence in the desert.

Genotypes of Rhizophora Propagules From a Non-mangrove Beach Provide Evidence of Recent Long-Distance Dispersal

Dispersal plays a crucial role in the connectivity of established mangrove populations and in species range dynamics. As species ranges shift in response to climate change, range expansions can occur from incremental short-distance dispersal events and from stochastic long-distance dispersal events. Most population genetic research dealt with historically accumulated events though evidence of actual propagule dispersal allows to estimate genotypic features and origin of founders. In this study, we aim to disentangle a contemporary dispersal event. Using microsatellite markers, we genotyped 60 Rhizophora racemosa drift propagules obtained on a bare unforested coastal area in southern Cameroon, estimated their relationship to 109 adult trees from most proximate sites (which were 3–85 km away), and assessed their relative difference with 873 trees of major mangrove areas (> 300 km) along the Cameroonian coastline. Proximate mangrove populations were considered as potential source populations in assignment tests. However, drift propagules could not be assigned to any of the Cameroonian mangrove sites and were genetically isolated from Cameroonian populations. Drift propagules showed higher levels of genetic diversity and private alleles giving a higher relatedness to each other than to any putative source population. Chloroplast sequences were used to confirm the identity of drift propagules as R. racemosa. We postulate that a complex interaction of ocean currents, estuarine geomorphology, and tidal patterns explain drift propagule dispersal to an area. Most likely the investigated cohort of propagules originated from more southern mangrove areas of the West African range beyond the Cameroonian border. This study unraveled the allelic, genetic, and genotypic features of stranded propagules following a stochastic long-distance dispersal. Transboundary dispersal of these propagules highlights the need for intergovernmental efforts in the management of biodiversity.

Local environment-driven adaptive evolution in a marine invasive ascidian (Molgula manhattensis)

Elucidating molecular mechanisms of environment-driven adaptive evolution in marine invaders is crucial for understanding invasion success and further predicting their future invasions. Although increasing evidence suggests that adaptive evolution could contribute to organisms' adaptation to varied environments, there remain knowledge gaps regarding how environments influence genomic variation in invaded habitats and genetic bases underlying local adaptation for most marine invaders. Here, we performed restriction-site-associated DNA sequencing (RADseq) to assess population genetic diversity and further investigate genomic signatures of local adaptation in the marine invasive ascidian, Molgula manhattensis. We revealed that most invasive populations exhibited significant genetic differentiation, low recent gene flow, and no signal of significant population bottleneck. Based on three genome scan approaches, we identified 109 candidate loci potentially under environmental selection. Redundancy analysis and variance partitioning analysis suggest that local environmental factors, particularly the salinity-related variables, represent crucial evolutionary forces in driving adaptive divergence. Using the newly developed transcriptome as a reference, 14 functional genes were finally obtained with potential roles in salinity adaptation, including SLC5A1 and SLC9C1 genes from the solute carrier gene (SLC) superfamily. Our findings confirm that differed local environments could rapidly drive adaptive divergence among invasive populations and leave detectable genomic signatures in marine invaders.

Latitudinal variation in climate-associated genes imperils range edge populations

The ecological impacts of increasing global temperatures are evident in most ecosystems on Earth, but our understanding of how climatic variation influences natural selection and adaptive resilience across latitudes remains largely unknown. Latitudinal gradients allow testing general ecosystem-level theories relevant to climatic adaptation. We assessed differences in adaptive diversity of populations along a latitudinal region spanning highly variable temperate to subtropical climates. We generated and integrated information from environmental mapping, phenotypic variation and genome-wide data from across the geographical range of the rainbowfish Melanotaenia duboulayi, an emerging aquatic system for studies of climate change. We detected, after controlling for spatial population structure, strong interactions between genotypes and environment associated with variation in stream flow and temperature. Some of these hydroclimate-associated genes were found to interact within functional protein networks that contain genes of adaptive significance for projected future climates in rainbowfish. Hydroclimatic selection was also associated with variation in phenotypic traits, including traits known to affect fitness of rainbowfish exposed to different flow environments. Consistent with predictions from the "climatic variability hypothesis," populations exposed to extremes of important environmental variables showed stronger adaptive divergence and less variation in climate-associated genes compared to populations at the centre of the environmental gradient. Our findings suggest that populations that evolved at environmental range margins and at geographical range edges may be more vulnerable to changing climates, a finding with implications for predicting adaptive resilience and managing biodiversity under climate change.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Stress responses in fish: From molecular to evolutionary processes

In the context of global changes, fish are increasingly exposed to multiple stressors that have cascading effects from molecules to the whole individual, thereby affecting wild fish populations through selective processes. In this review, we synthetize recent advances in molecular biology and evolutionary biology to outline some potentially important effects of stressors on fish across biological levels. Given the burgeoning literature, we highlight four promising avenues of research. First, (1) the exposure to multiple stressors can lead to unexpected synergistic or antagonistic effects, which should be better taken into account to improve our predictions of the effects of actual and future human activities on aquatic organisms. Second, (2) we argue that such interactive effects might be due to switches in energy metabolism leading to threshold effects. Under multiple stress exposure, fish could switch from a "compensation" strategy, i.e. a reallocation of energy to defenses and repair to a "conservation" strategy, i.e. blocking of stress responses leading to strong deleterious effects and high mortality. Third, (3) this could have cascading effects on fish survival and population persistence but multiscale studies are still rare. We propose emerging tools merging different levels of biological organization to better predict population resilience under multiple stressors. Fourth (4), there are strong variations in sensitivity among populations, which might arise from transgenerational effects of stressors through plastic, genetic, and epigenetic mechanisms. This can lead to local adaptation or maladaptation, with strong impacts on the evolutionary trajectories of wild fish populations. With this review, we hope to encourage future research to bridge the gap between molecular ecology, ecotoxicology and evolutionary biology to better understand the evolution of responses of fishes to current and future multiple stressors in the context of global changes.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

Shared genomic outliers across two divergent population clusters of a highly threatened seagrass

The seagrass, Zostera capensis, occurs across a broad stretch of coastline and wide environmental gradients in estuaries and sheltered bays in southern and eastern Africa. Throughout its distribution, habitats are highly threatened and poorly protected, increasing the urgency of assessing the genomic variability of this keystone species. A pooled genomic approach was employed to obtain SNP data and examine neutral genomic variation and to identify potential outlier loci to assess differentiation across 12 populations across the ∼9,600 km distribution of Z. capensis. Results indicate high clonality and low genomic diversity within meadows, which combined with poor protection throughout its range, increases the vulnerability of this seagrass to further declines or local extinction. Shared variation at outlier loci potentially indicates local adaptation to temperature and precipitation gradients, with Isolation-by-Environment significantly contributing towards shaping spatial variation in Z. capensis. Our results indicate the presence of two population clusters, broadly corresponding to populations on the west and east coasts, with the two lineages shaped only by frequency differences of outlier loci. Notably, ensemble modelling of suitable seagrass habitat provides evidence that the clusters are linked to historical climate refugia around the Last Glacial Maxi-mum. Our work suggests a complex evolutionary history of Z. capensis in southern and eastern Africa that will require more effective protection in order to safeguard this important ecosystem engineer into the future.

From conservation genetics to conservation genomics: a genome-wide assessment of blue whales (Balaenoptera musculus) in Australian feeding aggregations

Genetic datasets of tens of markers have been superseded through next-generation sequencing technology with genome-wide datasets of thousands of markers. Genomic datasets improve our power to detect low population structure and identify adaptive divergence. The increased population-level knowledge can inform the conservation management of endangered species, such as the blue whale (Balaenoptera musculus). In Australia, there are two known feeding aggregations of the pygmy blue whale (B. m. brevicauda) which have shown no evidence of genetic structure based on a small dataset of 10 microsatellites and mtDNA. Here, we develop and implement a high-resolution dataset of 8294 genome-wide filtered single nucleotide polymorphisms, the first of its kind for blue whales. We use these data to assess whether the Australian feeding aggregations constitute one population and to test for the first time whether there is adaptive divergence between the feeding aggregations. We found no evidence of neutral population structure and negligible evidence of adaptive divergence. We propose that individuals likely travel widely between feeding areas and to breeding areas, which would require them to be adapted to a wide range of environmental conditions. This has important implications for their conservation as this blue whale population is likely vulnerable to a range of anthropogenic threats both off Australia and elsewhere.