Search

Abstract Amazonian blackwaters are extremely biodiverse systems containing some of Earth’s most naturally acidic, dissolved organic carbon -rich and ion‐poor waters. Physiological adaptations of fish facing these ionoregulatory challenges are unresolved but could involve microbially-mediated processes. Here, we characterize the physiological response of 964 fish-microbe systems from four blackwater Teleost species along a natural hydrochemical gradient, using dual RNA-Seq and 16 S rRNA of gill samples. We find that host transcriptional responses to blackwaters are species-specific, but occasionally include the overexpression of Toll-receptors and integrins associated to interkingdom communication. Blackwater gill microbiomes are characterized by a transcriptionally-active betaproteobacterial cluster potentially interfering with epithelial permeability. We explore further blackwater fish-microbe interactions by analyzing transcriptomes of axenic zebrafish larvae exposed to sterile, non-sterile and inverted (non-native bacterioplankton) blackwater. We find that axenic zebrafish survive poorly when exposed to sterile/inverted blackwater. Overall, our results suggest a critical role for endogenous symbionts in blackwater fish physiology.

The black piranha ( Serrasalmus rhombeus ), a widely spread species in the rivers of the Amazon basin, plays a vital role as both key predator and important prey. Despite its essential contribution to ecosystem stability, there is a lack of information regarding its genetic diversity and population dynamics in the central Amazon region. As the Amazon continues to undergo environmental changes in the context of growing anthropogenic threats, such knowledge is fundamental for assist in the conservation of this species. This study is the first to analyze the genetic diversity and population structure of S. rhombeus in the central Amazon region using high-resolution genomic data. We employed a Genotyping-by-Sequencing approach with 248 samples across 14 study sites from various tributaries, encompassing diverse water types (black, white, and clear water) and characterized by 34 physiochemical parameters. The data reveals low diversity accompanied by pronounced signs of inbreeding in half of the sites and robust genetic differentiation and variation among sites and within-sites. Surprisingly, we also found evidence of higher dispersal capacity than previously recognized. Our analysis exposed a complex and high population structure with genetic groups exclusive to some sites. Gene flow was low and some groups presented ambiguous genealogical divergence index ( gdi ) signals, suggesting the occurrence of potential cryptic species. Moreover, our results suggest that the population structure of black piranha appears more influenced by historical events than contemporary factors. These results underscore the need to give greater attention to this keystone species, for which no regulatory framework or conservation strategies is presently in effect., Competing Interests: The authors declare no conflicts of interest., (© 2025 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Winter flounder (Pseudopleuronectes americanus) is a North Atlantic flatfish that inhabits cold-water environments already affected by global warming. Flatfishes are particularly sensitive during their juvenile stage to a phenomenon known as temperature-dependent sex determination (TSD). In this study, we hypothesized that many genes involved in the steroidogenesis pathway are already expressed at the larval stage in winter flounder and that temperature conditions may influence this pathway prior to the juvenile stage, which is usually considered the TSD-sensitive period. We also hypothesized that temperature effects on the steroidogenesis pathway may carry over from the larval to the juvenile stage. We surveyed three cohorts issued from three different spawning events during spring. Rearing temperature and salinity conditions followed seasonal conditions in the St. Lawrence Estuary (Québec, Canada). As a result, cohorts were exposed to different temperature regimes during egg and larval development. Once metamorphosis occurred, all juveniles were held at a stable temperature of 10.1 ± 0.8 °C. Larvae were sampled at 20, 30, and, 40 days post-hatching and juveniles at one and two months post-metamorphosis. RNA was extracted from these individuals and then sequenced. Transcriptome analysis showed clear differentiation between larvae and juveniles. Gene expression patterns showed that the steroidogenesis pathway, including genes involved in gonadal differentiation, was influenced by temperature during larval development, which indicates that sexual determination at the cellular level is an early process in winter flounder development. Carry-over effects of larval rearing temperature conditions were present at the juvenile stage, with an increased number of genes involved in the steroidogenesis pathway being affected. Altogether, our study highlights the important role of larval development and temperature exposure on the transcriptome of winter flounder., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Single nucleotide polymorphism (SNP) panels are powerful tools for assessing the genetic population structure and dispersal of fishes and can enhance management practices for commercial, recreational and subsistence mixed-stock fisheries. Arctic Char (Salvelinus alpinus), Brook Trout (Salvelinus fontinalis) and Lake Whitefish (Coregonus clupeaformis) are among the most harvested and consumed fish species in Northern Indigenous communities in Canada, contributing significantly to food security, culture, tradition and economy. However, genetic resources supporting Indigenous fisheries have not been widely accessible to northern communities (e.g. Inuit, Cree, Dene). Here, we developed Genotyping-in-Thousands by sequencing (GT-seq) panels for population assignment and mixed-stock analyses of three salmonids, to support fisheries stewardship or co-management in Northern Canada. Using low-coverage Whole Genome Sequencing data from 418 individuals across source populations in Cambridge Bay (Nunavut), Great Slave Lake (Northwest Territories), James Bay (Québec) and Mistassini Lake (Québec), we developed a bioinformatic SNP filtering workflow to select informative SNP markers from genotype likelihoods. These markers were then used to design GT-seq panels, thus enabling high-throughput genotyping for these species. The three GT-seq panels yielded an average of 413 autosomal loci and were validated using 525 individuals with an average assignment accuracy of 83%. Thus, these GT-seq panels are powerful tools for assessing population structure and quantifying the relative contributions of populations/stocks in mixed-stock fisheries across multiple regions. Interweaving genomic data derived from these tools with Traditional Ecological Knowledge will ensure the sustainable harvest of three culturally important salmonids in Indigenous communities, contributing to food security programmes and the economy in Northern Canada., (© 2024 The Author(s). Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex ("ciscoes") radiated rapidly prior to the Last Glacial Maximum (80-90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr (Salvelinus namaycush), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation., (© 2024. The Authors. Parts of this work were authored by US Federal Government authors and are not under copyright protection in the US; foreign copyright protection may apply.)

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects., (© 2023 John Wiley & Sons Ltd.)

Plasmodiophora brassicae (Woronin, 1877), a biotrophic, obligate parasite, is the causal agent of clubroot disease in brassicas. The clubroot pathogen has been reported in more than 80 countries worldwide, causing economic losses of hundreds of millions every year. Despite its widespread impact, very little is known about the molecular strategies it employs to induce the characteristic clubs in the roots of susceptible hosts during infection, nor about the mechanisms it uses to overcome genetic resistance. Here, we provide the first telomere-to-telomere complete genome of P. brassicae. We generated ∼27 Gb of Illumina, Oxford Nanopore, and PacBio HiFi data from resting spores of strain Pb3A and produced a 25.3 Mb assembly comprising 20 chromosomes, with an N50 of 1.37 Mb. The BUSCO score, the highest reported for any member of the group Rhizaria (Eukaryota: 88.2%), highlights the limitations within the Eukaryota database for members of this lineage. Using available transcriptomic data and protein evidence, we annotated the Pb3A genome, identifying 10,521 protein-coding gene models. This high-quality, complete genome of P. brassicae will serve as a crucial resource for the plant pathology community to advance the much-needed understanding of the evolution of the clubroot pathogen., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The primary focus of medicinal cannabis research is to ensure the stability of cannabis lines for consistent administration of chemically uniform products to patients. In recent years, tissue culture has emerged as a valuable technique for genetic preservation and rapid multiplication of cannabis clones. However, there is concern that the physical and chemical conditions of the growing media can induce somaclonal variation, potentially impacting the viability and uniformity of clones. To address this concern, we developed Comparative Restriction Enzyme Analysis of Methylation (CREAM), a novel method to assess DNA methylation patterns and used it to study a population of 78 cannabis clones maintained in tissue culture. Through bioinformatics analysis of the methylome, we successfully detected 2,272 polymorphic methylated regions among the clones. Remarkably, our results demonstrated that DNA methylation patterns were preserved across subcultures within the clonal population, allowing us to distinguish between two subsets of clonal lines used in this study. These findings significantly contribute to our understanding of the epigenetic variability within clonal lines in medicinal cannabis produced through tissue culture techniques. This knowledge is crucial for understanding the effects of tissue culture on DNA methylation and ensuring the consistency and reliability of medicinal cannabis products with therapeutic properties. Additionally, the CREAM method is a fast and affordable technology to get a first glimpse at methylation in a biological system. It offers a valuable tool for studying epigenetic variation in other plant species, thereby facilitating broader applications in plant biotechnology and crop improvement., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Boissinot, Adamek, Jones, Normandeau, Boyle and Torkamaneh.)

Brook charr is a cold-water species which is highly sensitive to increased water temperatures, such as those associated with climate change. Environmental variation can potentially induce phenotypic changes that are inherited across generations, for instance, via epigenetic mechanisms. Here, we tested whether parental thermal regimes (intergenerational plasticity) and offspring-rearing temperatures (within-generational plasticity) modify the brain transcriptome of Brook charr progeny (fry stage). Parents were exposed to either cold or warm temperatures during final gonad maturation and their progeny were reared at 5 or 8 °C during the first stages of development. Illumina Novaseq6000 was used to sequence the brain transcriptome at the yolk sac resorption stage. The number of differentially expressed genes was very low when comparing fry reared at different temperatures (79 differentially expressed genes). In contrast, 9,050 differentially expressed genes were significantly differentially expressed between fry issued from parents exposed to either cold or warm temperatures. There was a significant downregulation of processes related to neural and synaptic activity in fry originating from the warm parental group vs fry from the cold parental one. We also observed significant upregulation of DNA methylation genes and of the most salient processes associated with compensation to warming, such as metabolism, cellular response to stress, and adaptive immunity., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

While Atlantic salmon ( Salmo salar ) of the northernmost American populations is alimentary, economically, and culturally important for Ungava Inuit communities (Nunavik, Canada) and might play a key role in the persistence of the species in a global warming context, many mysteries remain about those remote and atypical populations. Thus, our first aim was to document the genomic structure of the Nunavik populations. The second objective was to determine whether salmon only migrating to the estuary without reaching the sea, apparently unique to those populations, represent distinct populations from the typical anadromous salmons and subsequently explore the genetic basis of migratory life-history tactics in the species. Finally, the third goal was to quantify the contribution of each genetically distinct population and life-history tactic in the mixed-stock subsistence fishery of the Koksoak R. estuary. We used Genotyping-by-Sequencing to genotype 14,061 single nucleotide polymorphisms in the genome of 248 individuals from 8 source populations and 280 individuals from the Koksoak estuary mixed-stock fishery. Life-history tactics were identified by a visual assessment of scales. Results show a hierarchical structure mainly influenced by isolation-by-distance with 7 populations out of the 8 studied rivers. While no obvious structure was detected between marine and estuarine salmon within the population, we have identified genomic regions putatively associated with those migration tactics. Finally, all salmon captured in the Koksoak estuary originated from the Koksoak drainage and mostly from 2 tributaries, but no inter-annual variation in the contribution of these tributaries was found. Our results indicate, however, that both marine and estuarine salmon contribute substantially to estuarine fisheries and that there is inter-annual variation in this contribution. These findings provide crucial information for the conservation of salmon populations in a rapidly changing ecosystem, as well as for fishery management to improve the food security of Inuit communities., Competing Interests: We declare that there are no conflicts of interest., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The American eel ( Anguilla rostrata ) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species., Competing Interests: The authors report no conflict of interest., (© 2023 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Differences between sexes in trait fitness optima can generate intralocus sexual conflicts that have the potential to maintain genetic diversity through balancing selection. However, these differences are unlikely to be associated with strong selective coefficients and are challenging to detect. Additionally, recent studies have highlighted that duplications on sexual chromosomes can create artifactual signals of intralocus sexual conflicts. Thus, testing the relationship between intralocus sexual conflicts and balancing selection requires stringent filtering of duplicated regions, and dedicated methods to detect loci with low levels of intersex differentiation. In this study, we investigated intralocus sexual conflicts in the three-spined stickleback using whole-genome sequencing (mean coverage = 12×) of 50 females and 49 males from an anadromous population in the St. Lawrence River, Québec, Canada. After stringent filtering of duplications from the sex chromosomes, we compared three methods to detect intralocus sexual conflicts. We found only two genomic regions under potential intralocus sexual conflict that also showed signals of balancing selection. Overall, our results suggest that most intralocus sexual conflicts do not drive long-term balancing selection and are most likely transient., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

Dense single nucleotide polymorphism (SNP) arrays are essential tools for rapid high-throughput genotyping for many genetic analyses, including genomic selection and high-resolution population genomic assessments. We present a high-density (200 K) SNP array developed for the Eastern oyster (Crassostrea virginica), which is a species of significant aquaculture production and restoration efforts throughout its native range. SNP discovery was performed using low-coverage whole-genome sequencing of 435 F1 oysters from families from 11 founder populations in New Brunswick, Canada. An Affymetrix Axiom Custom array was created with 219,447 SNPs meeting stringent selection criteria and validated by genotyping more than 4,000 oysters across 2 generations. In total, 144,570 SNPs had a call rate >90%, most of which (96%) were polymorphic and were distributed across the Eastern oyster reference genome, with similar levels of genetic diversity observed in both generations. Linkage disequilibrium was low (maximum r2 ∼0.32) and decayed moderately with increasing distance between SNP pairs. Taking advantage of our intergenerational data set, we quantified Mendelian inheritance errors to validate SNP selection. Although most of SNPs exhibited low Mendelian inheritance error rates overall, with 72% of called SNPs having an error rate of <1%, many loci had elevated Mendelian inheritance error rates, potentially indicating the presence of null alleles. This SNP panel provides a necessary tool to enable routine application of genomic approaches, including genomic selection, in C. virginica selective breeding programs. As demand for production increases, this resource will be essential for accelerating production and sustaining the Canadian oyster aquaculture industry., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

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., Competing Interests: Conflicts of interest statement The authors declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Genetics Society of America.)

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Captive rearing in salmon hatcheries can have considerable impacts on both fish phenotype and fitness within a single generation, even in the absence of genetic change. Evidence for hatchery-induced changes in DNA methylation is becoming abundant, though questions remain on the sex-specificity of these effects, their persistence until spawning and potential for transmission to future generations. Here we performed whole genome methylation sequencing of fin tissue for 16 hatchery and 16 wild Atlantic salmon (Salmo salar) returning to spawn in the Rimouski River, Québec, Canada. We identified two cohorts of hatchery-reared salmon through methylation analysis, one of which was epigenetically similar to wild fish, suggesting that supplementation efforts may be able to minimize the epigenetic effects of hatchery rearing. We found considerable sex-specific effects of hatchery rearing, with few genomic regions being affected in both males and females. We also analysed the methylome of 32 F 1 offspring from four groups (pure wild, pure hatchery origin and reciprocal hybrids). We found that few epigenetic changes due to parental hatchery rearing persisted in the F 1 offspring though the patterns of inheritance appear to be complex, involving nonadditive effects. Our results suggest that the epigenetic effects of hatchery rearing can be minimal in F 0 . There may also be minimal epigenetic inheritance and rapid loss of epigenetic changes associated with hatchery rearing. However, due to sex-specificity and nonadditive patterns of inheritance, methylation changes due to captive rearing are rather complex and the field would benefit from further research on minimizing the epigenetic effects of captive rearing in conservation efforts., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications., (© 2022 John Wiley & Sons Ltd.)

Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation., (© The Author(s) 2022. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon ( Oncorhynchus kisutch ) from 149 sampling locations and a total of 24,542 high-quality SNPs obtained using genotyping-by-sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype-environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well-suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined., Competing Interests: We declare no conflict of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Associations between host genotype and host-associated microbiomes have been shown in a variety of animal clades, but studies on teleosts mostly show weak associations. Our study aimed to explore these relationships in four sympatric Serrasalmidae (i.e., piranha) teleosts from an Amazonian lake, using data sets from the hosts genomes (single nucleotide polymorphisms from genotyping by sequencing), skin and gut microbiomes (16S rRNA gene metataxonomics) and diets (COI metabarcoding) from the same fish individuals. First, we investigated whether there were significant covariations of microbiome and fish genotypes at the inter- and intraspecific levels. We also assessed the extent of covariation between Serrasalmidae diet and microbiome, to isolate genotypic from dietary effects on community structure. We observed a significant covariation of skin microbiomes and host genotypes at interspecific (R 2  = 24.4%) and intraspecific (R 2  = 6.2%) levels, whereas gut microbiomes correlated poorly with host genotypes. Serrasalmidae diet composition was significantly correlated to fish genotype only at the interspecific level (R 2  = 5.4%), but did not covary with gut microbiome composition (Mantel R = -.04). Second, we investigated whether the study of interspecific differentiation could benefit from considering host-associated microbial communities in addition to host genotypes. By using a nonmetric multidimensional scaling (NMDS) ordination-based approach, we observed that ordinations from skin- and gut species-specific bacterial biomarkers identified through a random forest algorithm could significantly increase the average interspecific differentiation detected through host genotype data alone. Although future studies encompassing additional species and environments are needed, our results suggest Serrasalmidae microbiomes could constitute an insightful trait to be considered when studying the interspecific differences between members of this clade., (© 2022 John Wiley & Sons Ltd.)

Epigenetic inheritance can result in plastic responses to changing environments being faithfully transmitted to offspring. However, it remains unclear how epigenetic mechanisms such as DNA methylation can contribute to multigenerational acclimation and adaptation to environmental stressors. Brook charr ( Salvelinus fontinalis ), an economically important salmonid, is highly sensitive to thermal stress and is of conservation concern in the context of climate change. We studied the effects of temperature during parental sexual maturation and offspring rearing on whole-genome DNA methylation in brook charr juveniles (fry). Parents were split between warm and cold temperatures during sexual maturation, mated in controlled breeding designs, then offspring from each family were split between warm (8°C) and cold (5°C) rearing environments. Using whole-genome bisulfite sequencing, we found 188 differentially methylated regions (DMRs) due to parental maturation temperature after controlling for family structure. By contrast, offspring rearing temperature had a negligible effect on offspring methylation. Stable intergenerational inheritance of DNA methylation and minimal plasticity in progeny could result in the transmission of acclimatory epigenetic states to offspring, priming them for a warming environment. Our findings have implications pertaining to the role of intergenerational epigenetic inheritance in response to ongoing climate change.

The island syndrome hypothesis (ISH) stipulates that, as a result of local selection pressures and restricted gene flow, individuals from island populations should differ from individuals within mainland populations. Specifically, island populations are predicted to contain individuals that are larger, less aggressive, more sociable, and that invest more in their offspring. To date, tests of the ISH have mainly compared oceanic islands to continental sites, and rarely smaller spatial scales such as inland watersheds. Here, using a novel set of genome-wide SNP markers in wild deer mice (Peromyscus maniculatus) we conducted a genomic assessment of predictions underlying the ISH in an inland riverine island system: analysing island-mainland population structure, and quantifying heritability of phenotypes thought to underlie the ISH. We found clear genomic differentiation between the island and mainland populations and moderate to high marker-based heritability estimates for overall variation in traits previously found to differ in line with the ISH between mainland and island locations. F ST outlier analyses highlighted 12 loci associated with differentiation between mainland and island populations. Together these results suggest that the island populations examined are on independent evolutionary trajectories, the traits considered have a genetic basis (rather than phenotypic variation being solely due to phenotypic plasticity). Coupled with the previous results showing significant phenotypic differentiation between the island and mainland groups in this system, this study suggests that the ISH can hold even on a small spatial scale., (© 2021. The Author(s), under exclusive licence to The Genetics Society.)

Here, we present an annotated, chromosome-anchored, genome assembly for Lake Trout (Salvelinus namaycush) - a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long-read sequencing, paired-end Illumina sequencing, proximity ligation (Hi-C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome-sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k-mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self-versus-self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex-averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed., (© 2021 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution., (© Crown copyright 2021.)

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F-measure, root-mean-square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F-measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study., (© 2021 John Wiley & Sons Ltd.)

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations., (Copyright © 2021. Published by Elsevier Inc.)

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon ( Oncorhynchus kisutch ) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Increasing evidence shows that structural variants represent an overlooked aspect of genetic variation with consequential evolutionary roles. Among those, copy number variants (CNVs), including duplicated genomic regions and transposable elements (TEs), may contribute to local adaptation and/or reproductive isolation among divergent populations. Those mechanisms suppose that CNVs could be used to infer neutral and/or adaptive population genetic structure, whose study has been restricted to microsatellites, mitochondrial DNA and Amplified fragment length polymorphism markers in the past and more recently the use of single nucleotide polymorphisms (SNPs). Taking advantage of recent developments allowing CNV analysis from RAD-seq data, we investigated how variation in fitness-related traits, local environmental conditions and demographic history are associated with CNVs, and how subsequent copy number variation drives population genetic structure in a marine fish, the capelin (Mallotus villosus). We collected 1538 DNA samples from 35 sampling sites in the north Atlantic Ocean and identified 6620 putative CNVs. We found associations between CNVs and the gonadosomatic index, suggesting that six duplicated regions could affect female fitness by modulating oocyte production. We also detected 105 CNV candidates associated with water temperature, among which 20% corresponded to genomic regions located within the sequence of protein-coding genes, suggesting local adaptation to cold water by means of gene sequence amplification. We also identified 175 CNVs associated with the divergence of three previously defined parapatric glacial lineages, of which 24% were located within protein-coding genes, making those loci potential candidates for reproductive isolation. Lastly, our analyses unveiled a hierarchical, complex CNV population structure determined by temperature and local geography, which was in stark contrast to that inferred based on SNPs in a previous study. Our findings underline the complementarity of those two types of genomic variation in population genomics studies., (© 2021 John Wiley & Sons Ltd.)

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus . We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus . Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter- and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units., (© 2020 Her Majesty the Queen in Right of Canada. Evolutionary Applications published by John Wiley & Sons Ltd. Reproduced with the permission of the Minister of Fishers and Ocean.)

Copy number variants (CNVs) are a major component of genotypic and phenotypic variation in genomes. To date, our knowledge of genotypic variation and evolution has largely been acquired by means of single nucleotide polymorphism (SNPs) analyses. Until recently, the adaptive role of structural variants (SVs) and particularly that of CNVs has been overlooked in wild populations, partly due to their challenging identification. Here, we document the usefulness of Rapture, a derived reduced-representation shotgun sequencing approach, to detect and investigate copy number variants (CNVs) alongside SNPs in American lobster (Homarus americanus) populations. We conducted a comparative study to examine the potential role of SNPs and CNVs in local adaptation by sequencing 1,141 lobsters from 21 sampling sites within the southern Gulf of St. Lawrence, which experiences the highest yearly thermal variance of the Canadian marine coastal waters. Our results demonstrated that CNVs account for higher genetic differentiation than SNP markers. Contrary to SNPs, for which no significant genetic-environment association was found, 48 CNV candidates were significantly associated with the annual variance of sea surface temperature, leading to the genetic clustering of sampling locations despite their geographic separation. Altogether, we provide a strong empirical case that CNVs putatively contribute to local adaptation in marine species and unveil stronger spatial signal of population structure than SNPs. Our study provides the means to study CNVs in nonmodel species and highlights the importance of considering structural variants alongside SNPs to enhance our understanding of ecological and evolutionary processes shaping adaptive population structure., (© 2020 John Wiley & Sons Ltd.)

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored., Competing Interests: The authors have declared that no competing interests exist.

Gene flow has tremendous importance for local adaptation, by influencing the fate of de novo mutations, maintaining standing genetic variation and driving adaptive introgression. Furthermore, structural variation as chromosomal rearrangements may facilitate adaptation despite high gene flow. However, our understanding of the evolutionary mechanisms impending or favouring local adaptation in the presence of gene flow is still limited to a restricted number of study systems. In this study, we examined how demographic history, shared ancestral polymorphism, and gene flow among glacial lineages contribute to local adaptation to sea conditions in a marine fish, the capelin (Mallotus villosus). We first assembled a 490-Mbp draft genome of M. villosus to map our RAD sequence reads. Then, we used a large data set of genome-wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,310 individuals collected from 31 spawning sites in the northwest Atlantic. We reconstructed the history of divergence among three glacial lineages and showed that they probably diverged from 3.8 to 1.8 million years ago and experienced secondary contacts. Within each lineage, our analyses provided evidence for large N e and high gene flow among spawning sites. Within the Northwest Atlantic lineage, we detected a polymorphic chromosomal rearrangement leading to the occurrence of three haplogroups. Genotype-environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Our study also suggests that both shared polymorphisms among lineages, resulting from standing genetic variation or introgression, and chromosomal rearrangements may contribute to local adaptation in the presence of high gene flow., (© 2020 John Wiley & Sons Ltd.)

Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next-generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut ( Reinhardtius hippoglossoides ) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex-linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland ( F ST  = 0.00146, p -value = .001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year-to-year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex-linked markers and temporal replicates in studies of population genomics., Competing Interests: None declared., (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

How natural diversity is maintained is an evolutionary puzzle. Genetic variation can be eroded by drift and directional selection but some polymorphisms persist for long time periods, implicating a role for balancing selection. Here, we investigate the maintenance of a chromosomal inversion polymorphism in the seaweed fly Coelopa frigida. Using experimental evolution and quantifying fitness, we show that the inversion underlies a life-history trade-off, whereby each haplotype has opposing effects on larval survival and adult reproduction. Numerical simulations confirm that such antagonistic pleiotropy can maintain polymorphism. Our results also highlight the importance of sex-specific effects, dominance and environmental heterogeneity, whose interaction enhances the maintenance of polymorphism through antagonistic pleiotropy. Overall, our findings directly demonstrate how overdominance and sexual antagonism can emerge from a life-history trade-off, inviting reconsideration of antagonistic pleiotropy as a key part of multi-headed balancing selection processes that enable the persistence of genetic variation.

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost-effective sequencing approaches, namely Pool-seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping-by-sequencing (GBS) approach. For this, we used the American lobster ( Homarus americanus ) as a case study. First, we found that GBS, Rapture, and Pool-seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool-seq data showed F ST estimates, indicating that individual-based approaches provided more congruent results than Pool-seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost-effective in project involving a relatively small numbers of genotyped individuals (e.g., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates.F ST estimates, indicating that individual-based approaches provided more congruent results than Pool-seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost-effective in project involving a relatively small numbers of genotyped individuals (e.g., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates., Competing Interests: None declared.

Because significant global changes are currently underway in the Arctic, creating a large-scale standardized database for Arctic marine biodiversity is particularly pressing. This study evaluates the potential of aquatic environmental DNA (eDNA) metabarcoding to detect Arctic coastal biodiversity changes and characterizes the local spatio-temporal distribution of eDNA in two locations. We extracted and amplified eDNA using two COI primer pairs from ~80 water samples that were collected across two Canadian Arctic ports, Churchill and Iqaluit, based on optimized sampling and preservation methods for remote regions surveys. Results demonstrate that aquatic eDNA surveys have the potential to document large-scale Arctic biodiversity change by providing a rapid overview of coastal metazoan biodiversity, detecting nonindigenous species, and allowing sampling in both open water and under the ice cover by local northern-based communities. We show that DNA sequences of ~50% of known Canadian Arctic species and potential invaders are currently present in public databases. A similar proportion of operational taxonomic units was identified at the species level with eDNA metabarcoding, for a total of 181 species identified at both sites. Despite the cold and well-mixed coastal environment, species composition was vertically heterogeneous, in part due to river inflow in the estuarine ecosystem, and differed between the water column and tide pools. Thus, COI-based eDNA metabarcoding may quickly improve large-scale Arctic biomonitoring using eDNA, but we caution that aquatic eDNA sampling needs to be standardized over space and time to accurately evaluate community structure changes.

Large chromosomal rearrangements are thought to facilitate adaptation to heterogeneous environments by limiting genomic recombination. Indeed, inversions have been implicated in adaptation along environmental clines and in ecotype specialization. Here, we combine classical ecological studies and population genetics to investigate an inversion polymorphism previously documented in Europe among natural populations of the seaweed fly Coelopa frigida along a latitudinal cline in North America. We test if the inversion is present in North America and polymorphic, assess which environmental conditions modulate the inversion karyotype frequencies, and document the relationship between inversion karyotype and adult size. We sampled nearly 2000 flies from 20 populations along several environmental gradients to quantify associations of inversion frequencies to heterogeneous environmental variables. Genotyping and phenotyping showed a widespread and conserved inversion polymorphism between Europe and America. Variation in inversion frequency was significantly associated with environmental factors, with parallel patterns between continents, indicating that the inversion may play a role in local adaptation. The three karyotypes of the inversion are differently favoured across micro-habitats and represent life-history strategies likely to be maintained by the collective action of several mechanisms of balancing selection. Our study adds to the mounting evidence that inversions are facilitators of adaptation and enhance within-species diversity., (© 2018 The Author(s).)