Marine populations are typically characterized by weak genetic differentiation due to the potential for long‐distance dispersal favouring high levels of gene flow. However, strong directional ...advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2 = .49) than geographic distance (IBD) (R2 = .18), and directional processes played an important role in shaping fine‐scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.
Freshwater eels (Anguilla sp.) have large economic, cultural, ecological and aesthetic importance worldwide, but they suffered more than 90% decline in global stocks over the past few decades. Proper ...genetic resources, such as sequenced, assembled and annotated genomes, are essential to help plan sustainable recoveries by identifying physiological, biochemical and genetic mechanisms that caused the declines or that may lead to recoveries. Here, we present the first sequenced genome of the American eel. This genome contained 305 043 contigs (N50 = 7397) and 79 209 scaffolds (N50 = 86 641) for a total size of 1.41 Gb, which is in the middle of the range of previous estimations for this species. In addition, protein‐coding regions, including introns and flanking regions, are very well represented in the genome, as 95.2% of the 458 core eukaryotic genes and 98.8% of the 248 ultra‐conserved subset were represented in the assembly and a total of 26 564 genes were annotated for future functional genomics studies. We performed a candidate gene analysis to compare three genes among all three freshwater eel species and, congruent with the phylogenetic relationships, Japanese eel (A. japanica) exhibited the most divergence. Overall, the sequenced genome presented in this study is a crucial addition to the presently available genetic tools to help guide future conservation efforts of freshwater eels.
•The transcriptional responses of Perca flavescens to both metal and non metal stressors were measured with a 1000 candidate-gene microarray.•475, 287 and 176 genes were differentially transcribed ...depending on temperature, Ni and Cd concentrations, respectively.•Genes involved in iron metabolism, transcriptional and translational processes, vitamin metabolism, blood coagulation, and calcium transport were impacted by metals.•The developed microarray contributes to a better characterization of the impact of different stressors on the transcriptome.
The molecular mechanisms underlying nickel (Ni) and cadmium (Cd) toxicity and their specific effects on fish are poorly understood. Documenting gene transcription profiles offers a powerful approach toward identifying the molecular mechanisms affected by these metals and to discover biomarkers of their toxicity. However, confounding environmental factors can complicate the interpretation of the results and the detection of biomarkers for fish captured in their natural environment. In the present study, a 1000 candidate-gene microarray, developed from a previous RNA-seq study on a subset of individual fish from contrasting level of metal contamination, was used to investigate the transcriptional response to metal (Ni and Cd) and non metal (temperature, oxygen, and diet) stressors in yellow perch (Perca flavescens). Specifically, we aimed at (1) identifying transcriptional signatures specific to Ni and Cd exposure, (2) investigating the mechanisms of their toxicity, and (3) developing a predictive tool to identify the sublethal effects of Ni and Cd contaminants in fish sampled from natural environments. A total of 475 genes displayed significantly different transcription levels when temperature varied while 287 and 176 genes were differentially transcribed at different concentrations of Ni and Cd, respectively. These metals were found to mainly affect the transcription level of genes involved in iron metabolism, transcriptional and translational processes, vitamin metabolism, blood coagulation, and calcium transport. In addition, a linear discriminant analysis (LDA) made using gene transcription levels yielded 94% correctly reassigned samples regarding their level of metal contamination, which indicates the potential of the microarray to detect perch response to Cd or Ni effects.
In Lar National Park (Caspian Sea basin, Iran), the Caspian trout (Salmo caspius) population faces different threats, including introduced fish species. Due to the harsh environmental conditions and ...limited accessibility, monitoring of fish species via conventional approaches proves difficult. Hence, environmental DNA metabarcoding may prove an appropriate tool for monitoring fishes within the park. Environmental DNA samples from eight stream sites in the National Park were sequenced via metabarcoding of the 12S rRNA gene, and the species identified via eDNA metabarcoding were compared to the results of electrofishing performed at the same localities on the same day. No significant difference in the number of Caspian Sea trout DNA sequence reads was detected among the collection sites (p > 0.05). The highest number of reads was detected in Dalichay Stream, but the highest population density determined via electrofishing was in Siahpalas Stream. The discrepancy between the eDNA read count and trout population density, as well as the limited sampling scheme within this study, limit our ability to provide a robust conclusion about the application of environmental DNA metabarcoding for assessment of fish density in Lar National Park. Environmental DNA metabarcoding detected more species than electrofishing, but no significant differences in the composition of the local fish community were observed. Introduced fish species were all observed or detected in Siahpalas Stream, which is characterized by high water temperature, muddy substrate, and lower flow rate. A significant effect of flow rate and total dissolved solids on the presence of introduced fish species (p = 0.02) and of flow rate alone on relative abundance of introduced fish species (p = 0.03) was detected. To standardize the application of eDNA metabarcoding as a biodiversity assessment tool in Lar National Park, future studies should characterize the parameters that affect eDNA persistence and detectability in the system.
In Lar National Park (Caspian Sea basin, Iran) Caspian Sea trout populations and non‐native fish diversity and distribution was assessed and compared using eDNA metabarcoding and electrofishing. The results revealed that non‐native fish species were distributed in limited habitats in accordance with distribution of suitable environmental factors, but the Caspian Sea trout existed all over the studied region. Overall, eDNA metabarcoding appeared to better represent fish diversity compared to electrofishing.
Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this ...study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site-associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.
Abstract
Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. ...Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.
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.
One of the most important steps in conservation of the subterranean life forms is to decipher their distribution and ecology, which is challenging using traditional approaches. Development of an ...environmental DNA (eDNA) assay provides an efficient means for discovering and monitoring subterranean life forms. In this study, the distribution of three Iranian blind cave fish species (blind Iran cave barb Garra typhlops, blind Lorestan cave barb Garra lorestanensis, and blind cave loach Eidinemacheilus smithi) was assessed using 12S rRNA eDNA metabarcoding performed using MiFish‐U PCR primers and preliminary species distribution modeling (SDM) using bioclimatic data. The majority of sampling localities with positive detection of cave barb eDNA fall within suitable habitats in the Zagros Mountains of Iran. Our results revealed that Lorestan and Iran cave barbs have differential distribution patterns, with some extent of habitat overlap in the vicinity of the originally discovered cave barb locality. According to the observed distribution patterns, the blind Lorestan cave barb and cave loach are mostly distributed in habitats close to the Sezar River (Dez River drainage, Iran), and the blind Iran cave barb is distributed towards the west and probably in a few springs in the Karkheh River drainage. Our data support the previously proposed distribution pattern for the cave barbs, in which the species show partial niche separation and reproductive isolation, with the Lorestan cave barb being a water flow‐dependent species and the Iran cave barb being a generalist species preferring variable flow rates. We showed eDNA metabarcoding to be a useful approach for ecological surveys of subterranean fish biodiversity with implications for conservation.
In this work a subterranean fish community consisting three species in the Zagros Mountains of Iran was studied using eDNA metabarcoding. The results showed differential distribution of the species in the subterranean habitats that is likely in line with previously developed hypotheses about distribution of Iran cave fish species
Large‐scale biomonitoring of Arctic coastal marine communities is essential to track temporal changes in ecosystems. Despite the potential of environmental DNA (eDNA) as an innovative coastal ...biomonitoring tool, important questions remain pertaining to its temporal and spatial variation and how this may affect the evaluation of ecosystem changes over time in hydrodynamic ecosystems. In this study, we used eDNA metabarcoding of coastal water samples in two Canadian Arctic ports to evaluate the potential of eDNA to detect temporal transition in marine coastal communities. We sequenced eDNA from approximately 20 surface water samples collected each month (N ≈ 150 samples) covering the transition period between summer and late fall using four different universal primer pairs (two pairs of COI mitochondrial genes and two pairs of 18S rRNA genes). Our results from both primer pairs highlighted a significant transition from the summer to the fall marine community. We also observed a putative link between eDNA peaks of read abundance and timing for different life stages (e.g., spawning and larvae) of several species with the most abundant sequence reads. As such, our results show that temporal variation must be considered in ensuring comprehensive coastal biomonitoring with eDNA. Although much remains to be investigated about the ecology of eDNA, our results contribute to fundamental knowledge on the origin of eDNA and highlight the importance of considering temporal variation in developing guidance for coastal biomonitoring with this approach.
Despite environmental DNA approach is an innovative biomonitoring tool, knowledge of temporal and spatial eDNA variation need to be improve to obtain reliable interpretation of the monitored communities. In this study, we investigated the temporal transition in eDNA of coastal marine communities between open‐water and ice‐covered periods in two Arctic ports. Our results suggest that ecological processes (i.e., life stage transitions and migration) are the main drivers of eDNA changes overtime, despite the influence of environmental factors on spatio‐temporal eDNA distribution.
Conservation of freshwater biodiversity requires being able to track the presence and abundance of entire fish communities. However, studying fish community composition within rivers remains a ...technical challenge because of high spatial and temporal physico‐chemical variability, anthropic activities and connections with other river catchments, which may all contribute to important variations in local ecology and communities. Here, we used environmental DNA metabarcoding to document spatial variation in fish communities at a small geographic scale in a large river system. The study was conducted in the Contrecoeur sector (5.5 km long and approximately 1–1.5 km wide) of the St. Lawrence River (Québec, Canada), where two water masses with different physico‐chemical properties, known as "brown waters" and "green waters," flow in parallel with limited admixing. Water samples were collected during two consecutive days at 53 stations located in both water masses. Using universal PCR MiFish 12S primers, Illumina MiSeq sequencing, and the Barque (www.github.com/enormandeau/barque) eDNA analysis software developed by our group, a total of 67 fish species were detected. PERMANOVA and redundancy analyses (RDA) performed on relative read abundance revealed that each water mass comprised distinct communities that depended on turbidity, depth, and to a lesser extent on the upstream versus downstream position along the study area. eDNA metabarcoding results were compared with those of traditional surveys conducted previously in the sector and up to 40 km upstream of it. As previously reported, higher species diversity was detected by eDNA and with substantially lower sampling effort. Our results represent one of the few studies documenting the potential of eDNA metabarcoding to investigate small‐scale variation in 2D spatial patterns of distribution of whole fish communities associated with habitat characteristics within a lotic system.
We conducted a eDNA metabarcoding study in the Contrecoeur sector of the St. Lawrence River (Québec, Canada), which is only 5.5 km long and is characterized by two water masses with very different physico‐chemical properties that flow in parallel with limited admixing. We collected water samples during two consecutive days at 53 stations localized in both water masses. Using universal PCR MiFish 12S primers and Illumina MiSeq sequencing, we found variation in fish community composition within each water mass that was significantly associated with turbidity, depth, and to a lesser extent with the upstream‐downstream position along the study area. Our results demonstrate that it is possible to detect lateral (from shore toward center) variation in eDNA, and consequently in fish community composition at a small geographic scale in rivers, and that eDNA metabarcoding can efficiently complement traditional capture methods in order to more accurately document the role of physico‐chemical factors in shaping local fish community composition in large fluvial systems.